站点CERES-Rice模型区域应用效果和误差来源

作物区域模拟是利用有限的空间数据,最大限度地反映出生育期、产量等作物性状的时空变化规律.由于目前的作物模型大多是田间尺度的站点模型,把它运用到区域水平的效果如何研究甚少.文章利用CERES-Rice模型,对作物模型在我国的区域应用效果进行了分析.首先利用田间观测数据在各实验点上对模型进行了详细的站点校准,以验证模型在我国的模拟能力;然后以我国水稻生态区(精确到亚区)为单位,运用均方根差(RMSE)法进行了区域校准和验证;最后利用区域校准后的CERES-Rice模型,模拟了1980～2000年的网格(50km×50km)水稻产量,并与同期农调队调查产量进行统计比较,以验证区域应用的效果,为区域模拟的推广和应用提供参考.结果表明:经过空间校准后的CERES-Rice模型,在水稻的主产区1～4区(占种植面积的95%)模拟的平均产量与调查产量相对均方根差在22%以内,两者的符合度也较好,个别区域(5、6) RMSE%在24%～30%之间;1980～2000年水稻各产区模拟的平均产量与调查产量随时间变化趋势也具有一定的一致性;全国1896个网格中,大部分网格(71.01%)模拟的21年水稻年产量与调查产量的RMSE%在30%之内,且大部分分布在水稻主产区,考虑到水稻种植面积的权重后,认为利用区域校准和验证后的CERES-Rice模型进行水稻区域模拟,可以反映出产量的时空分布特征,能够为宏观决策提供相应的信息.但目前区域模拟中还存在着一定的误差,有待今后进一步研究.     

不同气象数据空间内插对区域玉米生长模拟结果的影响

作物模型区域模拟已成为作物模型应用的一个新方向。运用作物模型进行区域研究时,遇到的问题之一就是输入模型的空间数据质量问题,研究不同空间内插法获得的气象数据对作物模型区域模拟结果的影响,可以为区域模拟对输入数据的敏感性研究提供一定的参考。利用区域校准的CERES-Maize模型,将3类内插方法(几何内插、统计内插、动力模型内插)产生的网格化天气数据分别输入到CERES-Maize模型中,模拟了50km×50km网格水平下1961—1990年我国玉米生产状况,并选取1980—1990年模拟的平均产量与同期农调队调查产量进行比较,以了解区域模拟中,不同空间内插方法所得的逐日气象数据对区域模拟结果的影响。结果表明:(1)作物模型区域应用时,所采用的3种内插方法都能满足作物模型区域模拟对网格化天气数据的要求,采用3种天气数据的区域模拟结果都能反映出玉米平均产量的空间变化特征,与网格调查平均产量之间具有极显著的相关关系,但采用不同内插天气数据对模拟结果造成了8%以内的偏差。(2)采用不同内插天气数据,在进行作物区域模拟时,各方法的模拟结果之间呈极显著的相关关系,但这些模拟结果之间,在全国大部分地区是差异显著。     

区域气候模式与作物模型联接的影响评估模拟实验及不确定性分析

利用英国Hadley中心开发的区域气候模式RCMPRECIS(网格分辨率50km×50km),与经过田间试验资料和历史气候资料验证和校准过的CERES系列作物模式相结合,就区域气候模式与作物模式联接的影响评估方法及其不确定性进行了评估。结果表明,相对于大气环流模型来说,区域气候模式与作物模型的结合省去了随机天气发生器的中间环节,减小了不确定性产生的因素。在站点模拟上,该方法在平原地区的模拟效果较好,而山区的模拟效果较差,但如果能用实测天气数据对模拟的天气数据进行验证,模拟效果明显提高。在区域模拟上,该方法可以较好地体现出产量变化的空间分布规律,但由于空间数据的限制,模拟产量与实际产量的偏差较站点水平要大。     

水稻模型ORYZA2000在湖南双季稻区的验证与适应性评价

校准与验证水稻生长模型ORYZA2000,为模型本地化、区域化研究应用提供依据。本文采用湖南双季稻区作物田间观测数据,结合栽培管理措施、土壤以及同期逐日气象数据等资料对ORYZA2000进行参数校正,调试确定了早稻、晚稻有代表性品种的作物参数;利用独立的数据资料,对双季稻生育期、叶面积指数、生物量、产量等指标的模拟结果进行了详细地验证与适应性评价。结果表明：模型对双季稻品种的生育期模拟较好,开花期和成熟期的相对模拟误差为1—2d;早稻和晚稻叶面积指数的归一化均方根误差（NRMSE）均为24%,地上部总生物量、绿叶生物量、茎生物量和穗生物量的NRMSE值分别为18%、22%、22%、24%和19%、24%、28%、28%,产量的NRMSE值分别为11%和16%。校验的作物参数反映了湖南早稻和晚稻品种的生物学特性,参数值合理、有效。通过校准作物参数,ORYZA2000可较为准确地模拟双季稻生长发育及其生物量的动态累积过程,适应性较强,能够应用于双季稻生产。     

基于更新和同化策略相结合的遥感信息与水稻生长模型耦合技术的研究

将遥感与作物模型耦合有利于提高作物模型在区域尺度应用时的精度。基于集合平方根滤波算法(Ensemble Square RootFilter,EnSRF)和粒子群优化算法(Particle Swarm Optimization,PSO),以叶面积指数(Leaf Area Index,LAI)和叶片氮积累量(Leaf Nitrogen Accumulation,LNA)共同作为同化耦合点和过程更新点,将同化与更新策略相结合,研究建立了基于遥感信息与水稻生长模型(RiceGrow)耦合的水稻生长与产量预测技术。结果表明,将更新和同化策略结合后,利用RiceGrow模型模拟的水稻生长指标和产量结果更接近于实测值。其中LAI、LNA和产量与实测值间的RMSE分别为0.94、0.47 g/m2和320.15 kg/hm2;RiceGrow模型直接模拟LAI、LNA和产量的RMSE为1.25、1.24 g/m2和516.83 kg/hm2;而单纯基于同化策略模拟LAI、LNA和产量的RMSE为1.01、0.59 g/m2和335.70 kg/hm2。此外,基于该技术的模型区域尺度预测结果能较好地描述水稻生长和产量的时空分布状况,生长指标及区域总产量的模拟相对误差均小于20%。显示基于更新和同化策略相结合的遥感与模型耦合技术具有较高的预测精度,从而为区域尺度作物生长和产量预测提供了技术支撑。     

利用遥感技术实现作物模拟模型区域应用的研究进展

作物模拟模型从单点发展到区域应用时,模型中一些宏观资料的获取和参数的区域化方面出现困难,利用遥感技术将实现作物模拟模型的区域应用．文中综述了近年来遥感反演作物模型所需的地表生物物理参数的方法、利用遥感信息直接获取生物量的途径和遥感信息与作物模拟模型之间时空匹配问题等方面的研究概况,重点介绍了利用遥感技术实现作物模拟模型区域应用的3种解决方案(强迫型、调控型和验证型)及其研究进展,并讨论了目前存在的问题和今后研究的方向．     

大麦叶面积指数模拟模型

准确模拟叶面积指数是作物生长模拟模型预测作物生长和产量的关键.本文通过系统分析扬州和武汉地区不同大麦品种高产群体叶面积指数变化动态,建立了大麦群体的叶面积指数模拟模型.大麦叶面积指数是品种叶面积指数扩展的遗传参数和气温日较差、日照时数、辐射量等气候因子及水肥丰缺因子的函数.孕穗抽穗期最大叶面积指数与该期最适叶面积指数是不同的概念,二者之间存在着极显著差异.利用扬州、南京和昆明地区不同品种的播期试验及氮肥试验资料对模型进行了检验,结果表明,模型对大麦叶面积指数的模拟效果较好,模拟值与观测值吻合度高,根均方差RMSE介于0.742～2.865,平均值为1.348.对模拟值与观测值进行y=x的线性回归分析,相关系数R²介于0.511～0.954,均呈极显著正相关.     

灌溉水稻生长发育和潜力产量的模拟模型

本文提出的HDRICE模型是灌溉水稻生长的生理生态模型,它由相互衍接的水稻形态发育、干物质积累和叶面积发育三模块组成。形态发育模块用以模拟逐日温度和日长对水稻发育的影响,其参数可反映水稻品种的基本营养性、感温性和感光性;干物质积累模块用以模拟冠层CO_2同化、作物的维持呼吸和生长呼吸及干物质分配等过程;叶面积发育模块用以模拟叶面积指数的动态。本文还讨论了模型的输入参数和模型检验。模型可应用于模拟水稻的生长发育,预测水稻品种潜在产量及为取得潜在产量所必需的群体数量指标。     

WOFOST模型在内蒙古河套灌区模拟玉米生长全程的适应性

在河套灌区引入成熟的作物模型并进行适应性验证,可为进一步开展玉米生长监测及估产提供依据和基础。本文利用河套灌区巴彦淖尔农业气象试验站2012年玉米观测数据,结合当地气象、土壤资料对荷兰瓦赫宁根大学开发的WOFOST模型进行参数校准,并利用2013年玉米观测数据和2001—2011年农业气象观测资料对模型的区域适用性进行验证,获得了玉米的基本作物参数,包括各发育阶段比叶面积、最大CO2同化率、单叶光能利用率等。结果表明:通过校准作物参数,WOFOST模型可以较好地模拟LAI扩展、生物量的动态积累过程,LAI、各器官生物量及最终产量的模拟值与实测值吻合较好;独立样本检验中,模型模拟LAI的绝对偏差平均值为0.75,叶生物量、茎生物量、贮存器官生物量、地上部总生物量、产量的归一化均方根误差分别为33%、26%、17%、18%和13%;模拟2001—2011年玉米产量的归一化均方根误差为7.5%。参数校准后的模型对LAI、各器官生物量、产量的模拟结果较为符合实际,WOFOST模型能够适用于河套地区玉米生产过程生理、生态因子诊断、评估等。     

作物生长模拟模型参数校正与有效化的理论和实践

以ＧＯＳＳＹＭ 模型为例系统阐述了作物生长模拟模型有关参数校正和模型有效化的一般原理和方法,同时用新疆棉区的试验具体校正了品种参数、土壤参数和修改了部分模块,并对校正结果进行了验证．结果表明,两个试点土壤２０ ～４０ 、６０ ～８０ｃｍ 两个土层生长季水分动态观测值与土壤参数校正后模型的模拟值吻合较好;系５ 品种试验生育期６ 项生物指标动态模拟结果与实测值拟合的相关系数都在０ ．９ 以上,并且不同栽培条件的３ 个处理的模拟产量的相对误差平均为７ ．５ ％ ,模拟结果较理想．     

Equivalent Porous Media (EPM) Simulation of Groundwater Hydraulics and Contaminant Transport in Karst Aquifers

Karst aquifers have a high degree of heterogeneity and anisotropy in their geologic and hydrogeologic properties which makes predicting their behavior difficult. This paper evaluates the application of the Equivalent Porous Media (EPM) approach to simulate groundwater hydraulics and contaminant transport in karst aquifers using an example from the North Coast limestone aquifer system in Puerto Rico. The goal is to evaluate if the EPM approach, which approximates the karst features with a conceptualized, equivalent continuous medium, is feasible for an actual project, based on available data and the study scale and purpose. Existing National Oceanic Atmospheric Administration (NOAA) data and previous hydrogeological U. S. Geological Survey (USGS) studies were used to define the model input parameters. Hydraulic conductivity and specific yield were estimated using measured groundwater heads over the study area and further calibrated against continuous water level data of three USGS observation wells. The water-table fluctuation results indicate that the model can practically reflect the steady-state groundwater hydraulics (normalized RMSE of 12.4%) and long-term variability (normalized RMSE of 3.0%) at regional and intermediate scales and can be applied to predict future water table behavior under different hydrogeological conditions. The application of the EPM approach to simulate transport is limited because it does not directly consider possible irregular conduit flow pathways. However, the results from the present study suggest that the EPM approach is capable to reproduce the spreading of a TCE plume at intermediate scales with sufficient accuracy (normalized RMSE of 8.45%) for groundwater resources management and the planning of contamination mitigation strategies.     

An efficient approach to automate the manual trial and error calibration of activated sludge models

An efficient approach is introduced to help automate the rather tedious manual trial and error way of model calibration currently used in activated sludge modeling practice. To this end, we have evaluated a Monte Carlo based calibration approach consisting of four steps: (i) parameter subset selection, (ii) defining parameter space, (iii) parameter sampling for Monte Carlo simulations and (iv) selecting the best Monte Carlo simulation thereby providing the calibrated parameter values. The approach was evaluated on a formerly calibrated full-scale ASM2d model for a domestic plant (located in The Netherlands), using in total 3 months of dynamic oxygen, ammonia and nitrate sensor data. The Monte Carlo calibrated model was validated successfully using ammonia, oxygen and nitrate data collected at high measurement frequency. Statistical analysis of the residuals using mean absolute error (MAE), root mean square error (RMSE) and Janus coefficient showed that the calibrated model was able to provide statistically accurate and valid predictions for ammonium, oxygen and nitrate. This shows that this pragmatic approach can perform the task of model calibration and therefore be used in practice to save the valuable time of modelers spent on this step of activated sludge modeling. The high computational demand is a downside of this approach but this can be overcome by using distributed computing. Overall we expect that the use of such systems analysis tools in the application of activated sludge models will improve the quality of model predictions and their use in decision making.     

ORCHIDEE‐STICS,a process‐based model of sugarcane biomass production: calibration of model parameters governing phenology

Agro‐Land Surface Models (agro‐LSM) combine detailed crop models and large‐scale vegetation models (DGVMs) to model the spatial and temporal distribution of energy, water, and carbon fluxes within the soil–vegetation–atmosphere continuum worldwide. In this study, we identify and optimize parameters controlling leaf area index (LAI) in the agro‐LSM ORCHIDEE‐STICS developed for sugarcane. Using the Morris method to identify the key parameters impacting LAI, at eight different sugarcane field trial sites, in Australia and La Reunion island, we determined that the three most important parameters for simulating LAI are (i) the maximum predefined rate of LAI increase during the early crop development phase, a parameter that defines a plant density threshold below which individual plants do not compete for growing their LAI, and a parameter defining a threshold for nitrogen stress on LAI. A multisite calibration of these three parameters is performed using three different scoring functions. The impact of the choice of a particular scoring function on the optimized parameter values is investigated by testing scoring functions defined from the model‐data RMSE, the figure of merit and a Bayesian quadratic model‐data misfit function. The robustness of the calibration is evaluated for each of the three scoring functions with a systematic cross‐validation method to find the most satisfactory one. Our results show that the figure of merit scoring function is the most robust metric for establishing the best parameter values controlling the LAI. The multisite average figure of merit scoring function is improved from 67% of agreement to 79%. The residual error in LAI simulation after the calibration is discussed.     

单木生物量模型估计区域尺度生物量的不确定性

采用系统抽样体系江西省固定样地杉木连续观测数据和生物量数据,通过Monte Carlo法反复模拟由单木生物量模型推算区域尺度地上生物量的过程,估计了江西省杉木地上总生物量。基于不同水平建模样本量n及不同决定系数R~2的设计,分别研究了单木生物量模型参数变异性及模型残差变异性对区域尺度生物量估计不确定性的影响。研究结果表明:2009年江西省杉木地上生物量估计值为(19.84±1.27)t/hm~2,不确定性占生物量估计值约6.41%。生物量估计值和不确定性值达到平稳状态所需的运算时间随建模样本量及决定系数R~2的增大而减小;相对于模型参数变异性,残差变异性对不确定性的影响更小。     

From GCM grid cell to agricultural plot: scale issues affecting modelling of climate impact

General circulation models (GCM) are increasingly capable of making relevant predictions of seasonal and long-term climate variability, thus improving prospects of predicting impact on crop yields. This is particularly important for semi-arid West Africa where climate variability and drought threaten food security. Translating GCM outputs into attainable crop yields is difficult because GCM grid boxes are of larger scale than the processes governing yield, involving partitioning of rain among runoff, evaporation, transpiration, drainage and storage at plot scale. This study analyses the bias introduced to crop simulation when climatic data is aggregated spatially or in time, resulting in loss of relevant variation. A detailed case study was conducted using historical weather data for Senegal, applied to the crop model SARRA-H (version for millet). The study was then extended to a 10 degrees N-17 degrees N climatic gradient and a 31 year climate sequence to evaluate yield sensitivity to the variability of solar radiation and rainfall. Finally, a down-scaling model called LGO (Lebel-Guillot-Onibon), generating local rain patterns from grid cell means, was used to restore the variability lost by aggregation. Results indicate that forcing the crop model with spatially aggregated rainfall causes yield overestimations of 10-50% in dry latitudes, but nearly none in humid zones, due to a biased fraction of rainfall available for crop transpiration. Aggregation of solar radiation data caused significant bias in wetter zones where radiation was limiting yield. Where climatic gradients are steep, these two situations can occur within the same GCM grid cell. Disaggregation of grid cell means into a pattern of virtual synoptic stations having high-resolution rainfall distribution removed much of the bias caused by aggregation and gave realistic simulations of yield. It is concluded that coupling of GCM outputs with plot level crop models can cause large systematic errors due to scale incompatibility. These errors can be avoided by transforming GCM outputs, especially rainfall, to simulate the variability found at plot level.     

Modeling hydrologic controls on denitrification: sensitivity to parameter uncertainty and landscape representation

This paper presents a general discussion of the interplay between model structure and hydrologic parameters in the context of denitrification estimation using coupled hydro-ecosystem models at a watershed scale. Given the key role played by hydrology in denitrification models, sensitivity analysis of hydrologic parameters is needed to determine both uncertainty in denitrification estimates and to suggest how measured data, such as streamflow, can be effectively used to reduce this uncertainty. This paper contributes to the broad goal of sensitivity analysis by examining the linkage between landscape tessellation, calibration, and the ability of models to capture hot-spot contributions to watershed scale denitrification across a range of N-loading. For a small mid-Atlantic forested watershed, denitrification estimates using RHESSys (regional hydro-ecologic simulation system) are compared across different strategies for calibration and landscape tessellation. Results demonstrate the utility of several potential approaches to account for hydrologically mediated hot-spots within landscapes.     

油菜绿色面积指数动态模拟模型

准确模拟绿色面积指数是作物生长模拟模型可靠预测作物生长和产量的关键。该研究的目的是以生理生态过程为基础,构建油菜(Brassica napus)叶面积指数和角果面积指数变化动态的模拟模型。油菜叶面积指数模型综合考虑了库或源限制下的叶面积增长模式,其中库限制下叶面积指数的增长呈指数方程,且受到温度、水分和氮素因子的影响;源限制下叶面积指数增长用比叶面积法来模拟。油菜角果面积指数由比角果面积和角果干物重来决定。比叶面积和比角果面积均为生理发育时间的函数。利用不同类型品种的播期试验及氮肥试验资料分别对模型进行了校正和检验,结果表明模型能较好地模拟不同条件下油菜叶面积指数和角果面积指数。