南部样流域:2010、2015、2020、2025 防风固沙空间条件试算 v1 本次已经用真实缓存和新增 MODIS 资料运行四期空间计算,并通过独立数值复核。范围是南部 HYBAS_ID 4121367960 样流域,边界交叠面积 201.5184334648 km²,250 m 计算网格为 74 列×91 行。本范围用于验证数据处理和计算流程;它不代表宁夏中北部风沙区,也不是宁夏全区正式生态服务成果。 方法是区域化 RWEQ 经验公式,不是 InVEST 内置模型。没有把缺月填成裸地、零损失或其他年份观测,没有依据本次结果调整参数。水源涵养、土壤保持、碳贮存、生境质量的既有运行与默认参数保持原版本。本地网站仍为 v0.8,本批附件尚未加入网站目录或发布到公网。 先看哪些文件 1. 植被资料完整性_四期.png:说明每月可用范围和完整名义合成时段的覆盖范围。图中“天数”是通过筛选的16日合成资料与月份的重叠时段,不是实际逐日清晰观测数。 2. 防风固沙条件年值_四期共同范围.png:只显示四期共同有效的744格、44.7894629777 km²,占样流域22.23%。四幅图使用相同对数色标;条件模型零值单独浅色,灰色区域不显示年值。不能把灰色解释为零服务量。 3. 2010/2015/2020/2025_区间差值_条件年值_t_ha.tif:四张年度栅格,单位t/(ha·年),NoData=-9999。它们保留各年自身的全部有效范围,不只保留共同范围;因此各年栅格有效面积不同。文件是已复核结果的逐字节副本。 4. 参数与适用范围.txt:本次实际使用的参数、单位、证据定位和未率定假设。 5. vegetation_summary_float64.json、model_display_summary.json:图示数值、有效面积、源文件与SHA256记录。完整验收与文件清单见manifest.json。 如何理解本次结果 主候选使用有当年资料支持的月代表FVC、显式日P/S异常处理预算及原物理雪厚月阈值分类。年值必须31个期全部有效。区间潜在—实际风蚀差值是本包的主要对照指标;局地潜在、局地实际、带符号局地差值及区间潜在、区间实际均另存,不混用两种定义。 年份 各年年值有效面积(km²) 共同范围均值[t/(ha·年)] 2010 198.936216 0.0141177509162 2015 196.790902 0.00000129585294019 2020 49.971507 0.000000118783824671 2025 178.465226 0.00513568608582 这些均值均来自相同44.7894629777 km²,而非各年自身不同有效范围。数值精度用于重现计算,不代表同等科学精度。未估算完整201.52 km²的年度总量,也不外推宁夏全区。 要求每月全部名义合成时段都通过植被QA的分支,年度有效面积依次仅0.25、0、0、1.807086 km²,四年共同范围为空。进一步要求严格日通量和非负原始物理雪深的分支,四年年度均无有效像元。相应年度统计为null,栅格为NoData,绝不能当作零。 年际差异为何这么大 独立检查从原始小时档案按字段名重算风速、风核、日降水/降雪处理和分期气候因子,没有发现数组轴、年份或单位接错。在共同范围按原生气候格所占面积加权,超2m起沙阈值5m/s的小时数四期约为55.260、12.726、4.332、111.854;这是空间加权小时数,不是观测站的整数记录。 风核同时受超阈值幅度影响,因此2025虽有更多超阈值小时,风核不一定高于2010。年分期Wf之和约为22.6524、0.461271、0.029009、15.2302;加入空气密度、重力、湿润与积雪门后,WF之和约为1.510108、0.00823376、0.00216546、1.389112。经验损失公式还会进一步放大小因子的差异。 共同范围的2010年结果约99.725349%集中于4月16—30日这一个计算期,说明个别强风时段具有很大的控制作用。以上解释的是该输入和参数下的算法响应,不是对真实风蚀变化的观测验证。尤其不能据此得出“2015或2020生态改善了几个数量级”。后续需核验强风尾部代表性、适用区域与当地风蚀观测。 数据与处理范围 • 气象:复用四年真实ERA5-Land全年逐时缓存,精确识别9个原生0.1°气候格后按目标中心映射。2025源经度分支的行列号转换已单独核对。250m网格不代表气候具有250m信息量。 • 植被:四年每年23景MOD13Q1.061,新增4个原始71波段Int32栅格;包括NDVI、SummaryQA、DetailedQA和参考原生像元坐标。共取回2,451,825字节原始TIFF,大小、MD5、前后源版本和格网均已检查。云端导出保留,未删除。原生约231.656m格距,产品名义分辨率250m;目标250m采用最近邻采样。 • QA:SummaryQA=0,排除邻云、混合云、可能雪冰与阴影,NDVI原值范围[-2000,10000]后乘0.0001。以合成名义时段重叠天数加权,最后一期截到次年1月1日;不补缺。 • 土壤:SoilGrids 0–5cm砂/粉/黏粒和SOC,HWSD2 D1/SEQ1 0–20cm碳酸钙仅作表层代理。原有8格、0.5km²土壤缺测保留;不是新增年度土壤观测。 • 地形:原生SRTM 1角秒DEM,显式本地Horn 3×3与UTM距离计算坡度,先cos再按投影后原生像元交叠面积聚合。K′是地形代理,非实测空气动力或耕作粗糙度;与旧GEE四邻域算法不同,未悄然替换旧结果。 • 土地覆盖:四年同版CLCD,逐年读取。静态土壤、碳酸钙和地形在四期间共享。 验证完成到什么程度 植被:所有原始和派生栅格、NPZ的网格、mask、时间权重、FVC与支持天数独立复算通过;另一次GEE原生像元抽查的828原值、828个mask及552个坐标全部一致。这些远端控制样本全部有效,没有抽到服务器masked样本,不把这次抽查扩大解释成对所有缺测状态的远端验证。 静态输入:所有土壤格重算、561750个原生坡度/cos重算和25个目标像元交叠加权重算通过。数值闭合不代表地理精度或实测粗糙度验证。 气象及模型:独立从原始小时量复算日处理和分期因子;另按直接Q/s公式复算约396万个有限期值、完整性、共同范围、汇总及72张年度TIFF,均通过。没有新增InVEST原生运行。本批实际没有负局地差样本;负值保留是公式及代码规则,不宣称已用真实负样本验证。 原始vegetation_qa.json中的面积采用Float32求和,总面积与Float64累计相差不足1m²;原记录保留。展示摘要与最终模型统计使用原fraction值的Float64累计,像元值和数量不受影响。 完整来源与文件入口(相对项目根目录) data/wind_spatial_trial_v1/vegetation/manifest.json data/wind_spatial_trial_v1/static/manifest.json data/wind_spatial_trial_v1/weather/actual/report.json data/wind_spatial_trial_v1/model/actual/report.json data/wind_spatial_trial_v1/weather_driver_diagnostic/report.json docs/wind-spatial-v1/final-acceptance.json 完整模型目录有4年份×3输入分支×6指标=72张年度栅格,以及31期、12月和年度数组。其中严格分支的空栅格也计入文件数,并不等于72份都有有效年度结果。生产报告保持运行时“待独审”字样以保存原字节,后续验收记录通过路径和SHA256绑定它。 本次水文独立推进 已明确Han盐池2015观测站同年同址ERA5对照合同,目标原生格r523/c2872,与既有盐池2020格和南部2015缓存不同,不能借用冒充。合同在docs/han2015-era5-v1/acquisition-and-comparison-contract.json;本轮尚未采集该目标格气象。观测时区及半小时区间含义未清,后续即使完成取数,也须保留时间对齐假设,不能直接率定Kc。 数据产品说明: MOD13Q1.061 https://developers.google.com/earth-engine/datasets/catalog/MODIS_061_MOD13Q1 ERA5-Land Hourly https://developers.google.com/earth-engine/datasets/catalog/ECMWF_ERA5_LAND_HOURLY 各原始数据许可、署名及修改说明按原清单保留;本包没有替第三方资料重新授予统一许可。