CoLM模式地表温度变分同化研究

本文采用变分方法对通用陆面模式 (CoLM) 中的地表温度进行同化.同化伴随约束条件采用CoLM模式中的地表及植被能量平衡方程,调节因子采用裸土及植被蒸发比.采用美国通量网 (AmeriFlux) 中的Bonville站数据对同化方法进行了单点验证,验证结果表明同化后地表温度以及蒸散结果更加接近于实测值.选取中国华北地区对同化方法进行区域验证,结果显示每天仅采用白天一次观测值对地表温度进行同化的方法是有效的.通过对同化前后地表温度误差直方图比较可以发现,在有MODIS观测值的区域,同化后白天地表温度误差大大降低,同时,同化后地表蒸散空间分布图也发生了变化.单点验证以及区域验证结果都表明了变分同化方法是可靠的.变分同化方法可以改进陆面模式模拟结果,对于地表过程研究中的植被生态、水文等研究具有重要意义,同时,陆面模式可以与数值预报模式进行耦合,改进数值预报结果.     

Coupling the Common Land Model to ECHAM5 Atmospheric General Circulation Model

The ECHAM5 model is coupled with the widely used Common Land Model(CoLM). ECHAM5 is a state-of-theart atmospheric general circulation model incorporated into the integrated weather and climate model of the Chinese Academy of Meteorological Sciences(CAMS-CSM). Land surface schemes in ECHAM5 are simple and do not provide an adequate representation of the vegetation canopy and snow/frozen soil processes. Two AMIP(Atmospheric Model Intercomparison Project)-type experiments using ECHAM5 and ECHAM5-CoLM are run over 30 yr and the results are compared with reanalysis and observational data. It is found that the pattern of land surface temperature simulated by ECHAM5-CoLM is significantly improved relative to ECHAM5. Specifically, the cold bias over Eurasia is removed and the root-mean-square error is reduced in most regions. The seasonal variation in the zonal mean land surface temperature and the in situ soil temperature at 20-and 80-cm depths are both better simulated by ECHAM5-CoLM. ECHAM5-CoLM produces a more reasonable spatial pattern in the soil moisture content, whereas ECHAM5 predicts much drier soils. The seasonal cycle of soil moisture content from ECHAM5-CoLM is a better match to the observational data in six specific regions. ECHAM5-CoLM reproduces the observed spatial patterns of both sensible and latent heat fluxes. The strong positive bias in precipitation over land is reduced in ECHAM5-CoLM, especially over the southern Tibetan Plateau and middle–lower reaches of the Yangtze River during the summer monsoon rainy season.     

基于陆面模式基准平台ILAMB对陆面模式CoLM的评估

模式评估是模式发展中不可或缺的重要一环。本文利用最新版陆面模式评估软件—ILAMB（International Land Model Benchmarking）对通用陆面模式 (The Common Land Model, CoLM)进行客观评估,并与NCAR陆面模式CLM5(Community Land Model version 5)结果进行比较。作为一个陆面模式评估软件,ILAMB能对参与评估的模式变量自动生成诊断图形并对模式性能进行评分。评估结果表明,CoLM总体性能良好,模拟结果与基准数据较为接近。与CLM5相比,CoLM在总初级生产力及水文方面的表现略微逊色,在辐射方面则表现相当,对部分变量如地表向上长波辐射、地表净辐射等的模拟甚至优于CLM5。通过对比CRUNCEPv7和GSWP3v1两种强迫资料发现,它们在气候平均态上具有一定的差异,模式在GSWP3v1强迫下的表现相对较好。CoLM和CLM5在CRUCNEPv7强迫下模拟的潜热通量在亚马逊平原、亚洲东部和南部地区以及北美东部一带正偏差显著,而在GSWP3v1强迫下的模拟则有明显改善。这两个模式对感热通量的模拟在非洲北部、亚洲中部一带均明显偏高。在辐射方面,CoLM模拟的地表向上短波辐射在全球以偏高为主,这在一定程度上造成了地表净辐射的模拟偏低。各组试验模拟的地表向上长波辐射得分相差不大,但在具体的空间分布上有一定的差别。     

陆面模式中土壤和植被经验参数随机误差的传播研究

数据质量问题和模式参数化方案的非完备性是陆面模拟中不确定性的主要来源。本文将高斯误差传播原理 (Gaussian Error Propagation, GEP) 应用于通用陆面模式 (the Common Land Model, CoLM), 研究关键的植被和土壤属性参数随机误差在模式中的传播, 确定由此类误差导致的CoLM模拟的不确定性。结果表明: (1) 基于本研究给定的土壤和植被参数的不确定性, CoLM模拟的表层土壤温度、 土壤湿度和植被蒸散通量 (植冠蒸腾+地表蒸发) 的相对误差分别为0.11%、 34.07%和5.58%; 砂土和稀疏森林上模拟效果最差。土壤参数随机误差对CoLM模拟的影响高于植被参数, 而土壤水文参数 (孔隙率、 饱和基质势、 气孔尺寸分布指数和饱和导水率) 对各模拟量不确定性的贡献率均远大于热力参数 (饱和反照率和热容)。对于本研究涉及的所有模拟变量而言, 最关键的参数均是气孔尺寸分布指数b, 这可能与描述基质势与体积水含量关系的函数有关, 其次重要的是砂土的孔隙度和粘土的饱和导水率。混交森林上的根深分布和苔原上的动力学粗糙度对蒸散通量贡献显著。本身相对误差大的经验参数对CoLM模拟不确定性的贡献不一定多。 (2) 干燥条件下 (表层液态水饱和度小于0.1) 土壤温度的不确定性大; 相变发生时刻附近 (表层土壤温度在0℃附近且表层液态水含量大于0) 土壤湿度不确定性显著; 蒸散通量的不确定性随本身绝对值的增大而增大, 在相对温暖干燥环境中 (表层土壤温度高于280 K且表层液态水饱和度小于0.3) 其不确定性最高。研究证实, GEP能够辨识CoLM中需优先提高观测精度的关键参数和关键参数化过程, 对陆面模拟的参数选定、 不确定性评估和模式完善具有重要意义。     

Testing and improving the performance of the Common Land Model: A case study for the Gobi landscape

Land surface processes take place on the interface between the earth and atmosphere, exerting significant influences on the weather and climate. Correct modeling of these processes is important to numerical weather forecast and climate prediction. In order to obtain a more thorough understanding of the land surface processes over the Gobi landscape, we evaluated the performance of the Common Land Model (CoLM) at Dunhuang station in Gansu Province of China to determine whether the model formulation, driven by observational data, is capable of simulating surface fluxes over the underlying desert surface. In comparison with the enhanced observation data collected at Dunhuang station over the period 22–28 August 2008, the results showed that the surface albedo simulated by CoLM was larger than that in the observation, and the simulated surface temperature was lower than the observed. After the measured values were used to correct the surface albedo, the solar radiation absorbed by the ground surface was more consistent with the measurements. A new empirical relationship of the surface thermal exchange coefficient r_ah was used to modify the thermal aerodynamic impedance. The simulated soil surface temperature became significantly closer to the observed value, and the simulated surface sensible heat as well as net radiative fluxes were also improved.     

Uncertainty of Climate Response to Natural and Anthropogenic Forcings Due to Different Land Use Scenarios

The A.M.Obukhov Institute of Atmospheric Physics,Russian Academy of Sciences (IAP RAS) climate model (CM) of intermediate complexity is extended by a spatially explicit terrestrial carbon cycle module.Numerical experiments with the IAP RAS CM are performed forced by the reconstructions of anthropogenic and natural forcings for the 16th to the 20th centuries and by combined SRES (Special Report on Emission Scenarios) A2-LUH (Land Use Harmonization) anthropogenic scenarios for the 21st century.Hereby,the impact of uncertainty in land-use scenarios on results of simulations with a coupled climate-carbon cycle model is tested.The simulations of the model realistically reproduced historical changes in carbon cycle characteristics.In the IAP RAS CM,climate warming reproduced in the 20th and 21st centuries enhanced terrestrial net primary production but terrestrial carbon uptake was suppressed due to an overcompensating increase in soil respiration.Around year 2100,the simulations the model forced by different land use scenarios diverged markedly,by about 70 Pg (C) in terms of biomass and soil carbon stock but they differed only by about 10 ppmv in terms of atmospheric carbon dioxide content.     

北亚热带落叶阔叶林土壤呼吸时间变化特征

利用LI-8100开路式土壤碳通量系统测定龙王山森林土壤呼吸速率,研究北亚热带落叶阔叶林土壤呼吸速率的日变化和季节性变化规律.结果表明:北亚热带落叶阔叶林土壤呼吸速率在12—14时达到最大,与土壤温度变化基本一致;森林土壤呼吸速率随土壤温度的季节性变化而变化,在夏季土壤呼吸速率较高,在冬季土壤呼吸速率较低;土壤呼吸速率与土壤温度间存在着明显的指数关系,土壤呼吸温度敏感系数Q₁₀为2.81.     

城市陆面模式设计及检验

针对城市冠层模式(Urban Canopy Model,UCM)的各种局限性,本文首次提出整体城市陆面模式(Bulk Urbanized Land Surface Model,BULSM)的概念,即在现有先进的陆面模式的基础上直接构建城市陆面模式。本文在通用陆面模式(Common Land Model,CoLM)基础上构建整体城市陆面模式,根据城市的特点对CoLM进行了发展和重新参数化。模式采用高分辨率地表分类数据;改进了模式的地表能量平衡方程与水分平衡方程;对反照率、城市地表粗糙度、城市大气廓线、城市人为热、不透水面水分蒸发及积水深度等进行了重新参数化。模式保留了CoLM在自然下垫面的全部特性,加强了CoLM在城市或人为下垫面的模拟和预报能力。区域验证以及与CoLM对比结果表明整体城市陆面模式模拟结果能够显示出北京城区各参数的空间精细化分布情况,能够反映出地表分类对模拟结果的影响。单站点验证以及与CoLM和高分辨率城市陆面同化系统(u-HRLDAS)对比结果表明整体城市陆面模式地表温度模拟结果与CoLM相比有很大提高;和u-HRLDAS相比大部分站点地表温度模拟结果也有一定的改进。     

The Spatial Variability of Energy and Carbon Dioxide Fluxes at the Floor of a Deciduous Forest

Fluxes of carbon dioxide, water and sensible heat were measured using three different eddy covariance systems above the forest floor of a closed deciduous forest (leaf area index approx 6). The primary objective was to examine the representativeness of a single eddy covariance system in estimating soil respiration for time scales ranging from one-half hour to more than one week. Experiments were conducted in which the eddy covariance sensors were in one of three configurations: i) collocated, ii) separated horizontally or iii) separated vertically. A measure of the variation between the three systems (CV',related to the coefficient of variation) for half-hour carbon dioxide fluxes was 0.14 (collocated systems), 0.34 (vertically separated systems at 1, 2 and 4 m above the surface), and 0.57 (systems horizontally separated by 30 m). A similar variation was found for other scalar fluxes (sensible and latent heat). Variability between systems decreased as the number of half-hour sampling periods used to obtain mean fluxes was increased. After forty-eight hours (means from ninety-six half-hour samples), CV' values for carbon dioxide fluxes were 0.07, 0.09 and 0.16 in the collocated, vertically separated and horizontally separated experiments, respectively. The time dependence of variability has implications on the appropriateness of using short-term measurements in modelling validation studies. There are also implications concerning the appropriate number of half-hour samples necessary to obtain reliable causal relationships between flux data and environmental parameters. Based on the longer-term measurements, we also discuss the representativeness of a single eddy covariance system in long-term monitoring of soil respiration and evaporation beneath forest canopies using the eddy covariance method.     

一个区域地气耦合柱状模式的发展及其对半干旱区的模拟

在现有的区域气候模式的基础上构建了用于研究地气耦合模拟的区域柱状模式,其中包括了详细的物理过程,可以作为研究半干旱区地气相互作用过程的工具。通过区域柱状模式与不同陆面模式的耦合,检验了不同陆面模式在非耦合和耦合两种情况下对半干旱区的模拟性能。结果表明,陆面模式在耦合后对感热的模拟偏差放大,而潜热的模拟效果得到了改进,不同的陆面模式对区域柱状模式的模拟性能影响较大,耦合了CoLM后的区域柱状模式对地表气象要素的模拟得到了较大的改进。     

基于FLUXNET数据集对陆面模式CoLM能量通量的单点评估

模式评估是模式发展中的重要一环。本文利用来自FLUXNET2015数据集的30个站点的涡动相关系统观测数据,重点关注能量通量,对通用陆面模式（Common Land Model version 2014,CoLM2014）在不同典型下垫面的模拟能力进行评估。结果表明,模式总体上能抓住感热、潜热和净辐射通量在日、季节和年平均等不同时间尺度上的变化特征,对感热、潜热和净辐射通量都有较好的模拟能力,净辐射的模拟效果最好,潜热通量次之。季节变化模拟中,感热、潜热通量在夏季不同植被型下站点的空间离散程度大于冬季,不同站点间模拟效果相差较大,净辐射多站点标准差变化幅度要小于感热、潜热,不同站点间模拟效果偏差较小。CoLM在常绿针叶林、稀树林地、草地、农田模拟感热、潜热通量的效果相对较好,在永久湿地、落叶阔叶林下模拟感热通量较差。本研究对CoLM2014在未来的改进和发展中提供了有用的参考。     

陆面过程模型CoLM与区域气候模式RegCM3的耦合及初步评估

陆面过程通过影响陆面和大气之间物质（如,水分）和能量的交换影响气候, 其参数化方案对数值天气预报、全球及区域气候模拟有重要影响。本研究利用对生物物理、生物化学过程考虑更全面的陆面模式Common Land Model(CoLM) 替代区域气候模式RegCM3原有的陆面模式BATS, 发展了耦合区域气候模式C－RegCM3; 将其应用于东亚地区典型洪涝年份夏季气候模拟以进行评估, 结果表明新耦合的模式C－RegCM3能合理模拟大尺度环流场、近地表气温和降水的分布特征, 对西北半干旱地区降水模拟比RegCM3有所改进。通过利用区域气候模式C－RegCM3及RegCM3对地表能量和水文过程模拟结果的比较, 发现在半干旱、半湿润过渡区C－RegCM3模拟的潜热增大、感热减小; 模拟的地表吸收太阳辐射差异较明显的地区位于模式模拟的主要雨区; C-RegCM3在上述过渡区模拟的夏季地表土壤湿度比RegCM3偏干, 这与它在过渡区降水模拟偏少、蒸散发模拟偏大相对应, 体现了该模式在半干旱、半湿润过渡带模拟出比RegCM3更明显的局地土壤湿度－降水－蒸散发之间的正反馈作用。     

Evaluation of the New Dynamic Global Vegetation Model in CAS-ESM

In the past several decades, dynamic global vegetation models(DGVMs) have been the most widely used and appropriate tool at the global scale to investigate vegetation-climate interactions. At the Institute of Atmospheric Physics, a new version of DGVM(IAP-DGVM) has been developed and coupled to the Common Land Model(CoLM) within the framework of the Chinese Academy of Sciences' Earth System Model(CAS-ESM). This work reports the performance of IAP-DGVM through comparisons with that of the default DGVM of CoLM(CoLM-DGVM) and observations. With respect to CoLMDGVM, IAP-DGVM simulated fewer tropical trees, more "needleleaf evergreen boreal tree" and "broadleaf deciduous boreal shrub", and a better representation of grasses. These contributed to a more realistic vegetation distribution in IAP-DGVM,including spatial patterns, total areas, and compositions. Moreover, IAP-DGVM also produced more accurate carbon fluxes than CoLM-DGVM when compared with observational estimates. Gross primary productivity and net primary production in IAP-DGVM were in better agreement with observations than those of CoLM-DGVM, and the tropical pattern of fire carbon emissions in IAP-DGVM was much more consistent with the observation than that in CoLM-DGVM. The leaf area index simulated by IAP-DGVM was closer to the observation than that of CoLM-DGVM; however, both simulated values about twice as large as in the observation. This evaluation provides valuable information for the application of CAS-ESM, as well as for other model communities in terms of a comparative benchmark.     

青藏高原地区不同下垫面陆面过程的数值模拟研究

利用陆面过程模式Common Land Model(CoLM),选取青藏高原上3个不同下垫面观测站(藏东南站、纳木错站和珠峰站)的观测资料,对这3个野外观测站进行了单点数值模拟试验。根据3个测站的试验数据,对模式中土壤孔隙度和饱和导水率进行了优化,针对青藏高原地区土壤层薄的特点,对模式中土壤分层方案进行了调整。结果表明,调整分层方案后的CoLM模式对3个测站土壤湿度的模拟性能较原分层方案有明显提高,平均偏差均减小0.014以上。但是与观测值相比,藏东南站土壤湿度的模拟整体偏低,纳木错站和珠峰站则整体偏高。对土壤温度而言,3个测站模拟与观测的相关系数都达到了0.9以上,珠峰站偏差较大,调整分层方案后模拟的偏差有一定的改进。模式较好地模拟了3个测站的净辐射、感热通量和潜热通量的日变化和季节变化情况,调整分层方案后潜热通量的改进最为明显。     

Improvements of a Dynamic Global Vegetation Model and Simulations of Carbon and Water at an Upland-Oak Forest

The interest in the development and improvement of dynamic global vegetation models （DGVMs）, which have the potential to simulate fluxes of carbon, water and nitrogen, along with changes in the vegetation dynamics, within an integrated system, has been increasing. In this paper, some numerical schemes and a higher resolution soil texture dataset were employed to improve the Sheffield Dynamic Global Vegetation Model （SDGVM）. Using eddy covariance-based measurements, we then tested the standard version of the SDGVM and the modified version of the SDGVM. Detailed observations of daily carbon and water fluxes made at the upland oak forest on the Walker Branch Watershed in Tennessee, USA offered a unique opportunity for these comparisons. The results revealed that the modified version of the SDGVM did a reasonable job of simulating the carbon and water flux and the variation of soil water content （SWC）. However, at the end of the growing season, it failed to simulate the effect of the limitations on the soil respiration dynamics and as a result underestimated this respiration. It was also noted that the modified version overestimated the increase in the SWC following summer rainfall, which was attributed to an inadequate representation of the ground water and thermal cycle.     

Earth System Model FGOALS-s2: Coupling a dynamic global vegetation and terrestrial carbon model with the physical climate system model

Earth System Models （ESMs） are fundamental tools for understanding climate-carbon feedback. An ESM version of the Flexible Global Ocean-Atmosphere-Land System model （FGOALS） was recently developed within the IPCC AR5 Coupled Model Intercomparison Project Phase 5 （CMIP5） modeling framework, and we describe the development of this model through the coupling of a dynamic global vegetation and terrestrial carbon model with FGOALS-s2. The performance of the coupled model is evaluated as follows. The simulated global total terrestrial gross primary production （GPP） is 124.4 PgC yr-I and net pri- mary production （NPP） is 50.9 PgC yr-1. The entire terrestrial carbon pools contain about 2009.9 PgC, comprising 628.2 PgC and 1381.6 PgC in vegetation and soil pools, respectively. Spatially, in the tropics, the seasonal cycle of NPP and net ecosystem production （NEP） exhibits a dipole mode across the equator due to migration of the monsoon rainbelt, while the seasonal cycle is not so significant in Leaf Area Index （LAI）. In the subtropics, especially in the East Asian monsoon region, the seasonal cycle is obvious due to changes in temperature and precipitation from boreal winter to summer. Vegetation productivity in the northern mid-high latitudes is too low, possibly due to low soil moisture there. On the interannual timescale, the terrestrial ecosystem shows a strong response to ENSO. The model- simulated Nifio3.4 index and total terrestrial NEP are both characterized by a broad spectral peak in the range of 2-7 years. Further analysis indicates their correlation coefficient reaches -0.7 when NEP lags the Nifio3.4 index for about 1-2 months.     

三个陆面模式对新疆地区陆面过程模拟的对比研究

利用实测的大气观测数据建立陆面模式驱动场, 分别驱动Biosphere-Atmosphere Transfer Scheme（BATS）, Land Surface Model（LSM）, Common Land Model（CoLM）三个陆面模式, 对新疆地区99个测站的陆面热力过程进行了off-line模拟研究, 对比分析了三个模式模拟结果的异同, 并利用观测的土壤温度检验了各个模式的模拟性能。结果表明： 三个模式对所吸收的太阳辐射、 地表感热通量、 地表潜热通量和Bowen比的模拟结果有明显差别, 其中CoLM模拟的太阳辐射、 感热通量和土壤温度比BATS和LSM大。用观测的土壤温度对三个模式模拟性能的检验表明 BATS模拟的土壤温度最接近观测值, 大约偏低1℃, 而LSM和CoLM的模拟结果偏高, CoLM模拟的年际变化趋势与观测值相关性达到0.877, 高于另外两个模式的模拟; BATS在天山以北模拟值偏低, CoLM在天山以南模拟结果偏高, LSM模拟的区域分布状况与实际观测最一致。这些结果为研究陆面模式在新疆地区的应用提供了一些参考, 对于深入了解新疆地区的陆面过程具有重要意义。     

Surface energy, water and carbon cycle in China simulated by the Australian community land surface model (CABLE)

Through an Australia-China climate change bilateral project, we analyzed results of 51-year global offline simulations over China using the Australian community atmosphere biosphere land exchange (CABLE) model, focusing on integrated studies of its surface energy, water and carbon cycle at seasonal, interannual and longer time-scales. In addition to the similar features in surface climatology between the CABLE simulation and those derived from the global land-surface data assimilation system, comparison of surface fluxes at a CEOP reference site in northeast China also suggested that the seasonal cycles of surface evaporation and CO₂ flux are reasonably simulated by the model. We further assessed temporal variations of model soil moisture with the observed variations at a number of locations in China. Observations show a soil moisture recharge–discharge mechanism on a seasonal time scale in central-east China, with soil moisture being recharged during its summer wet season, retained in its winter due to low evaporation demand, and depleted during early spring when the land warms up. Such a seasonal cycle is shown at both 50- and 100-cm soil depths in observations while the model only shows a similar feature in its lower soil layers with its upper layer soil moisture varying tightly with rainfall seasonal cycle. In the analysis of the model carbon cycle, the net primary productivity (NPP) has similar spatial patterns as the ones derived from an ecosystem model with remote sensing. The simulated interannual variations of NPP by CABLE are consistent with the results derived from remote sensing-based and process-based studies over the period of 1981–2000. Nevertheless an upward trend from observations is not presented in the model results. The model shows a downward trend primarily due to the constant CO₂ concentration used in the experiment and a large increase of autotrophic respiration caused by an upward trend in surface temperature forcing data. Furthermore, we have compared river discharge data from the model experiments with observations in the Yangtze and Yellow River basins in China. In the Yangtze River basin, while the observed interannual variability is reasonably captured, the model significantly underestimates its river discharge, which is consist with its overestimation of evaporation in the region. In the Yellow River basin, the magnitudes of the river discharge is similar between modeled and observed but its variations are less skillfully captured as seen in the Yangtze River region.     

A New Approach for Parameter Optimization in Land Surface Model

In this study,a new parameter optimization method was used to investigate the expansion of conditional nonlinear optimal perturbation (CNOP) in a land surface model (LSM) using long-term enhanced field observations at Tongyu station in Jilin Province,China,combined with a sophisticated LSM (common land model,CoLM).Tongyu station is a reference site of the international Coordinated Energy and Water Cycle Observations Project (CEOP) that has studied semiarid regions that have undergone desertification,salination,and degradation since late 1960s.In this study,three key land-surface parameters,namely,soil color,proportion of sand or clay in soil,and leaf-area index were chosen as parameters to be optimized.Our study comprised three experiments:First,a single-parameter optimization was performed,while the second and third experiments performed triple-and six-parameter optimizations,respectively.Notable improvements in simulating sensible heat flux (SH),latent heat flux (LH),soil temperature (TS),and moisture (MS) at shallow layers were achieved using the optimized parameters.The multiple-parameter optimization experiments performed better than the single-parameter experminent.All results demonstrate that the CNOP method can be used to optimize expanded parameters in an LSM.Moreover,clear mathematical meaning,simple design structure,and rapid computability give this method great potential for further application to parameter optimization in LSMs.