Scenarios of major terrestrial ecosystems in China

The spatial pattern and mean-center shift of major terrestrial ecosystems, termed Holdridge Life Zones (HLZ), during the periods from 1961 to 1990 (T1), from 2010 to 2039 (T2), from 2040 to 2069 (T3) and from 2070 to 2099 (T4) were analyzed by combining the zonal patterns of climatic change in China and the climatic change scenarios of HadCM2 and HadCM3. The results showed that nival area would decrease rapidly with temperature increase in the future. HadCM2 and HadCM3 predicted that the nival areas might disappear in 552 years and 204 years, respectively. Using both HadCM2 and HadCM3, the five HLZ types with the largest areal extent are nival zone, cool temperate moist forest, warm temperate moist forest, subtropical moist forest and boreal wet forest, which collectively account for more than 50% of China's land mass. Among these five HLZ types, nival zone, warm temperate moist forest and boreal wet forest would decrease continuously, whereas subtropical moist forest and cool temperate forest would increase continuously during the four periods. HLZ diversity and patch connectivity would increase continuously in the 21st century. The shift distances of mean centers of HLZ types simulated using HadCM3 were markedly greater than those simulated using HadCM2, in general. The results from both HadCM2 and HadCM3 showed that boreal wet forest, subtropical moist forest, tropical dry forest, warm temperate moist forest and subtropical wet forest had bigger shift ranges, indicating that these HLZ types are more sensitive to the climatic change scenarios of HadCM2 and HadCM3.     

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.     

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.     

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.     

Impact of climate changes on pollution levels in Denmark

The impact of climate changes on the pollution levels in Denmark is the major topic of this paper. Variations of the Danish air pollution levels that are caused by climatic changes are studied together with variations caused by other factors (emissions, inter-annual variability of meteorological conditions, etc.). The Unified Danish Eulerian Model (UNI-DEM) was run on a fine, 10 km × 10 km, grid over a space domain covering all of Europe to minimize the influence of the boundary conditions on the Danish pollution levels. This study is based on four categories of scenarios: (i) traditional scenarios, (ii) climatic scenarios, (iii) scenarios with variations of the human-made (anthropogenic) emissions and (iv) scenarios in which the biogenic emissions were varied. The total number of applied scenarios was 14, and a time-period of 16 years was used. The results show clearly that although the concentrations of the major pollutants do not depend too much on the climatic changes, some quantities, in particular quantities related to high ozone levels, might be increased significantly as a result of the warming trends in the future climate. The reason for this phenomenon is explained.     

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.     

Modeling the carbon dynamics of the La Grande hydroelectric complex in northern Quebec

A mechanistic semi-empirical carbon cycle model of the La Grande reservoir complex in northern Quebec, Canada was conceived in order to investigate the climate impact of such a large alteration of the continental water cycle. The model includes inputs from the drainage basin, organic matter release from flooded soils, CO₂ emissions across the water-atmosphere interface and sedimentation. Most input data stems from previous research by our group on those ecosystems. The model includes the seven reservoirs of the La Grande complex and was run for periods of 50 and 100 years. Terrigeneous dissolved, particulate and suspended soil carbon fluxes and concentrations were computed. Over 100 years, 31.3 × 10¹² g C are released from flooded soils, equivalent to 28-29% of inputs from the drainage basin. 40-74% of dissolved organic carbon is mineralized. CO₂ fluxes over 100 years are 50.5-79.8 × 10¹² g C, 46.4-67.9 × 10¹² g C more than in the absence of reservoirs. The increase in mineralization of organic matter and in CO₂ emissions is a result of the increase in cumulated water residence time due to the creation of the reservoirs. Changes in other carbon sinks and sources likely offset a part of this additional carbon flux to the atmosphere. In the first years following flooding of the reservoir, organic carbon release from flooded soils exceeds CO₂ emissions, implying the downstream export of large quantities of eroded soil organic carbon. After this initial period, CO₂ emissions are fuelled by organic carbon originating from the drainage basin.     

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.     

Effects of climatic data low-pass filtering on the ICBM temperature- and moisture-based soil biological activity factors in a cool and humid climate

The air temperature (T_air), total precipitation (TP) and potential evapotranspiration (PET) are standard input data for soil carbon dynamic models, i.e., for calculating temperature and moisture effects on soil biological activity. The resolution needed depends on objectives, the complexity of models and inbuilt pedotransfer functions. The Introductory Carbon Balance Model (ICBM) soil climate front end model calculates a multiplicative soil-temperature (r_{e_temp}) and -moisture (r_{e_wat}) factor with a daily time-step to estimate soil biological activity, i.e., r_{e_crop} = r_{e_temp} × r_{e_wat}. Our objective was to determine how much r_{e_temp}, r_{e_wat} and r_{e_crop} are affected by low-pass filtering of the climatic input data for a cool, humid temperate region. To achieve this we conducted spectral analysis on T_air, TP, PET and r_{e_crop} in the frequency domain. Thereafter we applied Fourier low-pass filters of 5, 15, 30, 60 and 180 days on T_air, TP, PET and tracked their effects through the soil climate model's state variables and outputs. This was done using a sandy and a heavy clay soil and an 89-year daily time-series from a meteorological station in Quebec (Canada). The Fourier spectra showed that the variance for T_air, PET and r_{e_crop} was dominated by an annual cycle, as could be expected. There was no yearly cycle for TP. The variation in r_{e_temp} explained most of the variance in r_{e_crop}. The soil climate module outputs were not sensitive to low-pass filtering of PET. A daily time-step was needed to avoid overestimating r_{e_crop} for the sandy soil. Using a weekly time-step for TP and T_air allowed us to explain about 80% of the variance in r_{e_crop} for the heavy clay soil. This study also indicates that the standard leaf (and green) area index functions for calculating transpiration should receive more attention, since they have significant effects on the state and output variables.     

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.     

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.     

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.     

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.     

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.     

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.     

Assessing spawning ground hydraulic suitability for Chinese sturgeon (Acipenser sinensis) from horizontal mean vorticity in Yangtze River

The Chinese sturgeon, Acipenser sinensis, is an anadromous protected species. Its migration pattern has been blocked since the construction of Gezhouba dam and the length of the natural spawning site reduced to less than 7 km. However, the fish eventually established an alternative spawning ground in the narrow downstream reach of Gezhouba dam. In this article, we applied Delft3D-Flow model to simulate the hydraulic suitability of the spawning ground downstream Gezhouba dam. Horizontal mean vorticity was used to assess the hydraulic environment of spawning ground. The flow field state was determined through model simulation and field-measured data used to validate the model. The computational method was improved by calculating absolute horizontal mean vorticity from estimates the literature. The final vorticity was determined from the simulation output and its distribution pattern retrieved. The horizontal mean vorticity range was 0.71–4.61 10⁻³ s⁻¹ for the entire spawning grounds, with egg mass field upper limit of 1.0 × 10⁻³ s⁻¹. Vorticity strength selection of Chinese sturgeon spawning can enhance our understanding of egg fertilization rate, hence the protection of fertilized eggs. Furthermore, the results of the study would add to existing scientific database for spawning ground hydraulic environment protection, especially in the ecological regulation drive of the Three Gorges Reservoir.     

An optimized policy for the reduction of CO2 emission in the Brazilian Legal Amazon

The Brazilian government has already acknowledged the importance of investing in the development and application of technologies to reduce or prevent CO₂ emissions resulting from human activities in the Legal Brazilian Amazon (BA). The BA corresponds to a total area of 5 × 10⁶ km² from which 4 × 10⁶ km² was originally covered by the rain forest. One way to interfere with the net balance of greenhouse gases (GHG) emissions is to increase the forest area to sequester CO₂ from the atmosphere. The single most important cause of depletion of the rain forest is cattle ranching. In this work, we present an effective policy to reduce the net balance of CO₂ emissions using optimal control theory to obtain a compromising partition of investments in reforestation and promotion of clear technology to achieve a CO₂ emission target for 2020. The simulation indicates that a CO₂ emission target for 2020 of 376 million tonnes requires an estimated forest area by 2020 of 3,708,000 km², demanding a reforestation of 454,037 km². Even though the regional economic growth can foster the necessary political environment for the commitment with optimal emission targets, the reduction of 38.9% of carbon emissions until 2020 proposed by Brazilian government seems too ambitious.     

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.     

An inventory of tree species in Europe—An essential data input for air pollution modelling

We present a detailed tree species inventory covering Europe, parts of Africa and parts of Asia. The inventory contains 39 groups of species that are important for biogenic VOCs or pollen emission calculations. For example: oak (Quercus), poplar (Populus), pines (Pinus), spruce (Picea), birch (Betula) and alder (Alnus). The inventory is based on national forest inventories and national statistics and gives tree species distribution in percentage within broadleaved as well as conifer forests. The inventory includes data from 799 regions and is redistributed into the 50 km × 50 km EMEP grid. The inventory is therefore prepared for easy implementation into atmospheric transport models by providing an extension to already applied land use data such as the Corine Land Cover (CLC2000) or Global Land Cover (GLC2000). The gridded version of the data set will be available on the webpage http://www.dmu.dk/International/Air/Models/Background/Trees/.