LNAPL Volume Calculation: Parameter Estimation by Nonlinear Regression of Saturation Profiles

An estimation of the volume of light nonaqueous phase liquids (LNAPL) is often required during site assessment, remedial design, or litigation. LNAPL volume can be estimated by a strictly empirical approach whereby core samples, distributed throughout the vertical and lateral extent of LNAPL, are analyzed for LNAPL content, and these data are then integrated to compute a volume. Alternatively, if the LNAPL has obtained vertical equilibrium, the thickness of LNAPL in monitoring wells can be used to calculate of LNAPL in monitoring wells can be used to calculate LNAPL volume at the well locations if appropriate soil and LNAPL properties can be estimated.
A method is described for estimating key soil and LNAPL properties by nonlinear regression of vertical profiles of LNAPL saturation. The methods is relatively fast, cost effective, and amenable to quantitative analysis of uncertainty. Optionally, the method allows statistical determination of best-fit values for the Van Genuchten capillary parameters (n, α_oil-water and α_oil-air), residual water saturation and ANAPL density. The sensitivity of the method was investigated by fitting field LNAPL saturation profiles and then determining the variation in misfit (mean square residual) as a function of parameter value for each parameter. Using field data from a sandy aquifer, the fitting statistics were found to be highly sensitive to LNAPL density, α_oil-water and α_oil-air moderately sensitive to the Van Genuchten n value, and weakly sensitive to residual water saturation. The regression analysis also provides information that can be used to estimate uncertainty in the estimated parameters, which can then be used to estimate uncertainty in calculated values of specific volume.     

Biological and Physical Attenuation of Endocrine Disruptors and Pharmaceuticals: Implications for Water Reuse

A select group of endocrine disrupters, pharmaceuticals, and personal care products was studied to determine the degree of biological attenuation in water reuse applications. Laboratory investigations involved both batch reactors using biologically active sand and continuous flow simulated aquifer storage and recovery experiments. All laboratory experiments were conducted using Colorado River water spiked with various target compounds at concentrations between 10 and 100 ng/L. Field studies were also conducted to determine the occurrence and attenuation of target compounds in water reuse applications. Two golf courses irrigated with reuse water were studied to determine if turf applications led to contamination of nearby ground water. A waste water treatment facility that uses rapid infiltration basins seasonally was also tested to determine the degree of attenuation of detectable target compounds along a subsurface flowpath. A qualitative structural activity relationship model was applied to the target compounds to predict the general rate of aerobic biological degradation. Significant attenuation of many target compounds was observed in both laboratory and field experiments. Conversely, several compounds displayed limited removal during these studies. Field experiments were limited to detectable compounds and various nonbiological removal or concentration effects that may impact data interpretations, which are discussed in this paper. The predictive model was found to be moderately accurate within the confines of the project scope.     

Fluoride-Selective Electrode as a Tool to Evaluate the Degradation of PFAS in Groundwater: A Bench-Scale Investigation

Fluoride (F⁻) concentration in groundwater can be used as an independent measure of the degradation of per- and polyfluoroalkyl substances (PFAS) at contaminated sites. This study assessed the impact of the groundwater matrix, oxidizing reagents, quenching agent, and sample handling procedures (filtration) on determinations of F⁻. F⁻ was quantified using a fluoride-selective electrode (FSE) using matrix spike recovery (acceptable range 75 to 125%) and electrode slope performance (ideal range −54 to −60 mV). For the unaltered aqueous matrices considered (ultrapure water, ultrapure water with sodium chloride, and simulated groundwater), matrix spike recoveries were greater than 93% and acceptable electrode slopes (−59.5 ± 0.2 mV) were obtained. Matrix spike recoveries and electrode slopes when persulfate or permanganate was added were 87 and 93%, and −58.4 ± 1.3 and −47.3 ± 2.2 mV, respectively. Therefore, quantifying F⁻ in matrices containing permanganate should involve use of either a matrix-matched calibration curve or the method of standard additions due to the low electrode slope. The addition of ascorbic acid as a quenching agent resulted in a decrease in matrix spike recoveries to <74% and suggests alternate quenching procedures should be used if an FSE is employed to measure F⁻. Removal of sediments by filtration prior to the addition of the fluoride spike led to improved matrix spike recoveries (>96%), which were otherwise underestimated using the FSE in the presence of sediments. While the FSE may be a quick and portable tool, its significant limitations need to be fully understood before it can be used to quantify the production of F⁻ resulting from the potential degradation of PFAS at field sites.     

Large grains matter: contrasting bed stability and morphodynamics during two nearly identical experiments

While the stabilizing function of large grains in step‐pool streams has long been recognized, the role they play in gravel‐bed streams is less clear. Most researchers have ignored the role of large grains in gravel‐bed streams, and have assumed that the median bed surface size controls the erodibility of alluvial boundaries. The experiments presented herein challenge this convention. Two experiments were conducted that demonstrate the significant morphodynamic implications of a slight change to the coarse tail of the bed material. The two distributions had the same range of particle sizes, and nearly identical bulk d₅₀ values (1.6 mm); however the d₉₀ of experiment GSD1 was slightly finer (3.7 mm) than that for experiment GSD2 (3.9 mm). Transport rates during GSD1 were nearly four times greater than during GSD2 (even though the dimensionless shear stress was slightly lower), and the channel developed a sinuous pattern with well‐developed riffles, pools and bars. During GSD2 the initial rectangular channel remained virtually unchanged for the duration of the experiment. The relative stability of GSD2 seems to be associated with a slightly larger proportion of stable (large) grains on the bed surface: at the beginning of GSD1, 3.5% of the bed was immobile, while almost twice as much of it (6.1%) was immobile at the beginning of GSD2. The results demonstrate that the largest grains (not the median size) exert first‐order control on channel stability. Copyright © 2017 John Wiley & Sons, Ltd.     

A decadal-scale numerical model for wandering,cobble-bedded rivers subject to disturbance

Alluvial rivers are composed of self-formed channels which are sensitive to disturbances in their flow and sediment-supply regimes. Regime changes commonly occur over decadal and longer timescales and can be caused by anthropogenic alterations such as dam construction and removal. Advances in numerical modeling have increased our ability to explore geomorphic adjustments over long timescales; however, many models designed to be run for decades or longer assume that banks are immovable or that channel width is constant. Since river channels often respond to disturbance by adjusting their geometry, this is a significant shortcoming. To investigate the impact of long-term sediment supply alterations on channel geometry and stability, we have adapted MAST-1D, a reach-scale bed evolution model, to incorporate functions for bank erosion, vegetation encroachment, and local avulsions. The model is designed for medium-large, coarse multithreaded rivers and can be run over long (decades–centuries) timescales. Bank erosion is a function of the mobility and transport capacity for structurally-important grains which protect the bank toe. Vegetation growth is proportional to point bar width and occurs during conditions of low shear stress. Local avulsions occur when aggradation causes channel depth to drop below a threshold. We apply the model to the Elwha River in Washington, USA with the goal of investigating if and when the river recovers from dam emplacement and removal. The Elwha was dammed for nearly 100 years, and then two dams were removed, releasing a large pulse of sediment. We have modeled the set of reaches between the two dams. Our simulations suggest that channel response to dam emplacement occurs gradually over several decades but that the channel recovers to near pre-dam conditions within about a decade following the removal. The dams leave a lasting legacy on the floodplain, which does not completely recover, even after two centuries. © 2019 John Wiley & Sons, Ltd.     

Architecture and preservation in the fluvial to marine transition zone of a mixed-process humid-tropical delta: Middle Miocene Lambir Formation,Baram Delta Province,north-west Borneo

The interaction of river and marine processes in the fluvial to marine transition zone fundamentally impacts delta plain morphology and sedimentary dynamics. This study aims to improve existing models of the facies distribution, stratigraphic architecture and preservation in the fluvial to marine transition zone of mixed-process deltas, using a comprehensive sedimentological and stratigraphic dataset from the Middle Miocene Lambir Formation, Baram Delta Province, north-west Borneo. Eleven facies associations are identified and interpreted to preserve the interaction of fluvial and marine processes in a mixed-energy delta, where fluvial, wave and tidal processes display spatially and temporally variable interactions. Stratigraphic successions in axial areas associated with active distributary channels are sandstone-rich, comprising fluvial-dominated and wave-dominated units. Successions in lateral areas, which lack active distributary channels, are mudstone-rich, comprising fluvial-dominated, tide-dominated and wave-dominated units, including mangrove swamps. Widespread mudstone preservation in axial and lateral areas suggests well-developed turbidity maximum zones, a consequence of high suspended-sediment concentrations resulting from tropical weathering of a mudstone-rich hinterland. Within the fluvial to marine transition zone of distributary channels, interpreted proximal–distal sedimentological and stratigraphic trends suggest: (i) a proximal fluvial-dominated, tide-influenced subzone; (ii) a distal fluvial-dominated to wave-dominated subzone; and (iii) a conspicuously absent tide-dominated subzone. Lateral areas preserve a more diverse spectrum of facies and stratigraphic elements reflecting combined storm, tidal and subordinate river processes. During coupled storm and river floods, fluvial processes dominated the fluvial to marine transition zone along major and minor distributary channels and channel mouths, causing significant overprinting of preceding interflood deposits. Despite interpreted fluvial–tidal channel units and mangrove influence implying tidal processes, there is a paucity of unequivocal tidal indicators (for example, cyclical heterolithic layering). This suggests that process preservation in the fluvial to marine transition zone preserved in the Lambir Formation primarily records episodic (flashy) river discharge, river flood and storm overprinting of tidal processes, and possible backwater dynamics.     

Symplectic integrators for long-term integrations in celestial mechanics

We discuss the use of symplectic integration algorithms in long-term integrations in the field of celestial mechanics. The methods' advantages and disadvantages (with respect to more common integration methods) are discussed. The numerical performance of the algorithms is evaluated using the 2-body and circular restricted 3-body problems. Symplectic integration methods have the advantages of linear phase error growth in the 2-body problem (unlike most other methods), good conservation of the integrals of the motion, good performance for moderately eccentric orbits, and ease of use. Its disadvantages include a relatively large number of force evaluations and an inability to continuously vary the step size.     

A final non-redundant catalogue for the 7C 151-MHz survey

We present a final unified catalogue for the 7C survey at 151 MHz with resolution  70 × 70 cosec(δ) arcsec²  . This has been constructed by amalgamating the existing catalogues derived from individual fields imaged at this resolution and eliminating redundancy in regions of mutual overlap. This is a non-trivial procedure because the flux in multiple-component sources may be fitted differently on alternative images, owing, for example, to differences in local noise and beam distortion. The final catalogue lists 43683 sources over an area of about 1.7 sr. Separate final catalogues have been published for the 7C Galactic Plane survey (7CG) and the lower-resolution survey of the low-declination strip  9^h < RA < 16^h, 20° < Dec < 35°  .     

Aircraft observations of Venus' near-infrared reflection spectrum: Implications for cloud composition

We have obtained measurements of Venus' reflection spectrum in the 1.2 to 4.1-μm spectral region from a NASA-Ames operated Lear jet. This was accomplished by observing both Venus and the sun with a spectrometer that contained a circular, variable interference filter, whose effective spectral resolution was 2%. The aircraft results were compared with computer generated spectra of a number of cloud candidates. The only substance which gave an acceptable match to the profile of Venus' strong 3-μm absorption feature, was a water solution of sulfuric acid, that had a concentration of 75% or more H₂SO₄ by weight. However, our spectra also show a modest decline in reflectivity from 2.3 μm towards 1.2-μm wavekength, which is inconsistent with the flat spectrum of sulfuric acid in this spectral region. We hypothesize that this decline is due to impurities in the sulfuric acid droplets.We also compared our list of cloud candidates with several other observed properties of the Venus clouds. While this comparison does not provide as unique an answer as did our analysis of the 3-μm band, we find that, in agreement with the results of Young (1973) and Sill (1973), concentrated sulfuric acid solutions are compatible with these additional observed properties of the Venus clouds. We conclude that the visible cloud layer of Venus is composed of sulfuric acid solution droplets, whose concentration is 75% H₂SO₄, or greater, by weight.