Calibration of an integrated land surface process and radiobrightness (LSP/R) model during summertime

In this study, a soil vegetation and atmosphere transfer (SVAT) model was linked with a microwave emission model to simulate microwave signatures for different terrain during summertime, when the energy and moisture fluxes at the land surface are strong. The integrated model, land surface process/radiobrightness (LSP/R), was forced with weather and initial conditions observed during a field experiment. It simulated the fluxes and brightness temperatures for bare soil and brome grass in the Northern Great Plains. The model estimates of soil temperature and moisture profiles and terrain brightness temperatures were compared with the observed values. Overall, the LSP model provides realistic estimates of soil moisture and temperature profiles to be used with a microwave model. The maximum mean differences and standard deviations between the modeled and the observed temperatures (canopy and soil) were 2.6 K and 6.8 K, respectively; those for the volumetric soil moisture were 0.9% and 1.5%, respectively. Brightness temperatures at 19 GHz matched well with the observations for bare soil, when a rough surface model was incorporated indicating reduced dielectric sensitivity to soil moisture by surface roughness. The brightness temperatures of the brome grass matched well with the observations indicating that a simple emission model was sufficient to simulate accurate brightness temperatures for grass typical of that region and surface roughness was not a significant issue for grass-covered soil at 19 GHz. Such integrated SVAT-microwave models allow for direct assimilation of microwave observations and can also be used to understand sensitivity of microwave signatures to changes in weather forcings and soil conditions for different terrain types.     

基于双曲线模型的冻土中锚杆荷载传递特性研究

为研究不同冻结温度及含水率对冻结粉土中锚杆抗拔性能的影响,根据冻土与混凝土接触面的力学特性和变形规律,采用双曲线模型描述锚杆-冻土接触面的剪切特性。基于荷载传递法,建立考虑温度和含水率影响的锚杆荷载传递方程,采用有限差分法进行求解,得到锚固段的剪应力、轴力及承载力计算式;通过ABAQUS数值模拟验证荷载传递方程的合理性,结合算例分析表明：冻土中锚杆抗拔承载力大于常温土,且温度越低承载力越高;在相同荷载作用下,冻土温度越低,轴力沿深度衰减越快,剪应力分布越不均匀;相同冻结温度下,锚杆承载力随含水率的增大呈先增大后减小的趋势;增加锚固体直径能有效地提高锚杆的极限承载力。     

Research on the method for retrieving soil moisture using thermal inertia model

1 Introduction Thermal inertia is a bulk property that shows the re- sistance of a material to an input or output of heat. This plays a very important role in certain geological and hydrological studies, and climate modeling. In the 1970s, a simple thermal inertia model was proposed by Watson et al.[1―3]. Pratt (1979)[4] improved the thermal inertia model based on application tests where more factors were considered such as solar ra- diance, thermal conductivity effect, average humidity of g…     

Prediction of hourly actual evapotranspiration using neural networks,genetic programming,and statistical models

The complexity of the evapotranspiration process and its variability in time and space have imposed some limitations on previously developed evapotranspiration models. In this study, two data‐driven models: genetic programming (GP) and artificial neural networks (ANNs), and statistical regression models were developed and compared for estimating the hourly eddy covariance (EC)‐measured actual evapotranspiration (AET) using meteorological variables. The utility of the investigated data‐driven models was also compared with that of HYDRUS‐1D model, which makes use of conventional Penman–Monteith (PM) model for the prediction of AET. The latent heat (LE), which is measured using the EC method, is modelled as a function of five climatic variables: net radiation, ground temperature, air temperature, relative humidity, and wind speed in a reconstructed landscape located in Northern Alberta, Canada. Several ANN models were evaluated using two training algorithms of Levenberg–Marquardt and Bayesian regularization. The GP technique was used to generate mathematical equations correlating AET to the five climatic variables. Furthermore, the climatic variables, as well as their two‐factor interactions, were statistically analysed to obtain a regression equation and to indicate the climatic factors having significant effect on the evapotranspiration process. HYDRUS‐1D model as an available physically based model was examined for estimating AET using climatic variables, leaf area index (LAI), and soil moisture information. The results indicated that all three proposed data‐driven models were able to approximate the AET reasonably well; however, GP and regression models had better generalization ability than the ANN model. The results of HYDRUS‐1D model exhibited that a physically based model, such as HYDRUS‐1D, might be comparable or even inferior to the data‐driven models in terms of the overall prediction accuracy. Based on the developed GP and regression models, net radiation and ground temperature had larger contribution to the AET process than other variables. Copyright © 2010 John Wiley & Sons, Ltd.     

Assessing temperature‐based PET equations under a changing climate in temperate,deciduous forests

Despite the availability of numerous approaches to estimate potential evapotranspiration (PET), temperature‐based PET equations such as the Hamon equation and Thornthwaite equation are still used to predict changes in hydrology in a changing climate as temperature is one of the primary reported outputs from general circulation models. To isolate the actual dependence of PET on temperature, we analysed meteorological and energy balance measurements from five AmeriFlux deciduous forest sites in the eastern United States during periods with minimal soil moisture control on transpiration. For all five sites, when PET measurements with similar net radiation are grouped, temperature does not correlate to PET within each group. Conversely when PET measurements with similar temperature are grouped, net radiation strongly correlates to PET within each group. In terms of assessing standard PET models, when dormant and growing season PET are separated, we found that the Priestley–Taylor equation (a model primarily dependent on net radiation) consistently explained more of the variation in PET than temperature‐based methods such as the Hamon equation (median R² of 0·88 vs 0·66). We illustrate that the moderate ability of temperature‐based equations to predict PET arises from the correlation between temperature and net radiation when the meteorological observations are averaged over at least several days. However, we suggest that because temperature is not the fundamental driver of PET and because the relationship between temperature and net radiation underlying temperature‐based equations will shift with climate change, temperature‐based equations in their current state will likely exaggerate PET in a changing climate. Copyright © 2010 John Wiley & Sons, Ltd.     

A new saturated/unsaturated model for stormwater infiltration systems

Infiltration systems are widely used as an effective urban stormwater control measure. Most design methods and models roughly approximate the complex physical flow processes in these systems using empirical equations and fixed infiltration rates to calculate emptying times from full. Sophisticated variably saturated flow models are available, but rarely applied owing to their complexity. This paper describes the development and testing of an integrated one‐dimensional model of flow through the porous storage of a typical infiltration system and surrounding soils. The model accounts for the depth in the storage, surrounding soil moisture conditions and the interaction between the storage and surrounding soil. It is a front‐tracking model that innovatively combines a soil‐moisture‐based solution of Richard's equation for unsaturated flow with piston flow through a saturated zone as well as a reservoir equation for flow through a porous storage. This allows the use of a simple non‐iterative numerical solution that can handle ponded infiltration into dry soils. The model is more rigorous than approximate stormwater infiltration system models and could therefore be valuable in everyday practice. A range of test cases commonly used to test soil water flow models for infiltration in unsaturated conditions, drainage from saturation and infiltration under ponded conditions were used to test the model along with an experiment with variable depth in a porous storage over saturated conditions. Results show that the model produces a good fit to the observed data, analytical solutions and Hydrus. Copyright © 2008 John Wiley & Sons, Ltd.     

An improved algorithm to estimate the surface soil heat flux over a heterogeneous surface: A case study in the Heihe River Basin

Surface soil heat flux(G0) is an indispensable component of the surface energy balance and plays an important role in the estimation of surface evapotranspiration(ET). This study calculated G0 in the Heihe River Basin based on the thermal diffusion equation, using the observed soil temperature and moisture profiles, with the aim to analyze the spatial-temporal variations of G0 over the heterogeneous area(with alpine grassland, farmland, and forest). The soil ice content was estimated by the difference in liquid soil water content before and after the melting of the frozen soil and its impact on the calculation of G0 was further analyzed. The results show that:(1) the diurnal variation of G0 is obvious under different underlying surfaces in the Heihe River Basin, and the time when the daily maximum value of G0 occurs is a few minutes to several hours earlier than that of the net radiation flux, which is related to the soil texture, soil moisture, soil thermal properties, and the vegetation coverage;(2) the net radiation flux varies with season and reaches the maximum in summer and the minimum in winter, whereas G0 reaches the maximum in spring rather than in summer, because more vegetation in summer hinders energy transfer into the soil;(3) the proportions of G0 to the net radiation flux are different with seasons and surface types, and the mean values in January are 25.6% at the Arou site, 22.9% at the Yingke site and 4.3% at the Guantan site, whereas the values in July are 2.3%, 1.6% and 0.3%, respectively; and(4) G0 increases when the soil ice content is included in thermal diffusion equation, which improves the surface energy balance closure by 4.3%.     

An improved methodology to measure evaporation from bare soil based on comparison of surface temperature with a dry soil surface

Accurate estimation of the resistances to water vapor movement is a major difficulty in evaluating evaporation from soil. By including the temperature of a dry soil surface (the temperature of the surface of a dry soil column buried in the field), a method for estimating evaporation from soil is proposed. The necessary input variables for the suggested method are temperature, net radiation, and soil heat flux. There are three advantages of the proposed method over the conventional methods. First, soil surface resistance and aerodynamic resistance are not required. Second, the variables included are fewer. Third, measurement and analysis of the parameters involved are relatively easy. Sensitivity analysis shows that the suggested method is sensitive to temperatures. Test experiments were conducted in a sandy field, where a weighing lysimeter was installed. Evaporation from soil, together with the variables specified above, were measured. For temperatures measured by thermocouples, experimental results showed that the mean absolute error (MAE) for the daily evaporation over 22 days was 0.17 mm day⁻¹. The regression between calculated and measured evaporation was highly significant (r²=0.89). Moreover, the intercept and slope of the regression equation were not significantly different from zero and unity, respectively, at the 0.05 probability level. Furthermore, by using the temperatures measured by infrared thermometers, the MAE between measured evaporation and estimated evaporation was 0.15 mm day⁻¹. The regression between them was highly significant (r²=0.94). In addition, the intercept and slope of the regression equation were not significantly different from zero and unity, respectively, at the 0.05 probability level. These results show that evaporation calculated using the proposed method is in good agreement with lysimeter measured values. By comparing with the temperature difference method, it was shown that the suggested method estimated soil evaporation more accurately than the temperature difference method. Therefore, it is concluded that the proposed method is not only a simple way for application, but also an accurate way to estimate soil evaporation.     

A study of the infiltration characteristics of undisturbed soil under simulated rainfall

Infiltration is the single most important parameter in deriving the net quick response rainfall which contributes to stream flood discharges. Rainfall simulation is used to study the infiltration characteristics in a typical catchment, the Six Mile Water in N. Ireland. The design of the simulator was such that it could be easily moved from one test area to another within the catchment to examine the effect of soil and slope variation. The simulator was first calibrated in controlled laboratory conditions and later the calibration was checked in the field. The simulator was mounted over an undisturbed plot of 37 m² and the surface runoff from the area measured by means of a collecting channel located along a lower edge of the plot. Soil moisture variations were monitored using a soil moisture neutron probe. Soil classification tests and gravimetric moisture contents were carried out on each plot. The field tests were carried out with variations in rainfall intensity, initial conditions, changing seasons, and for different plots within the catchment area. The results obtained are unique in that they present data obtained under field conditions for undisturbed soil within a natural catchment. The infiltration behaviour was found to depend upon rainfall intensity, initial conditions of the plot under consideration, seasonal temperature, and a slope of the plot. The data showed that while a classical Horton type equation for infiltration was suitable for the later stages of each test result when significant surface runoff was taking place, the model failed to represent early response adequately due to storage effects being omitted in the equation. A modified form of Horton equation is proposed, which models more accurately the infiltration characteristics of the full period of each test run.     

On the validation of a coupled water and energy balance model at small catchment scales

Catchment runoff is the most widely used catchment scale measurement in modelling studies, and we have a reasonable degree of confidence in its accuracy. The advent of satellites gives access to a new suite of measurements taken over a defined spatial range. These measurements, principally reflected or emitted radiation, provide hydrologists with new possibilities for quantifying the state of a catchment. Surface temperatures can be readily measured by a satellite on a scale comparable to the size of a small catchment.

In this paper we show that satellite sensed temperatures can provide an important measure of catchment status, which can complement runoff measurements in water balance studies. A one-dimensional model, which couples the land surface energy balance with the soil and surface water balance is tested by comparison with runoff and with remotely sensed surface temperature measurements. Simulations have been run over four years for two small catchments which have a fairly homogeneous vegetation, one being forest and its neighbour pasture. Satellite “surface” temperatures have been interpreted in terms of the energy balance, and used as a test of modelling accuracy. An “effective” surface temperature is calculated as a weighted mean of temperatures of the separate soil and leaf surfaces. This modelled “effective” temperature correlates well with Landsat TM surface temperatures.

When pasture replaces forest, the model predicts a reduction in evapotranspiration of around 30%, a three-fold increase in runoff, and an increase in mean soil moisture status. The change to pasture also results in a rise in mean effective surface temperature of about 4°C, and an increase in summer diurnal temperature range from 10 to 22°C. The winter diurnal temperature range is similar for both vegetation systems.

Inclusion of soil moisture variability in thermal properties results in an increase in mean daily maximum temperature of about 2°C in summer and winter, without much change in daily minima. The daily mean temperature is not significantly affected.     

Observations of urban snow properties in Calgary,Canada

Urban winter hydrology has garnered very little attention owing to the general notion that high‐intensity rainfalls are the major flood‐generating events in urban areas. As a result, few efforts have been made to research urban snow and its melt characteristics. This study investigates the characteristics of urban snow that differentiate it from rural snow, and makes recommendations for incorporating these characteristics into an urban snowmelt model. A field study was conducted from the fall of 2001 to the spring of 2002 in the city of Calgary, Canada. Snow depths and densities, soil moisture, soil temperature, snow albedo, net radiation, snow evaporation, and surface temperature were measured at several locations throughout the winter period. The combination of urban snow removal practices and the physical elements that exist in urban areas were found to influence the energy balance of the snowpack profoundly. Shortwave radiation was found to be the main source of energy for urban snow; as a consequence, the albedo of urban snow is a very important factor in urban snowmelt modelling. General observations lead to the classification of snow as one of four types: snow piles, snow on road shoulders, snow on sidewalk edges, and snow in open areas. This resulted in the development of four separate functions for the changing snow albedo values. A study of the frozen ground conditions revealed that antecedent soil moisture conditions had very little impact on frozen ground, and thus frozen ground very nearly always acts as a near impervious area. Improved flood forecasting for urban catchments in cold regions can only be achieved with accurate modelling of urban winter runoff that involves the energy balance method, incorporating snow redistribution and urban snow‐cover characteristics, and using small time steps. Copyright © 2004 John Wiley & Sons, Ltd.     

An analysis of mean transition time between flood and drought in the large river basin

 A soil moisture balance equation over large spatial regions is studied at seasonal and annual time scales for the Arkansas river basin. Interaction and feedback effects between land-surface and atmospheric moisture are studied in the parameterization for this basin. Due to the interaction between the land-surface and atmosphere at large scales, the surface hydrology of large land areas is susceptible to two distinct stable modes in the long-term probability density function: a dry and a wet state. In the soil moisture balance equation, stochastic fluctuations lead to separate preferred statistical stable states with transitions between these stable states induced by environmental fluctuations. On the basis of historical data, the soil moisture balance equation is calibrated for the Arkansas river basin. The transition times between the stable modes in the model are studied based on the stochastic representation of the physical processes and the calibrated model parameters. This study has implications for prediction of the transition times between stable modes or residence times, that is, the time the system spends in a given stable mode, since this would be equivalent to predicting the duration of droughts or wet conditions.     

How soil moisture mediates the influence of transpiration on streamflow at hourly to interannual scales in a forested catchment

The water balance equation dictates that streamflow may be reduced by transpiration. Yet temporal disequilibrium weakens the relationship between transpiration and streamflow in many cases where inputs and outputs are unbalanced. We address two critical knowledge barriers in ecohydrology with respect to time, scale dependence and lags. Study objectives were to correlate components of the water balance equation at hourly to annual scales, quantify time lags, and simplify critical components of the water budget during wet and dry conditions. We tested interrelationships among precipitation, vapour pressure deficit, transpiration, soil moisture, and streamflow within the confines of a 60‐hectare forested watershed in the western Cascades of Oregon. The Pacific Northwest is an ideal location to compare wet and dry seasons because of its Mediterranean climate. Soil moisture explained more than 80% of the variation in streamflow at all temporal scales investigated. Streamflow was most strongly coupled to soil moisture in the wet season because of gravitational drainage patterns; strong coupling of transpiration to vapour pressure deficit was dominant in the dry season and driven by low humidity. We observed progressively longer hourly time lags between soil moisture and streamflow in the dry season, which relates to an increasing soil moisture deficit that took an average of 48 days to refill after the onset of winter rains. We propose that transpiration drives seasonal patterns in soil moisture that relate to patterns in streamflow only after long time lags. In other words, soil moisture mediates the influence of transpiration on streamflow. Copyright © 2011 John Wiley & Sons, Ltd.     

Response of water and energy exchange to the environmental variable in a desert‐oasis wetland of Northwest China

A case study on a desert‐oasis wetland ecosystem in the arid region of Northwest China measured the seasonal and interannual variation in energy partitioning and evapotranspiration to analyse the response of water and energy exchange on soil moisture, groundwater, and environmental variables. Energy partitioning showed a clear seasonal and interannual variability, and the process of water and energy exchange differed significantly in the monthly and interannual scales. The net radiation was 7.31 MJ m^?2· day^?1, and sensible heat flux accounted for 50.42% of net radiation in energy fluxes, 40.56% for latent heat flux, and 9.02% for ground heat flux. The parameters in energy fluxes were best described by a unimodal curve, whereas sensible heat flux followed a bimodal curve. Variations in the daily evapotranspiration and crop evapotranspiration also exhibited a single peak curve with annual values of 569.84 and 644.47 mm, respectively. Canopy conductance averaged 20.77 ± 13.75 mm s^?1 and varied from 0.16 to 83.96 mm s^?1 during the two hydrological years. The variation in water and energy exchange reflected environmental conditions and depended primarily on vapour pressure deficit, net radiation, soil moisture, and water depth. Although the effects of precipitation on evapotranspiration showed that the response of this ecosystem to climate changes was not obvious, the variation of air temperatures had a strong influence on evapotranspiration, resulting in a significant increase in evapotranspiration (R = 0.730; P < 0.01). Copyright © 2013 John Wiley & Sons, Ltd.     

A temperature prediction-correction method for estimating surface soil heat flux from soil temperature and moisture data

Surface soil heat flux is a component of surface energy budget and its estimation is needed in land-atmosphere interaction studies. This paper develops a new simple method to estimate soil heat flux from soil temperature and moisture observations. It gives soil temperature profile with the thermal diffusion equation and, then, adjusts the temperature profile with differences between observed and computed soil temperatures. The soil flux is obtained through integrating the soil temperature profile. Compared with previous methods, the new method does not require accurate thermal conductivity. Case studies based on observations, synthetic data, and sensitivity analyses show that the new method is preferable and the results obtained with it are not sensitive to the availability of temperature data in the topsoil. In addition, we pointed out that the soil heat flux measured with a heat-plate can be quite erroneous in magnitude though its phase is accurate.     

Use of remote sensing for evapotranspiration monitoring over land surfaces

Abstract

Monitoring evapotranspiration (ET) at large scales is important for assessing climate and anthropogenic effects on natural and agricultural ecosystems. This paper describes techniques used in evaluating ET with remote sensing, which is the only technology that can efficiently and economically provide regional and global coverage. Some of the empirical/statistical techniques have been used operationally with satellite data for computing daily ET at regional scales. The more complex numerical simulation models require detailed input parameters that may limit their application to regions containing a large database of soils and vegetation properties. Current efforts are being directed towards simplifying the parameter requirements of these models. Essentially all energy balance models rely on an estimate of the available energy (net radiation less soil heat flux). Net radiation is not easily determined from space, although progress is being made. Simplified approaches for estimating soil heat flux appear promising for operational applications. In addition, most ET models utilize remote sensing data in the shortwave and thermal wavelengths to measure key boundary conditions. Differences between the radiometric surface temperature and aerodynamic temperature can be significant and progress in incorporating this effect is evident. Atmospheric effects on optical data are significant, and optical sensors cannot see through clouds. This has led some to use microwave observations as a surrogate for optical data to provide estimates of surface moisture and surface temperature; preliminary results are encouraging. The approaches that appear most promising use surface temperature and vegetation indices or a time rate of change in surface temperature coupled to an atmospheric boundary layer model. For many of these models, differences with ET observations can be as low as 20% from hourly to daily time scales, approaching the level of uncertainty in the measurement of ET and contradicting some recent pessimistic conclusions concerning the utility of remotely sensed radiometric surface temperature for determining the surface energy balance.     

Energy balance above a boreal coniferous forest: a difference in turbulent fluxes between snow‐covered and snow‐free canopies

To evaluate the interactive effects of snow and forest on turbulent fluxes between the forest surface and the atmosphere, the surface energy balance above a forest was measured by the eddy correlation method during the winter of 1995–1996. The forest was a young coniferous plantation comprised of spruce and fir. The study site, in Sapporo, northern Japan, had heavy and frequent snowfalls and the canopy was frequently covered with snow during the study period. A comparison of the observed energy balance above the forest for periods with and without a snow‐covered canopy and an analysis using a single‐source model gave the following results: during daytime when the canopy was covered with snow, the upward latent heat flux was large, about 80% of the net radiation, and the sensible heat flux was positive but small. On the other hand, during daytime when the canopy was dry and free from snow, the sensible heat flux was dominant and the latent heat flux was minor, about 10% of the net radiation. To explain this difference of energy partition between snow‐covered and snow‐free conditions, not only differences in temperature but also differences in the bulk transfer coefficients for latent heat flux were necessary in the model. Therefore, the high evaporation rate from the snow‐covered canopy can be attributed largely to the high moisture availability of the canopy surface. Evaporation from the forest during a 60‐day period in midwinter was estimated on a daily basis as net radiation minus sensible heat flux. The overall average evaporation during the 60‐day period was 0·6 mm day⁻¹, which is larger than that from open snow fields. Copyright © 1999 John Wiley & Sons, Ltd.     

Net ecosystem CO<Subscript>2</Subscript> exchange and controlling factors in a steppe—<Emphasis Type="Italic">Kobresia</Emphasis> meadow on the Tibetan Plateau

Knowledge of seasonal variation of net ecosystem CO₂ exchange (NEE) and its biotic and abiotic controllers will further our understanding of carbon cycling process, mechanism and large-scale modelling. Eddy covariance technique was used to measure NEE, biotic and abiotic factors for nearly 3 years in the hinterland alpine steppe—Korbresia meadow grassland on the Tibetan Plateau, the present highest fluxnet station in the world. The main objectives are to investigate dynamics of NEE and its components and to determine the major controlling factors. Maximum carbon assimilation took place in August and maximum carbon loss occurred in November. In June, rainfall amount due to monsoon climate played a great role in grass greening and consequently influenced interannual variation of ecosystem carbon gain. From July through September, monthly NEE presented net carbon assimilation. In other months, ecosystem exhibited carbon loss. In growing season, daytime NEE was mainly controlled by photosynthetically active radiation (PAR). In addition, leaf area index (LAI) interacted with PAR and together modulated NEE rates. Ecosystem respiration was controlled mainly by soil temperature and simultaneously by soil moisture. Q ₁₀ was negatively correlated with soil temperature but positively correlated with soil moisture. Large daily range of air temperature is not necessary to enhance carbon gain. Standard respiration rate at referenced 10°C (R ₁₀) was positively correlated with soil moisture, soil temperature, LAI and aboveground biomass. Rainfall patterns in growing season markedly influenced soil moisture and therefore soil moisture controlled seasonal change of ecosystem respiration. Pulse rainfall in the beginning and at the end of growing season induced great ecosystem respiration and consequently a great amount of carbon was lost. Short growing season and relative low temperature restrained alpine grass vegetation development. The results suggested that LAI be usually in a low level and carbon uptake be relatively low. Rainfall patterns in the growing season and pulse rainfall in the beginning and at end of growing season control ecosystem respiration and consequently influence carbon balance of ecosystem.     

北京大学陆面过程模式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能够较好地模拟西北半干旱区农田下垫面地气交换过程.     

Net ecosystem CO2 exchange and controlling factors in a steppe—Kobresia meadow on the Tibetan Plateau

Knowledge of seasonal variation of net ecosystem CO₂ exchange (NEE) and its biotic and abiotic controllers will further our understanding of carbon cycling process, mechanism and large-scale modelling. Eddy covariance technique was used to measure NEE, biotic and abiotic factors for nearly 3 years in the hinterland alpine steppe—Korbresia meadow grassland on the Tibetan Plateau, the present highest fluxnet station in the world. The main objectives are to investigate dynamics of NEE and its components and to determine the major controlling factors. Maximum carbon assimilation took place in August and maximum carbon loss occurred in November. In June, rainfall amount due to monsoon climate played a great role in grass greening and consequently influenced interannual variation of ecosystem carbon gain. From July through September, monthly NEE presented net carbon assimilation. In other months, ecosystem exhibited carbon loss. In growing season, daytime NEE was mainly controlled by photosynthetically active radiation (PAR). In addition, leaf area index (LAI) interacted with PAR and together modulated NEE rates. Ecosystem respiration was controlled mainly by soil temperature and simultaneously by soil moisture. Q ₁₀ was negatively correlated with soil temperature but positively correlated with soil moisture. Large daily range of air temperature is not necessary to enhance carbon gain. Standard respiration rate at referenced 10°C (R ₁₀) was positively correlated with soil moisture, soil temperature, LAI and aboveground biomass. Rainfall patterns in growing season markedly influenced soil moisture and therefore soil moisture controlled seasonal change of ecosystem respiration. Pulse rainfall in the beginning and at the end of growing season induced great ecosystem respiration and consequently a great amount of carbon was lost. Short growing season and relative low temperature restrained alpine grass vegetation development. The results suggested that LAI be usually in a low level and carbon uptake be relatively low. Rainfall patterns in the growing season and pulse rainfall in the beginning and at end of growing season control ecosystem respiration and consequently influence carbon balance of ecosystem.