Model prediction capacity of ephemeral gully evolution in conservation tillage systems

Ephemeral gully (EG) erosion has an important impact on agricultural soil losses and increases field surface hydrology connectivity and transport of pollutants to nearby water bodies. Watershed models including an EG component are scarce and not yet properly evaluated. The objective of this study is to evaluate the capacity of one such tool, AnnAGNPS, to simulate the evolution of two EG formed in a conservation tillage system. The dataset for model testing included runoff measurements and EG morphological characteristics during 3 years. Model evaluation focused on EG evolution of volume, width, and length model outputs, and included calibration and testing phases and a global sensitivity analysis (GSA). While the model did not fully reproduce width and length, the model efficiency to simulate EG volume was satisfactory for both calibration and testing phases, supporting the watershed management objectives of the model. GSA revealed that the most sensitive factors were EG depth, critical shear stress, headcut detachment exponent coefficient b, and headcut detachment leading coefficient a. For EG outputs the model was additive, showing low sensitivity to interactions between the inputs. Prediction of EG spatial evolution on conservation tillage systems requires improved development of gully erosion components, since many of the processes were developed originally for traditional tillage practices or larger channel systems. Our results identify the need for future research when EG form within conservation tillage systems, in particular to study gully headcut, soil erodibility, and width functions specific to these practices.     

Modeling the low-latitude ionospheric electron density and plasma turbulence in the November 2004 storm period

The storm period of 8–12 November 2004 offers an opportunity for insight into the phenomena of low-latitude ionospheric structure during geomagnetically disturbed times because of the strength of the disturbances, the timing of the storms, and the instrumentation that was operating during the interval. We will take advantage of these factors to model the ambient ionosphere and the plasma turbulence responsible for radio scintillation within it, using the AFRL low-latitude ambient/turbulent ionospheric model and the storm-time model features described in the companion paper [Retterer, J.M., Kelley, M.C., 2009. Solar-wind drivers for low-latitude ionospheric models during geomagnetic storms. J. Atmos. Solar-Terr. Phys., this issue]. The model plasma densities show very good agreement with the densities measured by the Jicamarca ISR as well as with the total electron content (TEC) measured by the Boston College South American chain of GPS receivers. The detection by the radar of coherent returns from plasma turbulence match well the times of predicted ionospheric instability. The predicted geographic extent of the occurrence of equatorial plasma bubbles was matched by DMSP satellite observations and our forecasts of scintillation strength were validated with measurements of S4 at Ancon and Antofagasta by stations of the AFRL SCINDA network.     

Numerical Modeling of Thermal Conductive Heating in Fractured Bedrock

Numerical modeling was employed to study the performance of thermal conductive heating (TCH) in fractured shale under a variety of hydrogeological conditions. Model results show that groundwater flow in fractures does not significantly affect the minimum treatment zone temperature, except near the beginning of heating or when groundwater influx is high. However, fracture and rock matrix properties can significantly influence the time necessary to remove all liquid water (i.e., reach superheated steam conditions) in the treatment area. Low matrix permeability, high matrix porosity, and wide fracture spacing can contribute to boiling point elevation in the rock matrix. Consequently, knowledge of these properties is important for the estimation of treatment times. Because of the variability in boiling point throughout a fractured rock treatment zone and the absence of a well-defined constant temperature boiling plateau in the rock matrix, it may be difficult to monitor the progress of thermal treatment using temperature measurements alone.     

Methodology for comparing source and plume remediation alternatives

It is often difficult at contaminated sites to decide whether remediation effort should be focused on the contaminant source, the dissolved plume, or on both zones. The decision process at these sites is hampered by a lack of quantitative tools for comparing remediation alternatives. A new screening-level mass balance approach is developed for simulating the transient effects of simultaneous ground water source and plume remediation. The contaminant source model is based on a power function relationship between source mass and source discharge, and it can consider partial source remediation at any time after the initial release. The source model serves as a time-dependent mass flux boundary condition to a new analytical plume model, where flow is assumed to be one dimensional, with three-dimensional dispersion. The plume model simulates first-order sequential decay and production of several species, and the decay rates and parent/daughter yield coefficients are variable functions of time and distance. This new method allows for flexible simulation of natural attenuation or remediation efforts that enhance plume degradation. The plume remediation effort may be temporary or delayed in time, limited in space, and it may have different chemical effects on different contaminant species in the decay chain.     

Modelling spatial variations in dissolved oxygen in fine‐grained streams under uncertainty

This paper presents a novel triple‐layer model, called VART DO‐3L, for simulation of spatial variations in dissolved oxygen (DO) in fine‐grained streams, characterized by a fluid mud (fluff or flocculent) layer (an advection‐dominated storage zone) as the interface between overlying stream water and relatively consolidated streambed sediment (a diffusion‐dominated storage zone). A global sensitivity analysis is conducted to investigate the sensitivity of VART DO‐3L model input parameters. Results of the sensitivity analysis indicate that the most sensitive parameter is the relative size of the advection‐dominated storage zones (A_s/A), followed by a lumped reaction term (R) for the flocculent layer, biological reaction rate (μ_o) in diffusive layer and biochemical oxygen demand concentration (L) in water column. In order to address uncertainty in model input parameters, Monte Carlo simulations are performed to sample parameter values and to produce various parameter combinations or cases. The VART DO‐3L model is applied to the Lower Amite River in Louisiana, USA, to simulate vertical and longitudinal variations in DO under the cases. In terms of longitudinal variation, the DO level decreases from 7.9 mg l at the Denham Springs station to about 2.89 mg l^?1 at the Port Vincent station. In terms of vertical variation, the DO level drops rapidly from the overlying water column to the advection‐dominated storage zone and further to the diffusive layer. The DO level (C_F) in the advective layer (flocculent layer) can reach as high as 40% of DO concentration (C) in the water column. The VART DO‐3L model may be applied to similar rivers for simulation of spatial variations in DO level. Copyright © 2014 John Wiley & Sons, Ltd.     

Experimental synthesis of seismic horizontal free-field motion of soil in finite-domain simulations with absorbing boundary

In this paper, a fundamental assessment of the method of physical wave-absorbing boundary and centrifuge modeling is presented in the context of experimental simulations of seismic free-field ground motion. Focusing on the characteristics of a sand stratum, a series of seismic tests on models of uniform density and a large width-to-depth ratio with Duxseal as the side boundary were performed using an in-box shake-table system. By means of the transfer function approach in the frequency domain, the complex three-dimensional nature of the dynamic response of the finite soil model with the boundary treatment is demonstrated in terms of its variable resonant frequency distribution at different g-levels. Apart from being helpful in quantifying the difficulty in interpreting the finite-domain response simulations using one-dimensional theories or homogenized representations, the measured data substantiates the need and usefulness of coupling the Duxseal boundary approach with a three-dimensional elastodynamic synthesis. With the aid of a corresponding boundary element implementation, the feasibility of identifying the soil's in-flight shear modulus variation, Poisson's ratio and horizontal-to-vertical earth pressure ratio from the centrifuge model's free-field measurements is also explored.     

A comparison of ARM and TRM in magnetite

Experiments comparing anhysteretic remanence (ARM) and thermoremanence (TRM) in samples containing natural and synthetic magnetite, whose mean particle sizes range from single domain to multidomain, show that ARM and TRM are very similar (but not identical) in their stabilities with respect to alternating field (AF) demagnetization, temperature cycles in zero field to below magnetite's isotropic temperature near 130°K, and stability with respect to spontaneous decay in zero field. Therefore, for magnetites, ARM can be used to model (with reasonable success) these stability properties of TRM. The field dependence of the acquisition of ARM and TRM shows that the low field susceptibility ratio, χ_ARM/χ_TRM, has a particle size dependence, increasing from 0.1 for certain submicron particles to 2.0 for large multidomain crystals. Even for samples whose remanence is predominantly carried by submicron particles χ_ARM/χ_TRM is highly variable, 0.11 ≤ χ_ARM/χ_TRM ≤ 0.50. Therefore, ARM paleointensity methods which do not take into account the large variability in and the particle size dependence of χ_ARM/χ_TRM are subject to order-of-magnitude uncertainties.