Intercomparison of soil pore water extraction methods for stable isotope analysis and interpretation of hillslope runoff sources

Intercomparison of soil pore water extraction methods for stable isotope analysis has been a focus of recent studies in relation to plant source waters, which found a wide isotopic variance depending on the extraction method. Few studies have yet explored extraction effects for mobile pore waters that relate to hillslope runoff. This is because it is extremely difficult in natural systems to control the boundary conditions in order to assess and compare impacts of pore water extraction on resulting hillslope flow. With our new semicontrolled experiments on outdoor mini‐hillslopes, we studied mixing and runoff processes by means of stable isotopes of water and quantified relations between pore water extraction methods. We tested the null hypothesis that nondestructive and destructive pore water sampling methods sample the same soil water pool. Three hillslopes were mounted on load cells, filled with loamy sand textured soils from the Landscape Evolution Observatoryat Biosphere 2, equipped with soil moisture and temperature sensors, a bottom outflow, and a surface runoff gauge for isotope sampling. We followed the precipitation isotopic composition over and through the soil profile. One hillslope was instrumented with suction cups, on the second we installed sampling ports for in‐situ soil water vapour measurements, and the third hillslope was sampled destructively for applying the centrifugation and vapour equilibrium methods. All hillslopes were sampled at four depths (0–10, 10–20, 20–30, and 30–40 cm) at three different downslope positions. ²H and ¹⁸O analyses were performed via laser spectroscopy. We found no isotopic differences between rainfall, surface runoff, and bottom outflow. The in situ vapour ports' soil isotope data showed the widest spread over all hillslope positions and depths. Centrifugation's and suction cups' isotope results plotted closest to the local meteoric water line and within the range of hillslope runoff and bottom outflow data. Hillslope position did not influence the soil isotope results. These results suggest caution be used in the field when selecting an extraction technique for matching soil waters to runoff waters. Soil suction lysimeters and centrifugation appeared to be the most appropriate tools in this regard.     

Dynamic domain kinematic modelling for predicting interflow over leaky impeding layers

Traditional Boussinesq or kinematic simulations of interflow (i.e., lateral subsurface flow) assume no leakage through the impeding layer and require a no-flow boundary condition at the ridge top. However, recent analyses of many interflow-producing landscapes indicate that leaky impeding layers are common, that most interflow percolates well before reaching the toe slope, and therefore, the downslope contributing length is shorter than the hillslope length. In watersheds characterised by perched interflow over a low conductivity layer through permeable topsoil, interflow with percolation may be modelled with a kinematic wave model using a mobile upslope boundary condition defining the hillslope portion contributing interflow to valleys. Here, we developed and applied a dynamic interflow model to simulate interflow using a downslope travel distance concept such that only the active contributing length is modelled at any time. The model defines a variable active area based on the depth of the perched layer, the topographic slope and the ratio of the hydraulic conductivity of topsoil to that of the impeding layer. It incorporates a two-layer soil moisture accounting water balance analysis, a pedo-transfer function, and percolation and evaporation routines to predict interflow rates in continuous and event-based scenarios. We tested the modelling concept on two sets of data (2-year dataset of rainfall observations for the continuous simulation and a multi-day irrigation experiment for the event simulation) from a 121-m-long open interflow collection trench on an experimental hillslope at the Savannah River Site, South Carolina. The continuous model simulation partially represented the observed interflow hydrograph and perched water depth in the experimental hillslope with correlation coefficients of 0.85 and 0.35, respectively. Model performance improved significantly at event-scale analysis. The modelling approach realistically represents interflow dynamics in hillslopes with leaky impeding layers and can be integrated into catchment-scale hydrology models for more detailed hillslope process modelling.     

White mica trace element and boron isotope evidence for distinctive infiltration events during exhumation of deeply subducted continental crust

ABSTRACT

Previous study of subducted continental crust within the Luliang Shan terrane in Northwest China has documented metasomatic formation of thick, hydrated phengite + garnet-rich selvages at the interface between mafic eclogite blocks and quartzofeldspathic host gneiss. Whole rock concentrations of Cs and Ba within the selvage are enriched by two orders of magnitude relative to the eclogite blocks and host gneiss. We performed in situ ion microprobe analyses of Li, Be, B, Rb, Sr, Cs and Ba and δ¹¹B of phengite within the Luliang Shane terrane to better constrain the source(s) of the infiltrating fluid. The phengite within the selvage are enriched in Li, Cs and Ba and yield δ¹¹B values between ?30‰ and ?9‰, values that are lower than mantle values. High Ba/Rb, Cs/Rb coupled with low B/Be, B/Li and highly negative δ¹¹B values indicate that the high-pressure fluid that formed the selvage was derived from highly devolatilized rocks within the subduction channel. In contrast, muscovite, which crystallized in the adjacent host gneiss during a subsequent lower pressure phase of fluid infiltration at approximately 0.9 GPa depths, has much lower Li, Cs and Ba relative to the high-pressure phengite. These retrograde muscovite have very high concentrations of B (up to 5500 ppm) and Be (up to 50 ppm) and high (?2 to +8‰) δ¹¹B values that are consistent with crystallization from a fluid derived from shallower and less devolatilized regions of the subduction zone. Additional host gneiss samples, regionally distributed and kilometres away from the studied area lack the B-rich signature and indicate that the late stage fluids were likely localized to the region near the studied traverse.     

Multiple pulses in dE/dt and the fine-structure of E during the onset of first return strokes in cloud-to-ocean lightning

We have analyzed the fine-structure of 131 electric field (E) waveforms that were radiated during the onset of first return strokes in cloud-to-ocean lightning. The dE/dt waveforms were recorded using an 8-bit waveform digitizer sampling at 100 MHz, and the E waveforms were sampled at 10 MHz using a 10-bit digitizer. 49 (or 37%) of the dE/dt waveforms contain one or more large pulses within ± 1 μs of the largest (or dominant) peak in dE/dt, i.e. within an interval from − 1 μs to + 1 μs, where t = 0 μs is the time of the dominant peak, and 37 (or 28%) have one or more large pulses in the interval from 4 μs before to 1 μs before the dominant peak, i.e. − 4 μs to − 1 μs, and only the dominant peak within ± 1 μs . We give statistics on the amplitude and timing of dE/dt pulses that are near the dominant peak, and we show how the presence of these pulses adds considerable fine-structure to the shape of E_int, the integrated dE/dt waveform, on a time-scale of tens to hundreds of nanoseconds. This fine-structure includes fast pulses near the beginning of the slow front, large pulses and shoulders within the slow front and during the fast-transition, and very narrow peaks in E_int. Our overall conclusion is that the electromagnetic environment near the point(s) where lightning leaders attach to the surface is often more complicated than what would be produced by a single current pulse propagating up a single channel at the time of onset.     

Direct and Indirect Interactions Between Lower Estuarine Mangrove and Saltmarsh Habitats and a Commercially Important Penaeid Shrimp

Vegetated habitats in estuaries may provide a structural refuge and food supply in the same place, but benefits are also derived where a productive food source and suitable habitat are adjacent to each other. Quantifying these relationships is fundamental to understanding the structural and functional characteristics of estuarine ecosystems and for informing management actions. Effective juvenile habitat (habitat that contributes greater-than-average numbers of recruits to the adult population), recruitment patterns and trophic relationships were studied for Eastern King Prawn (Penaeus plebejus) in the lower Clarence River estuary, New South Wales, between 2014 and 2016. Effective juvenile habitat was identified in both the north arm and main river channel of the estuary, and these areas also supported a higher abundance of juvenile prawns. There was minimal recruitment to the southern channels of the estuary, possibly due to reduced connectivity with the incoming tide arising from a rock wall. Trophic relationships in parts of the lower estuary were evaluated using stable isotopes, and saltmarsh grass (Sporobolus virginicus) was the dominant primary producer supporting juvenile Eastern King Prawn productivity across the area. Mangroves were of minimal importance, and seagrass cover was minimal in the area studied. The patterns observed indicate that nursery function of different areas within the lower estuary is a product of connectivity, recruitment and nutrition derived from primary productivity of vascular plants. Habitats within the lower Clarence River estuary have seen substantial degradation over decadal time scales, and the implications of our findings for targeting future habitat repair are discussed.     

Responding to flood risk in Louisiana: the roles of place attachment,emotions, and location

Natural Hazards - Drawing from protection motivation theory (PMT), we examined how place attachment and negative emotions, alongside threat and coping appraisals, personal experiences, and...     

An examination of midwestern US cities’ preparedness for climate change and extreme hazards

Disasters continue to have a dramatic impact on lives, livelihoods and environments communities depend on. In response to these losses, the global community has developed various theories, assessment methodologies and policies aimed at reducing global losses. A contemporary outcome of these interventions is to build the disaster resilience. However, despite the disaster resilience-building endeavours espoused by policies, theories and methodologies, very little progress is being made in reducing disaster losses. This paper argues that a possible reason behind the limitations of current resilience-building policies and methodologies could be that most of these policies are based a mechanistic scientific paradigm that places an emphasis on system components that are perceived to build resilience and not the function of systems as a whole. This often leads to resilience-building initiatives that are based on a ‘one-size-fits-all’ approach. This paper argues for the use of a complex adaptive systems approach to building resilience. This approach argues that contextual factors within different social systems will have a nonlinear affect on disaster resilience-building efforts. Therefore, it is crucial to move away from ‘one-size-fits-all’ approaches to more flexible approaches to building resilience. These hypotheses are tested by means of a correlation statistical analysis of agricultural communities in Southern Africa. Results of this analysis indicate that unique resilience profiles are evident in almost all of the communities studied. This indicates that resilience is not the same for everybody, and that resilience-building endeavours should be flexible enough to be adapted for different contexts.     

A Meta‐Analysis on the Return on Investment of Geospatial Data and Systems: A Multi‐Country Perspective

The availability, quality and accessibility of Geographic Information (GI) have significant socio‐economic and environmental benefits, but the collection and maintenance of GI require substantial investments. Cost‐benefit assessments (CBAs) attempt to justify the costs of geospatial data investments, applying different methodologies and focusing on diverse areas. Therefore, the Returns on Investment (ROI) vary considerably across studies, regions and sectors. The objective of this study is to explain some of the variation in the average ROI of GI by conducting a meta‐analysis of 82 cost‐benefit assessments between 1994 and 2013. In a first step, CBAs are systematically reviewed and relevant information is extracted. Particular emphasis is given to investment conditions and study characteristics. In a second step, multivariate regression methods are used to assess the size, significance and direction of individual effects. The results suggest that regional factors have the largest impact on the profitability of GI. Returns in Australia and New Zealand, for example, are four times larger than in Europe. In addition, small‐scale regional investments have a 2.5 times lower return than large‐scale international investments. Overall, the expected benefits of GI investments are approximately 3.2 times larger than the costs.     

Alternative measures of dispersion applied to flow in a convoluted channel

Steady flow in a convoluted channel is studied via Particle Tracking Velocimetry. The channel is constructed from a sequence of closed parallel cylindrical tubes welded together in plane which are then sliced down the lateral mid-plane and the lower complex is laterally shifted relative to the upper complex. Flow is induced in the lateral direction normal to the axis of the tubes. The a-time, T_a, finite-size Lyapunov exponent, λ_a, and the real-space self- and distinct-part of the intermediate scattering functions, G_s and G_d, and the pair density function, G_p, are computed from the data. Particle trajectories, velocity maps and streamlines show the channel has two prominent recirculation zones and a main flow region. The first passage time probability density function of tagged particles past a plane transverse to the mean flow illustrates how particles are delayed by recirculation zones. The delay caused by fluid element folding is manifested in single particle statistics such as the first passage time and the slowing increase in horizontal evolution of G_s. G_p describes initial particle distribution and allows areas in the flow domain trapping particles to be identified and visualized. G_d shows the evolution of the average separation of pairs of particles and when examined in a recirculation zone, it evolves little because of fluid element rotation. λ_a gives information on what transpires at a fixed scale and provides an estimate of the rate at which particles initially separated by a distance x separate to a distance ax as opposed to G_d which allows one to view changes over time. At small separations, λ_1.3 approaches a constant and for intermediate separations it scales as x^−0.8.