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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Modeling estuarine-shelf exchanges in a deltaic estuary: Implications for coastal carbon budgets and hypoxia

The export of wetland-derived materials to the coastal ocean (i.e., the “Outwelling” hypothesis) has received considerable attention over the past several decades. While a number of studies have shown that estuaries export appreciable amounts of nutrients and carbon, few studies have attempted to estimate the importance of estuarine sources for the coastal carbon budgets in river-dominated coastal ecosystems. A novel tidal prism model was developed to examine estuarine-shelf exchanges in the Barataria estuary, a deltaic estuary located in the north-central Gulf of Mexico. This estuary has been the site of a massive wetland loss, and it has been hypothesized that carbon export from the eroding coastal wetlands supports the development of a large hypoxic zone in the coastal Gulf of Mexico. The model results show that the Barataria estuary receives nitrogen through the tidal passes and releases carbon to the coastal ocean. The mean calculated tidal water discharge of 6930 m³ s⁻¹ is equivalent to about 43% of the lower Mississippi River discharge. The annual total organic carbon (TOC) export is 109 million kg, or 57 gC m² yr⁻¹ when prorated to the total water area of the estuary. This carbon export is equivalent to a loss of 0.5 m of wetland soil horizon over an area of 8.4 km², and accounts for about 34% of the observed annual wetland loss in the estuary between 1978 and 2000. Compared to the lower Mississippi River, the Barataria estuary appears to be a very small source of TOC for the northern Gulf of Mexico (2.7% of riverine TOC), and is unlikely to have a significant influence on the development of the Gulf's hypoxia.     

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

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

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.     

Modelling forest management within a global vegetation model—Part 1: Model structure and general behaviour

This article describes a new forest management module (FMM) that explicitly simulates forest stand growth and management within a process-based global vegetation model (GVM) called ORCHIDEE. The net primary productivity simulated by ORCHIDEE is used as an input to the FMM. The FMM then calculates stand and management characteristics such as stand density, tree size distribution, tree growth, the timing and intensity of thinnings and clear-cuts, wood extraction and litter generated after thinning. Some of these variables are then fed back to ORCHIDEE. These computations are made possible with a distribution-based modelling of individual tree size. The model derives natural mortality from the relative density index (rdi), a competition index based on tree size and stand density. Based on the common forestry management principle of avoiding natural mortality, a set of rules is defined to calculate the recurrent intensity and frequency of forestry operations during the stand lifetime. The new-coupled model is called ORCHIDEE-FM (forest management).The general behaviour of ORCHIDEE-FM is analysed for a broadleaf forest in north-eastern France. Flux simulation throughout a forest rotation compare well with the literature values, both in absolute values and dynamics.Results from ORCHIDEE-FM highlight the impact of forest management on ecosystem C-cycling, both in terms of carbon fluxes and stocks. In particular, the average net ecosystem productivity (NEP) of 225 gC m⁻² year⁻¹ is close to the biome average of 311 gC m⁻² year⁻¹. The NEP of the “unmanaged” case is 40% lower, leading us to conclude that management explains 40% of the cumulated carbon sink over 150 years. A sensitivity analysis reveals 4 major avenues for improvement: a better determination of initial conditions, an improved allocation scheme to explain age-related decline in productivity, and an increased specificity of both the self-thinning curve and the biomass-diameter allometry.     

Restoration of ecosystems for biodiversity and carbon sequestration: Simulating growth dynamics of brigalow vegetation communities in Australia

Restoration of abandoned and degraded ecosystems through enhanced management of mature remnant patches and naturally regenerating (regrowth) forests is currently being used in the recovery of ecosystems for biodiversity protection and carbon sequestration. Knowledge of long-term dynamics of these ecosystems is often very limited. Vegetation models that examine long-term forest growth and succession of uneven aged, mixed-species forest ecosystems are integral to the planning and assessment of the recovery process of biodiversity values and biomass accumulation. This paper examined the use of the Ecosystem Dynamics Simulator (EDS) in projecting growth dynamics of mature remnant brigalow forest communities and recovery process of regrowth brigalow thickets. We used data from 188 long-term monitored plots of remnant and regrowth forests measured between 1963 and 2010. In this study the model was parameterised for 34 tree and shrub species and tested with independent long-term measurements. The model closely approximated actual development trajectories of mature forests and regrowth thickets but some inaccuracies in estimating regeneration through asexual reproduction and mortality were noted as reflected in stem density projections of remnant plots that had a mean of absolute relative bias of 46.2 (±12.4)%. Changes in species composition in remnant forests were projected with a 10% error. Basal area values observed in all remnant plots ranged from 6 to 29 m² ha⁻¹ and EDS projections between 1966 and 2005 (39 years) were 68.2 (±10.9)% of the observed basal area. Projected live aboveground biomass of remnant plots had a mean of 93.5 (±5.9) t ha⁻¹ compared to a mean of 91.3 (±8.0) t ha⁻¹ observed in the plots. In regrowth thicket, the model produced satisfactory projections of tree density (91%), basal area (89%), height (87%) and aboveground biomass (84%) compared to the observed attributes. Basal area and biomass accumulation in 45-year-old regrowth plots was approximately similar to that in remnant forests but recovery of woody understorey was very slow. The model projected that it would take 95 years for the regrowth to thin down to similar densities observed in original or remnant brigalow forests. These results indicated that EDS can produce relatively accurate projections of growth dynamics of brigalow regrowth forests necessary for informing restoration planning and projecting biomass accumulation.     

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.     

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 long-term hydrological impacts of regional urbanisation in Hanyang,China, using a GIS model and remote sensing

China has experienced a rapid urbanisation, especially since the 1980s; however, the environmental impacts of this process are not fully investigated. Hanyang (Hubei Province, south China) was selected as a typical case to investigate runoff and non-point source (NPS) pollution impacts of urbanisation. A water quality simulation model (L-THIA) was applied to determine the long-term implications of different degrees of regional urbanisation impacts on NPS pollutants. Land-use patterns in 1987, 1998 and 2003 were analysed to evaluate the temporal variation of urbanisation, and the precipitation dataset from 1975 to 2003 was used to estimate the mean annual runoff and NPS pollutants. The contributions of different land-use categories to average annual runoff and NPS pollutant production were assessed by the means of a regression model. Results show that urban/impervious lands increased by 30.4% between 1987 and 2003, with the most increase occurring in 1998–2003. Industrial and forestlands have the most and least impact, respectively, on mean annual runoff and NPS pollutants. A combination of L-THIA and the regression model was found to be useful as a decision support tool for regional and urban planning from the perspectives of water quality control.     

Economic and environmental performance of electricity production in Finland: A multicriteria assessment framework

Meeting environmental, economic, and societal targets in energy policy is complex and requires a multicriteria assessment framework capable of exploring trade-offs among alternative energy options. In this study, we integrated economic analysis and biophysical accounting methods to investigate the performance of electricity production in Finland at plant and national level. Economic and environmental costs of electricity generation technologies were assessed by evaluating economic features (direct monetary production cost), direct and indirect use of fossil fuels (GER cost), environmental impact (CO₂ emissions), and global environmental support (emergy cost). Three scenarios for Finland's energy future in 2025 and 2050 were also drawn and compared with the reference year 2008. Accounting for an emission permit of 25 €/t CO₂, the production costs calculated for CHP, gas, coal, and peat power plants resulted in 42, 67, 68, and 74 €/MWh, respectively. For wind and nuclear power a production cost of 63 and 35 €/MWh were calculated. The sensitivity analysis confirmed wind power's competitiveness when the price of emission permits overcomes 20 €/t CO₂. Hydro, wind, and nuclear power were characterized by a minor dependence on fossil fuels, showing a GER cost of 0.04, 0.13, and 0.26 J/J_e, and a value of direct and indirect CO₂ emissions of 0.01, 0.04, and 0.07 t CO₂/MWh. Instead, peat, coal, gas, and CHP plants showed a GER cost of 4.18, 4.00, 2.78, and 2.33 J/J_e. At national level, a major economic and environmental load was given by CHP and nuclear power while hydro power showed a minor load in spite of its large production. The scenario analysis raised technological and environmental concerns due to the massive increase of nuclear power and wood biomass exploitation. In conclusion, we addressed the need to further develop an energy policy for Finland's energy future based on a diversified energy mix oriented to the sustainable exploitation of local, renewable, and environmentally friendly energy sources.     

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

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