Individual-based modeling of flooding and salinity effects on a coastal swamp forest

Coastal swamps are among the rapidly vanishing wetland habitats in Louisiana. Increased flooding, nutrient and sediment deprivation, and salt-water intrusion have been implicated as probable causes of the decline of coastal swamps. We developed a two-species individual-based forest succession model to compare the growth and composition of a cypress-tupelo swamp under various combinations of flooding intensity and salinity levels, using historical time-series of stage and salinity data as inputs. Our model simulates forest succession over 500 years by representing the growth, mortality, and reproduction of individual Taxodium distichum (baldcypress) and Nyssa aquatica (water tupelo) trees in a 1-km² spatial grid of 10 m × 10 m cells that vary in water levels and salinity through differences in elevation. We independently adjusted the elevations of each cell to obtain different grid-wide mean elevations and standard deviations of elevation; this affected the temporal and spatial pattern of flooding. We calibrated the model by adjusting selected parameters until averaged basal area, stem density and wood production rates under two different mean elevations (partially versus highly flooded) were qualitatively similar to comparable values reported for swamps in the literature. Corroboration involved comparing model predictions to four well-monitored contrasting habitat sites within the Maurepas Basin, Louisiana, USA. Model predictions of both species combined showed the same patterns among sites as the data, but the model overestimated wood production and the dominance of T. distichum. Exploratory simulations predicted that increased flooding leads to swamps with reduced basal areas and stem densities, while increased salinity resulted in lower basal areas at low salinity concentration (∼1-3 psu) and complete tree mortality at higher salinity concentrations (∼2-6 psu). Our model can provide insight into the succession dynamics of coastal swamps and information for the effective design of restoration actions.     

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

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

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

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

Development of TRIPLEX-Management model for simulating the response of forest growth to pre-commercial thinning

In order to simulate forest growth response to pre-commercial thinning (PCT), TRIPLEX1.0 - a process-based model designed to predict forest growth as well as carbon (C) and nitrogen (N) dynamics - was modified and improved to also simulate managed forest ecosystem thinning practices. A three-parameter Weibull distribution model was integrated to simulate thinning treatments within the newly developed TRIPLEX-Management model. The thinning intensity component within the model allows users to simulate thinning treatments by applying basal area, stand density and volume to quantify thinning intensity. Natural mortality decreased following thinning due to an increase in growing space for residual stems. Predicted litterfall pools also increased after thinning events took place. The TRIPLEX-Management model was tested against published observational data for Jack Pine (Pinus banksiana Lamb.) stands subjected to PCT in Northwestern Ontario, Canada. The coefficients of determination (R²) between the predicted and observed variables including stand density, mean DBH (diameter at breast height), the quadratic mean DBH, total volume and merchantable volume as well as belowground, aboveground, and total biomass ranged from 0.50 to 0.88 (n = 20, P < 0.001) with the exception of mean tree height (R² = 0.25, n = 20, P < 0.05). Overall, the Willmott index of agreement between predicted and observed variables ranged from 0.97 to 1.00. Results show that the TRIPLEX-Management model is generally capable of simulating growth response to PCT for Jack Pine stands.     

Mortality of large trees and lianas following experimental drought in an Amazon forest

Severe drought episodes such as those associated with El Ni?o Southern Oscillation (ENSO) events influence large areas of tropical forest and may become more frequent in the future. One of the most important forest responses to severe drought is tree mortality, which alters forest structure, composition, carbon content, and flammability, and which varies widely. This study tests the hypothesis that tree mortality increases abruptly during drought episodes when plant-available soil water (PAW) declines below a critical minimum threshold. It also examines the effect of tree size, plant life form (palm, liana, tree) and potential canopy position (understory, midcanopy, overstory) on drought-induced plant mortality. A severe, four-year drought episode was simulated by excluding 60% of incoming throughfall during each wet season using plastic panels installed in the understory of a 1-ha forest treatment plot, while a 1-ha control plot received normal rainfall. After 3.2 years, the treatment resulted in a 38% increase in mortality rates across all stems >2 cm dbh. Mortality rates increased 4.5-fold among large trees (>30 cm dbh) and twofold among medium trees (10-30 cm dbh) in response to the treatment, whereas the smallest stems were less responsive. Recruitment rates did not compensate for the elevated mortality of larger-diameter stems in the treatment plot. Overall, lianas proved more susceptible to drought-induced mortality than trees or palms, and potential overstory tree species were more vulnerable than midcanopy and understory species. Large stems contributed to 90% of the pretreatment live aboveground biomass in both plots. Large-tree mortality resulting from the treatment generated 3.4 times more dead biomass than the control plot. The dramatic mortality response suggests significant, adverse impacts on the global carbon cycle if climatic changes follow current trends.     

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.     

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.     

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.     

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.     

Impact of litter removal on soil characteristics under a Pinus caribaea stand

ABSTRACT

Forest productivity is dependent on soil quality, which in turn is related to litter; yet there is limited understanding of the relationships between litter and soil quality. The effects of litter removal treatment on tree growth and soil properties were examined in a Pinus caribaea stand with the aim of providing an understanding of consequences of litter removal on soil patterns. This knowledge is important for planning appropriate long-term forest management. Three pairs of 30 × 30 m² plots (each pair a control and treatment) were established in the center of a P. caribaea stand in April 2002. The controls were left undisturbed with an intact litter layer, while litter was removed monthly from treated plots. The diameter and height of the P. caribaea decreased and soil quality was degraded over the experimental period in litter-removed (treated) plots. Litter removal also produced a significant increase in soil bulk density at 0–20 and 20–40 cm depths. The capillary porosity, capillary moisture capacity, and natural water content in controls at 0–20 and 20–40 cm depths were significantly greater than treated plots. The non-capillary porosity in controls at the 0–20 cm depth was also significantly higher than treated plots. The organic carbon, total nitrogen (N), total potassium (K) and alkalized N content in each layer, and available P and exchangeable K at 0–20 and 40–60 cm depths in control plots were significantly greater than litter-removed plots. The numbers of bacteria at 0–20 and 40–60 cm depths and of fungi at 20–40 and 40–60 cm depths were higher in control than treated plots. The number of actinomyces and urease, catalase, and acid phosphatase activities in controls at each depth were significantly greater than litter-removed plots.     

Simulation modeling to understand how selective foraging by beaver can drive the structure and function of a willow community

Beaver–willow (Castor-Salix) communities are a unique and vital component of healthy wetlands throughout the Holarctic region. Beaver selectively forage willow to provide fresh food, stored winter food, and construction material. The effects of this complex foraging behavior on the structure and function of willow communities is poorly understood. Simulation modeling may help ecologists understand these complex interactions. In this study, a modified version of the SAVANNA ecosystem model was developed to better understand how beaver foraging affects the structure and function of a willow community in a simulated riparian ecosystem in Rocky Mountain National Park, Colorado (RMNP). The model represents willow in terms of plant and stem dynamics and beaver foraging in terms of the quantity and quality of stems cut to meet the energetic and life history requirements of beaver. Given a site where all stems were equally available, the model suggested a simulated beaver family of 2 adults, 2 yearlings, and 2 kits required a minimum of 4 ha of willow (containing about10 stems m⁻²) to persist in a steady-state condition. Beaver created a willow community where the annual net primary productivity (ANPP) was 2 times higher and plant architecture was more diverse than the willow community without beaver. Beaver foraging created a plant architecture dominated by medium size willow plants, which likely explains how beaver can increase ANPP. Long-term simulations suggested that woody biomass stabilized at similar values even though availability differed greatly at initial condition. Simulations also suggested that willow ANPP increased across a range of beaver densities until beaver became food limited. Thus, selective foraging by beaver increased productivity, decreased biomass, and increased structural heterogeneity in a simulated willow community.     

Multi-metric evaluation of the models WARM,CropSyst, and WOFOST for rice

WARM (Water Accounting Rice Model) simulates paddy rice (Oryza sativa L.), based on temperature-driven development and radiation-driven crop growth. It also simulates: biomass partitioning, floodwater effect on temperature, spikelet sterility, floodwater and chemicals management, and soil hydrology. Biomass estimates from WARM were evaluated and compared with the ones from two generic crop models (CropSyst, WOFOST). The test-area was the Po Valley (Italy). Data collected at six sites from 1989 to 2004 from rice crops grown under flooded and non-limiting conditions were split into a calibration (to estimate some model parameters) and a validation set. For model evaluation, a fuzzy-logic based multiple-metrics indicator (MQI) was used: 0 (best) ≤ MQI ≤ 1 (worst). WARM estimates compared well with the actual data (mean MQI = 0.037 against 0.167 and 0.173 with CropSyst and WOFOST, respectively). On an average, the three models performed similarly for individual validation metrics such as modelling efficiency (EF > 0.90) and correlation coefficient (R > 0.98). WARM performed best in a weighed measure of the Akaike Information Criterion: (worst) 0<wk<1$0<w_k<1$ (best), considering estimation accuracy and number of parameters required to achieve it (mean wk=0.983$w_k=0.983$ against 0.007 and ∼0.000 with CropSyst and WOFOST, respectively). WARM results were sensitive to 30% of the model parameters (ratio being lower with both CropSyst, <10%, and WOFOST, <20%), but appeared the easiest model to use because of the lowest number of crop parameters required (10 against 15 and 34 with CropSyst and WOFOST, respectively). This study provides a concrete example of the possibilities offered using a range of assessment metrics to evaluate model estimates, predictive capabilities, and complexity.     

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.     

Influences of canopy structure and physiological traits on flux partitioning between understory and overstory in an eastern Siberian boreal larch forest

Boreal forests play an important role in the global balance of energy and CO₂. Our previous study of elaborate eddy covariance observations in a Siberian boreal larch forest, conducted both above the forest canopy and at the forest floor, revealed a significant contribution of latent heat flux (LE) from the cowberry understory to the whole ecosystem LE. Thus, in the present study, we examined what factors control the partitioning of whole ecosystem LE and CO₂ flux into the understory and overstory vegetation, using detailed leaf-level physiology (for both understory and overstory vegetation) and soil respiration property measurements as well as a multilayer soil-vegetation-atmosphere transfer (SVAT) model. The modeling results showed that the larch overstory's leaf area index (LAI) and vertical profile of leaf photosynthetic capacity were major factors determining the flux partitioning in this boreal forest ecosystem. This is unlike other forest ecosystems that tend to have dense LAI. We concluded that control of the larch overstory's LAI had a relationship with both the coexistence of the larch with the cowberry understory and with the water resources available to the total forest ecosystem.     

Climate, canopy conductance and leaf area development controls on evapotranspiration in a boreal coniferous forest over a 10-year period: A united model assessment

The aim of this work was to test a process-based model (hydrological model combined with forest growth model) on the simulation of seasonal variability of evapotranspiration (ET) in an even-aged boreal Scots pine (Pinus sylvestris L.) stand over a 10 year period (1999-2008). The water flux components (including canopy transpiration (E_t) and evaporation from canopy (E_c) and ground surface (E_g) were estimated in order to output the long-term stand water budget considering the interaction between climate variations and stand development. For validation, half-hourly data on eddy water vapor fluxes were measured during the 10 growing seasons (May-September). The model predicted well the seasonal course of ET compared to the measured values, but slightly underestimated the water fluxes both in non-drought and drought (2000, 2003 and 2006) years. The prediction accuracy was, on average, higher in drought years. The simulated ET over the 10 years explained, on average, 58% of the daily variations and 84% of the monthly amount of ET. Water amount from E_t contributed most to the ET, with the fractions of E_t, E_c and E_g being, on average, 67, 11 and 23% over the 10-year period, respectively. Regardless of weather conditions, the daily ET was strongly dependent on air temperature (T_a) and vapor pressure deficit (D_a), but less dependent on soil moisture (W_s). On cloudy and rainy days, there was a non-linear relationship between the ET and solar radiation (R_o). During drought years, the model predicted lower daily canopy stomatal conductance (g_cs) compared with non-drought years, leading to a lower level of E_t. The modeled daily g_cs responded well to D_a and W_s. In the model simulation, the annual LAI increased by 35% between 1999 and 2008. The ratio of E_c: ET correlated strongly with LAI. Furthermore, LAI reduced the proportion of E_g as a result of the increased share of E_c and E_t and radiation interception. Although the increase of LAI affected positively E_t, the contribution of E_t in ET was not significantly correlated with LAI. To conclude, although the model predicted reasonably well the seasonal course of ET, the calculation time steps of different processes in the model should be homogenized in the future to increase the prediction accuracy.     

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

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

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.     

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.     

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