Towards closing the vertical water balance in Canadian atmospheric models: Coupling of the land surface scheme class with the distributed hydrological model watflood

Abstract

Second generation land surface schemes are the subject of much development activity among atmospheric modellers. This work is aimed at, among other things, improving the representation of the soil water balance in order to simulate, more properly, exchanges with the atmosphere and to permit the use of model output to generate streamflow for model validation. The Canadian development program is centred on CLASS, the Canadian Land Surface Scheme, developed at Environment Canada. This paper focuses on the improvement of hydrology in CLASS. This was accomplished by designing a two‐way interface to WATFLOOD, a distributed hydrologic model developed at the University of Waterloo. The two models share many features, which facilitated the coupling procedure.

The interface retains the three‐layer vertical moisture budget representation in CLASS but adds three horizontal runoff possibilities. Runoff from the surface water follows Manning's equation for overland flow. Interflow is generated from the near‐surface soil layer using a parametrization of Richard's equation and base flow is produced by Darcian flow from the bottom of layer 3. An approximation of the internal topography of grid elements is used to supply horizontal gradients for the runoff components.

Tests are in progress in four Canadian study areas. Initial results are presented for the summer of 1993 for the Saugeen River in southwestern Ontario. The new scheme produces realistic hydrographs, whereas the old scheme did not. Bare ground evaporation is reduced by about 17% as a consequence of reduced water availability in layer 1. Evapotranspiration is not affected because the rooting depth extends into layer 3, in which soil moisture does not change appreciably with the new scheme. These results suggest that the new scheme improves the representation of streamflow in WATFLOOD/CLASS and of the soil moisture budget in CLASS. Work is in progress to validate this result over basins, such as the BOREAS study watersheds, where both runoff and evapotranspiration measurements are available.     

Parametrization schemes of incident radiation in the North Water polynya

Abstract

An evaluation of the Canadian Land Surface Scheme (CLASS) 3.1 snow cover simulations at four sites included in the Snow Model Intercomparison Project (SnowMIP) revealed that CLASS was able to provide realistic representations of snow cover accumulation, melt and physical properties over a range of snow cover climates. The modified snow aging parametrization in CLASS 3.1 provided improved simulations of snowpack density which resulted in a marked reduction in the root‐mean‐square (rms) error for daily snow depth, and slight improvements in snow surface temperature. CLASS 3.1 still exhibited a tendency to overestimate snow cover duration which is attributed to the way shallow snow ablation is treated. CLASS provided generally realistic simulations of daily and seasonal variation in snow albedo although cold snow albedo was underpredicted by 0.10 to 0.15 at a site with a deep (> 2 m) cold snowpack. CLASS also exhibited a tendency to overpredict late spring snow albedo which was reduced by the addition of a snow layer subroutine that kept track of snow albedo by precipitation event. CLASS had a noticeable cold bias averaging 3°–4°C at two mountain sites included in the comparison. The bias was closely linked to atmospheric stability and could exceed 10°C under conditions of strong radiative cooling and low wind speeds. The CLASS energy deficit under these conditions was determined to be ～20–40 W m^?2 and was mostly accounted for by introducing a windless exchange coefficient into the calculation of sensible heat fluxes following the approach used in a number of other physical snowpack models. CLASS provided realistic simulations of daily snowmelt runoff with the exception of the Weissfluhjoch site which was characterized by a deep cold snowpack. A preliminary assessment of snow water equivalent (SWE) rms error for the 23 models participating in SnowMIP showed that CLASS was one of the better single layer snow models included in the comparison. CLASS performance was comparable to the multi‐layer CROCUS snowpack model in the evaluations carried out in this study.     

Validation of the Canadian land surface scheme (class) for a subarctic open woodland

Abstract

This study reports on testing of the peatland version of the Canadian Land Surface Scheme (CLASS) for simulating the energy balance of subarctic open woodland terrain. Model results are compared against several years of measured data from a site near Churchill, Manitoba. In contrast to most forest environments, the floor of the open forest plays a large role in total ecosystem energy exchange. This behaviour presents a significant challenge for land surface models like CLASS and their simplified treatment of vegetation canopies.

Simulations of summer energy balance for seven years encompassing a wide range of meteorological conditions produced consistent results. Root mean square errors for sensible and latent heat fluxes fell between 11 and 28 W m^?2. CLASS consistently underestimated slightly the daily latent heat flux and overestimated the sensible heat flux, average mean bias errors being ‐7.6 and 9.1 W m^?2, respectively. The soil heat flux was less well represented. In general, CLASS was able to capture the diurnal and seasonal behaviour of the measured fluxes under a range of conditions with reasonable accuracy.

In a full year simulation, CLASS reproduced the annual variations in energy balance with some discrepancies associated with snow accumulation and ablation periods. The model performance was sensitive to both snow density and specification of the surface cover. Recommendations for improving the model for subarctic woodlands and terrain types with similar features are discussed.     

Local climates simulated by two generations of Canadian GCM land surface schemes

Abstract

The performance of two Canadian land surface schemes of widely differing complexity is compared and contrasted in a pair of year‐long simulations using the GCM developed at Atmospheric Environment Service, Canada. The old land surface model incorporates the force‐restore method for soil temperatures and the bucket approximation for soil moisture; the new model, CLASS (Canadian Land Surface Scheme) features three soil layers, an explicitly modelled snow layer, a thermally separate vegetation canopy, and physically‐based calculations of heat and moisture transfers between all of the land surface components and the atmosphere.

It was reported in previous papers that compared with observations, the old scheme tends to generate a climate which is characterized by anomalously high precipitation rates and cold temperatures over land. In this paper, by reference to field measurements and to the energy fluxes and temperatures generated by the two models at local scales, the hypotheses earlier postulated as to the underlying reasons for this are validated. The main factor contributing to the climate anomalies observed with the old scheme is found to be its generation of excessive evaporation rates; this is caused by the fact that the evaporation rate is never directly energy‐limited, the fact that the scaling of the evaporation rale with decreasing soil moisture content underestimates the effect of vegetation stomatal resistance, and the fact that the evaporation rate over bare soil depends not on the surface soil moisture, but on the moisture content of whole modelled soil column. The cold surface temperatures are additionally attributed to systematic errors incurred by the forward‐stepping temperature scheme, and to the fact that soils subjected to subzero temperature forcing in the winter are modelling as freezing completely. Finally, the inability of the old scheme to simulate partially frozen soils means that it proves unable to handle either shallow frost penetration at temperature latitudes, or the development of an active layer in permafrost.     

土壤热传导方程解析解和那曲地区土壤热扩散率研究

文中用Laplace变换推导了土壤热传导方程的解析解和包含热对流项的土壤热传导方程的解析解。用青藏高原 8个土壤湿度、温度廓线观测站 1998年 9月 4日到 10日实测资料基础上 ,根据谐波方法和Laplace变换方法得到了土壤热传导方程的解析解 ,计算了这些站的总体土壤热扩散率 ;用包含热对流项的土壤热传导方程的解析解计算了土壤热扩散率。结果表明 :对于一个深度从 0 .0 4～ 0 .2 0m的浅薄土壤层 ,总体土壤热扩散率的值为 0 .30×10 -6～ 0 .98× 10 -6m2 /s,土壤热扩散率的值为 0 .15× 10 -6～ 0 .72× 10 -6m2 /s。由谐波方法得到的总体土壤热扩散率比由Laplace变换的值稍大 ;总体土壤热扩散率总是比土壤热扩散率大     

Sensitivity of the hydrological cycle to the parametrization of soil hydrology in a GCM

 The sensitivity of the hydrological cycle to soil hydrology is investigated with the LMD GCM. The reference simulation includes the land-surface scheme SECHIBA, with a two-reservoir scheme for soil water storage and runoff at saturation. We studied a non-linear drainage parametrization, and a distributed surface runoff parametrization, accounting for the subgrid scale variability (SSV) of soil moisture capacity, through a distribution where the shape parameter was b. GCM results show that the drainage parametrization induces significant reductions in soil moisture and evaporation rate compared to the reference simulation. They are related to changes in moisture convergence in the tropics, and to a precipitation decrease in the extratropics. When drainage is implemented, the effect of the SSV parametrization (b=0.2) is also to reduce soil moisture and evaporation rates compared to the simulation with drainage only. These changes are much smaller than the former, but the sensitivity of the hydrological cycle to the SSV parametrization is shown to be larger in dry periods, and to be enhanced by an increase of the shape parameter b. The comparison of simulated total runoffs with observed data shows that the soil hydrological parametrizations does not reduce the GCM systematic errors in the annual water balance, but that they can improve the representation of the total runoff’s annual cycle.     

Investigating Propagation of Short-Period Ocean Waves Into the Periphery of Arctic Pack Ice Using High-Resolution Upward-Looking Sonar

Springtime fetch in the Cape Bathurst Polynya System may present opportunities for winds to generate waves capable of propagating into the thick pack ice formed over the winter. A waves-in-ice event at a study site located on the Canadian Shelf in the southern Beaufort Sea that occurred 22–23 May 2011 is presented and analyzed for wave energy attenuation and dissipation characteristics. The event was monitored near the ice edge and, therefore, presents information on attenuation of waves from the ice edge into the pack. Waves of T?=?5?s, λ?=?37.5?m were observed up to approximately 143?m and approximately 77?m away from the ice edge during two separate observation periods of ice edge wave propagation. We estimated reflection coefficients of 53% and 52% and wave attenuation coefficients of α?=?2.4?×?10^?2?m^?1 and α?=?5.4?×?10^?2?m^?1, respectively, for the two periods. Estimated attenuation rates are an order of magnitude greater than in comparable studies and are inconsistent with previous findings of a “rollover” effect in attenuation rates for short-period waves.     

On the mean and tidal currents in Hudson strait

Abstract

The vertical structures of the mean and tidal flows in Hudson Strait are described from moored current‐meter data collected during an 8‐week period in August to October of 1982. The residual flow in the strongly stratified waters off Quebec is directed along the Strait to the southeast, is highly baroclinic and is concentrated near shore (within an offshore length scale of approximately an internal Rossby radius). Maximum mean speeds of 0.3 m s^?1 were observed near‐surface (30 m). In the weakly stratified waters on the northern side of the Strait along Baffin Island the mean flow is northwestward. The maximum speeds are 0.1 m s^?1 near‐surface (30 m) and the current amplitudes decrease to 0.05 m s^?1 at 100 m. The mean southeastward transport is estimated to be 0.93 ±0.23 × 10⁶ m³ s^?1 with a northwestward transport of 0.82 ± 0.24 × 10⁶ m³ s^?1. Over most of the Strait the across‐channel residual currents are directed towards the Quebec shore with velocities ranging from 0.02 to 0.1 ms^?1. Current variability is dominated by the tides, the M₂ being the major tidal constituent. In the vicinity of the mooring the M₂ tide is primarily barotropic, progressive in nature, and has along‐channel current amplitudes varying across the Strait from 0.20 to 0.45 m s^?1. Observed differences in tidal sea‐level elevations across the Strait can be accounted for by the cross‐channel variations characteristic of Kelvin waves.     

Annual and seasonal climate and climatic changes in the Canadian prairies simulated by the CCC GCM

Abstract

As part of a study on the effects of climatic variability and change on the sustainability of agriculture in Alberto, the modelling performance of the second‐generation Canadian Climate Centre GCM (general circulation model) is examined. For the region in general, the simulation of 1 × CO₂ mean temperature is generally better than that for mean precipitation, and summer is the season best modelled for each variable. At the scale of individual grid squares, DJF (December, January, February) (temperature) and JJA (June, July, August) (precipitation) are the seasons best modelled. The GCM‐simulated increases in mean annual temperature resulting from a doubling of CO₂ are of the order of 5 to 6°C in the Prairie region, with much of this increase resulting from substantial warming in the winter and spring. Increases in mean annual precipitation are of the order of 50 to 150 mm (changes of +5 to +15%), with the greatest changes again occurring in winter and spring. As far as the limited GCM run durations allow, temperature and precipitation variance generally show no significant changes from a 1 × CO₂ to a 2 × CO₂ climate. Increased precipitation in winter and spring does not result in greater snow accumulations owing to the magnitude of warming; and significant decreases in soil moisture content occur in summer and fall. The resulting effects on the growing season and moisture regime have the potential to affect agricultural practices in the area.     

Regional‐scale surface hydrologic simulations from global climate models: A case study

Abstract

Current understanding of the regional nature of global changes in the climato‐logical regime of the earth is limited. General circulation climate models (GCMs) cannot provide consistent and detailed information on the regional patterns of precipitation, soil moisture and runoff that are required by water resource planners. A case study is presented that couples a limited area model, with high spatial resolution and realistic land‐surface parametrization, to a global climate model. Results for July are presented for the continent of Australia.     

On improving cold region hydrological processes in the Canadian Land Surface Scheme

Regional and global climate model simulated streamflows for high-latitude regions show systematic biases, particularly in the timing and magnitude of spring peak flows. Though these biases could be related to the snow water equivalent and spring temperature biases in models, a good part of these biases is due to the unaccounted effects of non-uniform infiltration capacity of the frozen ground and other related processes. In this paper, the treatment of frozen water in the Canadian Land Surface Scheme (CLASS), which is used in the Canadian regional and global climate models, is modified to include fractional permeable area, supercooled liquid water and a new formulation for hydraulic conductivity. The impact of these modifications on the regional hydrology, particularly streamflow, is assessed by comparing three simulations performed with the original and two modified versions of CLASS, driven by atmospheric forcing data from the European Centre for Medium-Range Weather Forecast (ECMWF) reanalysis (ERA-Interim) for the 1990–2001 period over a northeast Canadian domain. The two modified versions of CLASS differ in the soil hydraulic conductivity and matric potential formulations, with one version being based on formulations from a previous study and the other one is newly proposed. Results suggest statistically significant decreases in infiltration and therefore soil moisture during the snowmelt season for the simulation with the new hydraulic conductivity and matric potential formulations and fractional permeable area concept compared to the original version of CLASS, which is also reflected in the increased spring surface runoff and streamflows in this simulation with modified CLASS over most of the study domain. The simulated spring peaks and their timing in this simulation are also in better agreement to those observed. This study thus demonstrates the importance of treatment of frozen water for realistic simulation of streamflows.     

Preliminary results of modeling studies on washout of stable and radioactive analogs of <Superscript>90</Superscript>Sr and <Superscript>137</Superscript>Cs from snow-covered areas

Experimental field and laboratory studies on washout of radionuclides from the snow cover during snow melting were carried out in the winter of 2005/06. In the field studies, a specially equipped runoff site was used. In the laboratory conditions, the experiments were conducted using prepared soil monoliths. In the winter of 2006, 25 g/m² of water-free cesium chloride (CsCl) and 25 g/m³ of strontium chloride (SrCl₂) were put onto the snow cover surface of the runoff site. The snow surface of the soil monolith was coated with a ¹³⁷Cs-bearing solution, then with SrCl₂. Under experimental conditions, practically no surface runoff from the runoff site was recorded. The experiments with the soil monoliths demonstrated that the coefficient of the liquid washout of ¹³⁷Cs normalized to the runoff layer was within 0.9 × 10^?6–1.2 × 10^?4 mm^?1, and that of ⁹⁰Sr normalized to the runoff layer was within 2 × 10^?–1.6 × 10^?4 mm^?1.     

Precambrian Shield Wetlands: Hydrologic Control of the Sources and Export of Dissolved Organic Matter

Most Precambrian Shield forested catchments have some wetland component. Even small riparian wetlands are important modifiers of stream chemistry. Dissolved organic matter (DOM) is one of the most important products exported by wetlands in streams. Stratigraphic control of hydraulic conductivity generally leads to decreasing conductivity with depth. Thus important flowpaths occur in the uppermost organic rich layers and are reflected in chemical profiles of dissolved organic carbon (DOC). Accumulation of DOC in peat porewaters is the net effect of production, consumption and transport. DOC profiles vary with degree of interaction with the surrounding upland catchment and distance from the edge of the wetland as well as internal processes within the wetland. In wetlands, DOM production is offset by flushing resulting in decreasing DOC concentrations with increasing flows. Despite old carbon (2,000 to 3,000 years) at relatively shallow depths, 14C activity in DOC exported from wetlands is mostly modern (recent carbon), consistent with shallow flowpaths and export of DOM from shallow organic rich horizons. In contrast, the source area for DOM in upland catchments with developed B horizon soils increases with antecedent soil moisture conditions resulting in increasing DOC concentrations with higher stream flows. Activity of 14C in stream DOC from upland catchments span a range from low activities (older carbon) similar to B horizon soil water during dry moisture conditions to values slightly less than modern (more recent carbon) during high moisture conditions. The more modern carbon activities reflect the increased contribution of the organic rich litter and A horizon soil layers in the area immediately bordering the stream under wet antecedent moisture conditions. Reduced hydrologic export or loss of wetlands under drier climatic conditions may result in in larger fluctuations in stream DOC concentrations and reduced DOM loads to lakes.     

Hydrologic Response of a Wetland to Changing Moisture Conditions: Modeling Effects of Soil Heterogeneity

Prediction of the effects of external influences such as climate change on wetland systems requires the prediction of hydrologic effects. Because wetland soils are typically heterogeneous, it is particularly important to understand the extent and connectedness of hydraulically conductive soil units, since water flow may be concentrated in such units while bypassing others of lower conductivity. However, subsurface hydrologic models typically do not represent heterogeneity adequately, being limited by sparse parameterization of soil properties. Conventional techniques for mapping units of soil within wetlands are highly laborious, requiring soil coring and laboratory testing. As an alternative, we developed a portable piezocone driver and highly sensitive piezocone designed to map wetland soil units with centimeter-scale resolution in the vertical and meter-scale resolution in the horizontal dimension. This system successfully delineated several different layers of peat, sand, and limnetic sediments, and their degree of interconnectedness in an eight-meter-thick peat deposit. Monitoring of wetland response to precipitation, changes in stream stage, and overbank flooding was then used in conjunction with the piezocone data and a two-dimensional flow model to constrain the hydraulic properties of the soil units. Thus parameterized, a standard subsurface flow model was able to realistically simulate a variety of hydrologic processes relevant to climate change, including wetland-stream water exchange, the movement of wetland porewaters to the root zone of plants, and wetland desaturation under dry conditions.     

Simulating the Carbon Cycling of Northern Peatlands Using a Land Surface Scheme Coupled to a Wetland Carbon Model (CLASS3W-MWM)

Northern peatlands store approximately one-third of the terrestrial soil carbon (C), although they cover only 3% of the global land mass. Northern peatlands can be subdivided into bogs and fens based on their hydrology and biogeochemistry. Peatland hydrology and biogeochemistry are tightly coupled to climate and, therefore, may be very sensitive to climate variability and change. To address the fate of the large peatland soil C storage under a future changed climate, a peatland C model, the McGill Wetland Model (MWM), was coupled to a land surface climate model (the wetland version of the Canadian Land Surface Scheme, CLASS3W), referred as CLASS3W-MWM. We evaluated the CLASS3W-MWM for a bog (Mer Bleue, located at 45.41°N, 75.48°W, in eastern Canada) and a poor fen (Degerö Stormyr, located at 64°11′N, 19°33′E, in northern Sweden).

CLASS3W-MWM captured the magnitude and direction of the present day C cycling very well for both bogs and fens. Moreover, the seasonal and interannual variability were reproduced reasonably well. Root mean square errors (RMSE) were <0.65 and the degree of agreements (d*) were >0.8 for the components of net ecosystem production (NEP) for both the Mer Bleue bog and the Degerö Stormyr fen. The performance of the coupled model for both bog and fen is similar to that of the stand-alone MWM driven by observed weather rather than simulated surface and soil climate. This modelling study suggests that northern peatlands are hydrologically and thermally conservative ecosystems. It was also shown that C cycling for bogs and fens was more sensitive to changes in air temperature than precipitation. Changes in temperature within the Intergovernmental Panel on Climate Change (IPCC) projected range switch the peatlands from a present-day C sink to a source, but projected changes in precipitation still maintain the peatlands as a C sink, although to a somewhat lesser degree. Increase in atmospheric CO₂ concentration enhances C sequestration for both bogs and fens. Our sensitivity analysis suggests that northern peatlands respond to changes in temperature, precipitation and doubled CO₂ concentration in a highly non-linear way. The sensitivity of C cycling in northern peatlands with respect to changes in air temperature, precipitation and the concentration of atmospheric CO₂ together is not a simple addition or subtraction of the sensitivity of the individual changes. Therefore, the sensitivity of a combination of changes in temperature, precipitation and doubled CO₂ concentration is very different from the sensitivity of peatlands to each environmental variable on their own. Our sensitivity analysis suggests that fens have a narrower tolerance to climate changes than bogs.

RÉSUMÉ [Traduit par la rédaction] Les tourbières du Nord renferment approximativement le tiers du carbone se trouvant dans le sol terrestre, même si elles ne couvrent que 3% des terres du globe. On peut subdiviser les tourbières du Nord en tourbières hautes et en tourbières basses selon leur hydrologie et leur biogéochimie. L'hydrologie et la biogéochimie des tourbières sont intimement liées au climat et peuvent donc être très sensibles à la variabilité et au changement climatique. Pour étudier comment évoluera le stockage du carbone dans les grands terrains tourbeux sous un climat futur modifié, nous avons couplé un modèle de carbone de tourbière, le McGill Wetland Model (MWM), à un modèle climatique de surface terrestre (la version terres humides du CLASS3W canadien), c'est-à-dire le CLASS3W–MWM. Nous avons évalué le CLASS3W–MWM pour une tourbière haute (Mer Bleue, situé à 45,41°N, 75,48°O, dans l'est du Canada) et pour une tourbière basse ombrotrophe (Degerö Stormyr, situé à 64°11′N, 19°33′E, dans le nord de la Suède).

Le CLASS3W–MWM a très bien capturé la grandeur et la direction du recyclage actuel du carbone, tant pour les tourbières hautes que pour les tourbières basses. De plus, la variabilité saisonnière et interannuelle a été raisonnablement bien reproduire. Les écarts-types étaient <0,65 et les degrés de concordance (d*) étaient >0,8 pour les composantes de la production nette de l’écosystème tant pour la tourbière haute Mer Bleue que pour la tourbière basse Degerö Stormyr. La performance du modèle couplé pour la tourbière haute et la tourbière basse est semblable à celle du MWM autonome piloté par des conditions observées plutôt que par un climat simulé de la surface et du sol. Cette étude par modèle suggère que les tourbières du Nord sont des écosystèmes hydrologiquement et thermiquement conservatifs. Il a aussi été démontré que le recyclage du carbone pour les tourbières hautes et basses était plus sensible aux changements dans la température de l'air que dans les précipitations. Des changements de température de l'ordre de ceux projetés par le Groupe d'experts intergouvernemental sur l’évolution du climat (GIEC) font que les actuels puits de carbone que constituent les tourbières se transforment en sources, mais les changements projetés dans les précipitations maintiennent encore les tourbières comme des puits de carbone, quoique dans une moindre mesure. L'accroissement de la concentration du CO₂ atmosphérique améliore la séquestration du carbone à la fois pour les tourbières hautes et les tourbières basses. Notre analyse de sensibilité suggère que les tourbières du Nord réagissent aux changements dans la température et les précipitations et à une concentration doublée de CO₂ d'une façon fort peu linéaire. La sensibilité du recyclage du carbone dans les tourbières du Nord par rapport aux changements dans la température de l'air, les précipitations et la concentration du CO₂ atmosphérique ensemble n'est pas une simple addition ou soustraction de la sensibilité aux changements individuels. Par conséquent, la sensibilité à une combinaison de changements dans la température et les précipitations et à une concentration doublée de CO₂ est très différente de la sensibilité des tourbières à chaque variable environnementale prise seule. Notre analyse de sensibilité suggère que les tourbières basses ont une plus faible tolérance aux changements climatiques que les tourbières hautes.     

Northern Canadian Wetlands: Net Ecosystem CO2 Exchange and Climatic Change

Northern Canadian peatlands represent a long term sink for atmospheric carbon dioxide (CO2), however there is concern they may become a net source of CO2 due to climatic change. Climatic change is expected to result in significant changes in regional hydrology in boreal and subarctic regions of Canada. A hydrologic model predicted a summer water table drop of 0.14 m in northern Canadian fens given an increase in summer temperature and rainfall of 3°C and 1 mm d-1, respectively. Moreover, surface peat temperature increased by 2.3°C. Net ecosystem exchange of CO2 was modelled using these modelled hydrologic and thermal changes with respiration:peat temperature and water table:net ecosystem production relationships developed from measurements at wetlands in northern Sweden and near Churchill, Manitoba. Model results indicate that the net atmospheric CO2 sink function of fens may be enhanced under future 2 × CO2 scenarios, while bogs may become a net source of atmospheric CO2. If the net ecosystem productivity response to the new hydrologic conditions was ignored then the model predicts a decrease in summer carbon storage for all peatland types.     

Comparative analysis of different underlying surfaces using a high-resolution assimilation dataset in semi-arid areas in China

To study the land surface and atmospheric meteorological characteristics of non-uniform underlying surfaces in the semi-arid area of Northeast China, we use a “High-Resolution Assimilation Dataset of the water-energy cycle in China (HRADC)”. The grid points of three different underlying surfaces were selected, and their meteorological elements were averaged for each type (i.e., mixed forest, grassland, and cropland). For 2009, we compared and analyzed the different components of leaf area index (LAI), soil temperature and moisture, surface albedo, precipitation, and surface energy for various underlying surfaces in Northeast China. The results indicated that the LAI of mixed forest and cropland during the summer is greater than 5 m² m^?2 and below 2.5 m² m^?2 for grassland; in the winter and spring seasons, the Green Vegetation Fraction (GVF) is below 30%. The soil temperature and moisture both vary greatly. Throughout the year, the mixed forest is dominated by latent heat evaporation; in grasslands and croplands, the sensible heat flux and the latent heat flux are approximately equal, and the GVF contributed more to latent heat flux than sensible heat flux in the summer. This study compares meteorological characteristics between three different underlying surfaces of the semi-arid area of Northeast China and makes up for the insufficiency of purely using observations for the study. This research is important for understanding the water-energy cycle and transport in the semi-arid area.     

Effect of surface layer thickness on simultaneous transport of heat and water in a bare soil and its implications for land surface schemes

Abstract

A physically‐based multi‐layer numerical model is developed to determine the coupled transport of heat and water in the soil and in the soil‐atmosphere boundary layer. Using inputs of standard weather data and initial soil conditions the model is capable of predicting the surface energy balance components as well as water content and temperature profiles in the soil. It is used to predict these variables for a bare silt loam soil under two tillage treatments, viz. culti‐packed and left loose after disc‐harrowing, and the predicted results are compared with measurements. Very good agreement between the model predictions and measured evaporation and heat fluxes and soil water and temperatures for a ten‐day period shows that the model is capable of simulating the coupled transport of soil heat and soil water and their transfer across the soil surface‐atmosphere interface adequately.

Model predictions were compared with those of CLASS (Canadian Land Surface Scheme). It is shown that CLASS, version 2.6, provides good estimates of evaporation and hence the latent heat flux density, Q_E, under wetter soil conditions, but overestimates Q_E at moderately wet soil conditions and underestimates it under dry soil conditions. Under dry to moderately wet soil conditions the calculation of evaporation from bare soil is very sensitive to the thickness of the top layer particularly as the thickness approaches 10 cm.     

青藏高原湿地土壤冻结、融化期间的陆面过程特征

利用青藏高原中部玉树隆宝湿地2015年7月-2016年7月的观测资料,分析了土壤冻结、融化前后土壤温、湿度和地表能量收支特征,结果表明：冻土持续时期为12月至次年4月,深层土壤的冻结较浅层土壤滞后,融化过程快于冻结过程,5-40 cm土壤全部冻结历时51 d,全部融化历时19 d。土壤体积含水量年变化幅度达0.6 m³/m³。冻结过程5-40 cm土壤体积含水量下降,融化过程5-10 cm土壤体积含水量升高。土壤冻结之后,感热通量白天的值升高,潜热通量白天的值降低,净辐射和土壤热通量均降低,土壤热通量日变化幅度增大。土壤融化之后,潜热通量、净辐射和土壤热通量白天的值升高。地表反照率、鲍恩比、土壤热导率和土壤热扩散率冻结后增大融化后减小,土壤热容量冻结后减小融化后增大。     

Soil moisture influence on summertime surface air temperature over East Asia

Soil moisture influence on surface air temperature in summer is statistically quantified across East Asia using the Global Land Data Assimilation System soil moisture and observational temperature. The analysis uses a soil moisture feedback parameter computed based on lagged covariance ratios. It is found that significant negative soil moisture feedbacks on temperature mainly appear over the transition zones between dry and wet climates of northern China and Mongolia. Over these areas, the feedbacks account for typically 5–20% of the total temperature variance, with the feedback parameter of ?0.2°C to ?0.5°C (standardized soil moisture)^?1. Meanwhile, positive feedbacks may exist over some areas of Northeast Asia but are much less significant. These findings emphasize the importance of soil moisture-temperature feedbacks in influencing summer climate variability and have implications for seasonal temperature forecasting.