The influence of increasing atmospheric CO2 , temperature, and vapor pressure deficit on seawater-induced tree mortality.

Increasing seawater exposure is killing coastal trees globally, with expectations of accelerating mortality with rising sea levels. However, the impact of concomitant changes in atmospheric CO2 concentration, temperature, and vapor pressure deficit (VPD) on seawater-induced tree mortality is uncertain. We examined the mechanisms of seawater-induced mortality under varying climate scenarios using a photosynthetic gain and hydraulic cost optimization model validated against observations in a mature stand of Sitka spruce (Picea sitchensis) trees in the Pacific Northwest, USA, that were dying from recent seawater exposure. The simulations matched well with observations of photosynthesis, transpiration, nonstructural carbohydrates concentrations, leaf water potential, the percentage loss of xylem conductivity, and stand-level mortality rates. The simulations suggest that seawater-induced mortality could decrease by c. 16.7% with increasing atmospheric CO2 levels due to reduced risk of carbon starvation. Conversely, rising VPD could increase mortality by c. 5.6% because of increasing risk of hydraulic failure. Across all scenarios, seawater-induced mortality was driven by hydraulic failure in the first 2 yr after seawater exposure began, with carbon starvation becoming more important in subsequent years. Changing CO2 and climate appear unlikely to have a significant impact on coastal tree mortality under rising sea levels.

Seawater exposure causes hydraulic damage in dying Sitka-spruce trees.

Sea-level rise is one of the most critical challenges facing coastal ecosystems under climate change. Observations of elevated tree mortality in global coastal forests are increasing, but important knowledge gaps persist concerning the mechanism of salinity stress-induced nonhalophytic tree mortality. We monitored progressive mortality and associated gas exchange and hydraulic shifts in Sitka-spruce (Picea sitchensis) trees located within a salinity gradient under an ecosystem-scale change of seawater exposure in Washington State, USA. Percentage of live foliated crown (PLFC) decreased and tree mortality increased with increasing soil salinity during the study period. A strong reduction in gas exchange and xylem hydraulic conductivity (Ks) occurred during tree death, with an increase in the percentage loss of conductivity (PLC) and turgor loss point (πtlp). Hydraulic and osmotic shifts reflected that hydraulic function declined from seawater exposure, and dying trees were unable to support osmotic adjustment. Constrained gas exchange was strongly related to hydraulic damage at both stem and leaf levels. Significant correlations between foliar sodium (Na+) concentration and gas exchange and key hydraulic parameters (Ks, PLC, and πtlp) suggest that cellular injury related to the toxic effects of ion accumulation impacted the physiology of these dying trees. This study provides evidence of toxic effects on the cellular function that manifests in all aspects of plant functioning, leading to unfavourable osmotic and hydraulic conditions.

Changes in carbon and nitrogen metabolism during seawater-induced mortality of Picea sitchensis trees.

Increasing seawater exposure is causing mortality of coastal forests, yet the physiological response associated with seawater-induced tree mortality, particularly in non-halophytes, is poorly understood. We investigated the shifts in carbon and nitrogen (N) metabolism of mature Sitka-spruce trees that were dying after an ecosystem-scale manipulation of tidal seawater exposure. Soil porewater salinity and foliar ion concentrations increased after seawater exposure and were strongly correlated with the percentage of live foliated crown (PLFC; e.g., crown 'greenness', a measure of progression to death). Co-occurring with decreasing PLFC was decreasing photosynthetic capacity, N-investment into photosynthesis, N-resorption efficiency and non-structural carbohydrate (soluble sugars and starch) concentrations, with the starch reserves depleted to near zero when PLFC dropped below 5%. Combined with declining PLFC, these changes subsequently decreased total carbon gain and thus exacerbated the carbon starvation process. This study suggests that an impairment in carbon and N metabolism during the mortality process after seawater exposure is associated with the process of carbon starvation, and provides critical knowledge necessary to predict sea-level rise impacts on coastal forests.

Declining carbohydrate content of Sitka-spruce treesdying from seawater exposure.

Increasing sea levels associated with climate change threaten the survival of coastal forests, yet the mechanisms by which seawater exposure causes tree death remain poorly understood. Despite the potentially crucial role of nonstructural carbohydrate (NSC) reserves in tree survival, their dynamics in the process of death under seawater exposure are unknown. Here we monitored progressive tree mortality and associated NSC storage in Sitka-spruce (Picea sitchensis) trees dying under ecosystem-scale increases in seawater exposure in western Washington, USA. All trees exposed to seawater, because of monthly tidal intrusion, experienced declining crown foliage during the sampling period, and individuals with a lower percentage of live foliated crown (PLFC) died faster. Tree PLFC was strongly correlated with subsurface salinity and needle ion contents. Total NSC concentrations in trees declined remarkably with crown decline, and reached extremely low levels at tree death (2.4% and 1.6% in leaves and branches, respectively, and 0.4% in stems and roots). Starch in all tissues was almost completely consumed, while sugars remained at a homeostatic level in foliage. The decreasing NSC with closer proximity to death and near zero starch at death are evidences that carbon starvation occurred during Sitka-spruce mortality during seawater exposure. Our results highlight the importance of carbon storage as an indicator of tree mortality risks under seawater exposure.

Water Table Dynamics and Biogeochemical Cycling in a Shallow, Variably-Saturated Floodplain.

Three-dimensional variably saturated flow and multicomponent biogeochemical reactive transport modeling, based on published and newly generated data, is used to better understand the interplay of hydrology, geochemistry, and biology controlling the cycling of carbon, nitrogen, oxygen, iron, sulfur, and uranium in a shallow floodplain. In this system, aerobic respiration generally maintains anoxic groundwater below an oxic vadose zone until seasonal snowmelt-driven water table peaking transports dissolved oxygen (DO) and nitrate from the vadose zone into the alluvial aquifer. The response to this perturbation is localized due to distinct physico-biogeochemical environments and relatively long time scales for transport through the floodplain aquifer and vadose zone. Naturally reduced zones (NRZs) containing sediments higher in organic matter, iron sulfides, and non-crystalline U(IV) rapidly consume DO and nitrate to maintain anoxic conditions, yielding Fe(II) from FeS oxidative dissolution, nitrite from denitrification, and U(VI) from nitrite-promoted U(IV) oxidation. Redox cycling is a key factor for sustaining the observed aquifer behaviors despite continuous oxygen influx and the annual hydrologically induced oxidation event. Depth-dependent activity of fermenters, aerobes, nitrate reducers, sulfate reducers, and chemolithoautotrophs (e.g., oxidizing Fe(II), S compounds, and ammonium) is linked to the presence of DO, which has higher concentrations near the water table.

Variably saturated flow and multicomponent biogeochemical reactive transport modeling of a uranium bioremediation field experiment.

Three-dimensional, coupled variably saturated flow and biogeochemical reactive transport modeling of a 2008 in situ uranium bioremediation field experiment is used to better understand the interplay of transport and biogeochemical reactions controlling uranium behavior under pulsed acetate amendment, seasonal water table variation, spatially variable physical (hydraulic conductivity, porosity) and geochemical (reactive surface area) material properties. While the simulation of the 2008 Big Rusty acetate biostimulation field experiment in Rifle, Colorado was generally consistent with behaviors identified in previous field experiments at the Rifle IFRC site, the additional process and property detail provided several new insights. A principal conclusion from this work is that uranium bioreduction is most effective when acetate, in excess of the sulfate-reducing bacteria demand, is available to the metal-reducing bacteria. The inclusion of an initially small population of slow growing sulfate-reducing bacteria identified in proteomic analyses led to an additional source of Fe(II) from the dissolution of Fe(III) minerals promoted by biogenic sulfide. The falling water table during the experiment significantly reduced the saturated thickness of the aquifer and resulted in reactants and products, as well as unmitigated uranium, in the newly unsaturated vadose zone. High permeability sandy gravel structures resulted in locally high flow rates in the vicinity of injection wells that increased acetate dilution. In downgradient locations, these structures created preferential flow paths for acetate delivery that enhanced local zones of TEAP reactivity and subsidiary reactions. Conversely, smaller transport rates associated with the lower permeability lithofacies (e.g., fine) and vadose zone were shown to limit acetate access and reaction. Once accessed by acetate, however, these same zones limited subsequent acetate dilution and provided longer residence times that resulted in higher concentrations of TEAP reaction products when terminal electron donors and acceptors were not limiting. Finally, facies-based porosity and reactive surface area variations were shown to affect aqueous uranium concentration distributions with localized effects of the fine lithofacies having the largest impact on U(VI) surface complexation. The ability to model the comprehensive biogeochemical reaction network, and spatially and temporally variable processes, properties, and conditions controlling uranium behavior during engineered bioremediation in the naturally complex Rifle IFRC subsurface system required a subsurface simulator that could use the large memory and computational performance of a massively parallel computer. In this case, the eSTOMP simulator, operating on 128 processor cores for 12h, was used to simulate the 110-day field experiment and 50 days of post-biostimulation behavior.

Identifying the potential loss of monitoring wells using an uncertainty analysis.

From the mid-1940s through the 1980s, large volumes of waste water were discharged at the Hanford Site in southeastern Washington State, causing a large-scale rise (>20 m) in the water table. When waste water discharges ceased in 1988, ground water mounds began to dissipate. This caused a large number of wells to go dry and has made it difficult to monitor contaminant plume migration. To identify monitoring wells that will need replacement, a methodology has been developed using a first-order uncertainty analysis with UCODE, a nonlinear parameter estimation code. Using a three-dimensional, finite-element ground water flow code, key parameters were identified by calibrating to historical hydraulic head data. Results from the calibration period were then used to check model predictions by comparing monitoring wells' wet/dry status with field data. This status was analyzed using a methodology that incorporated the 0.3 cumulative probability derived from the confidence and prediction intervals. For comparison, a nonphysically based trend model was also used as a predictor of wells' wet/dry status. Although the numerical model outperformed the trend model, for both models, the central value of the intervals was a better predictor of a wet well status. The prediction interval, however, was more successful at identifying dry wells. Predictions made through the year 2048 indicated that 46% of the wells in the monitoring well network are likely to go dry in areas near the river and where the ground water mound is dissipating.