TBIA 生物開放資料數據松 決賽簡報
白鷺鷥去哪了？
臺灣原生鷺的時空分布＆族群變遷
觀測紀錄與數量是否必然相關
老祖宗保佑我 ｜ 梁凱皓、陳彥均、盧姸諠
2026.08.27 ｜ 國立臺灣大學
資料與方法
以 TBIA 四種原生鷺共 1,390,243 筆有效紀錄建立「物種 × 5 公里網格 × 年 × 月」表；
核心分析期間定為 1995–2025，因 1995 年前 91 年僅佔全部紀錄 1.1%。
資料來源 A1、A12
早期較為零散
1904–1994
15,857 筆
佔比 1.1%
核心長期分析期間
1995–2025
1,334,170 筆
佔比 96.0%
尚未完整
2026（僅 1–5 月）
40,216 筆
佔比 2.9%
資料清理原則
依資料內物種名稱重新辨識，不依檔名日期整理為 year 與 month排除日期或座標無法解析、超出臺灣範圍者保留 individualCount 作輔助指標保留 recordedBy 供後續觀測活動診斷
研究物種與最早紀錄
大白鷺、中白鷺、小白鷺、夜鷺夜鷺與小白鷺最早可追溯至 1904 年大白鷺 1968 年、中白鷺 1974 年空間單位為 5 公里網格，時間保留到月
季節節奏
四種原生鷺月別振幅由 2.6 倍至 1.3 倍不等且峰值月份不同，證實季節本身即造成出現水準差異，故後續跨期比較必須控制月份組成。
資料來源 A1
指數＝各月 Occupied Grid 平均 ÷ 該物種 12 個月平均（1974–2025）
振幅與峰值月份
大白鷺
2.6×
峰值 1 月
中白鷺
2.2×
峰值 1 月
小白鷺
1.5×
峰值 3 月
夜鷺
1.3×
峰值 5 月
四種共同低點落在 6–8 月
前後期月份組成若不同，季節差異會被誤判為長期變化；此為後續要求三階段涵蓋完整 12 個月的直接理由
長期空間變化
以出現年份比例比較前後各 15 年，增加訊號網格為大白鷺 63、中白鷺 45、夜鷺 34、小白鷺 0；四種加總的下降訊號僅 3 格，單向性本身即指向紀錄擴張而非族群變動。
資料來源 A1、A2
期間與指標定義
前期 1995–2009、後期 2011–2025，各 15 年2010 年為過渡年，不納入比較指標為該物種有紀錄年份 ÷ 該網格有有效資料年份用比例而非筆數，以抵銷各地被觀察年數不同
網格層級變化訊號
物種 | 增加 | 穩定 | 下降 |
大白鷺 | 63 | 418 | 1 |
中白鷺 | 45 | 436 | 1 |
夜鷺 | 34 | 448 | 0 |
小白鷺 | 0 | 481 | 1 |
增加訊號分散在局部，相鄰網格分析未發現大型連續區塊
真實族群變化不可能四種、六千餘格幾乎全無下降，故此處僅能視為「出現訊號變化」
觀測偏差
年度觀察者中位數由 26 人增至 2,674 人，且與 Occupied Grid 的 Spearman 相關達 0.72–0.82；前頁增加無法與觀測擴張分離，而 TBIA 缺乏調查努力量資料亦無法正式校正。
資料來源 A1、A3
年度觀察者中位數 1995–1999
26 人
年度觀察者中位數 2020–2025
2,674 人
訊號與觀測活動 Spearman 相關
0.72–0.82
重要界線
此診斷用於提醒偏差，而非宣稱已校正偏差。TBIA 為 presence only 資料，缺完整 checklist、調查時數與未偵測紀錄，做不到 effort correction。
資料空缺
1,551 個 5 公里網格中僅 482 格（31.08%）具備前後期可比較性，其餘 68.92% 屬資料空缺；換算為物種 × 網格則為 6,204 組合中 1,928 個可比較。
資料來源 A2
68.92%
1,069 個網格（＝4,276 個物種×網格組合）無法判讀
判定門檻
前後任一期有效年份少於 5 年，即判定為無法進行長期比較。判準以網格的有效觀測年份計算，四種原生鷺共用同一組判定。
空間結構
可比較區集中於西部平原、沿海、河口與都會濕地；東部、內陸山區與部分離島特別缺乏。
灰色不代表沒有鳥，而是代表目前資料數不足以判斷
埃及聖鹮與固定樣本
以 2011–2017 聖鹮歷史強度為暴露變數，四種原生鷺後期變化初步皆正向顯著；但收緊為 95 格固定樣本後，中白鷺與夜鷺部分不再顯著，小白鷺係數縮至三分之一。
資料來源 A4、A5、A6
2011–2017 聖鹮歷史紀錄強度
強度定義與三階段
強度＝有聖鹮紀錄月份 ÷ 同網格有原生鷺紀錄月份連續值不切組；分母使用原生鷺觀測月份，使強度不因當地觀察者多而虛高移除前 2011–2017、密集移除期 2018–2020、殘存監測期 2021–2025
固定樣本篩選
1,053
→
368
→
197
→
124
→
95
條件為每年皆有觀測，且三階段各涵蓋完整 12 個月；最終 95 格中 61 格有聖鹮歷史紀錄。固定樣本後中白鷺密集移除期 p=0.140、夜鷺殘存監測期 p=0.411 皆不再顯著，小白鷺係數由 0.055 降至 0.018。
原本的顯著性有一部分來自觀測覆蓋差異，故需再做一道更嚴格的檢驗
事件前
以同一批 95 格、僅用移除前資料比較 2011–2014 與 2015–2017 之假事件，四種原生鷺仍全部正向顯著，證明關聯早於事件本身。
資料來源 A7
檢驗邏輯 若正向關聯真由移除造成，則在移除發生之前，不應看到同樣的關聯
歷史聖鹮強度與原生鷺偵測率變化的迴歸係數
安慰劑結果 95 格固定樣本
物種 | 係數 | p 值 |
大白鷺 | 0.101 | 0.004 |
中白鷺 | 0.135 | 0.028 |
小白鷺 | 0.059 | <0.001 |
夜鷺 | 0.109 | 0.007 |
聖鹮多的地方本來就是濕地與河口，本來就是原生鷺棲地，也本來就是賞鳥熱點；三者共享同一組地理條件。
兩者為長期空間共域而非因果。開放資料分析應主動證明關聯可能不成立，而非找到好看的關聯就宣布因果
淡水河個案與結論
19 個穩定樣點顯示小白鷺六格全部同向，而夜鷺冬季都市河岸反增 130%；資料設計足夠時即可看見物種特異的棲地利用。
資料來源 A9、A10、A12
樣點設計與指標
保留 2017–2025 每年皆有調查的 19 個穩定樣點自然河口 4、人工營造 9、都市河岸 6排除 2023 年後新增的 6 個樣點，確保前後期對等指標為每次有效樣線調查平均個體數
限制
棲地為依樣點環境背景的固定描述性分類，非土地利用圖資，屬探索性結果，不可外推全臺。
關鍵結果
小白鷺
春冬兩季在自然河口與人工濕地皆上升、都市河岸皆下降，六格全部同向
夜鷺（反例）
冬季都市河岸增加 130%，四種鷺並非同進同退
全臺僅 95 個網格能通過最基本的可比較性門檻
需要填補的不只是地圖上的空白，而是缺少固定、跨季節、可比較的調查設計所造成的理解空白
附錄一 資料引用
TBIA 開放資料下載（ARK 永久識別碼）
ARK | 內容 |
ark:/35232/d176h8xl | 原生鷺物種下載 |
ark:/35232/d17jigas | 原生鷺物種下載 |
ark:/35232/d1wm4325 | 原生鷺物種下載 |
ark:/35232/d1t3tqay | 原生鷺物種下載 |
ark:/35232/d1a5sxi0 | 埃及聖鹮（44,327 筆，涵蓋 129 個資料集） |
解析網址 https://pid.tbiadata.tw/
非 TBIA 來源 淡水河流域同步鳥類調查（Tamsui River Bird Count）
提供者 Feathered Forest Ecological Work Co. Ltd.；資料集 UUID 0854e938-5718-4fbe-8a26-26b817fdc713，版本 v1.12；來源 TaiBIF IPT（https://ipt.taibif.tw）。
附錄二 主張與來源對照
編號 | 主張 | 來源檔案 |
A1 | 1,390,243 筆有效紀錄、三段期間佔比、月別季節指數 | grid_month_species.csv |
A2 | 6,204 組合、1,928 可比較、增加訊號 63/45/34/0、下降共 3 | native_heron_change_signals.csv |
A3 | 觀察者 26→2,674、Spearman 0.72–0.82 | native_heron_analysis.ipynb |
A4 | 聖鹮歷史強度、三階段強度 0.290/0.418/0.144 | ibis_historical_intensity_by_grid.csv |
A5 | 主分析四物種係數與 p 值 | ibis_native_heron_association_models.csv |
A6 | 固定樣本 1,053→95 格、中白鷺 p=0.140、夜鷺 p=0.411 | ibis_native_heron_balanced_panel_models.csv |
A7 | 安慰劑係數 0.101/0.135/0.059/0.109 | ibis_native_heron_strict_panel_placebo.csv |
A9 | 淡水河修正版前後期數值（19 穩定樣點，含夜鷺） | trbc_egret_absolute_change_summary_fixed.csv |
A10 | 19 個穩定樣點名單與排除的 6 個樣點 | trbc_transect_habitat_classification_fixed.csv |