Groundwater control of mangrove surface elevation: Shrink and swell varies with soil depth

We measured monthly soil surface elevation change and determined its relationship to groundwater changes at a mangrove forest site along Shark River, Everglades National Park, Florida. We combined the use of an original design, surface elevation table with new rod-surface elevation tables to separately track changes in the mid zone (0–4 m), the shallow root zone (0–0.35 m), and the full sediment profile (0–6 m) in response to site hydrology (daily river stage and daily groundwater piezometric pressure). We calculated expansion and contraction for each of the four constituent soil zones (surface [accretion and erosion; above 0 m], shallow zone [0–0.35 m], middle zone [0.35–4 m], and bottom zone [4–6]) that comprise the entire soil column. Changes in groundwater pressure correlated strongly, with changes in soil elevation for the entire profile (Adjusted R² = 0.90); this relationship was not proportional to the depth of the soil profile sampled. The change in thickness of the bottom soil zone accounted for the majority (R² = 0.63) of the entire soil profile expansion and contraction. The influence of hydrology on specific soil zones and absolute elevation change must be considered when evaluating the effect of disturbances, sea level rise, and water management decisions on coastal wetland systems.     

Effects of pressure on aqueous chemical equilibria at subzero temperatures with applications to Europa

Pressure plays a critical role in controlling aqueous geochemical processes in deep oceans and deep ice. The putative ocean of Europa could have pressures of 1200 bars or higher on the seafloor, a pressure not dissimilar to the deepest ocean basin on Earth (the Mariana Trench at 1100 bars of pressure). At such high pressures, chemical thermodynamic relations need to explicitly consider pressure. A number of papers have addressed the role of pressure on equilibrium constants, activity coefficients, and the activity of water. None of these models deal, however, with processes at subzero temperatures, which may be important in cold environments on Earth and other planetary bodies. The objectives of this work were to (1) incorporate a pressure dependence into an existing geochemical model parameterized for subzero temperatures (FREZCHEM), (2) validate the model, and (3) simulate pressure-dependent processes on Europa. As part of objective 1, we examined two models for quantifying the volumetric properties of liquid water at subzero temperatures: one model is based on the measured properties of supercooled water, and the other model is based on the properties of liquid water in equilibrium with ice.The relative effect of pressure on solution properties falls in the order: equilibrium constants(K) > activity coefficients (γ) > activity of water (a_w). The errors (%) in our model associated with these properties, however, fall in the order: γ > K > a_w. The transposition between K and γ is due to a more accurate model for estimating K than for estimating γ. Only activity coefficients are likely to be significantly in error. However, even in this case, the errors are likely to be only in the range of 2 to 5% up to 1000 bars of pressure. Evidence based on the pressure/temperature melting of ice and salt solution densities argue in favor of the equilibrium water model, which depends on extrapolations, for characterizing the properties of liquid water in electrolyte solutions at subzero temperatures, rather than the supercooled water model. Model-derived estimates of mixed salt solution densities and chemical equilibria as a function of pressure are in reasonably good agreement with experimental measurements.To demonstrate the usefulness of this low-temperature, high-pressure model, we examined two hypothetical cases for Europa. Case 1 dealt with the ice cover of Europa, where we asked the question: How far above the putative ocean in the ice layer could we expect to find thermodynamically stable brine pockets that could serve as habitats for life? For a hypothetical nonconvecting 20 km icy shell, this potential life zone only extends 2.8 km into the icy shell before the eutectic is reached. For the case of a nonconvecting icy shell, the cold surface of Europa precludes stable aqueous phases (habitats for life) anywhere near the surface. Case 2 compared chemical equilibria at 1 bar (based on previous work) with a more realistic 1460 bars of pressure at the base of a 100 km Europan ocean. A pressure of 1460 bars, compared to 1 bar, caused a 12 K decrease in the temperature at which ice first formed and a 11 K increase in the temperature at which MgSO₄·12H₂O first formed. Remarkably, there was only a 1.2 K decrease in the eutectic temperatures between 1 and 1460 bars of pressure. Chemical systems and their response to pressure depend, ultimately, on the volumetric properties of individual constituents, which makes every system response highly individualistic.     

The influence of edge sites on the development of surface charge on goethite nanoparticles: A molecular dynamics investigation

Large-scale molecular simulation of proton accumulations were carried out on (i) (110) and (021) slabs immersed in aqueous solution and (ii) a series of model goethite nanoparticles of dimension 2 to 8 nm with systematically varying acicularity and (110)/(021) surface areas. In the slab systems, the (021) surface exhibits 15% more proton charge per unit area than the (110) surface. In the particulate systems, the acicular particles having the highest (110)/(021) ratio accumulate the most charge, opposite to the trend expected from the slab simulations, indicating that, at length scales on the order of 10 nm, the slab results are not a good indicator of the overall charging behavior of the particles. The primary reason for the discrepancy between the particulate systems and slab systems is the preferential accumulation of protons at acute (110)-(110) intersections. Charge accumulates preferentially in this region because excess proton charge at an asperity is more effectively solvated than at a flat interface.     

Environmental significance of C/C and O/O ratios of modern land-snail shells from the southern great plains of North America

¹³C/¹²C and ¹⁸O/¹⁶O ratios of aragonite shells of modern land snails from the southern Great Plains of North America were measured for samples from twelve localities in a narrow east-west corridor that extended from the Flint Hills in North Central Oklahoma to the foothills of the Sangre de Cristo Mountains in Northern New Mexico, USA. Across the study area, shell δ¹⁸O values (PDB scale) ranged from −4.1‰ to 1.2‰, while δ¹³C values ranged from −13.2‰ to 0.0‰. δ¹⁸O values of the shell aragonite were predicted with a published, steady state, evaporative flux balance model. The predicted values differed (with one exception) by less than 1‰ from locality averages of measured δ¹⁸O values. This similarity suggests that relative humidity at the time of snail activity is an important control on the δ¹⁸O values of the aragonite and emphasizes the seasonal nature of the climatic information preserved in the shells. Correlated δ¹³C values of coexisting Vallonia and Gastrocopta suggest similar feeding habits and imply that these genera can provide information on variations in southern Great Plains plant ecology. Although there is considerable scatter, multispecies, transect average δ¹³C values of the modern aragonite shells are related to variations in the type of photosynthesis (i.e., C₃, C₄) in the local plant communities. The results of this study emphasize the desirability of obtaining isotope ratios representing averages of many shells in a locale to reduce possible biases associated with local variations among individuals, species, etc., and thus better represent the “neighborhood” scale temporal and/or spatial environmental variations of interest in studies of modern and ancient systems.     

Critical rainfall statistics for predicting watershed flood responses: rethinking the design storm concept

Recent advances have been made to modernize estimates of probable precipitation scenarios; however, researchers and engineers often continue to assume that rainfall events can be described by a small set of event statistics, typically average intensity and event duration. Given the easy availability of precipitation data and advances in desk‐top computational tools, we suggest that it is time to rethink the ‘design storm’ concept. Design storms should include more holistic characteristics of flood‐inducing rain events, which, in addition to describing specific hydrologic responses, may also be watershed or regionally specific. We present a sensitivity analysis of nine precipitation event statistics from observed precipitation events within a 60‐year record for Tompkins County, NY, USA. We perform a two‐sample Kolmogorov–Smirnov (KS) test to objectively identify precipitation event statistics of importance for two related hydrologic responses: (1) peak outflow from the Six Mile Creek watershed and (2) peak depth within the reservoir behind the Six Mile Creek Dam. We identify the total precipitation depth, peak hourly intensity, average intensity, event duration, interevent duration, and several statistics defining the temporal distribution of precipitation events to be important rainfall statistics to consider for predicting the watershed flood responses. We found that the two hydrologic responses had different sets of statistically significant parameters. We demonstrate through a stochastic precipitation generation analysis the effects of starting from a constrained parameter set (intensity and duration) when predicting hydrologic responses as opposed to utilizing an expanded suite of rainfall statistics. In particular, we note that the reduced precipitation parameter set may underestimate the probability of high stream flows and therefore underestimate flood hazard. Copyright © 2016 John Wiley & Sons, Ltd.     

Age of Gimli beach of Lake Agassiz based on new OSL dating

Gimli beach in Manitoba is one of the lowest elevation beaches in the southern Lake Agassiz basin, and is a distinct ridge composed of bedded sand and gravel that rises above the lake plain and extends for more than 40 km. Ten new optically stimulated luminescence (OSL) ages from Gimli beach yield ages mostly ranging from 9.7 ± 0.7 to 10.5 ± 0.8 ka (average 10.3 ± 0.5 ka), which is older by 0.6 to >1.0 ka than age estimates of previous researchers. Two of our new OSL ages are notably older than the others, dating to ~11.3 ± 0.8 and 13.9 ± 1.0 ka, which we attribute to poorly bleached sands. We ascribe an age of about 10 ka to Gimli beach, which is several centuries before overflow from Lake Agassiz and its vast drainage basin shifted from the western Great Lakes to glacial Lake Ojibway and the St. Lawrence Valley.     

Landscape controls on long‐term runoff in subhumid heterogeneous Boreal Plains catchments

We compared median runoff (R) and precipitation (P) relationships over 25 years from 20 mesoscale (50 to 5,000 km²) catchments on the Boreal Plains, Alberta, Canada, to understand controls on water sink and source dynamics in water‐limited, low‐relief northern environments. Long‐term catchment R and runoff efficiency (RP^?1) were low and varied spatially by over an order of magnitude (3 to 119 mm/year, 1 to 27%). Intercatchment differences were not associated with small variations in climate. The partitioning of P into evapotranspiration (ET) and R instead reflected the interplay between underlying glacial deposit texture, overlying soil‐vegetation land cover, and regional slope. Correlation and principal component analyses results show that peatland‐swamp wetlands were the major source areas of water. The lowest estimates of median annual catchment ET (321 to 395 mm) and greatest R (60 to 119 mm, 13 to 27% of P) were observed in low‐relief, peatland‐swamp dominated catchments, within both fine‐textured clay‐plain and coarse‐textured glacial deposits. In contrast, open‐water wetlands and deciduous‐mixedwood forest land covers acted as water sinks, and less catchment R was observed with increases in proportional coverage of these land covers. In catchments dominated by hummocky moraines, long‐term runoff was restricted to 10 mm/year, or 2% of P. This reflects the poor surface‐drainage networks and slightly greater regional slope of the fine‐textured glacial deposit, coupled with the large soil‐water and depression storage and higher actual ET of associated shallow open‐water marsh wetland and deciduous‐forest land covers. This intercatchment study enhances current conceptual frameworks for predicting water yield in the Boreal Plains based on the sink and source functions of glacial landforms and soil‐vegetation land covers. It offers the capability within this hydro‐geoclimatic region to design reclaimed catchments with desired hydrological functionality and associated tolerances to climate or land‐use changes and inform land management decisions based on effective catchment‐scale conceptual understanding.     

Dynamic storage: a potential metric of inter‐basin differences in storage properties

The potential for dynamic storage to serve as a metric of basin behaviour was assessed using data from five drainage basins with headwaters on the thick sand and gravel deposits of the Oak Ridges Moraine in southern Ontario, Canada. Dynamic storage was directly correlated with the ratio of variability of δ²H in streamflow relative to that in precipitation. This ratio has previously been shown to be inversely related to basin mean transit time (MTT), suggesting an inverse relationship between dynamic storage and MTT for the study basins. Dynamic storage was also directly correlated with interannual variability in stream runoff, baseflow and baseflow:runoff ratio, implying that basins with smaller dynamic storage have less interannual variability in their streamflow regimes. These preliminary results suggest that dynamic storage may serve as a readily derived and useful metric of basin behaviour for inter‐basin comparisons. Copyright © 2016 John Wiley & Sons, Ltd.