Simulating the Energy and Water Fluxes from Two Alkaline Desert Ecosystems over Central Asia

The Central Asia region is covered by vast desert ecosystems, where the characteristic of energy and water fluxes is different from other humid ecosystems. The application of land surface models (LSMs) in arid and semiarid ecosystems was largely limited. This paper presents a detailed evaluation of Common Land Model (CoLM) at two eddy covariance (EC) sites in alkaline desert ecosystems over Central Asia. Simulations of the net radiation (Rnet), latent heat flux (Qle), sensible heat flux (Qh), and soil temperature showed that refined estimate of roughness length (Z0m) significantly improved the performance of CoLM in simulating turbulent heat fluxes. Qle was increased but Qh was decreased, which were in better agreement with the observations from EC system.The results indicated that accurate parameterization of Z0m is of crucial importance for predicting energy and water fluxes in LSM when applied in Central Asia desert ecosystems. Sensitivity analysis regarding leaf area index (LAI), Z0m, and albedo (α) showed that Rnet is very sensitive to α butQle,Qh, and soil temperature (Tsoil) are sensitively varying with the estimate of Z0m at the two EC sites over Central Asia.


Introduction
Terrestrial ecosystem is an important part of the earth system as it controls land-atmosphere interactions.Land surface model, which described plant physiological behavior in relation to soil and atmospheric processes, was widely used to quantify the land-atmosphere interactions for improving the predictability of Earth System Models (ESMs) [1][2][3][4].LSM consists of different biogeophysical and biogeochemical processes, which described energy, momentum, water, and CO 2 exchange between the atmosphere and the land surface.LSMs became increasingly complicated since Manabe introduced the first bucket model in 1969 [5].Increases in observations obtained from in situ instruments and satellites provided more opportunities to test, validate, and evaluate the performance of LSMs in specific terrestrial ecosystem [6,7].Despite the significant improvement of LSMs during the past few decades, large errors and uncertainties still existed [8].Errors and uncertainties underlying the LSMs can result from model structure, model initial conditions, model parameters, and model forcing data, and evaluating the performance of LSMs among different ecosystems is helpful to identify and diagnose the errors of the model in structure and parameters and hence to reduce the uncertainties [9].
Evaluating LSMs would provide a reliable indicator of performance of models under various climates and among different biomes.Comparing the model results with measurements is a routine approach in the assessment of a particular LSM [10].Application of LSMs was intensively implemented in humid and dense vegetation regions, but LSM was rarely investigated in dryland or sparse-vegetated areas (e.g., deserts ecosystems in Central Asia) [11,12].Absence of the exquisite description on specially ecophysiological characteristics of dryland ecosystems further enlarged the errors of the LSMs in simulating turbulent fluxes [13,14].Under humid areas, LSMs were found to perform well [12].However, around onethird of the world's area was belonging to dryland ecosystems [15], where the land surface fluxes, such as  ℎ ,  net , and 2 ground surface temperature are less dependent on vegetation parameters and soil hydraulic parameters, but more on energy-related soil parameters and surface parameters (surface albedo, the ground surface emissivity, aerodynamic roughness length, and thermal roughness length) [16].Unfortunately, application of LSMs in arid and semiarid ecosystems was largely limited.For example, recent studies have found that CoLM extremely underpredicted latent heat fluxes in dryland areas [13,17].Prediction errors can be derived from the forcing data, physics processes, and the parameterizations of land characteristic (which mainly included vegetation and soil effects) [2,18].Vegetation effects consist of different plant types that differ in leaf areas, root profile, height, optical properties, stomatal conductance, and roughness length.Soil properties can be soil texture (percent of sand and clay) and soil thermal and hydraulic parameters.The errors resulting from parameterizations of these processes will cause large errors for LSMs in modelling the turbulent fluxes.

Advances in Meteorology
The Central Asia dryland area is characterized by a typical arid continental climate, complex terrain, sparse vegetation, and high heterogeneity in land surface, which inevitably proposes high challenges on the application of LSMs [19].Investigating the energy and water exchange between land surface and atmosphere was of great importance to hydrological and ecological research [15,20].Few recent studies have evaluated the performance of the CoLM to reproduce energy and water vapor fluxes in China's desert ecosystem [13,21].They have found that root distribution and root water uptake function have significant effects on the performance of CoLM in estimating energy and water fluxes.However, ,  0 , and LAI are also expected to significantly impact the performance of LSMs, which largely influence energy budget, thermodynamic resistance, and turbulence fluxes.The errors resulting from these critical parameterizations will cause large error and uncertainties for LSMs in the estimates of sensible and latent heat fluxes.Unfortunately, a comprehensive evaluation of LSMs and their uncertainties caused by ,  0 , and LAI was never investigated in Central Asia dryland ecosystems.
Therefore, the main objectives of this research are (1) to evaluate the performance of CoLM of two desert ecosystems in Kazakhstan using EC observations during the growing seasons, for the purpose of investigating how well the CoLM is able to simulate the energy and water fluxes over Central Asia land surface, and (2) to analyze model sensitivities to albedo, roughness length, and leaf area index.

Site Description.
Two sets of EC equipment were established in representative dryland ecosystem in Kazakhstan in 2012 (Figure 1).One site is close to Aral Sea and the other is close to Balkhash Lake.The Aral Sea site (KZ-Ara, 61.08 ∘ E, 45.96 ∘ N) is located northeast of the Aral Sea and at the edge of Aralkum Desert.Within radius of 5 km of the KZ-Ara site, the dominant plant species are weedgrass and reed vegetation in combination with xerophytic and halophytic desert semishrubs and shrubs.The average fraction of vegetation coverage was about 30%.The dominant soil type was solonchak.Historical climatic records with long term observations showed that mean annual precipitation is 140 mm and mean annual air temperature is 8.3 ∘ C (Table 1).The Balkhash Lake site (KZ-Bal, 76.63 ∘ E, 44.57∘ N) is located between the Balkhash Lake and the Kapchagay Reservoir and between original deserts and oasis croplands.Within radius of 5 km of the KZ-Bal site, the dominant plant species are irrigated crops, grasses, and desert shrubs.The average fraction of vegetation coverage was about 40%.The soil in the KZ-Bal site is saline solonchak [22].Historical climatic records with long term observations showed that mean annual precipitation is 140 mm and mean annual air temperature is 5.7 ∘ C (Table 1).For detailed information about the two sites, refer to Li et al. [23].

Eddy Covariance and Ancillary Measurements.
At each site, a set of EC systems was used to measure energy, water vapor, and CO 2 fluxes between the land surface and the atmosphere.The eddy covariance consists of an open path gas analyzer (LI-7500, LICOR) and a three-dimensional sonic anemometer thermometer (Wind Master Pro, Gill Instruments, Lymington, UK).The EC system also measures air temperature and humidity (HMP45C, Campbell), precipitation (TE525MM, Texas Electronics, Dallas, TX, USA), short-wave and long-wave radiation (CNR-1, Kipp & Zonen, Delft, Netherlands), soil temperature (TCAV, Campbell), soil moisture content (CS616, Campbell Sci.), and soil heat flux (HFP01, Hukseflux, Delft, Netherlands).The open path gas analyzer and the three-dimensional sonic anemometer thermometer are mounted at 2.0 m above the ground.Soil temperature, soil moisture content, and soil heat flux () were measured at 20, 40, 60, and 80 cm depth below the ground.Therefore, this study chose the interpolated soil temperature in CoLM on the corresponding depth for comparisons.These data were logged every 30 min by the data logger at 10 Hz, and all variables were carried out with data processing and gap filling.The data was used for the analysis in this study during growing seasons at the two sites and covered the period from 30 April to 18 August 2012 at the KZ-Ara site and from 23 May to 6 September 2012 at the KZ-Bal site.[18] is one of the widely used land surface models, which was originally proposed to provide a framework for a truly community developed land component of the National Centre for Atmospheric Research (NCAR) Community Climate System Model (CCSM) [24,25], and it was subsequently developed by an open collaboration of many scientists.A variety of comprehensive multiyear point observational data over different regions of the world have been used [26][27][28].These data were included in the project for the Intercomparison of Land Parameterization Schemes [29], and the results from these extensive tests were useful for CoLM improvement.

Common Land Model. CoLM
CoLM simulates the exchange of energy, mass, and momentum between the atmosphere and terrestrial biosphere through a wide range of ground and canopy biogeophysical processes and allows for the simulation of explicit biophysical mechanisms, including the representation of thermodynamic, hydrologic, and physiological processes [18,30].In CoLM, soil temperatures are predicted using a heat diffusion equation in 10 soil layers, and the soil thermal conductivity depends on soil water density.Surface evapotranspiration consists of evaporation of precipitation intercepted by leaves, canopy transpiration, and bare soil evaporation.The sensible heat and water vapor fluxes between the surface and the reference height can be written in the following form: where  ℎ is sensible heat flux and  is water vapor flux.  is air density and   is specific heat for dry air (1004.67J kg −1 K −1 );   and   are air potential temperature at reference height and surface potential temperature, respectively;   and   are water vapor specific humidity at reference height and surface air humidity, respectively. ℎ and   are the aerodynamic resistances for heat and water vapor, respectively, which are crucial parameters in determining  ℎ and , calculated by Monin-Obukhov similarity theory [31], and mainly dependent on the thermal roughness length ( 0ℎ ) and  0 [16], as well as the meteorological forcing variables in CoLM, such as air temperature, wind speed, and precipitation.

Model Simulations and Sensitivity
Analysis.The version of CoLM, which integrated an optimal root water uptake function for dryland ecosystem [21], was used (termed S0) in this study.To investigate the effect of roughness length on the performance of CoLM, a simulation driven by an empirical estimate of roughness length (S1) was conducted (Table 2).In S0 simulation, momentum roughness length, thermal roughness length, and the water vapor roughness length ( 0 ) were all defined the same in CoLM for bare soil ( 0 =  0ℎ =  0 = 0.05).However, some research in dryland ecosystems indicated that the abovementioned roughness length values have different magnitude [4,[32][33][34].
As in the previous studies, Chen et al. [16] found that the  surface flux and temperature have different sensitivities to  0ℎ and  0 and the sensible heat flux is very sensitive to the parameterization schemes of  0ℎ in arid regions [16].
A common method to calculate  0 and  0ℎ is where / is a stability parameter ( is the Monin-Obukhov length and  is the observational height) and   (/) is the stability function of wind profile and becomes 0 at the neutral condition;  is the average wind speed and  * is the surface friction velocity;  ℎ (/) is the stability function of the temperature profile and becomes 0 at the neutral condition,  is von Karman constant and equals 0.4, generally,  is the air temperature, and   is the surface temperature.The relationship between  0 and  0ℎ can be described as where  −1 can be obtained from the bulk transfer equation as where  is the observed sensible heat flux,  is the air density,   is the specific heat for dry air, and  0 is the zero plane displacement.Thus, the relationship among  0 ,  0ℎ , and  0 is described as where the excess resistance to heat transfer  −1 is important to the sensible heat exchange between land surface and atmosphere, and there are linear correlations between  −1 and surface temperature. 0 is physically related to the geometric roughness of surface elements and can be derived from the wind speed and temperature profiles.Bao et al. [4] and Yang et al. [34] argued that this scheme overestimated  0ℎ and would misestimate the energy and water fluxes.Momentum transport is more efficient than heat transport, due to the influence of pressure fluctuation, because individual roughness elements may enhance the momentum flux through form drag with little contribution to the area-averaged heat flux [35].
An appropriate method suggested by Zhang et al. [36] was employed to estimate  0 at desert ecosystems.In desert,  0 was estimated as 0.0019 ± 0.0071 m and  0ℎ and  0 were at same order but almost one magnitude lower than  0 [36].In this study,  0 was set as 0.002 m, and  0ℎ =  0 = 0.0002 m for bare soil at the two sites (Table 2).
This research further investigated the sensitivities of the CoLM to LAI,  0 , and .Therefore, other six simulations, by increasing or decreasing the value of each parameter, were conducted.The specifications of all simulations are listed in Table 3.In order to avoid the cross-influence of each parameter, all sensitivity tests took the simulation S0 as a reference.[35] was used to give an overall evaluation of energy balance closure by averaging over random errors in the half-hour measurements at two flux tower sites, and it was calculated by

Statistical Analysis. Energy balance ratio (EBR)
where  is the number of half hours of data.The values of EBR close to 1 indicate the best degree of energy balance closure.Additionally, coefficient of determination ( 2 ), root mean square error (RMSE), slope (  ), and intercept ( 0 ) are used to justify the performance of the model.The Taylor diagram [37] was used to quantify the degree of the sensitivities of the model to management parameters.RMSE,  2 , and standard error (STD) are used in the Taylor diagram.The output of the model simulation is specified by a single point, with the STD being the polar axis and  the polar angle.The "reference" point represents observations, and the other points refer to the model results from the simulations of sensitivity tests.The distances from the reference point to the other points representing the consequence of the relationship indicate the RMSE.The higher  and the smaller the STD and RMSE, the better the agreement between model and data.When comparing two simulations with different parameter values, the longer the distance between the two simulation points, the greater the sensitivity to that parameter.

Energy Balance
Closure.The slopes of the linear regression between the observed  le +  ℎ and  net −  were 0.76 and 0.95 at KZ-Ara and KZ-Bal, respectively.The coefficient of determination ( 2 ) of the observed  le +  ℎ and  net −  was 0.91 and 0.97, and the root mean square error (RMSE) was 46.45 and 35.19 W m −2 , respectively (Figure 2).Energy balance ratio (EBR) at KZ-Ara and KZ-Bal was 1.11 and 1.06, respectively.

3.2.
Modelled   ,   ,  ℎ , and   .Figure 3 shows the comparisons between the measurements and the simulations of the reference model (S0) and the model with refined roughness length (S1) for  net ,  le , and  ℎ at the two Kazakhstan sites.The reference model (S0) significantly underestimated the latent heat flux and overestimated the sensible heat flux at both sites.However, the performance of the refined roughness length (S1) was largely improved in simulating turbulent heat fluxes.The latent heat flux was increased and sensible heat flux was decreased obviously.At the KZ-Ara site, RMSE for  net decreased from 63.6 to 39.5 W m −2 . 2 values for  le given by the two simulations (S0 and S1) were 0.36 and 0.61, respectively, and RMSE decreased from 41.15 in S0 to 33.3 W m −2 in S1 (Table 4).For  ℎ simulations, the RMSE for S0 was 121.24 W m −2 , as compared to 58.47 W m −2 for S1 (Table 4).The results indicated that the simulation with refined roughness length (S1) significantly improved the performance of the model for both  net and  le and particularly for  ℎ .At the KZ-Bal site, the performance of the simulation with refined roughness (S1) was also greatly improved,  2 values for  le given by S0 and S1 were 0.9 and 0.92, respectively, and RMSE decreased from 46.42 in S0 to 43.11 W m −2 in S1.  2 values for  ℎ given by the two simulations were 0.7 and 0.67, respectively, and RMSE decreased from 77.38 in S0 to 37.97 W m −2 in S1 (Table 4, Figure 3).To further investigate the effects of refined roughness on the energy fluxes simulations, Figure 4 shows the mean diurnal turbulent fluxes during growing seasons at the two sites.Diurnal variations of the three components of energy fluxes showed typical characteristics at the KZ-Ara site; that is,  net >  ℎ >  le (Figure 4(a)).The KZ-Bal site is located between oasis croplands and original deserts.Although the amount of the average annual precipitation is similar to KZ-Ara site, the characteristic of the energy fluxes allocation showed higher latent heat flux and lower sensible heat flux compared to the KZ-Ara site (Figure 4(b)).At the KZ-Ara site, S0 overestimated  net with the peak value of 520 W m −2 , as compared to the observed peak value 450 W m −2 .However, the simulation with refined roughness (S1) produced a good agreement for  net between the simulation and the measurements.In addition, S0 overestimated  ℎ with the peak value of 400 W m −2 at noontime, as compared to 200 W m −2 for the measurements.S1 improved the simulation for  ℎ as well.At KZ-Bal site, both S0 and S1 agreed better with observations for  net and  le .S0 significantly overestimated  ℎ at the daytime.However, the simulation with refined roughness (S1) produced better agreement with the observations.Figure 5 showed the comparisons of soil temperature between the observation and the simulations of the reference model (S0) and the model with refined roughness length (S1) at each site.Although there are significant challenges when it comes to validating soil temperature from LSM because of the high sensitivities of simulated and measured soil temperature to soil texture, moisture conditions, and the limitations of the measurement, the results from the reference model (S0) and the model with refined roughness length (S1) appeared to be realistic and generally captured the seasonal variation for soil temperature.It was found that the modelled  soil in May was underestimated, and significantly improved in June and July, but overestimated in August.However, S1 produced higher values of soil temperature than S0.Additionally, the CoLM simulated a smaller vertical soil temperature gradient as shown that the contour for the simulations was more sparse compared with the observation.

Sensitivity of CoLM to LAI, 𝑍 0𝑚
, and .In order to investigate the sensitivities of CoLM, the values of the parameters were adjusted in large ranges.Key parameters consist of LAI,  0 , and  with regard to the simulation of  net ,  le ,  ℎ , and  soil at the two EC sites.A reference model (S0) and six independent sensitivity tests (S2-S7, see Table 3) were conducted at each site.This research only chose the modelled soil temperatures at the depth of 20 cm for sensitivity tests in this study.
The modelled  le and  ℎ in CoLM were divided into two parts, the fluxes on vegetation leaves and the fluxes on the ground.Thus,  le showed a more complicated sensitivity to LAI.Since the leaf temperature increased and ground temperature decreased with the increase in LAI, the evapotranspiration from the leaves improved and the ground evaporation weakened, but the total latent heat flux was increased and exhibited a strong sensitivity for  le at KZ-Ara.When it came to KZ-Bal, the vegetation coverage was larger than KZ-Ara, while continuing to increase LAI has slightly further improved the performance of  le .Similarly,  ℎ and  soil showed a strong sensitivity to LAI.The more the vegetation, the more the solar radiation intercepted.When it came to  0 , the values of  net and  ℎ were increased but  le and  soil were decreased with the increase in  0 .Figure 6 indicated that  le ,  ℎ , and  soil were highly sensitive to  0 .Taking the KZ-Ara site as an example,  2 values for the modelled  le improved from 0.4 to 0.7 when  0 decreased from 0.5 (S4) to 0.005 (S5).Figure 6 also showed a strong sensitivity of albedo to  net ,  ℎ , and  soil .The imprecise settings of albedo may enlarge the errors for sensible heat flux in CoLM.Ground sensible heat flux was negatively correlated to albedo.The increase in the surface albedo decreased the solar radiation absorbed by soil, and soil temperature was decreased.However,  le was slightly sensitive to albedo (Figure 6).
These sensitivity analysis results demonstrated that the improvement in model performance observed in S2-S7 was significantly affected by the values of these three parameters in CoLM, further justifying the significance of these key parameters (LAI,  0 , and albedo) to the Common Land Model.

Discussion
Arid and semiarid areas cover approximately one-third of the global terrestrial land surfaces [15].Central Asia has    Hour Hour  vast proportion of dryland ecosystems, where climate was featured as hot and dry during growing season [19] and the vegetation is sparse.The characteristic of dryland ecosystems is significantly different from those in humid areas in terms of surface energy budget.Many researchers have focused on ecosystem functioning and structure in response to climate change.Kariyeva et al. [38] examined spatiotemporal change patterns and interactions between terrestrial phenology and climate change in Central Asia during the period 1981-2008.Lioubimtseva et al. [19] have found that local and regional human impacts in arid zones can significantly modify land surface albedo, as well as water exchange and nutrient cycles that could have essential impacts on the climate at both the regional and global scales.Such kind of analyses advanced our understanding on the response of dryland ecosystem behaviour and functioning to climate change.However, applying LSM into dryland ecosystems was far more limited.Recently, the CoLM has been validated at one desert shrub ecosystem in Xinjiang, China, and the results found that root fraction and root water uptake have important effects on the performance of CoLM in simulating energy and water fluxes [13,21].In this study, a refined parameterization of CoLM has been evaluated at two newly built EC sites in Kazakhstan.
The most commonly used technique to obtain land surface turbulent fluxes is measurement of eddy covariance system, and the analyzer was based on flux footprint models.The footprint concept is the probability that a scalar coming from a given elemental source reaches the measurement point.Footprint models describe the relationship between the spatial distribution of surface sources and the measured signal using footprint functions.Several flux footprint models have been designed [39][40][41][42].But most of them cannot account for inhomogeneous turbulence or require larger computational resources.Göckede et al. [43] improved an Eulerian footprint model use of satellite maps for explicit assignment of surface type.Göckede et al. [44] and Rebmann et al. [45] applied this new model at the EC sites and obtained satisfactory results.At present, footprint models are used to estimate the source areas contributing to the flux observations.In addition, they provide a tool for quality control of the flux measurements and provide guidance in designing experiments [46].Thus, the footprint models have considerable potential in microclimatology investigations, especially in studies which include nonhomogeneous surfaces.
The momentum roughness length ( 0 ), thermal roughness length ( 0ℎ ), and the water vapor roughness length ( 0 ) are crucial parameters for calculating momentum and heat fluxes in bulk transfer equations, which is one of the essential components in LSMs.It has been widely observed that  0 differs from  0ℎ and  0 [4,47].Unfortunately,  0 ,  0ℎ , and  0 , up to date, are still treated as constants in most LSMs.Inaccurate estimates of roughness length would enlarge the bias of simulated energy and water fluxes in LSMs.Many researchers have found that roughness length strongly depended on surface heterogeneity, vegetation height, and coverage [48,49].Therefore, the values of roughness length vary considerably in different geographical context or vegetation types [50][51][52][53][54]. Dryland ecosystems were sparsely Advances in Meteorology   [9,29,55,56].In Central Asia dryland ecosystems, leaf area index, albedo, and roughness length should be much more important, due to low fraction of vegetation and high heterogeneity of the land surface.Compared with recent findings at a Chinese desert shrub ecosystem [21], which is quite similar to the two sites used in this study, that simulation of latent heat flux in CoLM was strongly dependent on the choice of different root water uptake functions [21].Land surface and vegetation parameters such as albedo, roughness length, and LAI also showed significant sensitivities in simulating the energy and water fluxes in Central Asia dryland ecosystems.These researches implied that CoLM can be very sensitive to both belowground ecological process (such as root water uptake) and aboveground soil and vegetation properties (LAI,  0 , and ) in Central Asia regions, which is different from the sensitivities of LSMs in humid regions [8,26,57].
Vegetation and soil parameters, such as leaf area index, roughness length, and albedo, are closely related to vegetation coverage on the land surface.With recent findings at a Chinese desert shrub ecosystem [21], root water uptake process is more relevant to root distribution and root water uptake efficiency.However, many of these crucial parameters are difficult to observe in large areas.Biases in land surface forcing data and parameterizations in representing soil moisture, soil temperature, and other variables in numerical weather forecast and climate models will enlarge the model errors of water and energy fluxes.Reinitialization of land surface states would mollify this problem if the land surface variables and parameters were reliable and available in large areas and with high spatial resolution.Fortunately, the improvement of remote sensing technologies made these operable.Remote sensing technologies also could obtain LAI, canopy and ground albedo, vegetation height, and other parameters in large regions.To make best use of the satellite-based and ground-based observational data in land surface modelling for investigating global climate change issue in regional area, many land data assimilation systems have been developed [58,59].Such treatment is to generate optimal fields for parameterizing and forcing LSMs [59].Broad use of land data assimilation systems' results is valuable for predicting climate change, weather, and biological and agricultural productivity and for performing a wide array of studies in the broader biogeosciences.

Conclusions
In this study, the CoLM, for the first time, has been evaluated at two Central Asia desert ecosystems.Additionally, sensitivities of the model to LAI,  0 , and  were conducted.Evaluation of the CoLM and their sensitivities against the observed energy fluxes using eddy covariance system and the sensitivity tests resulted in the following conclusions: (1) The reference simulations (S0) significantly underestimated the latent heat flux and overestimated the sensible heat flux at two sites, especially at KZ-Ara.However, refined estimate of roughness length (S1) significantly improved the performance in simulating turbulent heat fluxes.The latent heat flux was increased but sensible heat flux was decreased, which were in better agreement with the observations from eddy covariance system.
(2) Sensitivity analysis regarding leaf area index, roughness length, and albedo showed that net radiation is very sensitive to albedo but latent and sensible heat fluxes and soil temperature are sensitively varying with the estimate of  0 at two EC sites over Central Asia.

Figure 1 :
Figure 1: Study area and location of the two eddy covariance (EC) sites in Kazakhstan.

Figure 2 :
Figure 2: Energy balance closure at the KZ-Ara site and the KZ-Bal site.The slope of the fitted line represents energy closure ratio and  2 is the coefficient of determination.RMSE (W m −2 ) is the root mean square error.The energy fluxes include sensible heat flux ( ℎ ), latent heat flux ( le ), net radiation ( net ), and ground heat flux ().

Figure 3 :
Figure 3: Comparison between the measured half-hourly net radiation ( net ), latent heat flux ( le ), and sensible heat flux ( ℎ ) and the values simulated by the reference model (S0) and the model with refined roughness length (S1) at the KZ-Ara and KZ-Bal sites.The solid red line represents the linear regression between the simulation and the observed data, and the dashed line represents a 1 : 1 relationship between the datasets.

Figure 4 :
Figure 4: Comparison between the measured net radiation ( net ), latent heat flux ( le ), and sensible heat flux ( ℎ ) and the values simulated by the reference model (S0) and the model with refined roughness length (S1) on a diurnal course at the KZ-Ara and KZ-Bal sites.The diurnal flux values were calculated as the mean values of all data at same measurement time in a day for the entire time period.

Figure 5 :
Figure 5: Panels (a) and (b) are the isograms changed with time and depth variation for the reference model (S0) and the model with refined roughness length (S1), and (c) shows the corresponding observed data at the depth of 20-80 cm below ground at the KZ-Ara site; panels (d) and (e) are the isograms for S0 and S1 and (f) shows the observations at the KZ-Bal site.

Figure 6 :
Figure 6: Taylor diagram plot of the simulated  net (a),  le (b),  ℎ (c), and  soil (d) computed from a reference model (S0) and six independent sensitivity tests (S2-S7, see Table3) from the CoLM against EC observations at the two EC sites.Standard deviation (STD, W m −2 ) is calculated as the simulated variables divided by the observed data."Obs" refers to observed data points.Root mean square error (RMSE, W m −2 ) is represented by green lines. is the correlation coefficient.The higher the  and the smaller the STD and RMSE, the better the agreement between model and data.When comparing two simulations with different parameter values, the longer the distance between the two simulation points, the greater the sensitivity to that parameter.

Table 1 :
Location, soil, vegetation, climate factors, and flux time series at KZ-Ara site and KZ-Bal site.

Table 2 :
The parameterizations of the reference model (S0) and the model with refined roughness length (S1) in CLM at KZ-Ara site and KZ-Bal site.RWUF is an optimal root water uptake function for dryland ecosystem(Jing et al. 2014 [21]), and   ,   , and  are empirical constants in RWUF.

Table 3 :
A baseline of reference exprement (S0) and six different configurations of CLM at each site (S2-S7), as used in this study for sensitivity testing.Two sets of leaf area index (LAI), average roughness length ( 0 ), and average albedo () were used in CLM."" refers to the model default parameter setting.

Table 4 :
Model performance for simulating  net ,  le , and  ℎ , indicated by coefficient of determination ( 2 ), slope (  ), intercept ( 0 ), and root mean square error (RMSE, W m −2 ) of linear regressions between model and observed data at the KZ-Ara site and the KZ-Bal site.

Table 3
, with vegetation fraction of 10%-40% at Central Asia desert ecosystems.Bare soil has larger proportion in the Central Asia desert ecosystems, and its characteristic is entirely different from the high-vegetated land surface, particularly for the roughness length.The values of roughness length were replaced by empirical formula based on experimental observations in this study, and the results found that such treatment significantly improved the performance of the model.The simulated turbulent heat fluxes with default roughness length showed very large variability during a day, especially at daytime.The most possible reason is that the roughness length was extremely overestimated and hence aerodynamic resistance was overestimated in the default version of CoLM.As a consequence, the turbulent exchange was strong and of high fluctuations.However, observed diurnal dynamics of energy fluxes are quite harmonious.Using an empirical approach to estimate roughness length into CoLM significantly improved the performance in simulating turbulent heat fluxes, which indicated that accurate parameterization of roughness length is of crucial importance for predicting energy and water fluxes in LSM when applied in Central Asia desert ecosystems.Many researches have indicated that LSMs are very sensitive to model physics, land characteristics (including leaf area index, roughness length, albedo, and soil properties), and forcing ) from the CoLM against EC observations at the two EC sites.Standard deviation (STD, W m −2 ) is calculated as the simulated variables divided by the observed data."Obs" refers to observed data points.Root mean square error (RMSE, W m −2 ) is represented by green lines. is the correlation coefficient.The higher the  and the smaller the STD and RMSE, the better the agreement between model and data.When comparing two simulations with different parameter values, the longer the distance between the two simulation points, the greater the sensitivity to that parameter.vegetated