Evaluating Common Land Model Energy Fluxes Using FLUXNET Data

Given the crucial role of land surface processes in global and regional climates, there is a pressing need to test and verify the performance of land surface models via comparisons to observations. In this study, the eddy covariance measurements from 20 FLUXNET sites spanning more than 100 site-years were utilized to evaluate the performance of the Common Land Model （CoLM） over different vegetation types in various climate zones. A decomposition method was employed to separate both the observed and simulated energy fluxes, i.e., the sensible heat flux, latent heat flux, net radiation, and ground heat flux, at three timescales ranging from stepwise （30 rain） to monthly. A comparison between the simulations and observations indicated that CoLM produced satisfactory simulations of all four energy fluxes, although the different indexes did not exhibit consistent results among the different fluxes, A strong agreement between the simulations and observations was found for the seasonal cycles at the 20 sites, whereas CoLM underestimated the latent heat flux at the sites with distinct dry and wet seasons, which might be associated with its weakness in simulating soil water during the dry season. CoLM cannot explicitly simulate the midday depression of leaf gas exchange, which may explain why CoLM also has a maximum diurnal bias at noon in the summer. Of the eight selected vegetation types analyzed, CoLM performs best for evergreen broadleaf forests and worst for croplands and wetlands.

北京大学陆面过程模式PKULM(Peking University Land Model)介绍及检验

陆面过程模式是气候模式和天气模式的核心组成部分之一.在土壤—植被—大气耦合模式(Soil-PlantAtmosphere Model,SPAM)的基础上,发展了新一代北京大学陆面过程模式PKULM(Peking University Land Model).本文首先介绍了PKULM的辐射传输、湍流输送、光合作用、土壤水热输送等过程的参数化方案;采用隐式迭代计算框架,发展并应用了一个快速的线性方程组求解算法,提高了模式计算稳定性;提出并使用了二分搜索算法计算气孔阻抗,避免了CLM(Community Land Model)等使用的迭代方法在干旱区不稳定的情况,提高了模式的适用性;采用水势为基础的土壤水分扩散方程,使模式能够模拟土壤饱和区的水分输送过程,为进一步与水文过程模式耦合奠定了基础;还发展了一个地表积水与径流过程的机理模型,提高了模式对地表水分平衡过程的模拟能力;最后,使用"中国西北干旱区陆—气相互作用观测试验"平凉站的资料对模式进行了检验并与NOAH(National Center for Environmental Prediction,Oregon State University,Air Force,and Hydrology Lab model)陆面过程模式的模拟结果进行了比较,结果表明PKULM能够较好地模拟西北半干旱区农田下垫面地气交换过程.

Evaluation of the Community Microwave Emission Model Coupled with the Community Land Model over East Asia

The Community Microwave Emission Model (CMEM) developed by the European Centre for Me-dium-Range Weather Forecasts (ECMWF) can provide a link between surface states and satellite observations and simulate the passive microwave brightness temperature of the surface at low frequencies (from 1 GHz to 20 GHz).This study evaluated the performance of the CMEM cou-pled with the Community Land Model (CLM) (CMEM-CLM) using C-band (6.9 GHz) microwave brightness temperatures from the Advanced Microwave Scanning Radiometer on Earth Observing System (AMSR-E) over East Asia.Preliminary results support the argument that the simulated brightness temperatures of CMEM-CLM from July 2005 to June 2010 are comparable to AMSR-E observational data.CMEM-CLM performed better for vertical polarization,for which the root mean square error was approximately 15 K,compared to over 30 K for horizontal polarization.An evaluation performed over seven sub-regions in China indicated that CMEM-CLM was able to capture the temporal evolution of C-band brightness temperatures well,and the best correlation with AMSR-E appeared over western Northwest China (over 0.9 for vertical polarization).However,larger biases were found over southern North China and the middle and lower reaches of the Yangtze River.

Land Response to Atmosphere at Different Resolutions in the Common Land Model over East Asia

Towards a better understanding of hydrological interactions between the land surface and atmosphere, land surface mod- els are routinely used to simulate hydro-meteorological fluxes. However, there is a lack of observations available for model forcing, to estimate the hydro-meteorological fluxes in East Asia. In this study, Common Land Model （CLM） was used in offline-mode during the summer monsoon period of 2006 in East Asia, with different forcings from Asiaflux, Korea Land Data Assimilation System （KLDAS）, and Global Land Data Assimilation System （GLDAS）, at point and regional scales, separately. The CLM results were compared with observations from Asiaflux sites. The estimated net radiation showed good agreement, with r = 0.99 for the point scale and 0.85 for the regional scale. The estimated sensible and latent heat fluxes using Asiaflux and KLDAS data indicated reasonable agreement, with r = 0.70. The estimated soil moisture and soil temperature showed similar patterns to observations, although the estimated water fluxes using KLDAS showed larger discrepancies than those of Asiaflux because of scale mismatch. The spatial distribution of hydro-meteorological fluxes according to KLDAS for East Asia were compared to the CLM results with GLDAS, and the GLDAS provided online. The spatial distributions of CLM with KLDAS were analogous to CLM with GLDAS, and the standalone GLDAS data. The results indicate that KLDAS is a good potential source of high spatial resolution forcing data. Therefore, the KLDAS is a promising alternative product, capable of compensating for the lack of observations and low resolution grid data for East Asia.

Elucidating Dominant Factors Affecting Land Surface Hydrological Simulations of the Community Land Model over China

In order to compare the impacts of the choice of land surface model(LSM)parameterization schemes,meteorological forcing,and land surface parameters on land surface hydrological simulations,and explore to what extent the quality can be improved,a series of experiments with different LSMs,forcing datasets,and parameter datasets concerning soil texture and land cover were conducted.Six simulations are run for the Chinese mainland on 0.1°×0.1°grids from 1979 to 2008,and the simulated monthly soil moisture(SM),evapotranspiration(ET),and snow depth(SD)are then compared and assessed against observations.The results show that the meteorological forcing is the most important factor governing output.Beyond that,SM seems to be also very sensitive to soil texture information;SD is also very sensitive to snow parameterization scheme in the LSM.The Community Land Model version 4.5(CLM4.5),driven by newly developed observation-based regional meteorological forcing and land surface parameters(referred to as CMFD_CLM4.5_NEW),significantly improved the simulations in most cases over the Chinese mainland and its eight basins.It increased the correlation coefficient values from 0.46 to 0.54 for the SM modeling and from 0.54 to 0.67 for the SD simulations,and it decreased the root-mean-square error(RMSE)from 0.093 to 0.085 for the SM simulation and reduced the normalized RMSE from 1.277 to 0.201 for the SD simulations.This study indicates that the offline LSM simulation using a refined LSM driven by newly developed observation-based regional meteorological forcing and land surface parameters can better model reginal land surface hydrological processes.

Preliminary Assessment of the Common Land Model Coupled with the IAP Dynamic Global Vegetation Model

The Common Land Model(CoLM) was coupled with the IAP Dynamic Global Vegetation Model(IAPDGVM), and the performance of this combined CoLMIAP model was evaluated. Offline simulations using both the original Common Land Model(CoLM-LPJ) and CoLM-IAP were conducted. The CoLM-IAP coupled model showed a significant improvement over CoLMLPJ, as the deciduous tree distribution decreased over temperate and boreal regions, while the distribution of evergreen trees increased over the tropics. Some biases in CoLM-LPJ were preserved, including the overestimation of evergreen trees in tropical savanna, the underestimation of boreal evergreen trees, and the absence of boreal shrubs. However, most of these biases did not exist in a further coupled simulation of IAP-DGVM with the Community Land Model(CLM), for which the parameters of IAP-DGVM were optimized. This implies that further improvement is needed to deal with the differences between CoLM and CLM in parameterizations of landbased physical and biochemical processes.

Origin and advances in implementing blowing-snow effects in the Community Land Model

Snow cover on the Tibetan Plateau(TP)is closely related to regional and continental biological and hydrological processes.The vast snow cover,special climatic conditions,and sparse vegetative cover over the TP facilitate the occurrence of blowing snow,leading to substantial heterogeneities in the snow cover and great promotion in the moisture supply from the land surface to the overlying atmospheric boundary layer.However,blowing-snow processes are significantly misrepresented or even neglected in current models,which causes considerable uncertainties of numerical model simulations and leads to erroneous estimates of snow-related processes in mountainous terrain.We present in this paper a brief review of our work in the past 5 years to serve as a basis for further development and improvement of the land-surface model.These studies can be divided into three parts:detection of the problems,development of the land-surface model,and application of the coupled model over the TP(the logical framework is presented in Figure 1).The origin and advances in the development of a land-surface model with consideration of blowing-snow effects are described herein;and the importance of blowing-snow processes in the land-surface model,especially over the TP,is highlighted.We expect that the blowingsnow studies over the TP will play a key role in documenting and understanding the land-surface processes(LSPs)and the cryospheric changes over the TP.

Estimates of net primary productivity and actual evapotranspiration over the Tibetan Plateau from the Community Land Model version 4.5 with four atmospheric forcing datasets

The accuracy of the simulation of carbon and water processes largely relies on the selection of atmospheric forcing datasets when driving land surface models(LSM).Particularly in high-altitude regions,choosing appropriate atmospheric forcing datasets can effectively reduce uncertainties in the LSM simulations.Therefore,this study conducted four offline LSM simulations over the Tibetan Plateau(TP)using the Community Land Model version 4.5(CLM4.5)driven by four state-of-the-art atmospheric forcing datasets.The performances of CRUNCEP(CLM4.5 model default)and three other reanalysis-based atmospheric forcing datasets(i.e.ITPCAS,GSWP3 and WFDEI)in simulating the net primary productivity(NPP)and actual evapotranspiration(ET)were evaluated based on in situ and gridded reference datasets.Compared with in situ observations,simulated results exhibited determination coefficients(R2)ranging from 0.58 to 0.84 and 0.59 to 0.87 for observed NPP and ET,respectively,among which GSWP3 and ITPCAS showed superior performance.At the plateau level,CRUNCEP-based simulations displayed the largest bias compared with the reference NPP and ET.GSWP3-based simulations demonstrated the best performance when comprehensively considering both the magnitudes and change trends of TP-averaged NPP and ET.The simulated ET increase over the TP during 1982-2010 based on ITPCAS was significantly greater than in the other three simulations and reference ET,suggesting that ITPCAS may not be appropriate for studying long-term ET changes over the TP.These results suggest that GSWP3 is recommended for driving CLM4.5 in conducting long-term carbon and water processes simulations over the TP.This study contributes to enhancing the accuracy of LSM in water-carbon simulations over alpine regions.

Seasonal constraint of dynamic water temperature on riverine dissolved inorganic nitrogen transport in land surface modeling

水体温度变化对河流可溶性无机氮(DIN)输送有着强烈控制作用.然而,在全球尺度上河流DIN输送量对水温度变化的响应尚不清楚.因此,本文基于陆面过程模式,耦合河流水温估算和DIN传输方案,设定有,无动态水温情景,对比研究陆面模拟中水温变化对河流DIN通量变化的影响.结果表明:在考虑水温动态变化后,在30°N和30°S之间, DIN通量年振幅减小5%–25%.在中国东部地区,水温动态变化使河流DIN通量在夏季减少1%–3%,在冬季增加1%–5%,对DIN通量具有明显的季节性约束作用,表明动态水温的表达在河流DIN输送模拟中的重要性.

Modelling Land Use/Land Cover Change of River Rwizi Catchment, South-Western Uganda Using GIS and Markov Chain Model

Analysis of catchment Land use/Land cover (LULC) change is a vital tool in ensuring sustainable catchment management. The study analyzed land use/land cover changes in the Rwizi catchment, south western Uganda from 1989-2019 and projected the trend by 2040. Landsat images, field observations, key informant interviews and focus group discussions were used to collect data. Changes in cropland, forestland, built up area, grazing land, wetland and open water bodies were analyzed in ArcGIS version 10.2.2 and ERDAS IMAGINE 14 software and a Markov chain model. All the LULC classes increased in area except grazing land. Forest land and builtup area between 2009-2019 increased by 370.03% and 229.53% respectively. Projections revealed an increase in forest land and builtup area by 2030 and only built up area by 2040. LULCC in the catchment results from population pressure, reduced soil fertility and high value of agricultural products.

Parameter Sensitivities of the Community Land Model at Two Alpine Sites in the Three-River Source Region

The three-river source region plays an important role on China’s ecological security and Asia’s water supply. Historically, the region has experienced severe ecological degradation due to climate change and human activities. Reasonable simulations of the energy and water cycles are essential to predict the responses of land surface processes to future climate change. Current land surface models involve empirical functions that are associated with many parameters. These parameter uncertainties will largely affect the simulation when applied to a new domain. The Community Land Model(CLM) is a widely used land surface model, and version 5.0 is the newest version. Compared to the prior version CLM4.5, CLM5.0 has largely updated plant hydraulic and stomatal conductance schemes. How these changes affect parameter sensitivities is unknown. In our work, we tested 17 key parameters in CLM4.5 and 19 parameters in CLM5.0 at two eddy flux sites in the three-river source region: the Maqu and Maduo sites. We adopted the simplest one-at-a-time changes on each parameter and quantified their sensitivities by the parameter effect(PE).We found that the Maqu site was more sensitive to vegetation parameters, while the Maduo site was more sensitive to the initial soil water content in both CLM4.5 and CLM5.0. This is because Maduo grid cell has wetland that does not respond to vegetation parameters in CLM, which may not reflect the reality. Further model development on wetland vegetation parameterization is important. Our validation on the default simulation showed CLM5.0 did not always improve the simulations. The largest difference between CLM5.0 and CLM4.5 was that soil moisture(SM) showed a much stronger decrease in response to a higher leaf area index(LAI) in CLM5.0 than in CLM4.5, suggesting that SM is more sensitive to vegetation changes in CLM5.0.

A Land Surface Model(IAP94) for Climate Studies PartI:Formulation and Validation in Off-line Experiments

The IAP (Institute of Atmospheric Physics) land-surface model (IAP94) is described. This model is a comprehensive one with detailed description for the processes of vegetation, snow and soil. Particular attention has been paid to the cases with three water phases in the surface media. On the basis of the mixture theory and the theory of fluid dynamics of porous media, the system of universal conservational equations for water and heat of soil, snow and vegetation canopy has been constructed. On this background, all important factors that may affect the water and heat balance in media can be considered naturally, and each factor and term possess distinct physical meaning. In the computation of water content and temperature, the water phase change and the heat transportation by water flow are taken into account. Moreover, particular attention has been given to the water vapor diffusion in soil for arid or semi-arid cases, and snow compaction. In the treatment of surface turbulent fluxes, the difference between aerodynamic and thermal roughness is taken into account. The aerodynamic roughness of vegetation is calculated as a function of canopy density, height and zero-plane displacement. An extrapolation of log linear and exponential relationship is used when calculating the wind profile within canopy. The model has been validated against field measurements in off-line simulations. The desirable model′s performance leads to the conclusion that the IAP94 is able to reproduce the main physical mechanisms governing the energy and water balances in the global land surface. Part II of the present study will concern the validation in a 3-D experiment coupled with the IAP Two-Level AGCM.

Impact of Spin-up Forcing on Vegetation States Simulated by a Dynamic Global Vegetation Model Coupled with a Land Surface Model

A dynamic global vegetation model （DGVM） coupled with a land surface model （LSM） is generally initialized using a spin-up process to derive a physically-consistent initial condition. Spin-up forcing, which is the atmospheric forcing used to drive the coupled model to equilibrium solutions in the spin-up process, varies across earlier studies. In the present study, the impact of the spin-up forcing in the initialization stage on the fractional coverages （FCs） of plant functional type （PFT） in the subsequent simulation stage are assessed in seven classic climate regions by a modified Community Land Model’s Dynamic Global Vegetation Model （CLM-DGVM）. Results show that the impact of spin-up forcing is considerable in all regions except the tropical rainforest climate region （TR） and the wet temperate climate region （WM）. In the tropical monsoon climate region （TM）, the TR and TM transition region （TR-TM）, the dry temperate climate region （DM）, the highland climate region （H）, and the boreal forest climate region （BF）, where FCs are affected by climate non-negligibly, the discrepancies in initial FCs, which represent long-term cumulative response of vegetation to different climate anomalies, are large. Moreover, the large discrepancies in initial FCs usually decay slowly because there are trees or shrubs in the five regions. The intrinsic growth timescales of FCs for tree PFTs and shrub PFTs are long, and the variation of FCs of tree PFTs or shrub PFTs can affect that of grass PFTs.

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 influ- ences 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 simulat- ing 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~ 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 im- proved,

Improving the Vegetation Dynamic Simulation in a Land Surface Model by Using a Statistical-dynamic Canopy Interception Scheme

Canopy interception of incident precipitation, as a critical component of a forest＇s water budget, can affect the amount of water available to the soil, and ultimately vegetation distribution and function. In this paper, a statistical-dynamic approach based on leaf area index and statistical canopy interception is used to parameterize the canopy interception process. The statistical-dynamic canopy interception scheme is implemented into the Community Land Model with dynamic global vegetation model （CLM-DGVM） to improve its dynamic vegetation simulation. The simulation for continental China by the land surface model with the new canopy interception scheme shows that the new one reasonably represents the precipitation intercepted by the canopy. Moreover, the new scheme enhances the water availability in the root zone for vegetation growth, especially in the densely vegetated and semi-arid areas, and improves the model＇s performance of potential vegetation simulation.

Case Analyses and Numerical Simulation of Soil ThermalImpacts on Land Surface Energy Budget Based on anOff-Line Land Surface Model

The statistical relationship between soil thermal anomaly and short-term climate change is presented based on a typical case study. Furthermore, possible physical mechanisms behind the relationship are re-vealed through using an off-line land surface model with a reasonable soil thermal forcing at the bottom of the soil layer. In the first experiment, the given heat flux is 5 W m2 at the bottom of the soil layer (in depth of 6.3 m) for 3 months, while only a positive ground temperature anomaly of 0.06°C can be found compared to the control run. The anomaly, however, could reach 0.65°C if the soil thermal conductivity was one order of magnitude larger. It could be even as large as 0.81°C assuming the heat flux at bottom is 10 W m-2. Mean-while, an increase of about 10 W m−2 was detected both for heat flux in soil and sensible heat on land sur-face, which is not neglectable to the short-term climate change. The results show that considerable response in land surface energy budget could be expected when the soil thermal forcing reaches a certain spatial-tem-poral scale. Therefore, land surface models should not ignore the upward heat flux from the bottom of the soil layer. Moreover, integration for a longer period of time and coupled land-atmosphere model are also necessary for the better understanding of this issue.

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.

Evaluation of the WRF Model with Different Land Surface Schemes: A Drought Event Simulation in Southwest China during 2009–10

The authors examined the performance of version 3.4.1 of the Weather Research and Forecasting Model(WRF) with various land surface schemes in simulating a severe drought event in Southwest China. Five numerical experiments were completed using the Noah land surface scheme, the Pleim-Xiu land surface scheme, the Noah-MP land surface schemes, the Noah- MP scheme with dynamic vegetation, and the Noah-MP scheme with dynamic vegetation and groundwater processes. In general, all the simulations reasonably reproduced the spatial and temporal variations in precipitation, but significant bias was also found, especially for the spatial pattern of simulated precipitation. The WRF simulations with the Noah-MP series land surface schemes performed slightly better than the WRF simulation with the Noah and Pleim-Xiu land surface schemes in reproducing the severe drought events in Southwest China. The leaf area index(LAI) simulated by the different land surface schemes showed significant deviations in Southwest China. The Pleim-Xiu scheme overestimated the value of LAI by a factor of two. The Noah-MP scheme with dynamical vegetation overestimated the magnitude of the annual cycle of the LAI, although the annual mean LAI was close to observations. The simulated LAI showed a long-term lower value from autumn 2009 to spring 2010 relative to normal years. This indicates that the LAI is a potential indictor to monitor drought events.

Simulation of the Asian Monsoon by IAP AGCM Coupled with an Advanced Land Surface Model(IAP94)

In this paper, the global and regional features of the seasonal variation of general circulation, and especially the Asian monsoon simulated by the Institute of Atmospheric Physics Two-level AGCM coupled with a sophisticated land-surface model (IAP94-GCM) are presented and compared with the observation. The comparison is made by using the equilibrium multiyear seasonal cycle climate from a 100-year integration. In the integration sea surface temperature (SST) and sea ice are taken from the observed climatological data (with seasonal variation) because our purpose is to see the improvement of simulation due to the coupling with an advanced land surface model. Overall, the IAP94-GCM provides a reasonably realistic simulation of the interseasonal and intraseasonal climatology of the Asian monsoon and yields an important information that sheds light on the thermal underpinning and the thermodynamics of the seasonal and even multiscale variabilities associated with the Asian summer monsoon.

A Land Surface Model(IAP94)for Climate Studies Part II: Implementation and Preliminary Results of Coupled Model with IAP GCM

The Institute of Atmospheric Physics Land Surface Model (IAP94) has been incorporated into the IAP two-level atmospheric general circulation model (IAP GCM). Global and regional climatology averaged over the last 25 years of 100 year integrations from the IAP GCM with and without IAP94 (“bucket” scheme) is compared. The simulated results are also compared with the reanalysis data. Major findings are: \ \ (1) The IAP GCM simulation without IAP94 has extensive regions of warmer than observed surface air temperatures, while the simulation with IAP94 very much improves the surface air temperature. \ \ (2) The IAP GCM simulation with IAP94 gives improvement of the simulated precipitation pattern and intensity, especially the precipitation of East Asian summer monsoon and its intraseasonal migration of the rainbelts. \ \ (3) In five selected typical regions, for most of the surface variables such as surface air temperature, precipitation, precipitation minus evaporation, net radiation, latent heat flux and sensible heat flux, the IAP GCM with IAP94 provides better simulations.