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.     

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.     

Identification of important factors for water vapor flux and CO2 exchange in a cropland

Water vapor flux and carbon dioxide (CO₂) exchange in croplands are crucial to water and carbon cycle research as well as to global warming evaluation. In this study, a standard three-layer feed-forward back propagation neural network technique associated with the Bayesian technique of automatic relevance determination (ARD) was employed to investigate water vapor and CO₂ exchange between the canopy of summer maize and atmosphere in responses to variations of environmental and physiological factors. These factors, namely the photosynthetically active radiation (PAR), air temperature (T), vapor pressure deficient (VPD), leaf-area index (LAI), soil water content in root zone (W), and friction velocity (U*), were used as inputs in neural network analysis. Results showed that PAR, VPD, T and LAI were the primary factors regulating both water vapor and CO₂ fluxes with VPD and W more critical to water vapor flux and PAR and T more crucial to CO₂ exchange. Furthermore, two time variables “day of the year (DOY)” and “time of the day (TOD)” could also improve the simulation results of neural network analysis. The important factors identified by the neural network technique used in this study were in the order of PAR > T > VPD > LAI > U* > TOD for water vapor flux and in the order of VPD > W > LAI > T > PAR > DOY for CO₂ exchange. This study suggests that neural network technique associated with ARD could be a useful tool for identifying important factors regulating water vapor and CO₂ fluxes in terrestrial ecosystem.     

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.     

Surface modelling of global terrestrial ecosystems under three climate change scenarios

A high accuracy and speed method (HASM) of surface modelling is developed to find a solution for error problem and to improve computation speed. A digital elevation model (DEM) is established on spatial resolution of 13.5 km × 13.5 km. Regression formulations among temperature, elevation and latitude are simulated in terms of data from 2766 weather observation stations scattered over the world by using the 13.5 km × 13.5 km DEM as auxiliary data. Three climate scenarios of HadCM3 are refined from spatial resolution of 405 km × 270 km to 13.5 km × 13.5 km in terms of the regression formulations. HASM is employed to simulate surfaces of mean annual bio-temperature, mean annual precipitation and potential evapotranspiration ratio during the periods from 1961 to 1990 (T₁), from 2010 to 2039 (T₂), from 2040 to 2069 (T₃), and from 2070 to 2099 (T₄) on spatial resolution of 13.5 km × 13.5 km. Three scenarios of terrestrial ecosystems on global level are finally developed on the basis of the simulated climate surfaces. The scenarios show that all polar/nival, subpolar/alpine and cold ecosystem types would continuously shrink and all tropical types, except tropical rain forest in scenario A1Fi, would expand because of the climate warming. Especially at least 80% of moist tundra and 22% of nival area might disappear in period T₄ comparing with the ones in the period T₁. Tropical thorn woodland might increase by more than 97%. Subpolar/alpine moist tundra would be the most sensitive ecosystem type because its area would have the rapidest decreasing rate and its mean center would shift the longest distance towards west. Subpolar/alpine moist tundra might be able to serve as an indicator of climatic change. In general, climate change would lead to a continuous reduction of ecological diversity.     

Modeling amphibian energetics, habitat suitability, and movements of western toads, Anaxyrus (=Bufo) boreas, across present and future landscapes

Effective conservation of amphibian populations requires the prediction of how amphibians use and move through a landscape. Amphibians are closely coupled to their physical environment. Thus an approach that uses the physiological attributes of amphibians, together with knowledge of their natural history, should be helpful. We used Niche Mapper™ to model the known movements and habitat use patterns of a population of Western toads (Anaxyrus (=Bufo) boreas) occupying forested habitats in southeastern Idaho. Niche Mapper uses first principles of environmental biophysics to combine features of topography, climate, land cover, and animal features to model microclimates and animal physiology and behavior across landscapes. Niche Mapper reproduced core body temperatures (T_c) and evaporation rates of live toads with average errors of 1.6 ± 0.4 °C and 0.8 ± 0.2 g/h, respectively. For four different habitat types, it reproduced similar mid-summer daily temperature patterns as those measured in the field and calculated evaporation rates (g/h) with an average error rate of 7.2 ± 5.5%. Sensitivity analyses indicate these errors do not significantly affect estimates of food consumption or activity. Using Niche Mapper we predicted the daily habitats used by free-ranging toads; our accuracy for female toads was greater than for male toads (74.2 ± 6.8% and 53.6 ± 15.8%, respectively), reflecting the stronger patterns of habitat selection among females. Using these changing to construct a cost surface, we also reconstructed movement paths that were consistent with field observations. The effect of climate warming on toads depends on the interaction of temperature and atmospheric moisture. If climate change occurs as predicted, results from Niche Mapper suggests that climate warming will increase the physiological cost of landscapes thereby limiting the activity for toads in different habitats.     

Exploring the influence of lake water chemistry on chlorophyll a: A multivariate statistical model analysis

A multivariate statistical approach integrating the absolute principal components score (APCS) and multivariate linear regression (APCS-MLR), along with structural equation modeling (SEM), was used to model the influence of water chemistry variables on chlorophyll a (Chl a) in Lake Qilu, a severely polluted lake in southwestern China. Water quality was surveyed monthly from 2000 to 2005. APCS-MLR was used to identify key water chemistry variables, mine data for SEM, and predict Chl a. Seven principal components (PCs) were determined as eigenvalues >1, which explained 68.67% of the original variance. Four PCs were selected to predict Chl a using APCS-MLR. The results showed a good fit between the observed data and modeled values, with R² = 0.80. For SEM, Chl a and eight variables were used: NH₄-N (ammonia-nitrogen), total phosphorus (TP), Secchi disc depth (SD), cyanide (CN), arsenic (As), cadmium (Cd), fluoride (F), and temperature (T). A conceptual model was established to describe the relationships among the water chemistry variables and Chl a. Four latent variables were also introduced: physical factors, nutrients, toxic substances, and phytoplankton. In general, the SEM demonstrated good agreement between the sample covariance matrix of observed variables and the model-implied covariance matrix. Among the water chemistry factors, T and TP had the greatest positive influence on Chl a, whereas SD had the largest negative influence. These results will help researchers and decision-makers to better understand the influence of water chemistry on phytoplankton and to manage eutrophication adaptively in Lake Qilu.     

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.     

The MONICA model: Testing predictability for crop growth, soil moisture and nitrogen dynamics

A fundamentally revised version of the HERMES agro-ecosystem model, released under the name of MONICA, was calibrated and tested to predict crop growth, soil moisture and nitrogen dynamics for various experimental crop rotations across Germany, including major cereals, sugar beet and maize. The calibration procedure also included crops grown experimentally under elevated atmospheric CO₂ concentration. The calibrated MONICA simulations yielded a median normalised mean absolute error (nMAE) of 0.20 across all observed target variables (n = 42) and a median Willmott's Index of Agreement (d) of 0.91 (median modelling efficiency (ME): 0.75). Although the crop biomass, habitus and soil moisture variables were all within an acceptable range, the model often underperformed for variables related to nitrogen. Uncalibrated MONICA simulations yielded a median nMAE of 0.27 across all observed target variables (n = 85) and a median d of 0.76 (median ME: 0.30), also showing predominantly acceptable results for the crop biomass, habitus and soil moisture variables. Based on the convincing performance of the model under uncalibrated conditions, MONICA can be regarded as a suitable simulation model for use in regional applications. Furthermore, its ability to reproduce the observed crop growth results in free-air carbon enrichment experiments makes it suited to predict agro-ecosystem behaviour under expected future climate conditions.     

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.     

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.     

Modelling and mapping topographic variations in forest soils at high resolution: A case study

This article examines the utility of a digitally derived cartographic depth-to-water (DTW) index to model and map variations in drainage, vegetation and soil type and select soil properties within a forested area (40 ha) of the Swan Hills, Alberta, Canada. This index was derived from a LiDAR (Light Detection and Ranging) derived digital elevation model (DEM), with at least 1 ground return per m². The resulting DTW pattern was set to be zero along all DEM-derived flow channels, each with a 4 ha flow-initiation threshold. Soil type (luvisol, gleysol, mesisol), drainage type (very poor to excessive), vegetation type (hydric to xeric) and forest floor depth were determined along hillslope transects. These determinations conformed more closely to the DEM-derived log₁₀(DTW) variations (R² > 60%) than to the corresponding variations of the widely used topographic wetness index (TWI) (R² < 25%). Setting log₁₀(DTW) thresholds to represent the wet to moist to dry transitions between vegetation, drainage and soil type enabled a high-resolution mapping of these types across the study area. Also determined were soil moisture content, coarse fragment and soil particle composition (sand, silt, clay), pH, total C, N, S, P, Ca, Mg, K, Fe, Al, Mn, Zn, and available Ca, Mg, K, P and NH₄, by soil layer type and depth. Most of these variables were also more correlated with log₁₀(DTW) than with TWI, with and without soil layer depth as an additional regression variable. These variables are, therefore, subject to topographic controls to at least some extent, and can be modelled and mapped accordingly, as illustrated for soil moisture, forest floor depth and pH across the study area, from ridge tops to depressions. Further examinations revealed that the DEM-produced DTW and TWI patterns complemented one another, with DTW delineating soils in relation to local water-table influences, and with TWI delineating where the water would flow and accumulate.     

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.     

Structural changes in lake functioning induced from nutrient loading and climate variability

Climate variability is increasingly recognized as an important regulatory factor, capable of influencing the structural properties of aquatic ecosystems. Lakes appear to be particularly sensitive to the ecological impacts of climate variability, and several long time series have shown a close coupling between climate, lake thermal properties and individual organism physiology, population abundance, community structure, and food web dynamics. Thus, understanding the complex interplay among meteorological forcing, hydrological variability, and ecosystem functioning is essential for improving the credibility of model-based water resources/fisheries management. Our objective herein is to examine the relative importance of the ecological mechanisms underlying plankton seasonal variability in Lake Washington, Washington State (USA), over a 35-year period (1964–1998). Our analysis is founded upon an intermediate complexity plankton model that is used to reproduce the limiting nutrient (phosphate)–phytoplankton–zooplankton–detritus (particulate phosphorus) dynamics in the lake. Model parameterization is based on a Bayesian calibration scheme that offers insights into the degree of information the data contain about model inputs and allows obtaining predictions along with uncertainty bounds for modeled output variables. The model accurately reproduces the key seasonal planktonic patterns in Lake Washington and provides realistic estimates of predictive uncertainty for water quality variables of environmental management interest. A principal component analysis of the annual estimates of the underlying ecological processes highlighted the significant role of the phosphorus recycling stemming from the zooplankton excretion on the planktonic food web variability. We also identified a moderately significant signature of the local climatic conditions (air temperature) on phytoplankton growth (r = 0.41), herbivorous grazing (r = 0.38), and detritus mineralization (r = 0.39). Our study seeks linkages with the conceptual food web model proposed by Hampton et al. [Hampton, S.E., Scheuerell, M.D., Schindler, D.E., 2006b. Coalescence in the Lake Washington story: interaction strengths in a planktonic food web. Limnol. Oceanogr. 51, 2042–2051.] to emphasize the “bottom-up” control of the Lake Washington plankton phenology. The posterior predictive distributions of the plankton model are also used to assess the exceedance frequency and confidence of compliance with total phosphorus (15 μg L⁻¹) and chlorophyll a (4 μg L⁻¹) threshold levels during the summer-stratified period in Lake Washington. Finally, we conclude by underscoring the importance of explicitly acknowledging the uncertainty in ecological forecasts to the management of freshwater ecosystems under a changing global environment.     

Evaluating weather effects on interannual variation in net ecosystem productivity of a coastal temperate forest landscape: A model intercomparison

Forest productivity is strongly affected by seasonal weather patterns and by natural or anthropogenic disturbances. However weather effects on forest productivity are not currently represented in inventory-based models such as CBM-CFS3 used in national forest C accounting programs. To evaluate different approaches to modelling these effects, a model intercomparison was conducted among CBM-CFS3 and four process models (ecosys, CN-CLASS, Can-IBIS and 3PG) over a 2500 ha landscape in the Oyster River (OR) area of British Columbia, Canada. The process models used local weather data to simulate net primary productivity (NPP), net ecosystem productivity (NEP) and net biome productivity (NBP) from 1920 to 2005. Other inputs used by the process and inventory models were generated from soil, land cover and disturbance records. During a period of intense disturbance from 1928 to 1943, simulated NBP diverged considerably among the models. This divergence was attributed to differences among models in the sizes of detrital and humus C stocks in different soil layers to which a uniform set of soil C transformation coefficients was applied during disturbances. After the disturbance period, divergence in modelled NBP among models was much smaller, and attributed mainly to differences in simulated NPP caused by different approaches to modelling weather effects on productivity. In spite of these differences, age-detrended variation in annual NPP and NEP of closed canopy forest stands was negatively correlated with mean daily maximum air temperature during July-September (T_amax) in all process models (R² = 0.4-0.6), indicating that these correlations were robust. The negative correlation between T_amax and NEP was attributed to different processes in different models, which were tested by comparing CO₂ fluxes from these models with those measured by eddy covariance (EC) under contrasting air temperatures (T_a). The general agreement in sensitivity of annual NPP to T_amax among the process models led to the development of a generalized algorithm for weather effects on NPP of coastal temperate coniferous forests for use in inventory-based models such as CBM-CFS3: NPP′ = NPP − 57.1 (T_amax − 18.6), where NPP and NPP′ are the current and temperature-adjusted annual NPP estimates from the inventory-based model, 18.6 is the long-term mean daily maximum air temperature during July-September, and T_amax is the mean value for the current year. Our analysis indicated that the sensitivity of NPP to T_amax was nonlinear, so that this algorithm should not be extrapolated beyond the conditions of this study. However the process-based methodology to estimate weather effects on NPP and NEP developed in this study is widely applicable to other forest types and may be adopted for other inventory based forest carbon cycle models.     

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 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.     

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.     

The derivation of species response curves with Gaussian logistic regression is sensitive to sampling intensity and curve characteristics

We investigated quantitatively the sensitivity of plant species response curves to sampling characteristics (number of plots, occurrence and frequency of species), along a simulated pH gradient. We defined 54 theoretical unimodal response curves, issued from combinations of six values for optimum (opt = 3, 4, …, 8), three values for tolerance (tol = 0.5, 1.0, and 1.5, sensu ter Braak and Looman [ter Braak, C.J.F., Looman, C.W.N., 1986. Weighted averaging, logistic regression and the Gaussian response model. Vegetatio 65, 3–11]), and three values for maximum probability of presence (p_max = 0.05, 0.20, and 0.50). For each of these 54 theoretical response curves, we built artificial binary data sets (presence/absence) to test the influence of species occurrence, frequency, or number of available plots. With real data extracted from EcoPlant, a phytoecological database for French forests [Gégout, J.-C., Coudun, Ch., Bailly, G., Jabiol, B., 2005. EcoPlant: a forest sites database linking floristic data with soil characteristics and climatic conditions. J. Veg. Sci. 16, 257–260], we compared the ecological response of 50 plant species to soil pH, based first on a small data set (100 randomly sampled plots), and then based on the whole data set available (3810 plots).     

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.