基于SWAT模型的森林分布不连续流域水源涵养量多时间尺度分析

为解决森林分布不连续流域森林水源涵养功能及其多时间尺度特征的定量评价问题,根据分布式水文模型（SWAT）的特点,提出了反映森林斑块空间分布的水文响应单元划分方法,以及基于水量平衡法的森林不连续分布流域森林水源涵养量计算公式。以东南沿海的晋江流域为例,构建了2006年土地利用条件下的日时间步长SWAT模型,统计分析了2002—2010年降水条件下森林水源涵养量的时空变化规律。结果表明：① 构建的晋江流域SWAT模型精度较高,面积阈值为零生成的水文响应单元比较准确地反映流域森林斑块分布,提出的森林水源涵养量计算公式适用于森林空间分布不连续流域森林水源涵养量的多时间尺度分析,为流域森林水源涵养功能评价提供了一个新的方法。② 晋江流域森林水源年涵养量271.41~565.25 mm;月涵养量-29.15~154.59 mm;日尺度的极端降水期皆为正值,极端枯水期都为负值。表明年际之间不存在森林水源涵养的蓄丰补枯调节作用,但在年内的部分月份得到体现,而日尺度的森林蓄丰补枯功能充分发挥。从而揭示了不同时间尺度森林水源涵养量及其蓄丰补枯功能的差异。     

Estimating fractional cover of crop,crop residue,and soil in cropland using broadband remote sensing data and machine learning

The fractional vegetation cover (FVC), crop residue cover (CRC), and bare soil (BS) are three important parameters in vegetation–soil ecosystems, and their correct and timely estimation can improve crop monitoring and environmental monitoring. The triangular space method uses one CRC index and one vegetation index to create a triangular space in which the three vertices represent pure vegetation, crop residue, and bare soil. Subsequently, the CRC, FVC, and BS of mixed remote sensing pixels can be distinguished by their spatial locations in the triangular space. However, soil moisture and crop-residue moisture (SM-CRM) significantly reduce the performance of broadband remote sensing CRC indices and can thus decrease the accuracy of the remote estimation and mapping of CRC, FVC, and BS. This study evaluated the use of broadband remote sensing, the triangular space method, and the random forest (RF) technique to estimate and map the FVC, CRC, and BS of cropland in which SM-CRM changes dramatically. A spectral dataset was obtained using: (1) from a field-based experiment with a field spectrometer; and (2) from a laboratory-based simulation that included four distinct soil types, three types of crop residue (winter-wheat, maize, and rice), one crop (winter wheat), and varying SM-CRM. We trained an RF model [designated the broadband crop-residue index from random forest (CRRF)] that can magnify spectral features of crop residue and soil by using the broadband remote sensing angle indices as input, and uses a moisture-resistant hyperspectral index as the target. The effects of moisture on crop residue and soil were minimized by using the broadband CRRF. Then, the CRRF-NDVI triangular space method was used to estimate and map CRC, FVC, and BS. Our method was validated by using both laboratory- and field-based experiments and Sentinel-2 broadband remote-sensing images. Our results indicate that the CRRF-NDVI triangular space method can reduce the effect of moisture on the broadband remote-sensing of CRC, and may also help to obtain laboratory and field CRC, FVC, and BS. Thus, the proposed method has great potential for application to croplands in which the SM-CRM content changes dramatically.     

Mapping wetland characteristics using temporally dense Sentinel-1 and Sentinel-2 data: A case study in the St. Lucia wetlands,South Africa

Wetlands have been determined as one of the most valuable ecosystems on Earth and are currently being lost at alarming rates. Large-scale monitoring of wetlands is of high importance, but also challenging. The Sentinel-1 and -2 satellite missions for the first time provide radar and optical data at high spatial and temporal detail, and with this a unique opportunity for more accurate wetland mapping from space arises. Recent studies already used Sentinel-1 and -2 data to map specific wetland types or characteristics, but for comprehensive wetland characterisations the potential of the data has not been researched yet. The aim of our research was to study the use of the high-resolution and temporally dense Sentinel-1 and -2 data for wetland mapping in multiple levels of characterisation. The use of the data was assessed by applying Random Forests for multiple classification levels including general wetland delineation, wetland vegetation types and surface water dynamics. The results for the St. Lucia wetlands in South Africa showed that combining Sentinel-1 and -2 led to significantly higher classification accuracies than for using the systems separately. Accuracies were relatively poor for classifications in high-vegetated wetlands, as subcanopy flooding could not be detected with Sentinel-1’s C-band sensors operating in VV/VH mode. When excluding high-vegetated areas, overall accuracies were reached of 88.5% for general wetland delineation, 90.7% for mapping wetland vegetation types and 87.1% for mapping surface water dynamics. Sentinel-2 was particularly of value for general wetland delineation, while Sentinel-1 showed more value for mapping wetland vegetation types. Overlaid maps of all classification levels obtained overall accuracies of 69.1% and 76.4% for classifying ten and seven wetland classes respectively.     

Estimating aboveground biomass of the mangrove forests on northeast Hainan Island in China using an upscaling method from field plots,UAV-LiDAR data and Sentinel-2 imagery

The mangrove forests of northeast Hainan Island are the most species diverse forests in China and consist of the Dongzhai National Nature Reserve and the Qinglan Provincial Nature Reserve. The former reserve is the first Chinese national nature reserve for mangroves and the latter has the most abundant mangrove species in China. However, to date the aboveground ground biomass (AGB) of this mangrove region has not been quantified due to the high species diversity and the difficulty of extensive field sampling in mangrove habitat. Although three-dimensional point clouds can capture the forest vertical structure, their application to large areas is hindered by the logistics, costs and data volumes involved. To fill the gap and address this issue, this study proposed a novel upscaling method for mangrove AGB estimation using field plots, UAV-LiDAR strip data and Sentinel-2 imagery (named G∼LiDAR∼S2 model) based on a point-line-polygon framework. In this model, the partial-coverage UAV-LiDAR data were used as a linear bridge to link ground measurements to the wall-to-wall coverage Sentinel-2 data. The results showed that northeast Hainan Island has a total mangrove AGB of 312,806.29 Mg with a mean AGB of 119.26 Mg ha⁻¹. The results also indicated that at the regional scale, the proposed UAV-LiDAR linear bridge method (i.e., G∼LiDAR∼S2 model) performed better than the traditional approach, which directly relates field plots to Sentinel-2 data (named the G∼S2 model) (R² = 0.62 > 0.52, RMSE = 50.36 Mg ha⁻¹<56.63 Mg ha⁻¹). Through a trend extrapolation method, this study inferred that the G∼LiDAR∼S2 model could decrease the number of field samples required by approximately 37% in comparison with those required by the G∼S2 model in the study area. Regarding the UAV-LiDAR sampling intensity, compared with the original number of LiDAR plots, 20% of original linear bridges could produce an acceptable accuracy (R² = 0.62, RMSE = 51.03 Mg ha⁻¹). Consequently, this study presents the first investigation of AGB for the mangrove forests on northeast Hainan Island in China and verifies the feasibility of using this mangrove AGB upscaling method for diverse mangrove forests.     

In-season crop classification using elements of the Kennaugh matrix derived from polarimetric RADARSAT-2 SAR data

Accurate spatio-temporal classification of crops is of prime importance for in-season crop monitoring. Synthetic Aperture Radar (SAR) data provides diverse physical information about crop morphology. In the present work, we propose a day-wise and a time-series approach for crop classification using full-polarimetric SAR data. In this context, the 4 × 4 real Kennaugh matrix representation of a full-polarimetric SAR data is utilized, which can provide valuable information about various morphological and dielectric attributes of a scatterer. The elements of the Kennaugh matrix are used as the parameters for the classification of crop types using the random forest and the extreme gradient boosting classifiers.The time-series approach uses data patterns throughout the whole growth period, while the day-wise approach analyzes the PolSAR data from each acquisition into a single data stack for training and validation. The main advantage of this approach is the possibility of generating an intermediate crop map, whenever a SAR acquisition is available for any particular day. Besides, the day-wise approach has the least climatic influence as compared to the time series approach. However, as time-series data retains the crop growth signature in the entire growth cycle, the classification accuracy is usually higher than the day-wise data.Within the Joint Experiment for Crop Assessment and Monitoring (JECAM) initiative, in situ measurements collected over the Canadian and Indian test sites and C-band full-polarimetric RADARSAT-2 data are used for the training and validation of the classifiers. Besides, the sensitivity of the Kennaugh matrix elements to crop morphology is apparent in this study. The overall classification accuracies of 87.75% and 80.41% are achieved for the time-series data over the Indian and Canadian test sites, respectively. However, for the day-wise data, a ∼6% decrease in the overall accuracy is observed for both the classifiers.     

基于机器学习的海洋环境预报订正方法研究

结合中尺度数值模式 WRF 预报数据和 ERA5 再分析资料,利用机器学习方法对 WRF 预报场的风场、温度、气压进行预报订正。采用 ERA5 作为真值,与原始 WRF 预报相比,利用随机森林模型可以将预报结果整体均方根误差降低 44%以 上,利用深度神经网络模型可以将预报结果整体均方根误差降低 34%以上。通过随机森林模型实验得到不同输入特征对预报要素的影响程度,分析了关键的预报订正因子。     

Changes in air pollutants during the COVID-19 lockdown in Beijing:Insights from a machine-learning technique and implications for future control policy

基于2015-2020年北京35个环境空气站和20个气象站观测资料,应用机器学习方法(随机森林算法)分离了气象条件和源排放对大气污染物浓度的影响.结果发现,为应对疫情采取的隔离措施使北京2020年春节期间大气污染物浓度降低了35.1％-51.8％;其中,背景站氮氧化物和一氧化碳浓度的降幅最大,超过了以往报道较多的交通站...     

中国碳强度关键影响因子的机器学习识别及其演进

碳强度影响因子数量众多,通过在众多因子中评估其重要性以识别出关键影响因子进而解析碳强度关键因子的变化规律,是中国2030年碳强度能否实现比2005年下降60%~65%目标的科学基础。传统的回归分析方法对于评估众多因子的重要性存在多重共线性等问题,而机器学习处理海量数据则具有较好的稳健性等优点。本文从能源结构、产业结构、技术进步和居民消费等方面选取了56个中国碳强度影响因子指标,采用随机森林算法基于信息熵评估了1980-2014年逐年各项因子的重要性,通过指标数量与信息熵的对应关系统一筛选出每年重要性最大的前22个指标作为相应年度关键影响因子,最终依据关键影响因子的变化趋势划分了3个阶段作了演进分析。结果发现：1980-1991年,碳强度的关键因子主要以高耗能产业规模及占比、化石能源占比和技术进步为主;1992-2007年,中国经济进入快车道增长时期,服务业占比和化石能源价格对碳强度的影响作用开始显现,居民传统消费的影响作用在增大;2008年全球金融危机后,中国进入经济结构深化调整时期,节能减排力度大大增强,新能源占比和居民新兴消费的影响作用迅速显现。为实现2030年碳强度下降60%~65%目标,优化能源结构和产业结构,促进技术进步,提倡绿色消费,强化政策调控是未来需要采取的主要措施。     

石羊河下游人工梭梭林土壤呼吸变化特征及其与水热因子的关系

土壤呼吸不仅是反映土壤生物活性的重要指标,也是全球碳循环研究中备受关注的热点问题。在地处典型干旱区的石羊河下游,以流动沙丘和去除土壤结皮人工梭梭林为对照,采用LI-8100土壤碳通量监测系统研究了栽植约40 a、30 a、10 a和5 a的人工梭梭林生长季和非生长季的土壤呼吸日变化,并分析了土壤水分和温度对土壤呼吸的影响。结果表明：（1）不同林龄梭梭林生长季和非生长季土壤呼吸速率的日变化均为明显的单峰曲线,且呈现出一定的波动性,日最大排放速率出现在12:00～14:00时,最小值出现在8:00时左右。（2）梭梭林营造和去结皮处理显著提高了沙漠土壤呼吸速率,而且不同林龄土壤呼吸速率大体上随着种植年限的增加而递增,表现为MC >40 a>30 a>10 a>MS >5 a,非生长季表现为MC >40 a>10 a>5 a>30 a>MS。（3）不同林龄梭梭林土壤呼吸速率均具有明显的季节变化特征,生长季（8 月）的土壤呼吸作用明显强于非生长季（1月）。（4）相关性分析表明,生长季和非生长季土壤呼吸均与0～5 cm土壤水分显著相关,且均呈二次曲线关系,分别为Y =-0.205 8X ²+0.946 5X-0.316 6（R ² =0.506 ²,P= 0.041 7）和Y= 0.118 7 X ²+0.156 3X+0.118 8（R ²=0.675 7,P =0.001 1）;但与10 cm土壤温度的相关性不显著,土壤水分是影响人工梭梭林土壤呼吸的关键因素。该研究进一步证明了人工梭梭林的营造有效改善了沙漠土壤的生物活性,提高了土壤碳通量水平,以土壤结皮破坏为基本特征的人工梭梭林退化和沙漠化必然在短期内加剧碳排放。因此,需要在沙漠地区合理营造人工林,并在造林和林业管理过程中注意保护土壤结皮,以减少CO²排放。     

天山雪岭云杉个体生物量分配及其变化规律的研究

为了探讨天山雪岭云杉林生物量在个体组织中的分配情况及其变化规律,在研究区进行了大量的野外测量,利用已有的雪岭云杉林估算方程,分析了天山雪岭云杉林生物量在各器官（干、枝、叶、皮、根）中的分配及其变化规律。结果表明：（1） 研究区雪岭云杉林的平均生物量为388.74 t·hm^-2,树木各器官中,干、枝、根、叶和皮分别占生物量的43.65%、28.60%、13.49%、11.08%和3.18%。（2）各径级生物量所占百分比为：33.53%（40~50 cm）、20.13%（20~30 cm）、19.59%（30~40 cm）、18.19%（50~60 cm）和2.05%（10~20 cm）;树木生物量在不同树高中的分配表现为：48.78%（20~30 m）＞35.27%（10~20 m）＞14.70%（30~40 m）＞1.25%（0~10 m）;地上和地下生物量的分配比例为：87.54%和12.46%,分别为340.30 t·hm^-2和48.44 t·hm^-2。（3） 随海拔升高,天山雪岭云杉林生物量呈“单峰”变化,在海拔2 100~2 400 m处达到最大值611.58 t·hm^-2;干、皮生物量所占比例随海拔升高而减小,枝生物量逐渐增加,叶、根生物量呈先减小后增加的趋势;径级20~30 cm、30~40 cm和50~60 cm的生物量随海拔升高均呈“单峰型”变化趋势,都在海拔2 100~2 400 m处达到最大;雪岭云杉林不同树高生物量随海拔的升高呈现的趋势不同。天山雪岭云杉林生物量和年均降水量随经纬度的升高均呈降低变化,研究区林分生物量自西向东总体呈现逐渐降低的趋势;林分密度、海拔和降水共同决定了森林生物量的大小及其变化规律,海拔2 100~2 400 m是本研究区雪岭云杉林生长的最适宜场所。结果可为雪岭云杉林生态系统的恢复和重建提供基础资料,对研究区进行综合管理与生态健康分析具有重要意义。