Analytic Element Modeling of Steady Interface Flow in Multilayer Aquifers Using AnAqSim

This paper presents the analytic element modeling approach implemented in the software AnAqSim for simulating steady groundwater flow with a sharp fresh‐salt interface in multilayer (three‐dimensional) aquifer systems. Compared with numerical methods for variable‐density interface modeling, this approach allows quick model construction and can yield useful guidance about the three‐dimensional configuration of an interface even at a large scale. The approach employs subdomains and multiple layers as outlined by Fitts (2010) with the addition of discharge potentials for shallow interface flow (Strack 1989). The following simplifying assumptions are made: steady flow, a sharp interface between fresh‐ and salt water, static salt water, and no resistance to vertical flow and hydrostatic heads within each fresh water layer. A key component of this approach is a transition to a thin fixed minimum fresh water thickness mode when the fresh water thickness approaches zero. This allows the solution to converge and determine the steady interface position without a long transient simulation. The approach is checked against the widely used numerical codes SEAWAT and SWI/MODFLOW and a hypothetical application of the method to a coastal wellfield is presented.     

Flow‐Through Stream Modeling with MODFLOW and MT3D: Certainties and Limitations

This paper aims to assess MODFLOW and MT3D capabilities for simulating the spread of contaminants from a river exhibiting an unusual relationship with an alluvial aquifer, with the groundwater head higher than the river head on one side and lower on the other (flow‐through stream). A series of simulation tests is conducted using a simple hypothetical model so as to characterize and quantify these limitations. Simulation results show that the expected contaminant spread could be achieved with a specific configuration composed of two sets of parameters: (1) modeled object parameters (hydraulic groundwater gradient, hydraulic conductivity values of aquifer and streambed), and (2) modeling parameters (vertical discretization of aquifer, horizontal refinement of stream modeled with River [RIV] package). The influence of these various parameters on simulation results is investigated, and potential complications and errors are identified. Contaminant spread from stream to aquifer is not always reproduced by MT3D due to the RIV package's inability to simulate lateral exchange fluxes between stream and aquifer. This paper identifies the need for a MODFLOW streamflow package allowing lateral stream‐aquifer interactions and streamflow routine calculations. Such developments could be of particular interest for modeling contaminated flow‐through streams.     

The New Potential for Understanding Groundwater Contaminant Transport

The groundwater remediation field has been changing constantly since it first emerged in the 1970s. The remediation field has evolved from a dissolved‐phase centric conceptual model to a DNAPL‐dominated one, which is now being questioned due to a renewed appreciation of matrix diffusion effects on remediation. Detailed observations about contaminant transport have emerged from the remediation field, and challenge the validity of one of the mainstays of the groundwater solute transport modeling world: the concept of mechanical dispersion (Payne et al. 2008). We review and discuss how a new conceptual model of contaminant transport based on diffusion (the usurper) may topple the well‐established position of mechanical dispersion (the status quo) that is commonly used in almost every groundwater contaminant transport model, and evaluate the status of existing models and modeling studies that were conducted using advection‐dispersion models.     

Uncertainty assessment of ephemeral gully identification,characteristics and topographic threshold when using aerial photographs in agricultural settings

Manual digitizing on aerial photographs is still commonly used for characterizing gully erosion over large areas. Even when automated detection procedures are implemented, manual digitizing is frequently being resorted to in order to constitute reference datasets used for training and validation. In both cases, manual digitizing entails some subjective decisions on behalf of the operator, which introduces uncertainty into the resulting datasets. To assess the magnitude of this uncertainty, 11 experienced operators were asked to digitize and classify ephemeral gullies (EGs) on cropland following a standardized methodology. The resulting 11 datasets were compared in terms of number, type and location of EGs. Furthermore, for EGs located on a well‐defined runoff flow concentration axis, the slope versus contributing area topographic thresholds required for initiating gully channels were assessed using four thresholding methods, and compared across the 11 datasets. The operators identified 259 different EGs. However, the number (52–139) and sum total length (8.9–23.7 km) of EGs varied widely across operators. Only 34% of the EGs were digitized by more than half of the operators, and 7% were identified by all. Identification of EGs located on a well‐defined flow concentration axis proved least subjective. The longer the EG and the more fields the EG crossed, the larger the number of operators that were able to identify it. EGs were also most easily identified when located in sugar beet fields as compared to other crops. EG classification and topographic threshold lines were also found to be strongly operator‐dependent. Quantile regression appeared to be one of the most robust thresholding methods. Operator subjectivity when digitizing EGs on orthophotographs introduces uncertainty that should be taken into account in future remote sensing‐based studies of EG erosion whenever they rely, in part or in full, on manual photograph interpretation. Copyright © 2014 John Wiley & Sons, Ltd.     

Scoria cone formation through a violent Strombolian eruption: Irao Volcano, SW Japan

Scoria cones are common volcanic features and are thought to most commonly develop through the deposition of ballistics produced by gentle Strombolian eruptions and the outward sliding of talus. However, some historic scoria cones have been observed to form with phases of more energetic violent Strombolian eruptions (e.g., the 1943–1952 eruption of Parícutin, central Mexico; the 1975 eruption of Tolbachik, Kamchatka), maintaining volcanic plumes several kilometers in height, sometimes simultaneous with active effusive lava flows. Geologic evidence shows that violent Strombolian eruptions during cone formation may be more common than is generally perceived, and therefore it is important to obtain additional insights about such eruptions to better assess volcanic hazards. We studied Irao Volcano, the largest basaltic monogenetic volcano in the Abu Monogenetic Volcano Group, SW Japan. The geologic features of this volcano are consistent with a violent Strombolian eruption, including voluminous ash and fine lapilli beds (on order of 10^?1 km³ DRE) with simultaneous scoria cone formation and lava effusion from the base of the cone. The characteristics of the volcanic products suggest that the rate of magma ascent decreased gradually throughout the eruption and that less explosive Strombolian eruptions increased in frequency during the later stages of activity. During the eruption sequence, the chemical composition of the magma became more differentiated. A new K–Ar age determination for phlogopite crystallized within basalt dates the formation of Irao Volcano at 0.4?±?0.05 Ma.     

On‐Site Vapor‐Phase Analysis as a Novel Approach for Monitoring Groundwater Wells

The results of comprehensive field testing of on‐site vapor‐phase‐based groundwater monitoring methods are presented to demonstrate their utility as a robust and cost‐effective approach for rapidly obtaining volatile organic compounds (VOCs) concentration data from a monitoring well. These methods—which rely on sensitive, commercially available field equipment to analyze vapor in equilibrium with groundwater—proved easy to implement and can be tailored to site‐specific needs, including multilevel sampling. During field testing, low‐flow groundwater concentrations could be reasonably estimated using submerged passive vapor diffusion samplers or field equilibration of collected groundwater (R² = 0.85 to 0.96). These two methods are not as reliant on in‐well mixing to overcome vertical stratification within wells as simpler headspace methods. The importance of well and aquifer‐specific factors on concentration data (and therefore method selection) is highlighted, including the effect of changing in‐well patterns due to seasonal temperature gradients. Results indicated that vertical stratification was relatively limited within the set of wells included in these studies, resulting in similar performance for short depth‐discrete passive vapor diffusion samplers (constructed from 40‐mL vials) and longer samplers (2.5 to 5 feet in length) designed to cover a larger portion of the screened interval. A year‐long, multi‐event evaluation demonstrated that vapor‐phase‐based monitoring methods are no more variable than conventional groundwater monitoring methods, with both types subject to similar spatial and temporal variability that can be difficult to reduce. Vapor sampling methods represent a promising approach for estimation of groundwater concentrations by reducing the cost liabilities associated with monitoring while providing a more sustainable approach.     

Assessing wave climate trends in the Bay of Biscay through an intercomparison of wave hindcasts and reanalyses

The Bay of Biscay, located in the Northeast Atlantic Ocean, is exposed to energetic waves coming from the open ocean that have crucial effects on the coast. Knowledge of the wave climate and trends in this region are critical to better understand the last decade’s evolution of coastal hazards and morphology and to anticipate their potential future changes. This study aims to characterize the long-term trends of the present wave climate over the second half of the twentieth century in the Bay of Biscay through a robust and homogeneous intercomparison of five-wave datasets (Corrected ERA-40 (C-ERA-40), ECMWF Reanalysis Interim (ERA-Interim), Bay Of Biscay Wave Atlas (BOBWA-10kH), ANEMOC, and Bertin and Dodet 2010)). The comparison of the quality of the datasets against offshore and nearshore measurements reveals that at offshore locations, global reanalyses slightly underestimate wave heights, while regional hindcasts overestimate wave heights, especially for the highest quantiles. At coastal locations, BOBWA-10kH is the dataset that compares the best with observations. Concerning long time-scale features, the comparison highlights that the main significant trends are similarly present in the five datasets, especially during summer for which there is an increase of significant wave heights and mean wave periods (up to +15 cm and +0.6 s over the period 1970–2001) as well as a southerly shift of wave directions (around ?0.4° year^?1). Over the same period, an increase of high quantiles of wave heights during the autumn season (around 3 cm year^?1 for 90th quantile of significant wave heights (SWH₉₀)) is also apparent. During winter, significant trends are much lower than during summer and autumn despite a slight increase of wave heights and periods during 1958–2001. These trends can be related to modifications in the wave-type occurrence. Finally, the trends common to the five datasets are discussed by analyzing the similarities with centennial trends issued from longer time-scale studies and exploring the various factors that could explain them.     

A Snapshot of the Sun Near Solar Minimum: The Whole Heliosphere Interval

We present an overview of the data and models collected for the Whole Heliosphere Interval, an international campaign to study the three-dimensional solar?Cheliospheric?Cplanetary connected system near solar minimum. The data and models correspond to solar Carrington Rotation 2068 (20 March??C?16 April 2008) extending from below the solar photosphere, through interplanetary space, and down to Earth??s mesosphere. Nearly 200 people participated in aspects of WHI studies, analyzing and interpreting data from nearly 100 instruments and models in order to elucidate the physics of fundamental heliophysical processes. The solar and inner heliospheric data showed structure consistent with the declining phase of the solar cycle. A closely spaced cluster of low-latitude active regions was responsible for an increased level of magnetic activity, while a highly warped current sheet dominated heliospheric structure. The geospace data revealed an unusually high level of activity, driven primarily by the periodic impingement of high-speed streams. The WHI studies traced the solar activity and structure into the heliosphere and geospace, and provided new insight into the nature of the interconnected heliophysical system near solar minimum.     

Offshore sand sources for beach replenishment: Potential borrows on the continental shelf of the Eastern Gulf of Mexico

Erosion of sandy beaches is a worldwide problem that elicits innovative geoengineer‐ing techniques to reduce adverse impacts of shoreline retreat. Beach replenishment has emerged as the “soft”; shore‐stabilization technique of choice for mitigating beach erosion. This method of shore protection involves the addition of sand to the littoral sediment budget for sacrificial purposes. Because inland sand sources are often uneconomical or impractical to use, and known nearshore sources are limited, finding adequate quantities of suitable sand on the inner continental shelf is often vital to beach replenishment projects. The technical studies of survey and materials analysis that identify and delineate usable sand sources are sometimes almost as expensive as small‐project dredging, pumping, and placing the sand on the beach as fill. Inadequate quantity or substandard quality of shelf sand, as well as often‐prohibitive overhead expenses, thus compel shoreline managers to seek suitable sand sources offshore.

In the study area off the central‐west coast of Florida, offshore potential borrow areas (PBAs) were identified on the basis of studies conducted in reconnoitory and detailed phases. Sophisticated state‐of‐the‐art equipment used in this investigation provided more detailed subbottom mapping information than is normally obtained with conventional seismic equipment. An example of sand exploration studies was incorporated in a 215‐km² survey of offshore areas by conducting bathymetric surveys and subbottom seismic profiling, collecting jet probes, grab samples, and vibrocores, and analyzing sediment grading in subsamples from vibrocores. These combined analyses indicated that at least 8.8 ×10⁶ m³ of sand is available in potential borrow areas from 7.0 to 12 km offshore in water depths of 8.0 to 11.5 m. In the PBAs, mean grain size of sand falls into the range 0.13–0.53 mm, sorting averages 0.65–1.31ø, and the overall silt content varies from 3.9–8.5%. High silt contents (13–19%) mapped in some areas make these sedimentary deposits unsuitable as fill for artificial beach renourishment.