Joint isotopic mass balance: a novel approach to quantifying channel bed to channel margins sediment transfer during storm events

The important role of floodplains and the broader riparian zone in providing temporary storage for a large fraction of the annual sediment load of rivers is well established, but this understanding is largely based on observations of the long‐term average behavior of the catchment. Here we combine measurements of the fallout radionuclides ⁷Be and ²¹⁰Pb and the stable isotopes of hydrogen in water to quantify fine sediment mobilization and storage in a stream and its channel margins during individual intermediate‐sized storm events with recurrence intervals of a few months or less. We demonstrate this method using five storm events in a small (~15 km²), undeveloped, gravel‐bedded tributary of the Connecticut River (USA). We estimate that in each storm, the mass of sediment deposited onto the margins accounts for almost 90% of the sediment mobilized from the bed, with the remainder of the mobilized bed sediment transported downstream as suspended load. The result that the bed is a net source of sediment to the stream and the margins a net sink is robust, but estimates of the mass of material eroded from the bed and deposited on the margins are less certain. The source of sediment to the bed remains unclear as, consistent with earlier studies, we observe only limited deposition of sediment to the bed during the storm events. The suspended sediment is organic‐rich and thus its source may be associated with in‐channel organic decay between storm events. Understanding the coupled interactions between discharge magnitude and frequency and sediment resupply at the event time scale has important implications for stream restoration efforts seeking to connect the channel and the broader riparian zone, and for the development of accurate sediment budgets and predictions of sediment flux from a watershed. Copyright © 2015 John Wiley & Sons, Ltd.     

Improving Groundwater Model Calibration with Repeat Microgravity Measurements

Groundwater-flow models depend on hydraulic head and flux observations for evaluation and calibration. A different type of observation—change in storage measured using repeat microgravity—can also be used for parameter estimation by simulating the expected change in gravity from a groundwater model and including the observation misfit in the objective function. The method is demonstrated using new software linked to MODFLOW input and output files and field data from the vicinity of the All American Canal in southeast California, USA. Over a 10-year period following lining of the previously highly permeable canal with concrete, gravity decreased by over 100 μGal (equivalent to about 2.5 m of free-standing water) at some locations as seepage decreased and the remnant groundwater mound dissipated into the aquifer or was removed by groundwater pumping. Simulated gravity from a MODFLOW model closely matched observations, and repeat microgravity data proved useful for constraining both hydraulic conductivity and specific yield estimates. Specific yield estimated using the infinite-horizontal slab approximation agreed well with model-derived values, and the departure from the linear, flat-water-table approximation was small, less than 2%, despite relatively large and dynamic water-table slope. First-order second-moment parameter uncertainty analysis shows reduction in uncertainty for all hydraulic conductivity and specific yield parameter estimates with the addition of repeat microgravity data, as compared to drawdown data alone.     

Fishing dynamics associated with periodically harvested marine closures

Periodically-harvested fisheries closures are emerging as a socially acceptable and locally implementable way to balance concerns about conserving ecosystem function and sustaining livelihoods. Across the Indo-Pacific periodically-harvested closures are commonly employed, yet their contribution towards more sustainable fisheries remains largely untested in the social and ecological context of tropical small-scale fisheries. To address this, we use an interdisciplinary approach to examine harvesting dynamics that would affect sustainability, namely, fishing effort, yield, gear and method use, periodicity of harvesting, controls placed on harvesting and resource owners’ decisions to open and close four fishing grounds in Solomon Islands. We compare these fishing patterns with those on surrounding, continuously open fishing grounds. Our study shows that total effort and total catch from periodically-harvested reef closures are low to moderate compared to reefs open to continuous fishing. When periodically-harvested closures were opened, effort in the closures was relatively intense, however, in most cases yield did not exceed annual benchmarks of sustainability described by previous studies. In some cases, harvesting during openings was restricted to a single taxon and single fishing gear and method, while in others there was unrestricted multi-species and multi-method harvesting. The duration and frequency of openings were highly variable, with open periods ranging from a single night to one month in duration, and occurring between one and 15 times per year. Fishing during openings was permitted for entire fishing communities in some cases, and only for specific rights-holding families in others. Decisions to open periodically-harvested closures tended to be based on immediate social or economic needs, and the openings provided a small boost to fish catch landed in communities. While periodically-harvested closures may alleviate fishing pressure in a small area of fishing grounds by reducing the opportunity to fish, openings of long duration or high frequency, combined with heavy or destructive exploitation, may lead to unsustainable harvesting within the area.     

A simple carbon cycle representation for economic and policy analyses

Integrated Assessment Models (IAMs) that couple the climate system and the economy require a representation of ocean CO₂ uptake to translate human-produced emissions to atmospheric concentrations and in turn to climate change. The simple linear carbon cycle representations in most IAMs are not however physical at long timescales, since ocean carbonate chemistry makes CO₂ uptake highly nonlinear. No linearized representation can capture the ocean’s dual-mode behavior, with initial rapid uptake and then slow equilibration over ∽10,000 years. In a business-as-usual scenario followed by cessation of emissions, the carbon cycle in the 2007 version of the most widely used IAM, DICE (Dynamic Integrated model of Climate and the Economy), produces errors of ∽2^°C by the year 2300 and ∽6^°C by the year 3500. We suggest here a simple alternative representation that captures the relevant physics and show that it reproduces carbon uptake in several more complex models to within the inter-model spread. The scheme involves little additional complexity over the DICE model, making it a useful tool for economic and policy analyses.     

Calculating the effective permeability of sandstone with multiscale lattice Boltzmann/finite element simulations

The lattice Boltzmann (LB) method is an efficient technique for simulating fluid flow through individual pores of complex porous media. The ease with which the LB method handles complex boundary conditions, combined with the algorithm’s inherent parallelism, makes it an elegant approach to solving flow problems at the sub-continuum scale. However, the realities of current computational resources can limit the size and resolution of these simulations. A major research focus is developing methodologies for upscaling microscale techniques for use in macroscale problems of engineering interest. In this paper, we propose a hybrid, multiscale framework for simulating diffusion through porous media. We use the finite element (FE) method to solve the continuum boundary-value problem at the macroscale. Each finite element is treated as a sub-cell and assigned permeabilities calculated from subcontinuum simulations using the LB method. This framework allows us to efficiently find a macroscale solution while still maintaining information about microscale heterogeneities. As input to these simulations, we use synchrotron-computed 3D microtomographic images of a sandstone, with sample resolution of 3.34 μm. We discuss the predictive ability of these simulations, as well as implementation issues. We also quantify the lower limit of the continuum (Darcy) scale, as well as identify the optimal representative elementary volume for the hybrid LB–FE simulations.     

Continuum deformation and stability analyses of a steep hillside slope under rainfall infiltration

Rainfall weakens an earth slope in a number of ways. It increases the degree of saturation of the soil, thereby breaking the bonds created by surface tension between the soil particles. When the volume of infiltrating water is large enough to mobilize fluid flow inside the soil matrix, the fluid exerts a downhill frictional drag on the slope, creating a destabilizing effect. When excess fluid can no longer infiltrate the slope due to increased saturation in the soil, it is discharged as a surface runoff and erodes the slope. In this paper, we present a physics-based framework for continuum modeling of a hydrologically driven slope failure similar to what occurred in a steep experimental catchment CB1 near Coos Bay, Oregon. We quantify the rainfall-induced slope deformation and assess the failure potential of the slope using finite element modeling that couples solid deformation with fluid pressure in an unsaturated soil. Results of the studies suggest that for a steep hillside slope underlain by a shallow bedrock similar to the CB1 site, failure would occur by multiple slide blocks with the failure surfaces emerging on the slope face. These results suggest that an infinite slope mechanism would be insufficient to represent the failure kinematics for a slope similar to CB1.     

Eclipsing binary Trojan asteroid Patroclus: Thermal inertia from Spitzer observations

We present mid-infrared observations of the binary L5-Trojan system (617) Patroclus-Menoetius before, during, and after two shadowing events, using the Infrared Spectrograph (IRS) on board the Spitzer Space Telescope. For the first time, we effectively observe changes in asteroid surface temperature in real time, allowing the thermal inertia to be determined very directly. A new detailed binary thermophysical model is presented which accounts for the system’s known mutual orbit, arbitrary component shapes, and thermal conduction in the presence of eclipses.We obtain two local thermal-inertia values, representative of the respective shadowed areas: and . The average thermal inertia is estimated to be , potentially with significant surface heterogeneity. This first thermal-inertia measurement for a Trojan asteroid indicates a surface covered in fine regolith. Independently, we establish the presence of fine-grained (<a few μm) silicates on the surface, based on emissivity features near 10 and similar to those previously found on other Trojans.We also report V-band observations and report a lightcurve with complete rotational coverage. The lightcurve has a low amplitude of peak-to-peak, implying a roughly spherical shape for both components, and is single-periodic with a period equal to the period of the mutual orbit, indicating that the system is fully synchronized.The diameters of Patroclus and Menoetius are 106±11 and , respectively, in agreement with previous findings. Taken together with the system’s known total mass, this implies a bulk mass density of , significantly below the mass density of L4-Trojan asteroid (624) Hektor and suggesting a bulk composition dominated by water ice.All known physical properties of Patroclus, arguably the best studied Trojan asteroid, are consistent with those expected in icy objects with devolatilized surface (extinct comets), consistent with what might be implied by recent dynamical modeling in the framework of the Nice Model.     

Core Density Profile Survival in Dissipationless Mergers

We use N-body experiments to study the survival of power-law density cusps in mergers of spherical galaxy models. Steep cusps (e.g. ρ∝r ^-2) are preserved during mergers of equal-mass models; this may be viewed as a general consequence of the incomplete violent relaxation observed in such mergers. Progenitor profile has a strong effect on remnant shape; steep cusps produce oblate remnants. If bright ellipticals form by hierarchical mergers of faint ones, the lack of steep cusps in these galaxies is somewhat surprising. This revised version was published online in July 2006 with corrections to the Cover Date.     

The prospects for life on Mars: A pre-Viking assessment

Mariner 9 has provided a refutation or reinterpretation of several historical claims for Martian biology, and has permitted an important further characterization of the environmental constraints on possible Martian organisms. Four classes of conceivable Martian organisms are identified, depending on the environmental temperature, T, and water activity, a_w: Class I, high T, high a_w; Class II, low T, high a_w; Class III, high T, low a_w; and Class IV, low T, low a_w. The Viking lander biology experiments are essentially oriented toward Class I organisms, although arguments are given for the conceivable presence on Mars of organisms in any of the four classes. Organisms which extract their water requirements from hydrated minerals or from ice are considered possible on Mars, and the high ultraviolet flux and low oxygen partial pressure are considered to be negligible impediments to Martian biology. Large organisms, possibly detectable by the Viking lander cameras, are not only possible on Mars; they may be favored. The surface distribution of Martian organisms and future search strategies for life on Mars are discussed.