Climate change and coastal flooding in Metro Boston: impacts and adaptation strategies

Sea level rise (SLR) due to climate change will increase storm surge height along the 825 km long coastline of Metro Boston, USA. Land at risk consists of urban waterfront with piers and armoring, residential areas with and without seawalls and revetments, and undeveloped land with either rock coasts or gently sloping beachfront and low-lying coastal marshes. Risk-based analysis shows that the cumulative 100 year economic impacts on developed areas from increased storm surge flooding depend heavily upon the adaptation response, location, and estimated sea level rise. Generally it is found that it is advantageous to use expensive structural protection in areas that are highly developed and less structural approaches such as floodproofing and limiting or removing development in less developed or environmentally sensitive areas.     

A numerical study of the air flow within and around a single tree

A three-dimensional non-hydrostatic numerical model has been used to investigate the air flow and turbulence around a single tree.The results for velocity and turbulence distributions have been compared with available data from windtunnel experiments; the agreement is satisfactory.Simulations have been carried out for different meteorological conditions (wind speed, thermal stratification) as well as for different canopy characteristics (stem height, crown diameter, crown height, porosity).Dedicated to Prof. F. Wippermann on the occasion of his 65th birthday.     

Phase equilibria modelling of retrograde amphibole and clinozoisite in mafic eclogite from the Tso Morari massif,northwest India: constraining the P–T–M(H2O) conditions of exhumation

Phase equilibria modelling of post‐peak metamorphic mineral assemblages in (ultra)high‐P mafic eclogite from the Tso Morari massif, Ladakh Himalaya, northwest India, has provided new insights into the potential behaviour and source of metamorphic fluid during exhumation, and constrained the P–T conditions of hydration. A series of P–M(H₂O) pseudosections constructed in the Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O (NCKFMASHTO) system show that a number of petrographically distinct hydration episodes occurred during exhumation from peak P–T conditions (~640 °C, 27–28 kbar), resulting in the formation of abundant compositionally zoned amphibole and minor clinozoisite poikiloblasts at the expense of a peak assemblage dominated by garnet and omphacite. Initial hydration is interpreted to have occurred as a result of the destabilization of talc following isothermal decompression to ~23 kbar, which led to the formation of barroisite–winchite amphibole core domains. An episode of fluid infiltration from an external source at ~19 kbar, with or without syn‐decompressional cooling to ~560 °C, resulted in further barroisitic–winchitic amphibole growth, followed by the formation of clinozoisite poikiloblasts. Continued buoyancy‐driven exhumation to the base of the lower crust is constrained to have taken place with no additional fluid input. A final hydration event is characterized by the formation of magnesiohornblende rims on the barroisite–winchite cores, with the former interpreted to have formed during later prograde overprinting in the middle crust associated with the final stages of exhumation. Notably, the vast majority of externally sourced H₂O, comprising just over half of the current bulk rock fluid content, was added during this later hydration event. In a middle crustal setting, this is interpreted as the result of devolatilization reactions occurring in migmatitic host orthogneiss and/or metasedimentary units, or following the crystallization of partial melt.     

Simulation of less‐mobile porosity dynamics in contrasting sediment water interface porous media

Considering heterogeneity in porous media pore size and connectivity is essential to predicting reactive solute transport across interfaces. However, exchange with less‐mobile porosity is rarely considered in surface water/groundwater recharge studies. Previous research indicates that a combination of pore‐fluid sampling and geoelectrical measurements can be used to quantify less‐mobile porosity exchange dynamics using the time‐varying relation between fluid and bulk electrical conductivity. For this study, we use macro‐scale (10 s of cm) advection–dispersion solute transport models linked with electrical conduction in COMSOL Multiphysics to explore less‐mobile porosity dynamics in two different types of observed sediment water interface porous media. Modeled sediment textures contrast from strongly layered streambed deposits to poorly sorted lakebed sands and cobbles. During simulated ionic tracer perturbations, a lag between fluid and bulk electrical conductivity, and the resultant hysteresis, is observed for all simulations indicating differential loading of pore spaces with tracer. Less‐mobile exchange parameters are determined graphically from these tracer time series data without the need for inverse numerical model simulation. In both sediment types, effective less‐mobile porosity exchange parameters are variable in response to changes in flow direction and fluid flux. These observed flow‐dependent effects directly impact local less‐mobile residence times and associated contact time for biogeochemical reaction. The simulations indicate that for the sediment textures explored here, less‐mobile porosity exchange is dominated by variable rates of advection through the domain, rather than diffusion of solute, for typical low‐to‐moderate rate (approximately 3–40 cm/day) hyporheic fluid fluxes. Overall, our model‐based results show that less‐mobile porosity may be expected in a range of natural hyporheic sediments and that changes in flowpath orientation and magnitude will impact less‐mobile exchange parameters. These temporal dynamics can be assessed with the geoelectrical experimental tracer method applied at laboratory and field scales.     

Delimiting the eastern extent of the Altyn Tagh Fault: Insights from structural analyses of seismic reflection profiles

The Altyn Tagh Fault (ATF) serves as a key continental‐scale controlling structural element of the Tibetan Plateau. However, its eastward extent remains controversial. Here we use high‐resolution seismic reflection profiles to investigate the subsurface structures of the easternmost ATF and use these to delimit the easternmost extent of the fault. The structural analyses show an eastward geometric change from transpressional positive flower structures to compressional thrusts, with transpression‐induced shortening magnitudes decreasing eastwards from a maximum of ~5.3 km to being absent. Stratigraphic controls indicate that the deformation took place over the last ~<1.2 Ma. Our wider findings lead us to: (a) reject the suggestion that the ATF previously extended beyond the Kuantan Shan‐Hei Shan to link with the Alxa‐East Mongolia Fault; and (b) propose that the rigid block model used to describe the Tibetan Plateau crust is not consistent with the extent and structural details of the easternmost ATF.     

A luminescence‐based chronology for the Harletz loess sequence,Bulgaria

The Harletz loess‐palaeosol sequence is located in northwestern Bulgaria and represents an important link between well‐studied loess sequences in eastern Romania and further sites to the west of the Carpathians (e.g. Serbia and Hungary). The aim of this study was to establish a chronostratigraphy of the deposits, using various methods of luminescence dating, together with basic stratigraphical field observations as well as magnetic properties. Luminescence dating was carried out using the quartz fine grain fraction and a SAR protocol, and the feldspar coarse grain fraction, applying the MET‐pIRIR protocol. Due to underestimation of the quartz fine grain fraction in the lower parts of the sequence, the resulting chronology is mainly based on the feldspar ages, which are derived from the stimulation temperature at 150 °C. A comparison with nearby sequences from Serbia, Hungary and Romania, and interpretations obtained through the stratigraphical and sedimentological signature of the sequence, supports the established chronology. Our data suggest that the prominent palaeosol (soil complex) in the upper quarter of the sequence was formed during MIS 5. It would follow that large parts of the Last Glacial loess overlying this palaeosol were probably eroded, and that the thick loess accumulation underlying this soil complex can be allocated to the penultimate glacial (MIS 6). A prominent MIS 6 tephra, which has been reported from other sequences in the area, is also present at Harletz.     

Subduction metamorphism of serpentinite‐hosted carbonates beyond antigorite‐serpentinite dehydration (Nevado‐Filábride Complex,Spain)

At sub‐arc depths, the release of carbon from subducting slab lithologies is mostly controlled by fluid released by devolatilization reactions such as dehydration of antigorite (Atg‐) serpentinite to prograde peridotite. Here we investigate carbonate–silicate rocks hosted in Atg‐serpentinite and prograde chlorite (Chl‐) harzburgite in the Milagrosa and Almirez ultramafic massifs of the palaeo‐subducted Nevado‐Filábride Complex (NFC, Betic Cordillera, S. Spain). These massifs provide a unique opportunity to study the stability of carbonate during subduction metamorphism at P–T conditions before and after the dehydration of Atg‐serpentinite in a warm subduction setting. In the Milagrosa massif, carbonate–silicate rocks occur as lenses of Ti‐clinohumite–diopside–calcite marbles, diopside–dolomite marbles and antigorite–diopside–dolomite rocks hosted in clinopyroxene‐bearing Atg‐serpentinite. In Almirez, carbonate–silicate rocks are hosted in Chl‐harzburgite and show a high‐grade assemblage composed of olivine, Ti‐clinohumite, diopside, chlorite, dolomite, calcite, Cr‐bearing magnetite, pentlandite and rare aragonite inclusions. These NFC carbonate–silicate rocks have variable CaO and CO₂ contents at nearly constant Mg/Si ratio and high Ni and Cr contents, indicating that their protoliths were variable mixtures of serpentine and Ca‐carbonate (i.e., ophicarbonates). Thermodynamic modelling shows that the carbonate–silicate rocks attained peak metamorphic conditions similar to those of their host serpentinite (Milagrosa massif; 550–600°C and 1.0–1.4 GPa) and Chl‐harzburgite (Almirez massif; 1.7–1.9 GPa and 680°C). Microstructures, mineral chemistry and phase relations indicate that the hybrid carbonate–silicate bulk rock compositions formed before prograde metamorphism, likely during seawater hydrothermal alteration, and subsequently underwent subduction metamorphism. In the CaO–MgO–SiO₂ ternary, these processes resulted in a compositional variability of NFC serpentinite‐hosted carbonate–silicate rocks along the serpentine‐calcite mixing trend, similar to that observed in serpentinite‐hosted carbonate‐rocks in other palaeo‐subducted metamorphic terranes. Thermodynamic modelling using classical models of binary H₂O–CO₂ fluids shows that the compositional variability along this binary determines the temperature of the main devolatilization reactions, the fluid composition and the mineral assemblages of reaction products during prograde subduction metamorphism. Thermodynamic modelling considering electrolytic fluids reveals that H₂O and molecular CO₂ are the main fluid species and charged carbon‐bearing species occur only in minor amounts in equilibrium with carbonate–silicate rocks in warm subduction settings. Consequently, accounting for electrolytic fluids at these conditions slightly increases the solubility of carbon in the fluids compared with predictions by classical binary H₂O–CO₂ fluids, but does not affect the topology of phase relations in serpentinite‐hosted carbonate‐rocks. Phase relations, mineral composition and assemblages of Milagrosa and Almirez (meta)‐serpentinite‐hosted carbonate–silicate rocks are consistent with local equilibrium between an infiltrating fluid and the bulk rock composition and indicate a limited role of infiltration‐driven decarbonation. Our study shows natural evidence for the preservation of carbonates in serpentinite‐hosted carbonate–silicate rocks beyond the Atg‐serpentinite breakdown at sub‐arc depths, demonstrating that carbon can be recycled into the deep mantle.     

Testing a spatially distributed tracer‐aided runoff model in a snow‐influenced catchment: Effects of multicriteria calibration on streamwater ages

Integrating stable isotope tracers into rainfall‐runoff models allows investigation of water partitioning and direct estimation of travel times and water ages. Tracer data have valuable information content that can be used to constrain models and, in integration with hydrometric observations, test the conceptualization of catchment processes in model structure and parameterization. There is great potential in using tracer‐aided modelling in snow‐influenced catchments to improve understanding of these catchments' dynamics and sensitivity to environmental change. We used the spatially distributed tracer‐aided rainfall‐runoff (STARR) model to simulate the interactions between water storage, flux, and isotope dynamics in a snow‐influenced, long‐term monitored catchment in Ontario, Canada. Multiple realizations of the model were achieved using a combination of single and multiple objectives as calibration targets. Although good simulations of hydrometric targets such as discharge and snow water equivalent could be achieved by local calibration alone, adequate capture of the stream isotope dynamics was predicated on the inclusion of isotope data in the calibration. Parameter sensitivity was highest, and most local, for single calibration targets. With multiple calibration targets, key sensitive parameters were still identifiable in snow and runoff generation routines. Water ages derived from flux tracking subroutines in the model indicated a catchment where runoff is dominated by younger waters, particularly during spring snowmelt. However, resulting water ages were most sensitive to the partitioning of runoff sources from soil and groundwater sources, which was most realistically achieved when isotopes were included in the calibration. Given the paucity of studies where hydrological models explicitly incorporate tracers in snow‐influenced regions, this study using STARR is an important contribution to satisfactorily simulating snowpack dynamics and runoff generation processes, while simultaneously capturing stable isotope variability in snow‐influenced catchments.     

Storage in confined aquifer: Spectral analysis of groundwater in responses to Earth tides and barometric effect

Dilatation of aquifer and associated water level fluctuation in groundwater well is known to be driven periodically from lunar, solar, or other tidal forces. Time‐dependent variables in groundwater system, such as water level, can be converted to power spectra in the frequency domain using Fourier transform to evaluate significant fluctuation. The major innovation of this research is to develop spectral representation in frequency domain for the groundwater system that the storage in confined aquifer can be determined considering dilatation affected by Earth tides and barometric effect. In order to verify applicability of the evolved method, time series of Earth tides and barograph are collected; aquifer storage is then determined inversely by selecting significant semidiurnal and diurnal components in spectra computation. It suggests that to discover groundwater storage using groundwater level with barograph and tidal potential of Earth in frequency domain becomes accessible and feasible.