Regolith production and transport at the Susquehanna Shale Hills Critical Zone Observatory, Part 2: Insights from meteoric 10Be

Regolith-mantled hillslopes are ubiquitous features of most temperate landscapes, and their morphology reflects the climatically, biologically, and tectonically mediated interplay between regolith production and downslope transport. Despite intensive research, few studies have quantified both of these mass fluxes in the same field site. Here we present an analysis of 87 meteoric 10Be measurements from regolith and bedrock within the Susquehanna Shale Hills Critical Zone Observatory (SSHO), in central Pennsylvania. Meteoric 10Be concentrations in bulk regolith samples (n=73) decrease with regolith depth. Comparison of hillslope meteoric 10Be inventories with analyses of rock chip samples (n=14) from a 24 m bedrock core confirms that >80% of the total inventory is retained in the regolith. The systematic downslope increase of meteoric 10Be inventories observed at SSHO is consistent with 10Be accumulation in slowly creeping regolith (∼ 0.2 cm yr-1). Regolith flux inferred from meteoric 10Be varies linearly with topographic gradient (determined from high-resolution light detection and ranging-based topography) along the upper portions of hillslopes at SSHO. However, regolith flux appears to depend on the product of gradient and regolith depth where regolith is thick, near the base of hillslopes. Meteoric 10Be inventories at the north and south ridgetops indicate minimum regolith residence times of 10.5 ± 3.7 and 9.1 ± 2.9 ky, respectively, similar to residence times inferred from U-series isotopes in Ma et al. (2013). The combination of our results with U-series-derived regolith production rates implies that regolith production and erosion rates are similar to within a factor of two on SSHO hillcrests. ©2013. American Geophysical Union. All Rights Reserved

Asymptotic Analysis of Time-Dependent Neutron Transport Coupled with Isotopic Depletion and Radioactive Decay

We describe an asymptotic analysis of the coupled nonlinear system of equations describing time-dependent three-dimensional monoenergetic neutron transport and isotopic depletion and radioactive decay. The classic asymptotic diffusion scaling of Larsen and Keller [1], along with a consistent small scaling of the terms describing the radioactive decay of isotopes, is applied to this coupled nonlinear system of equations in a medium of specified initial isotopic composition. The analysis demonstrates that to leading order the neutron transport equation limits to the standard time-dependent neutron diffusion equation with macroscopic cross sections whose number densities are determined by the standard system of ordinary differential equations, the so-called Bateman equations, describing the temporal evolution of the nuclide number densities

Angularly Adaptive P1--Double P0 Diffusion Solutions of Non-Equilibrium Grey Radiative Transfer Problems in Planar Geometry

The double spherical harmonics angular approximation in the lowest order, i.e. double P{sub 0} (DP{sub 0}), is developed for the solution of time-dependent non-equilibrium grey radiative transfer problems in planar geometry. The standard P{sub 1} angular approximation represents the angular dependence of the radiation specific intensity using a linear function in the angular domain -1 {le} {mu} {le} 1. In contrast, the DP{sub 0} angular approximation represents the angular dependence as isotropic in each half angular range -1 {le} {mu} < 0 and 0 < {mu} {le} 1. Neglecting the time derivative of the radiation flux, both the P{sub 1} and DP{sub 0} equations can be written as a single diffusion equation for the radiation energy density. Although the DP{sub 0} diffusion approximation is expected to be less accurate than the P{sub 1} diffusion approximation at and near thermodynamic equilibrium, the DP{sub 0} angular approximation can more accurately capture the complicated angular dependence near the non-equilibrium wave front. We develop an adaptive angular technique that locally uses either the DP{sub 0} or the P{sub 1} diffusion approximation depending on the degree to which the radiation and material fields are in thermodynamic equilibrium. Numerical results are presented for a test problem due to Su and Olson for which a semi-analytic transport solution exists. The numerical results demonstrate that the adaptive P{sub 1}-DP{sub 0} diffusion approximation can yield improvements in accuracy over the standard P{sub 1} diffusion approximation for non-equilibrium grey radiative transfer

Angularly Adaptive P1 - Double P0 Flux-Limited Diffusion Solutions of Non-Equilibrium Grey Radiative Transfer Problems

The double spherical harmonics angular approximation in the lowest order, i.e. double P{sub 0} (DP{sub 0}), is developed for the solution of time-dependent non-equilibrium grey radiative transfer problems in planar geometry. Although the DP{sub 0} diffusion approximation is expected to be less accurate than the P{sub 1} diffusion approximation at and near thermodynamic equilibrium, the DP{sub 0} angular approximation can more accurately capture the complicated angular dependence near a non-equilibrium radiation wave front. In addition, the DP{sub 0} approximation should be more accurate in non-equilibrium optically thin regions where the positive and negative angular domains are largely decoupled. We develop an adaptive angular technique that locally uses either the DP{sub 0} or P{sub 1} flux-limited diffusion approximation depending on the degree to which the radiation and material fields are in thermodynamic equilibrium. Numerical results are presented for two test problems due to Su and Olson and to Ganapol and Pomraning for which semi-analytic transport solutions exist. These numerical results demonstrate that the adaptive P{sub 1}-DP{sub 0} diffusion approximation can yield improvements in accuracy over the standard P{sub 1} diffusion approximation, both without and with flux-limiting, for non-equilibrium grey radiative transfer

A Perspective on the Recent Progress in Solution-processed Methods for Highly Efficient Perovskite Solar Cells

Perovskite solar cells (PSCs) were developed in 2009 and have led to a number of significant improvements in clean energy technology. The power conversion efficiency (PCE) of PSCs has increased exponentially and currently stands at 22%. PSCs are transforming photovoltaic (PV) technology, outpacing many established PV technologies through their versatility and roll-to-roll manufacturing compatibility. The viability of low-temperature and solution-processed manufacturing has further improved their viability. This article provides a brief overview of the stoichiometry of perovskite materials, the engineering behind various modes of manufacturing by solution processing methods, and recommendations for future research to achieve large-scale manufacturing of high efficienc

Regolith production and transport at the Susquehanna Shale Hills Critical Zone Observatory, part 2: Insights from meteoric10Be

Regolith-mantled hillslopes are ubiquitous features of most temperate landscapes, and their morphology reflects the climatically, biologically, and tectonically mediated interplay between regolith production and downslope transport. Despite intensive research, few studies have quantified both of these mass fluxes in the same field site. Here we present an analysis of 87 meteoric 10Be measurements from regolith and bedrock within the Susquehanna Shale Hills Critical Zone Observatory (SSHO), in central Pennsylvania. Meteoric 10Be concentrations in bulk regolith samples (n=73) decrease with regolith depth. Comparison of hillslope meteoric 10Be inventories with analyses of rock chip samples (n=14) from a 24 m bedrock core confirms that >80% of the total inventory is retained in the regolith. The systematic downslope increase of meteoric 10Be inventories observed at SSHO is consistent with 10Be accumulation in slowly creeping regolith (∼ 0.2 cm yr-1). Regolith flux inferred from meteoric 10Be varies linearly with topographic gradient (determined from high-resolution light detection and ranging-based topography) along the upper portions of hillslopes at SSHO. However, regolith flux appears to depend on the product of gradient and regolith depth where regolith is thick, near the base of hillslopes. Meteoric 10Be inventories at the north and south ridgetops indicate minimum regolith residence times of 10.5 ± 3.7 and 9.1 ± 2.9 ky, respectively, similar to residence times inferred from U-series isotopes in Ma et al. (2013). The combination of our results with U-series-derived regolith production rates implies that regolith production and erosion rates are similar to within a factor of two on SSHO hillcrests

Variable Curvature Slab Molecular Dynamics as a Method to Determine Surface Stress

A thin plate or slab, prepared so that opposite faces have different surface stresses, will bend as a result of the stress difference. We have developed a classical molecular dynamics (MD) formulation where (similar in spirit to constant-pressure MD) the curvature of the slab enters as an additional dynamical degree of freedom. The equations of motion of the atoms have been modified according to a variable metric, and an additional equation of motion for the curvature is introduced. We demonstrate the method to Au surfaces, both clean and covered with Pb adsorbates, using many-body glue potentials. Applications to stepped surfaces, deconstruction and other surface phenomena are under study.Comment: 16 pages, 8 figures, REVTeX, submitted to Physical Review

Mars Sample Return: The Value of Depth Profiles

Sample return from Mars offers the promise of data from Martian materials that have previously only been available from meteorites. Return of carefully selected samples may yield more information about the history of water and possible habitability through Martian history. Here we propose that samples collected from Mars should include depth profiles of material across the interface between weathered material on the surface of Mars into unweathered parent rock material. Such profiles have the potential to yield chemical kinetic data that can be used to estimate the duration of water and information about potential habitats on Mars