Application of a three-dimensional model for assessing effects of small clear-cuttings on radiation and soil temperature

A three-dimensional model Mixfor-3D of soil–vegetation–atmosphere transfer (SVAT) was developed and applied to estimate possible effects of tree clear-cutting on radiation and soil temperature regimes of a forest ecosystem. The Mixfor-3D model consists of several closely coupled 3D sub-models describing: forest stand structure; radiative transfer in a forest canopy; turbulent transfer of sensible heat, H₂O and CO₂ between ground surface and the atmospheric surface layer; evapotranspiration of ground surface vegetation and soil; heat and moisture transfer in soil. The model operates with the horizontal grid resolution, 2 m × 2 m; vertical resolution, 1 m and primary time step, 1 h.     

A process-based model designed for filling of large data gaps in tower-based measurements of net ecosystem productivity

In this paper we present a simple hybrid gap-filling model (GFM) designed with a minimum number of parameters necessary to capture the ecological processes important for filling medium-to-large gaps in Flux data. As the model is process-based, the model has potential to be used in filling large gaps exhibiting a broad range of micro-meteorological and site conditions. The GFM performance was evaluated using “Punch hole” and extrapolation experiments based on data collected in west-central New Brunswick. These experiments indicated that the GFM is able to provide acceptable results (r² > 0.80) when >500 data points are used in model parameterization. The GFM was shown to address daytime evolution of NEP reasonably well for a wide range of weather and site conditions. An analysis of residuals indicated that for the most part no obvious trends were evident; although a slight bias was detected in NEP with soil temperature. To explore the portability of the GFM across ecosystem types, a transcontinental validation was conducted using NEP and ancillary data from seven ecosystems along a north-south transect (i.e., temperature–moisture gradient) from northern Europe (Finland) to the Middle East (Israel). The GFM was shown to explain over 75% of the variability in NEP measured at most ecosystems, which strongly suggests that the GFM maybe successfully applied to forest ecosystems outside Canada.     

Estimating the potential impact of vegetation on the water cycle requires accurate soil water parameter estimation

It is well known that vegetation dynamics at the catchment scale depends on the prevailing weather and soil moisture conditions. Soil moisture, however, is not equally distributed in space due to differences in topography, weather patterns, soil properties and the type and amount of vegetation cover. To elucidate the complex interaction between vegetation and soil moisture, the dynamic vegetation model LPJ-GUESS (Smith et al., 2001), which provides estimations of vegetation dynamics, but does not consider lateral water fluxes was coupled with the hydrological TOPMODEL (cf. Beven, 2001) in order to be able to evaluate the importance of these lateral fluxes. The new model LG-TM was calibrated and validated in two climatically different mountain catchments. The estimations of runoff were good, when monthly and weekly time scales were considered, although the low flow periods at winter time were somewhat underestimated. The uncertainty in the climate induced change vegetation carbon storage caused by the uncertainty in soil parameters was up to 3-5 kg C m⁻² (depending on elevation and catchment), compared to the total change in vegetation carbon storage of 5-9 kg C m⁻². Therefore accurate estimates of the parameters influencing the water holding capacity of the soil, for example depth and porosity, are necessary when estimating future changes in vegetation carbon storage. Similarly, changes in plant transpiration due to climatic changes could be almost double as high (88 mm m⁻²) in the not calibrated model compared to the new model version (ca 50 mm m⁻² transpiration change). The uncertainties in these soil properties were found to be more important than the lateral water exchange between grid cells, even in steep topography at least for the temporal and spatial resolution used here.     

Nitrogen controls on ecosystem carbon sequestration: a model implementation and application to Saskatchewan,Canada

A plant–soil nitrogen (N) cycling model was developed and incorporated into the Integrated BIosphere Simulator (IBIS) of Foley et al. [Foley, J.A., Prentice, I.C., Ramankutty, N., Levis, S., Pollard, D., Sitch, S., Haxeltine, A., 1996. An integrated biosphere model of land surface process, terrestrial carbon balance and vegetation dynamics. Global Biogeochem. Cycles 10, 603–628]. In the N-model, soil mineral N regulates ecosystem carbon (C) fluxes and ecosystem C:N ratios. Net primary productivity (NPP) is controlled by feedbacks from both leaf C:N and soil mineral N. Leaf C:N determines the foliar and canopy photosynthesis rates, while soil mineral N determines the N availability for plant growth and the efficiency of biomass construction. Nitrogen controls on the decomposition of soil organic matter (SOM) are implemented through N immobilization and mineralization separately. The model allows greater SOM mineralization at lower mineral N, and conversely, allows greater N immobilization at higher mineral N. The model's seasonal and inter-annual behaviours are demonstrated. A regional simulation for Saskatchewan, Canada, was performed for the period 1851–2000 at a 10 km × 10 km resolution. Simulated NPP was compared with high-resolution (1 km × 1 km) NPP estimated from remote sensing data using the boreal ecosystem productivity simulator (BEPS) [Liu, J., Chen, J.M., Cihlar, J., Park, W.M., 1997. A process-based boreal ecosystem productivity simulator using remote sensing inputs. Remote Sens. Environ. 44, 81–87]. The agreement between IBIS and BEPS, particularly in NPP spatial variation, was considerably improved when the N controls were introduced into IBIS.     

Are savannas patch-dynamic systems? A landscape model

Savannas are ecosystems characterized by the coexistence of woody species (trees and bushes) and grasses. Given that savanna characteristics are mainly formed from competition, herbivory, fire, woodcutting, and patchy soil and precipitation characteristics, we propose a spatially explicit model to examine the effects of the above-mentioned parameters on savanna vegetation dynamics in space and time. Furthermore, we investigate the effects of the above-mentioned parameters on tree–bush–grass ratios, as well as the degrees of aggregation of tree–bush–grass biomass. We parameterized our model for an arid savanna with shallow soil depth as well as a mesic one with generally deeper and more variable soil depths. Our model was able to reproduce savanna vegetation characteristics for periods of time over 2000 years with daily updated time steps. According to our results, tree biomass was higher than bush biomass in the arid savanna but bush biomass exceeded tree and grass biomass in the simulated mesic savanna. Woody biomass increased in our simulations when the soil's porosity values were increased (mesic savanna), in combination with higher precipitation. Savanna vegetation varied from open savanna to woodland and back to open savanna again. Vegetation cycles varied over ∼300-year cycles in the arid and ∼220-year cycles in the mesic-simulated savanna. Autocorrelation values indicated that there are both temporal and spatial vegetation cycles. Our model indicated cycling savanna vegetation at the landscape scale, cycles in cells, and patchiness, i.e. patch dynamics.     

Impact of climate change on soil moisture dynamics in Brandenburg with a focus on nature conservation areas

Global warming impacts the water cycle not only by changing regional precipitation levels and temporal variability, but also by affecting water flows and soil moisture dynamics. In Brandenburg, increasing average annual temperature and decreasing precipitation in summer have already been observed. For this study, past trends and future effects of climate change on soil moisture dynamics in Brandenburg were investigated, considering regional and specific spatial impacts. Special Areas of Conservation (SACs) were focused on in particular. A decreasing trend in soil water content was shown for the past by analyzing simulation results from 1951 to 2003 using the integrated ecohydrological model SWIM [Krysanova, V., Müller-Wohlfeil, D.-I., Becker, A., 1998. Development and test of a spatially distributed hydrological/water quality model for mesoscale watersheds. Ecol. Model. 106, 261–289]. The trend was statistically significant for some areas, but not for the entire region. Simulated soil water content was particularly low in the extremely dry year 2003. Comparisons of simulated trends in soil moisture dynamics with trends in the average annual Palmer Drought Severity Index for the region showed largely congruent patterns, though the modeled soil moisture trends are characterized by a much higher spatial resolution. Regionally downscaled climate change projections representing the range between wetter and drier realizations were used to evaluate future trends of available soil water. A further decrease of average available soil water ranging from −4% to −15% was projected for all climate realizations up to the middle of the 21st century. An average decrease of more than 25 mm was simulated for 34% of the total area in the dry realization. Available soil water contents in SACs were generally higher and trends in soil moisture dynamics were lower mainly due to their favorable edaphic conditions. Stronger absolute and relative changes in the simulated trends for the past and future were shown for SACs within Brandenburg than for the state as a whole, indicating a high level of risk for many wetland areas. Nonetheless, soil water content in SACs is expected to remain higher than average under climate change conditions as well, and SACs therefore have an important buffer function under the projected climate change. They are thus essential for local climate and water regulation and their status as protected areas in Brandenburg should be preserved.     

Development and optimization of an Agro-BGC ecosystem model for C4 perennial grasses

Extrapolating simulations of bioenergy crop agro-ecosystems beyond data-rich sites requires biophysically accurate ecosystem models and careful estimation of model parameters not available in the literature. To increase biophysical accuracy we added C₄ perennial grass functionality and agricultural practices to the Biome-BGC (BioGeochemical Cycles) ecosystem model. This new model, Agro-BGC, includes enzyme-driven C₄ photosynthesis, individual live and dead leaf, stem, and root carbon and nitrogen pools, separate senescence and litter fall processes, fruit growth, optional annual seeding, flood irrigation, a growing degree day phenology with a killing frost option, and a disturbance handler that simulates nitrogen fertilization, harvest, fire, and incremental irrigation. To obtain spatially generalizable vegetation parameters we used a numerical method to optimize five unavailable parameters for Panicum virgatum (switchgrass) using biomass yield data from three sites: Mead, Nebraska, Rockspring, Pennsylvania, and Mandan, North Dakota. We then verified simulated switchgrass yields at three independent sites in Illinois (IL). Agro-BGC is more accurate than Biome-BGC in representing the physiology and dynamics of C₄ grass and management practices associated with agro-ecosystems. The simulated two-year average mature yields with single-site Rockspring optimization have Root Mean Square Errors (RMSE) of 70, 152, and 162 and biases of 43, −87, 156 g carbon m⁻² for Shabbona, Urbana, and Simpson IL, respectively. The simulated annual yields in June, August, October, December, and February have RMSEs of 114, 390, and 185 and biases of −19, −258, and 147 g carbon m⁻² for Shabbona, Urbana, and Simpson IL, respectively. These RMSE and bias values are all within the largest 90% confidence interval around respective IL site measurements. Twenty-four of twenty-six simulated annual yields with Rockspring optimization are within 95% confidence intervals of Illinois site measurements during the mature fourth and fifth years of growth. Ten of eleven simulated two-year average mature yields with Rockspring optimization are within 65% confidence intervals of Illinois site measurements and the eleventh is within the 95% confidence interval. Rockspring optimized Agro-BGC achieves accuracies comparable to those of two previously published models: Agricultural Land Management Alternatives with Numerical Assessment Criteria (ALMANAC) and Integrated Farm System Model (IFSM). Agro-BGC suffers from static vegetation parameters that can change seasonally and as plants age. Using mature plant data for optimization mitigates this deficiency. Our results suggest that a multi-site optimization scheme using mature plant data from more sites would be adequate for generating spatially generalizable vegetation parameters for simulating mature bioenergy crop agro-ecosystems with Agro-BGC.     

Estimating soil carbon turnover using radiocarbon data: A case-study for European Russia

Turnover rates of soil carbon for 20 soil types typical for a 3.7 million km² area of European Russia were estimated based on ¹⁴C data. The rates are corrected for bomb radiocarbon which strongly affects the topsoil ¹⁴C balance. The approach is applied for carbon stored in the organic and mineral layers of the upper 1 m of the soil profile. The turnover rates of carbon in the upper 20 cm are relatively high for forest soils (0.16–0.78% year⁻¹), intermediate for tundra soils (0.25% year⁻¹), and low for grassland soils (0.02–0.08% year⁻¹) with the exception of southern Chernozems (0.32% year⁻¹). In the soil layer of 20–100 cm depth, the turnover rates were much lower for all soil types (0.01–0.06% year⁻¹) except for peat bog soils of the southern taiga (0.14% year⁻¹). Combined with a map of soil type distribution and a dataset of several hundred soil carbon profiles, the method provides annual fluxes for the slowest components of soil carbon assuming that the latter is in equilibrium with climate and vegetation cover. The estimated carbon flux from the soil is highest for forest soils (12–147 gC/(m² year)), intermediate for tundra soils (33 gC/(m² year)), and lowest for grassland soils (1–26 gC/(m² year)). The approach does not distinguish active and recalcitrant carbon fractions and this explains the low turnover rates in the top layer. Since changes in soil types will follow changes in climate and land cover, we suggest that pedogenesis is an important factor influencing the future dynamics of soil carbon fluxes. Up to now, both the effect of soil type changes and the clear evidence from ¹⁴C measurements that most soil organic carbon has a millennial time scale, are basically neglected in the global carbon cycle models used for projections of atmospheric CO₂ in 21st century and beyond.     

Dissolved organic carbon concentrations and fluxes in forest catchments and streams: DOC-3 model

Dissolved organic carbon (DOC) concentrations in south-western Nova Scotia streams, sampled at weekly to biweekly intervals, varied across streams from about 3 to 40 mg L⁻¹, being highest mid-summer to fall, and lowest during winter to spring. A 3-parameter model (DOC-3) was proposed to project daily stream DOC concentrations and fluxes from modelled estimates for daily soil temperature and moisture, year-round, and in relation to basin size and wetness. The parameters of this model refer to (i) a basin-specific DOC release parameter “k_DOC^”, related to the wet- and open-water area percentages per basin, (ii) the lag time “τ” between DOC production and subsequent stream DOC emergence, related to the catchment area above the stream sampling location; and (iii) the activation energy “Ea”, to deal with the temperature effect on DOC production. This model was calibrated with the 1988-2006 DOC concentration data from three streams (Pine Marten, Moosepit Brook, and the Mersey River sampled at or near Kejimkujik National Park, or KNP), and was used to interpret the biweekly 1999-2003 DOC concentrations data (stream, ground and lake water, soil lysimeters) of the Pockwock-Bowater Watershed Project near Halifax, Nova Scotia. The data and the model revealed that the DOC concentrations within the streams were not correlated to the DOC concentrations within the soil- and groundwater, but were predictable based on (i) the hydrologically inferred weather-induced changes in soil moisture and temperature next to each stream, and (ii) the topographically inferred basin area and wet- and open-water area percentages associated with each stream (R² = 0.53; RMSE = 3.5 mg L⁻¹). Model-predicted fluxes accounted 74% of the hydrometrically determined DOC exports at KNP.     

Simulating soil organic matter with CQESTR (v. 2.0): Model description and validation against long-term experiments across North America

Soil carbon (C) models are important tools for examining complex interactions between climate, crop and soil management practices, and to evaluate the long-term effects of management practices on C-storage potential in soils. CQESTR is a process-based carbon balance model that relates crop residue additions and crop and soil management to soil organic matter (SOM) accretion or loss. This model was developed for national use in U.S and calibrated initially in the Pacific Northwest. Our objectives were: (i) to revise the model, making it more applicable for wider geographic areas including potential international application, by modifying the thermal effect and incorporating soil texture and drainage effects, and (ii) to recalibrate and validate it for an extended range of soil properties and climate conditions. The current version of CQESTR (v. 2.0) is presented with the algorithms necessary to simulate SOM at field scale. Input data for SOM calculation include crop rotation, aboveground and belowground biomass additions, tillage, weather, and the nitrogen content of crop residues and any organic amendments. The model was validated with long-term data from across North America. Regression analysis of 306 pairs of predicted and measured SOM data under diverse climate, soil texture and drainage classes, and agronomic practices at 13 agricultural sites having a range of SOM (7.3–57.9 g SOM kg⁻¹), resulted in a linear relationship with an r² of 0.95 (P < 0.0001) and a 95% confidence interval of 4.3 g SOM kg⁻¹. Using the same data the version 1.0 of CQESTR had an r² of 0.71 with a 95% confidence interval of 5.5 g SOM kg⁻¹. The model can be used as a tool to predict and evaluate SOM changes from various management practices and offers the potential to estimate C accretion required for C credits.     

Response of tree growth to a changing climate in boreal central Canada: A comparison of empirical,process-based,and hybrid modelling approaches

The impact of 2 × CO₂ driven climate change on radial growth of boreal tree species Pinus banksiana Lamb., Populus tremuloides Michx. and Picea mariana (Mill.) BSP growing in the Duck Mountain Provincial Forest of Manitoba (DMPF), Canada, is simulated using empirical and process-based model approaches. First, empirical relationships between growth and climate are developed. Stepwise multiple-regression models are conducted between tree-ring growth increments (TRGI) and monthly drought, precipitation and temperature series. Predictive skills are tested using a calibration–verification scheme. The established relationships are then transferred to climates driven by 1× and 2 × CO₂ scenarios using outputs from the Canadian second-generation coupled global climate model. Second, empirical results are contrasted with process-based projections of net primary productivity allocated to stem development (NPP_s). At the finest scale, a leaf-level model of photosynthesis is used to simulate canopy properties per species and their interaction with the variability in radiation, temperature and vapour pressure deficit. Then, a top-down plot-level model of forest productivity is used to simulate landscape-level productivity by capturing the between-stand variability in forest cover. Results show that the predicted TRGI from the empirical models account for up to 56.3% of the variance in the observed TRGI over the period 1912–1999. Under a 2 × CO₂ scenario, the predicted impact of climate change is a radial growth decline for all three species under study. However, projections obtained from the process-based model suggest that an increasing growing season length in a changing climate could counteract and potentially overwhelm the negative influence of increased drought stress. The divergence between TRGI and NPP_s simulations likely resulted, among others, from assumptions about soil water holding capacity and from calibration of variables affecting gross primary productivity. An attempt was therefore made to bridge the gap between the two modelling approaches by using physiological variables as TRGI predictors. Results obtained in this manner are similar to those obtained using climate variables, and suggest that the positive effect of increasing growing season length would be counteracted by increasing summer temperatures. Notwithstanding uncertainties in these simulations (CO₂ fertilization effect, feedback from disturbance regimes, phenology of species, and uncertainties in future CO₂ emissions), a decrease in forest productivity with climate change should be considered as a plausible scenario in sustainable forest management planning of the DMPF.     

Mechanisms of moisture stress in a mid-latitude temperate forest: Implications for feedforward and feedback controls from an irrigation experiment

While it is well established that stomata close during moisture stress, strong correlations among environmental (e.g., vapor pressure deficit, soil moisture, air temperature, radiation) and internal (e.g., leaf water potential, sap flow, root-shoot signaling) variables obscure the identification of causal mechanisms from field experiments. Models of stomatal control fitted to field data therefore suffer from ambiguous parameter identification, with multiple acceptable (i.e., nearly optimal) model structures emphasizing different moisture status indicators and different processes. In an effort to minimize these correlations and improve parameter and process identification, we conducted an irrigation experiment on red maples (Acer rubrum L.) at Harvard Forest (summers of 2005 and 2006). Control and irrigated trees experienced similar radiative and boundary layer forcings, but different soil moisture status, and thus presumably different diurnal cycles of internal leaf water potential. Measured soil moisture and atmospheric forcing were used to drive a transient tree hydraulic model that incorporated a Jarvis-type leaf conductance in a Penman–Monteith framework with a Cowan-type (resistance and capacitance) tree hydraulic representation. The leaf conductance model included dependence on both leaf matric potential, Ψ_L (so-called feedback control) and on vapor pressure deficit, D (so-called feedforward control). Model parameters were estimated by minimizing the error between predicted and measured sap flow. The whole-tree irrigation treatment had the effect of elevating measured transpiration during summer dry-downs, demonstrating the limiting effect that subsurface resistance may have on transpiration during these times of moisture stress. From the best fitted model, we infer that during dry downs, moisture stress manifests itself in an increase of soil resistance with a resulting decrease in Ψ_L, leading to both feedforward and feedback controls in the control trees, but only feedforward control for the irrigated set. Increases in the sum-of-squares error when individual model components were disabled allow us to reject the following three null hypotheses: (1) the f(D) stress is statistically insignificant (p = 0.01); (2) the f(Ψ_L) stress is statistically insignificant (p = 0.07); and (3) plant storage capacitance is independent of moisture status (p = 0.07).     

Interpreting spatial heterogeneity of crop yield with a process model and remote sensing

A process-based crop growth model (Vegetation Interface Processes (VIP) model) is used to estimate crop yield with remote sensing over the North China Plain. Spatial pattern of the key parameter—maximum catalytic capacity of Rubisco (V_cmax) for assimilation is retrieved from Normalized Difference of Vegetation Index (NDVI) from Terra-MODIS and statistical yield records. The regional simulation shows that the agreements between the simulated winter wheat yields and census data at county-level are quite well with R² being 0.41-0.50 during 2001-2005. Spatial variability of photosynthetic capacity and yield in irrigated regions depend greatly on nitrogen input. Due to the heavy soil salinity, the photosynthetic capacity and yield in coastal region is less than 50 μmol C m⁻² s⁻¹ and 3000 kg ha⁻¹, respectively, which are much lower than that in non-salinized region, 84.5 μmol C m⁻² s⁻¹ and 5700 kg ha⁻¹. The predicted yield for irrigated wheat ranges from 4000 to 7800 kg ha⁻¹, which is significantly larger than that of rainfed, 1500-3000 kg ha⁻¹. According to the path coefficient analysis, nitrogen significantly affects yield, by which water exerts noticeably indirect influences on yield. The effect of water on yield is regulated, to a certain extent, by crop photosynthetic capacity and nitrogen application. It is believed that photosynthetic parameters retrieved from remote sensing are reliable for regional production prediction with a process-based model.     

Modelling the impact of defoliation and leaf damage on forest plantation function and production

After presenting a short review of process-based model requirements to capture the plant dynamic response to defoliation, this paper describes the development and testing of a model of crown damage and defoliation for Eucalyptus. A model that calculates light interception and photosynthetic production for canopies that vary spatially and temporally in leaf area and photosynthetic properties is linked to the forest growth model CABALA. The process of photosynthetic up-regulation following defoliation is modelled with a simple conditional switch that triggers up-regulation when foliar damage or removal causes the ratio of functional leaf area to living tissue in the tree to change.We show that the model predicts satisfactorily when validated with trees of Eucalyptus nitens and Eucalyptus globulus from a range of sites of different ages, subject to different types of stress and different types of defoliation events (R² = 0.96 across a range of sites). However, the complexity of particular situations can cause the model to fail (e.g. very heavy defoliation events where branch death occurs).It is concluded that while the model will not cope with all situations, an appropriate level of generality has been captured to represent many of the physiological processes and feedbacks that occur following defoliation or leaf damage. This makes the model useful for guiding management interventions following pest attack and allows the development of scenarios including climate change impact analyses and decision-making on the merits of post-defoliation fertilisation to expedite recovery.     

Spatially distributed modeling of the long-term carbon balance of a boreal landscape

Spatially and temporally distributed information on the sizes of biomass carbon (C) pools (BCPs) and soil C pools (SCPs) is vital for improving our understanding of biosphere-atmosphere C fluxes. Because the sizes of C pools result from the integrated effects of primary production, age-effects, changes in climate, atmospheric CO₂ concentration, N deposition, and disturbances, a modeling scheme that interactively considers these processes is important. We used the InTEC model, driven by various spatio-temporal datasets to simulate the long-term C-balance in a boreal landscape in eastern Canada. Our results suggested that in this boreal landscape, mature coniferous stands had stabilized their productivity and fluctuated as a weak C-sink or C-source depending on the interannual variations in hydrometeorological factors. Disturbed deciduous stands were larger C-sinks (NEP₂₀₀₄ = 150 gC m⁻² yr⁻¹) than undisturbed coniferous stands (e.g. NEP₂₀₀₄ = 8 gC m⁻² yr⁻¹). Wetlands had lower NPP but showed temporally consistent C accumulation patterns. The simulated spatio-temporal patterns of BCPs and SCPs were unique and reflected the integrated effects of climate, plant growth and atmospheric chemistry besides the inherent properties of the C pool themselves. The simulated BCPs and SCPs generally compared well with the biometric estimates (BCPs: r = 0.86, SCPs: r = 0.84). The largest BCP biases were found in recently disturbed stands and the largest SCP biases were seen in locations where moss necro-masses were abundant. Reconstructing C pools and C fluxes in the ecosystem in such a spatio-temporal manner could help reduce the uncertainties in our understanding of terrestrial C-cycle.     

库布齐沙漠自然植被与人工植被土壤水盐的空间异质性比较研究

干旱半干旱地区土壤水盐的空间分布对土地利用和生态恢复具有重要作用。运用传统统计学和地统计学对库布齐沙漠5种自然植被和5种人工植被0～10和10～20cm深度土壤水分和盐分的空间异质性进行了小尺度比较分析。结果表明：9种群落的土壤水盐平均值下层大于表层,且盐分的空间相关性较水分更高;人工植被土壤水分（CV=5.3%～22.7%）和盐分（CV=15.7%～51.2%）具有空间分布均匀、层间差异不明显等特征,而自然植被水分（CV=9.9%～32.6%）和盐分的（CV=26.9%～180.0%）却与之相反;小尺度上,人工植被土壤不同层间的水分关系、盐分关系以及水盐关系可能随建植时间的增加会越来越不明显,格局强度将不断减弱,这也极有可能改变大尺度上的水盐运移状况,进而影响研究区生态系统的稳定性和安全性。     

Modeling seasonal course of carbon fluxes and evapotranspiration in response to low temperature and moisture in a boreal Scots pine ecosystem

Environmental conditions act above and below ground, and regulate carbon fluxes and evapotranspiration. The productivity of boreal forest ecosystems is strongly governed by low temperature and moisture conditions, but the understanding of various feedbacks between vegetation and environmental conditions is still unclear. In order to quantify the seasonal responses of vegetation to environmental factors, the seasonality of carbon and heat fluxes and the corresponding responses for temperature and moisture in air and soil were simulated by merging a process-based model (CoupModel) with detailed measurements representing various components of a forest ecosystem in Hyytiälä, southern Finland. The uncertainties in parameters, model assumptions, and measurements were identified by generalized likelihood uncertainty estimation (GLUE). Seasonal and diurnal courses of sensible and latent heat fluxes and net ecosystem exchange (NEE) of CO₂ were successfully simulated for two contrasting years. Moreover, systematic increases in efficiency of photosynthesis, water uptake, and decomposition occurred from spring to summer, demonstrating the strong coupling between processes. Evapotranspiration and NEE flux both showed a strong response to soil temperature conditions via different direct and indirect ecosystem mechanisms. The rate of photosynthesis was strongly correlated with the corresponding water uptake response and the light use efficiency. With the present data and model assumptions, it was not possible to precisely distinguish the various regulating ecosystem mechanisms. Our approach proved robust for modeling the seasonal course of carbon fluxes and evapotranspiration by combining different independent measurements. It will be highly interesting to continue using long-term series data and to make additional tests of optional stomatal conductance models in order to improve our understanding of the boreal forest ecosystem in response to climate variability and environmental conditions.     

A general process-based mass-balance model for phosphorus/eutrophication as a tool to estimate historical reference values for key bioindicators,as exemplified using data for the Gulf of Riga

This work describes how a general, process-based mass-balance model (CoastMab) for phosphorus for coastal areas may be used as a tool to estimate realistic values of “natural” or preindustrial reference levels of key bioindicators in coastal science, including the Secchi depth, a standard measure of water clarity, the chlorophyll-a concentration, an operational measure of phytoplankton biomass and the concentration of cyanobacteria, a measure of the concentration of harmful algae. The CoastMab-model is an ecosystem model giving monthly predictions to achieve seasonal variations of basin-wide properties. The selected case-study area, the Gulf of Riga, is sensitive to nutrient loading because of its shallowness and low openness towards the Baltic Proper. The morphometry of any coastal area, as given by the size and form parameters, influences all internal processes, such as sedimentation, resuspension, diffusion in water and from sediments to water, biouptake and retention in biota, stratification, mixing and outflow. There has been no mass-balance modeling for nitrogen (N) in this work because empirical data (from the HELCOM database) clearly indicate that the monthly primary production in the Gulf of Riga is regulated by phosphorus (P) – the mean monthly total-N to total-P ratios are well over 7.2 (the Redfield-ratio) and generally higher than 15 for the data used in this study (from 1992 to 2005). At present anthropogenic loads, the average modeled monthly values for Secchi depth, chlorophyll (Chl), cyanobacteria (CB) and total-P (TP) are 3.2 m, 3.8 μg/l, 78 μg/l and 31.3 μg/l, respectively. If 50% of all anthropogenic sources to the Gulf of Riga via rivers, point sources and diffuse sources were to be removed, these values would be 3.6 m, 3.4 μg Chl/l, 63 μg CB/l and 29.1 μg TP/l. If 60% of the anthropogenic phosphorus fluxes to the Baltic Proper were to be omitted and as well as 75% of all direct anthropogenic sources to the Gulf of Riga, the values would be 4.6 m, 2.7 μg Chl/l, 45 μg CB/l and 25.4 μg TP/l. These values represent the “natural” reference levels and it is not realistic to expect that remedial measures would improve the conditions more than that. Using the CoastWeb-model, similar calculations can be made for any given coastal area and the data necessary for such calculations are discussed in this work.     

Modeling ecosystem responses to prescribed fires in a phosphorus-enriched Everglades wetland: I. Phosphorus dynamics and cattail recovery

We have developed and applied a process-based model, the Wetland Ecosystem Model (WEM), to evaluate the effects of a prescribed fire on the phosphorus (P) dynamics and cattail (Typha domingensis) growth in a P-enriched area in the Florida Everglades. The WEM couples major ecosystem processes including carbon (C), nitrogen (N) and P biogeochemical cycles, plant growth, hydrology, and fire disturbance. The model is used to assess the effects of a prescribed fire on P dynamics and cattail growth through dynamic interaction among four modules: fire, water chemistry, soil, and vegetation. The simulation results are in agreement with observed data including cattail above- and belowground biomass and dead mass, P concentration in surface-water, pore-water, and soil, and soil and water temperature. Cattail aboveground biomass reached the unburned level one year after burn; belowground biomass recovered to unburned level one and half years after the fire, however, dead mass did not completely reach unburned level two years after fires. The fire increased water and soil temperatures in the short term, while indirectly increasing the sensitivity of water and soil temperature post-fire response to air temperature by altering the energy exchange between air and water through a canopy gap created by fire. The fire also altered the P dynamics in surface-water and pore-water. A post-fire P pulse that lasted for less than one month was observed in surface-water. A similar P pulse, but in a small magnitude and a longer duration, was also observed in the pore-water total phosphorus (TP), and then came back to normal level after approximately three months. No significant changes in soil TP was observed during the study period. Meanwhile, no significant changes in water nutrients were observed downstream of the study plot. This finding indicated that the P-enriched wetlands in Everglades act as a buffer in regulating the P concentration in surface-water. Our study showed that the distance of fire effects on a 300 m × 300 m plot was less than 300 m downstream. Sensitivity analysis identified that the air temperature and hydrological conditions are two important driving factors which may alter the cattail community dynamics in response to prescribed fires. Similar to the filed studies, this study provided evidences that fire played an important role in managing plant growth and P dynamics in the Florida Everglades.     

Dealing with microtopography of an ombrotrophic bog for simulating ecosystem-level CO2 exchanges

Peatlands contain approximately 25% of the global soil carbon (C), despite covering only 3% of the earth's land surface. In order to evaluate the role of peatlands in global C cycling, models of ecosystem biogeochemistry are required, but peatland ecosystems present a number of unique challenges, particularly how to deal with the large variability that occurs at scales of one to several metres. In models, spatial variability is considered either explicitly for each individual unit and the outputs averaged, referred to as flux upscaling, or implicitly by weighting model parameters by the fractional occurrence of the individual units, referred to as parameter upscaling. The advantage of parameter upscaling is that it is much more computationally efficient: a requirement for hemispheric scale simulations. In this study we determined the differences between modelling a raised bog peatland with hummock-hollow microtopography using flux and parameter upscaling. We used the McGill Wetland Model (MWM), a process-based ecosystem C model for peatlands, configured for hummocks and hollows separately and then a weighted mixture of both. The simulated output based on flux and parameter upscaling was compared with eddy-covariance tower measurements. We found that net ecosystem production (NEP) for hollows was much larger than that for hummocks because total ecosystem respiration (TER) for hummocks was greater while gross primary production (GPP) did not differ significantly between the two topographic features. However, despite differences in components of NEP between hummocks and hollows, there was no statistically significant difference between the NEP based on flux and parameter upscaling using the MWM. Both flux and parameter upscaling show equivalent capability to capture the magnitude, direction, seasonality and inter-annual variability. The root-mean-square-errors (RMSE) are 0.66, 0.45, and 0.49 g C m⁻² day⁻¹, respectively for GPP, TER and NEP based on the flux upscaling, while 0.67, 0.44, and 0.48 g C m⁻² day⁻¹, respectively based on the parameter upscaling. The degree of agreement (d*) is 0.96, 0.97, and 0.88, respectively for GPP, TER and NEP based on the flux upscaling, while 0.96, 0.97, and 0.89, respectively based on the parameter upscaling. This result suggests that differences in processes caused by peatland microtopography scale linearly, which means an ecosystem-level model set-up (i.e. parameter upscaling scheme), is sufficient to simulate the C cycling.