Carbon sequestration in soils of SW-Germany as affected by agricultural management—Calibration of the EPIC model for regional simulations

Global emissions trading allows for agricultural measures to be accounted for the carbon sequestration in soils. The Environmental Policy Integrated Climate (EPIC) model was tested for central European site conditions by means of agricultural extensification scenarios. Results of soil and management analyses of different management systems (cultivation with mouldboard plough, reduced tillage, and grassland/fallow establishment) on 13 representative sites in the German State Baden-Württemberg were used to calibrate the EPIC model. Calibration results were compared to those of the Intergovernmental Panel on Climate Change (IPCC) prognosis tool. The first calibration step included adjustments in (a) N depositions, (b) N₂-fixation by bacteria during fallow, and (c) nutrient content of organic fertilisers according to regional values. The mixing efficiency of implements used for reduced tillage and four crop parameters were adapted to site conditions as a second step of the iterative calibration process, which should optimise the agreement between measured and simulated humus changes. Thus, general rules were obtained for the calibration of EPIC for different criteria and regions. EPIC simulated an average increase of +0.341 Mg humus-C ha⁻¹ a⁻¹ for on average 11.3 years of reduced tillage compared to land cultivated with mouldboard plough during the same time scale. Field measurements revealed an average increase of +0.343 Mg C ha⁻¹ a⁻¹ and the IPCC prognosis tool +0.345 Mg C ha⁻¹ a⁻¹. EPIC simulated an average increase of +1.253 Mg C ha⁻¹ a⁻¹ for on average 10.6 years of grassland/fallow establishment compared to an average increase of +1.342 Mg humus-C ha⁻¹ a⁻¹ measured by field measurements and +1.254 Mg C ha⁻¹ a⁻¹ according to the IPCC prognosis tool. The comparison of simulated and measured humus C stocks was r² ≥ 0.825 for all treatments. However, on some sites deviations between simulated and measured results were considerable. The result for the simulation of yields was similar. In 49% of the cases the simulated yields differed from the surveyed ones by more than 20%. Some explanations could be found by qualitative cause analyses. Yet, for quantitative analyses the available information from farmers was not sufficient. Altogether EPIC is able to represent the expected changes by reduced tillage or grassland/fallow establishment acceptably under central European site conditions of south-western Germany.     

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.     

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.     

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.     

Modelling soil carbon development in Swedish coniferous forest soils—An uncertainty analysis of parameters and model estimates using the GLUE method

Boreal forest soils such as those in Sweden contain a large active carbon stock. Hence, a relatively small change in this stock can have a major impact on the Swedish national CO₂ balance. Understanding of the uncertainties in the estimations of soil carbon pools is critical for accurately assessing changes in carbon stocks in the national reports to UNFCCC and the Kyoto Protocol. Our objective was to analyse the parameter uncertainties of simulated estimates of the soil organic carbon (SOC) development between 1994 and 2002 in Swedish coniferous forests with the Q model. Both the sensitivity of model parameters and the uncertainties in simulations were assessed. Data of forests with Norway spruce, Scots pine and Lodgepole pine, from the Swedish Forest Soil Inventory (SFSI) were used. Data of 12 Swedish counties were used to calibrate parameter settings; and data from another 11 counties to validate. The “limits of acceptability” within GLUE were set at the 95% confidence interval for the annual, mean measured SOC at county scale. The calibration procedure reduced the parameter uncertainties and reshaped the distributions of the parameters county-specific. The average measured and simulated SOC amounts varied from 60 t C ha⁻¹ in northern to 140 t C ha⁻¹ in the southern Sweden. The calibrated model simulated the soil carbon pool within the limits of acceptability for all calibration counties except for one county during one year. The efficiency of the calibrated model varied strongly; for five out of 12 counties the model estimates agreed well with measurements, for two counties agreement was moderate and for five counties the agreement was poor. The lack of agreement can be explained with the high inter-annual variability of the down-scaled measured SOC estimates and changes in forest areas over time. We conclude that, although we succeed in reducing the uncertainty in the model estimates, calibrating of a regional scale process-oriented model using a national scale dataset is a sensitive balance between introducing and reducing uncertainties. Parameter distributions showed to be scale sensitive and county specific. Further analysis of uncertainties in the methods used for reporting SOC changes to the UNFCCC and Kyoto protocol is recommended.     

Estimating California ecosystem carbon change using process model and land cover disturbance data: 1951-2000

Land use change, natural disturbance, and climate change directly alter ecosystem productivity and carbon stock level. The estimation of ecosystem carbon dynamics depends on the quality of land cover change data and the effectiveness of the ecosystem models that represent the vegetation growth processes and disturbance effects. We used the Integrated Biosphere Simulator (IBIS) and a set of 30- to 60-m resolution fire and land cover change data to examine the carbon changes of California's forests, shrublands, and grasslands. Simulation results indicate that during 1951-2000, the net primary productivity (NPP) increased by 7%, from 72.2 to 77.1 Tg C yr⁻¹ (1 teragram = 10¹² g), mainly due to CO₂ fertilization, since the climate hardly changed during this period. Similarly, heterotrophic respiration increased by 5%, from 69.4 to 73.1 Tg C yr⁻¹, mainly due to increased forest soil carbon and temperature. Net ecosystem production (NEP) was highly variable in the 50-year period but on average equalled 3.0 Tg C yr⁻¹ (total of 149 Tg C). As with NEP, the net biome production (NBP) was also highly variable but averaged −0.55 Tg C yr⁻¹ (total of -27.3 Tg C) because NBP in the 1980s was very low (-5.34 Tg C yr⁻¹). During the study period, a total of 126 Tg carbon were removed by logging and land use change, and 50 Tg carbon were directly removed by wildland fires. For carbon pools, the estimated total living upper canopy (tree) biomass decreased from 928 to 834 Tg C, and the understory (including shrub and grass) biomass increased from 59 to 63 Tg C. Soil carbon and dead biomass carbon increased from 1136 to 1197 Tg C.Our analyses suggest that both natural and human processes have significant influence on the carbon change in California. During 1951-2000, climate interannual variability was the key driving force for the large interannual changes of ecosystem carbon source and sink at the state level, while logging and fire were the dominant driving forces for carbon balances in several specific ecoregions. From a long-term perspective, CO₂ fertilization plays a key role in maintaining higher NPP. However, our study shows that the increase in C sequestration by CO₂ fertilization is largely offset by logging/land use change and wildland fires.     

Inter-annual variability of ecosystem production in boreal jack pine forests (1975-2004) estimated from tree-ring data using CBM-CFS3

We describe and apply a method of using tree-ring data and an ecosystem model to reconstruct past annual rates of ecosystem production. Annual data on merchantable wood volume increment and mortality obtained by dendrochronological stand reconstruction were used as input to the Carbon Budget Model of the Canadian Forest Sector (CBM-CFS3) to estimate net ecosystem production (NEP), net primary production (NPP), and heterotrophic respiration (Rh) annually from 1975 to 2004 at 10 boreal jack pine (Pinus banksiana Lamb.) stands in Saskatchewan and Manitoba, Canada. From 1975 (when sites aged 41-60 years) to 2004 (when they aged 70-89 years), all sites were moderate C sinks except during some warmer than average years where estimated Rh increased. Across all sites and years, estimated annual NEP averaged 57 g Cm⁻² yr⁻¹ (range −31 to 176 g Cm⁻² yr⁻¹), NPP 244 g Cm⁻² yr⁻¹ (147-376 g Cm⁻² yr⁻¹), and Rh 187 g Cm⁻² yr⁻¹ (124-270 g Cm⁻² yr⁻¹). Across all sites, NPP was related to stand age and density, which are proxies for successional changes in leaf area. Regionally, warm spring temperature increased NPP and defoliation by jack pine budworm 1 year previously reduced NPP. Our estimates of NPP, Rh, and NEP were plausible when compared to regional eddy covariance and carbon stock measurements. Inter-annual variability in ecosystem productivity contributes uncertainty to inventory-based assessments of regional forest C budgets that use yield curves predicting averaged growth over time. Our method could expand the spatial and temporal coverage of annual forest productivity estimates, providing additional data for the development of empirical models accounting for factors not presently considered by these models.     

Carbon, nitrogen, oxygen and sulfide budgets in the Black Sea: A biogeochemical model of the whole water column coupling the oxic and anoxic parts

Carbon, nitrogen, oxygen and sulfide budgets are derived for the Black Sea water column from a coupled physical-biogeochemical model. The model is applied in the deep part of the sea and simulates processes over the whole water column including the anoxic layer that extends from ?115 m to the bottom (?2000 m). The biogeochemical model involves a refined representation of the Black Sea foodweb from bacteria to gelatinous carnivores. It includes notably a series of biogeochemical processes typical for oxygen deficient conditions with, for instance, bacterial respiration using different types of oxidants (i.e denitrification, sulfate reduction), the lower efficiency of detritus degradation, the ANAMMOX (ANaerobic AMMonium OXidation) process and the occurrence of particular redox reactions. The model has been calibrated and validated against all available data gathered in the Black Sea TU Ocean Base and this exercise is described in Gregoire et al. (2008). In the present paper, we focus on the biogeochemical flows produced by the model and we compare model estimations with the measurements performed during the R.V. KNORR expedition conducted in the Black Sea from April to July 1988 (Murray and the Black Sea Knorr Expedition, 1991). Model estimations of hydrogen sulfide oxidation, metal sulfide precipitation, hydrogen sulfide formation in the sediments and water column, export flux to the anoxic layer and to the sediments, denitrification, primary and bacterial production are in the range of field observations.With a simulated Gross Primary Production (GPP) of 7.9 mol C m⁻² year⁻¹ and a Community Respiration (CR) of 6.3 mol C m⁻² year⁻¹, the system is net autotrophic with a Net Community Production (NCP) of 1.6 mol C m⁻² year⁻¹. This NCP corresponds to 20% of the GPP and is exported to the anoxic layer. In order to model Particulate Organic Matter (POM) fluxes to the bottom and hydrogen sulfide profiles in agreement with in situ observations, we have to consider that the degradation of POM in anoxic conditions is less efficient that in oxygenated waters as it has often been observed (see discussion in Hedges et al., 1999). The vertical POM profile produced by the model can be fitted to the classic power function describing the oceanic carbon rate (C_R=Z^−α) using an attenuation coefficient α of 0.36 which is the value proposed for another anoxic environment (i.e. the Mexico Margin) by Devol and Hartnett (2001). Due to the lower efficiency of detritus degradation in anoxic conditions and to the aggregation of particles that enhanced the sinking, an important part of the export to the anoxic layer (i.e. 33%, 0.52 mol C m⁻² year⁻¹) escapes remineralization in the water column and reaches the sediments. Therefore, sediments are active sites of sulfide production contributing to 26% of the total sulfide production.In the upper layer, the oxygen dynamics is mainly governed by photosynthesis and respiration processes as well as by air-sea exchanges. ?71% of the oxygen produced by phytoplankton (photosynthesis+nitrate reduction) is lost through respiration, ?21% by outgasing to the atmosphere, ?5% through nitrification and only ?2% in the oxidation of reduced components (e.g. Mn²⁺, Fe²⁺, H₂S).The model estimates the amount of nitrogen lost through denitrification at 307 mmol N m⁻² year⁻¹ that can be partitioned into a loss of ?55% through the use of nitrate for the oxidation of detritus in low oxygen conditions, ?40% in the ANAMMOX process and the remaining ?5% in the oxidation of reduced substances by nitrate.In agreement with data analysis performed on long time series collected since the 1960s (Konovalov and Murray, 2001), the sulfide and nitrogen budgets established for the anoxic layer are not balanced in response to the enhanced particle fluxes induced by eutrophication: the NH₄ and H₂S concentrations increase.     

Modeling carbon sequestration under zero-tillage at the regional scale. II. The influence of crop rotation and soil type

For policy decisions with respect to CO₂-mitigation measures in the agricultural sector, national and regional estimations of the efficiency of such measures are required. The conversion of ploughed cropland to zero-tillage is discussed as an option to reduce CO₂ emissions and promises at the same time effective soil and water conservation. Based on the upscaling of simulation results with the soil and land resources information system SLISYS-BW, estimations of CO₂-mitigation rates in relation to crop rotations and soil type have been made for the state of Baden-Württemberg (Germany). The results indicate considerable differences in the CO₂-mitigation rates between crop rotations ranging from 0.48 to 0.03 Mg C ha⁻¹ a⁻¹ for winter cereals–spring cereals–rape rotations and winter cereals–spring cereals–corn silage rotations, respectively. The efficiency of the crop rotations is strongly related to the total carbon input and in particular the amount of crop residues. Among the considered soil types, highest CO₂-mitigation rates are associated with Cumulic Anthrosols (0.62 Mg C ha⁻¹ a⁻¹) and the lowest with Gleysols (−0.01 Mg C ha⁻¹ a⁻¹). An agricultural extensification scenario with conventional plowing but conversion of the presently applied intensive crop rotations to a clover–clover–winter cereals rotation indicated a CO₂-mitigation potential of 466 Gg C a⁻¹. However, the present high market prices for cereals and increasing demand for energy production from biomass encourages an intensification of the agricultural production and an excessive removal of biomass which in future will seriously reduce the potential for carbon sequestration on cropland.     

Modelling the impact of thermal adaptation of soil microorganisms and crop system on the dynamics of organic matter in a tropical soil under a climate change scenario

No consensus currently exists about how climate change should affect the status of soil organic matter (SOM) in the tropics. In this study, we analyse the impact of climate change on the underlying mechanisms controlling SOM dynamics in a ferralsol under two contrasting tropical crops: maize (C4 plant) and banana (C3 plant). We model the effect of microbial thermal adaptation on carbon (C) mineralisation at the crop system scale and introduce it in the model STICS, which was previously calibrated for the soil-crop systems tested in this study. Microbial thermal adaptation modelling is based on a reported theory for thermal acclimation of plant and soil respiration. The climate is simulated from 1950 to 2099 for the tropical humid conditions of Guadeloupe (French Antilles), using the ARPEGE model and the IPCC emission scenario A1B. The model predicts increases of 3.4 °C for air temperature and 1100 mm yr⁻¹ for rainfall as a response to an increase of 375 ppm for atmospheric carbon dioxide concentration in the 2090-2099 decade compared with the 1950-1959 decade. The results of the STICS model indicate that the crop affects the response of SOM to climate change by controlling the change in several variables involved in C dynamics: C input, soil temperature and soil moisture. SOM content varies little until 2020, and then it decreases faster for maize than for banana. The decrease is weakened under the hypothesis of thermal adaptation, and this effect is greater for maize (−180 kg C ha⁻¹ yr⁻¹ without adaptation and −140 kg C ha⁻¹ yr⁻¹ with adaptation) than for banana (−60 kg C ha⁻¹ yr⁻¹ and −40 kg C ha⁻¹ yr⁻¹, respectively). The greater SOM loss in maize is mainly due to the negative effect of warming on maize growth decreasing C input from residues. Climate change has a small effect on banana growth, and SOM loss is linked to its effect on C mineralisation. For both crops, annual C mineralisation increases until 2040, and then it decreases continuously. Thermal adaptation reduces the initial increase in mineralisation, but its effect is lower on the final decrease, which is mainly controlled by substrate limitation. No stabilisation in SOM status is attained at the end of the analysed period because C mineralisation is always greater than C input. Model predictions indicate that microbial thermal adaptation modifies, but does not fundamentally change the temporal pattern of SOM dynamics. The vegetation type (C3 or C4) plays a major role in SOM dynamics in this tropical soil because of the different impact of climate change on crop growth and then on C inputs.     

The sensitivity of carbon sequestration to harvesting and climate conditions in a temperate cypress forest: Observations and modeling

The role of disturbance and climate factors in determining the forest carbon balance was investigated at a Japanese cypress forest in central Japan with eddy flux measurements, tree-ring analyses, and a terrestrial biosphere model. The forest was established as a plantation after intermittent harvesting and replanting between 1959 and 1977, and acted as a strong carbon sink of approximately 500 g C m⁻² year⁻¹ for the measurement years between 2001 and 2007. A terrestrial biosphere model, BIOME-BGC, was validated using the eddy flux data at daily to interannual timescales, and the tree-ring width data at interannual to decadal timescales. According to the model simulation, during the observation period 270 ± 55 g C m⁻² year⁻¹ was additionally sequestered due to the indirect effects of the harvesting and planting, whereas the increase of CO₂ concentration and the change in climate increased the sink of 110 ± 40 and 30 ± 80 g C m⁻² year⁻¹, respectively. The model simulation shows that the forest is now recovering from harvesting, and that harvesting is a more important determinant of the current carbon sink than either interannual climate anomalies or increased atmospheric CO₂ concentration. We found that harvesting with long rotation length could be effective management activity in order to increase carbon sequestration, if the harvested timber is converted into products with long lifecycles.     

Dynamics of a chemostat with three competitive hydrogen oxidizing denitrifying microbial populations and their efficiency for denitrification

In this work, competition for two nitrogen resources (nitrate-, nitrite-nitrogen) between three hydrogen oxidizing denitrifying populations (Acidovorax sp. strain Ic3 (X₁), Paracoccus sp. strain Ic1 (X₂), and Acinetobacter sp. strain Ic2 (X₃)) was examined. The dynamics of three systems of microbial populations (system I: X₁ − X₃, system II: X₂ − X₃, and system III: X₁ − X₂ − X₃), grown in a chemostat, was studied using bifurcation analysis. The chemostat is the most common type of biological reactor used for the study of microbial growth under controlled conditions. The effect of the operating parameters (i.e., dilution rate and feed nitrate nitrogen concentration) on the long-term behavior of the systems showed that X₃ was the predominant population for a wide range of combinations of dilution rate and feed nitrate nitrogen concentration. Also, coexistence of two populations (X₂X₃, X₁X₃) was observed. The results of the bifurcation analysis were also used to determine the denitrification rate and the nitrite nitrogen accumulation for each of the three systems as a function of the dilution rate (up to 0.17 h⁻¹) and the feed nitrate nitrogen concentration (up to 300 mg/L). The highest denitrification rate was achieved by system I (28 mg/Lh). A comparison between the three systems showed that the nitrite nitrogen concentration in system I was less than the one in system III, while the two systems gave similar denitrification rates. The second system had the greatest accumulation of nitrites with the lowest denitrification rate.     

Local calibration of coliforms parameters of water quality problem at Igapó I Lake, Londrina, Paraná, Brazil

This article presents results concerning the local calibration of the transport parameters (longitudinal and transversal diffusions and decay coefficient) for a two-dimensional problem of water quality at Igapó I Lake, located in Londrina, Paraná, Brazil, using fecal coliforms as an indicator of water quality. The simulation of fecal coliforms concentrations all over the water body is conducted by means of a structured discretization of the geometry of Igapó I Lake, together with the finite difference and finite element methods. By using the velocity field, modeled by the Navier-Stokes and Poisson equations, the flow of fecal coliforms is described by means of a transport model, which considers advective and diffusive processes, as well as a process of fecal coliforms decay. In the checkpoint, the longitudinal and transversal diffusion coefficients and the coliforms decay coefficient that best fitted the value of the fecal coliforms concentration were D_x = D_y = 0.001 m²/h and k = 0.5 d⁻¹ = 0.02083 h⁻¹. A qualitative and quantitative analysis of the numerical simulations conducted in function of the diffusion coefficients and of the coliforms decay parameter provided a better understanding of the local water quality at Igapó I Lake.     

Estimation of annual potential evapotranspiration at regional scale based on the effect of moisture on soil respiration

Potential evapotranspiration (PET) is an important component of water cycle. For traditional models derived from the principle of aerodynamics and the surface energy balance, its calculation always includes many parameters, such as net radiation, water vapor pressure, air temperature and wind speed. We found that it can be acquired in an easier way in specific regions. In this study, a new PET model (PETP model) derived from two empirical models of soil respiration was evaluated using the Penman-Monteith equation as a standard method. The results indicate that the PETP model estimation concur with the Penman-Monteith equation in sites where annual precipitation ranges from 717.71 mm to 1727.37 mm (R² = 0.68, p = 0.0002), but show large discrepancies in all sites (R² = 0.07, p = 0.1280). Then we applied our PETP model at the global scale to the regions with precipitation higher than 700 mm using 2.5° CMAP data to obtain the annual PET for 2006. As expected, the spatial pattern is satisfactory overall, with the highest PET values distributed in the lower latitudes or coastal regions, and with an average of 1292.60 ± 540.15 mm year⁻¹. This PETP model provides a convenient approach to estimate PET at regional scales.     

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.     

Simulation of thermodynamic transmission in green roof ecosystem

Green roofs entail the creation of vegetated space on the top of artificial structures. They can modify the thermal properties of buildings to bring cooling energy conservation and improve human comfort. This study evaluates the thermodynamic transmission in the green roof ecosystem under different vegetation treatments. Our model simulation is based on the traditional Bowen ratio energy balance model (BREBM) and a proposed solar radiation shield effectiveness model (SEM). The BREBM investigates energy absorption of different components of radiation, and the SEM evaluates the radiation shield effects. The proposed model is tested and validated to be efficient to simulate solar energy transmission in green roofs, with some major findings. Firstly, the solar radiation transmission processes might be considered as free vibration motion. Daytime positive heat storage of the green roof is 350-520 W·m⁻² on an hourly basis. Nighttime or afternoon negative value registers a rather constant magnitude of −60 W·m⁻². Daily net average is positive around 155-210 W·m⁻². Secondly, solar radiation vibration is highly correlated with plant structure. The canopy reflectance and transmittance are strongly correlated (R² = 0.87). The multi-layer shrub treatment has the highest shield effectiveness (0.34), followed by two-layer groundcover (0.27), and single-layer grass (0.16). Green roof vegetation absorbs and stores large amounts of heat to form an effective thermal buffer against daily temperature fluctuation. Vegetated roofs drastically depress air temperature in comparison with bare ground (control treatment). Finally, the thermodynamic model is relatively simple and efficient for investigating thermodynamic transmission in green roof ecosystem, and it could be developed into a broad solar radiant land cover model.     

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.     

Residual levels and health risk of polycyclic aromatic hydrocarbons in freshwater fishes from Lake Small Bai-Yang-Dian, Northern China

The residual levels of polycyclic aromatic hydrocarbons (PAHs) in the liver, brain, gill and muscle tissues of four common edible freshwater fish species including crucian carp, snakehead fish, grass carp and silver carp collected from Lake Small Bai-Yang-Dian in northern China were measured by GC-MS. The distribution and composition pattern of PAHs in the fish tissues, and the effects of lipid contents in fish tissues and the octanol-water partition coefficient (Kow) of PAHs congeners on them were analyzed. The human health risk of PAHs though fish consumptions was estimated. The following results were obtained: (1) The average residual levels of total PAHs (PAH₁₆) on wet weight base in the different tissues of each fish species ranged from 4.764 to 144.254 ng/g ww. The differences in the average residual levels on wet weight base for PAH₁₆ within four fish species were not statistically significant (P > 0.05); however, these within four fish tissues were statistically significant (P < 0.01). (2) There were very similar distribution patterns of PAH congeners among both the fish tissues and the fish species, as indicated by statistically significant positive interrelationships (R = 0.58-0.97, P < 0.01 or P < 0.05). Low molecular weight (LMW) PAHs predominated the distribution in the fish tissues, accounting for 89.97% of total PAHs. Phe was the most dominant component, according for 37.79% of total PAHs, followed by Ant (18.59%), Flo (12.59%), Nap (10.79%), Fla (9.82%) and Pyr (6.43%). (3) The PAHs residues and distribution in the fish tissues are dependent on both the Kow of PAH congeners and the lipid contents in the fish tissues. There was a significant positive relationship (R = 0.7116, P < 0.0001) between lipid contents and PAHs residual levels. The statistically significant negative relationships (P < 0.05) were found between LogKow and log-transformed PAHs contents on wet weight base for all fish tissues except for the muscle tissue of snakehead fish, the brain and liver tissues of crucian carp. (4) The risk levels of total PAHs were lower than 10⁻⁵ for the muscle tissues of four studied fish species and for the brain tissues of grass carp and snakehead fish; while these were higher than 10⁻⁵ for the brain tissues of crucian carp and silver carp. The risk levels of total PAHs in the liver tissues of four studied fish species except for snakehead fish exceeded 10⁻⁵ for 2-4.5 times. However, the potency equivalent concentration (PEC) of total PAHs in four studied fish tissues were still lower than the maximum permissible BaP limits for crops and baked meat and for plants in the national criterions. The distributions of PAH congeners in fish were well simulated by a level III fugacity model, especially for low molecule weight PAHs.     

Foliage profiles of individual trees determine competition,self-thinning,biomass and NPP of a Cryptomeria japonica forest stand: A simulation study based on a stand-scale process-based forest model

A simulation study was carried out to investigate simultaneously the effects of eco-physiological parameters on competitive asymmetry, self-thinning, stand biomass and NPP in a temperate forest using an atmosphere–vegetation dynamics interactive model (MINoSGI). In this study, we selected three eco-physiological relevant parameters as foliage profiles (i.e. vertical distribution of leaf area density) of individual trees (distribution pattern is described by the parameter η), biomass allocation pattern in individual tree growth (χ) and the maximum carboxylation velocity (V_max). The position of the maximal leaf area density shifts upward in the canopy with increasing η. For scenarios with η < 4 (foliage concentrated in the lowest canopy layer) or η > 12 (foliage concentrated in the uppermost canopy layer), a low degree of competitive asymmetry was produced. These scenarios resulted in the survival of subordinate trees due to a brighter lower canopy environment when η < 4 or the generation of spatially separated foliage profiles between dominant and subordinate trees when η > 12. In contrast, competition between trees was most asymmetric when 4 ≤ η ≤ 12 (vertically widespread foliage profile in the canopy), especially when η = 8. In such cases, vertically widespread foliage of dominant trees lowered the opportunity of light acquisition for subordinate trees and reduced their carbon gain. The resulting reduction in carbon gain of subordinate trees yielded a higher degree of competitive asymmetry and ultimately higher mortality of subordinate trees. It was also shown that 4 ≤ η ≤ 12 generated higher self-thinning speed, smaller accumulated NPP, litter-fall and potential stand biomass as compared with the scenarios with η < 4 or η > 12. In contrast, our simulation revealed small effects of χ or V_max on the above-mentioned variables as compared with those of η. In particular, it is notable that greater V_max would not produce greater potential stand biomass and accumulated NPP although it has been thought that physiological parameters relevant to photosynthesis such as V_max influence dynamic changes in forest stand biomass and NPP (e.g. the greater the V_max, the greater the NPP). Overall, it is suggested that foliage profiles rather than biomass allocation or maximum carboxylation velocity greatly govern forest dynamics, stand biomass, NPP and litter-fall.     

GIS-based modeling of the geographic distribution of Quercus emoryi Torr. (Fagaceae) in México and identification of significant environmental factors influencing the species’ distribution

The greatest concentration of oak species in the world is believed to be found in Mexico. These species are potentially useful for reforestation because of their capacity to adapt to diverse environments. Knowledge of their geographic distribution and of species–environment relations is essential for decision-making in the management and conservation of natural resources. The objectives of this study were to develop a model of the distribution of Quercus emoryi Torr. in Mexico, using geographic information systems and data layers of climatic and other variables, and to determine the variables that significantly influence the distribution of the species. The study consisted of the following steps: (A) selection of the target species from a botanical scientific collection, (B) characterization of the collecting sites using images with values or categories of the variables, (C) model building with the overlay of images that meet the habitat conditions determined from the characterization of sites, (D) model validation with independent data in order to determine the precision of the model, (E) model calibration through adjustment of the intervals of some variables, and (F) sensitivity analysis using precision and concordance non-parametric statistics applied to pairs of images. Results show that the intervals of the variables that best describe the species’ habitat are the following: altitude from 1650 to 2750 amsl, slope from 0 to 66°; average minimum temperature of January from −12 to −3 °C; mean temperature of June from 11 to 25 °C; mean annual precipitation from 218 to 1225 mm; soil units: lithosol, eutric cambisol, haplic phaeozem, chromic luvisol, rendzina, luvic xerosol, mollic planosol, pellic vertisol, eutric regosol; type of vegetation: oak forest, oak–pine forest, pine forest, pine–oak forest, juniperus forest, low open forest, natural grassland and chaparral. The resulting model of the geographic distribution of Quercus emoryi in Mexico had the following values for non-parametric statistics of precision and agreement: Kappa index of 0.613 and 0.788, overall accuracy of 0.806 and 0.894, sensitivity of 0.650 and 0.825, specificity of 0.963, positive predictive value of 0.945 and 0.957 and negative predictive value of 0.733 and 0.846. Results indicate that the variable average minimum temperature of January, with a maximum value of −3 °C, is an important factor in limiting the species’ distribution.