Assessing regional‐scale spatio‐temporal patterns of groundwater–surface water interactions using a coupled SWAT‐MODFLOW model

Interaction between groundwater and surface water in watersheds has significant impacts on water management and water rights, nutrient loading from aquifers to streams, and in‐stream flow requirements for aquatic species. Of particular importance are the spatial patterns of these interactions. This study explores the spatio‐temporal patterns of groundwater discharge to a river system in a semi‐arid region, with methods applied to the Sprague River Watershed (4100 km²) within the Upper Klamath Basin in Oregon, USA. Patterns of groundwater–surface water interaction are explored throughout the watershed during the 1970–2003 time period using a coupled SWAT‐MODFLOW model tested against streamflow, groundwater level and field‐estimated reach‐specific groundwater discharge rates. Daily time steps and coupling are used, with groundwater discharge rates calculated for each model computational point along the stream. Model results also are averaged by month and by year to determine seasonal and decadal trends in groundwater discharge rates. Results show high spatial variability in groundwater discharge, with several locations showing no groundwater/surface water interaction. Average annual groundwater discharge is 20.5 m³/s, with maximum and minimum rates occurring in September–October and March–April, respectively. Annual average rates increase by approximately 0.02 m³/s per year over the 34‐year period, negligible compared with the average annual rate, although 70% of the stream network experiences an increase in groundwater discharge rate between 1970 and 2003. Results can assist with water management, identifying potential locations of heavy nutrient mass loading from the aquifer to streams and ecological assessment and planning focused on locations of high groundwater discharge. Copyright © 2016 John Wiley & Sons, Ltd.     

Oxygen and Nutrient Exchanges at the Sediment-Water Interface: a Global Synthesis and Critique of Estuarine and Coastal Data

Estuarine and coastal marine sediment-water fluxes are considered to be important ecological features, but a global-scale assessment has yet to be developed. Goals of this work were to assemble a global-scale database of net sediment-water flux measurements, examine measurement techniques, characterize the geographic distribution and magnitude of sediment fluxes, explore the data for controls on sediment flux magnitude, and assess the importance of sediment fluxes in ecosystem-level metabolism and primary production. We examined 480 peer-reviewed sources and found sediment flux data for 167 estuarine and coastal systems. Most measurements were made in North America, Europe, and Australia. Fluxes varied widely among systems, some by several orders of magnitude. Inter-annual variability within sites was less than an order of magnitude but time series flux data to evaluate this were rare. However, limited time series data exhibited large and rapid responses to decreased external nutrient loading rates, climate change effects (possible temperature effects), and variability in trophic conditions. Comparative analyses indicated organic matter supply to sediments set the upper limits of flux magnitude, with other factors playing secondary roles. Two metrics were developed to assess ecosystem-level importance of sediment-water fluxes. Sediments represented 30% or more of depth-integrated rates of aerobic system respiration at depths of <10 m. An annual phytoplankton production data set was used to estimate N and P demand; sediments supplied an average of 15–32% of N and 17–100% of P demand and, in some cases, was as large or larger than external nutrient inputs. The percent of demand supplied by sediments was highest in temperate latitudes and lower in high and tropical latitudes.     

Dimethylsulfide production in Sargasso Sea eddies

Lagrangian time series of dimethylsulfide (DMS) concentrations from a cyclonic and an anticyclonic eddy in the Sargasso Sea were used in conjunction with measured DMS loss rates and a model of vertical mixing to estimate gross DMS production in the upper 60 m during summer 2004. Loss terms included biological consumption, photolysis, and ventilation to the atmosphere. The time- and depth (0–60 m)-averaged gross DMS production was estimated to be 0.73±0.09 nM d⁻¹ in the cyclonic eddy and 0.90±0.15 nM d⁻¹ in the anticyclonic eddy, with respective DMS replacement times of 5±1 and 6±1 d. The higher estimated rate of gross production and lower measured loss rate constants in the anticyclonic eddy were equally responsible for this eddy's 50% higher DMS inventory (0–60 m). When normalized to chlorophyll and total dimethylsulfoniopropionate (DMSP), estimated gross production in the anticyclonic eddy was about twice that in the cyclonic eddy, consistent with the greater fraction of phytoplankton that were DMSP producers in the anticyclonic eddy. Higher rates of gross production were estimated below the mixed layer, contributing to the subsurface DMS maximum found in both eddies. In both eddies, gas exchange, microbial consumption, and photolysis were roughly equal DMS loss terms in the surface mixed layer (0.2–0.4 nM d⁻¹). Vertical mixing was a substantial source of DMS to the surface mixed layer in both eddies (0.2–0.3 nM d⁻¹) owing to the relatively high DMS concentrations below the mixed layer. Estimated net biological DMS production rates (gross production minus microbial consumption) in the mixed layer were substantially lower (by almost a factor of 3) than those estimated in a previous study of the Sargasso Sea, which may explain the relatively low mixed-layer DMS concentrations found here during July 2004 (3 nM) compared to previous summers (4–6 nM).     

Chemical variations and genetic relationships between the Hess and Foy offset dikes at the Sudbury impact structure

Offset dikes are found concentrically around—and extending radially outward from—the Sudbury Igneous Complex (SIC), which represents an ~3 km thick differentiated impact melt sheet. The dikes are typically composed of an inclusion‐rich, so‐called quartz diorite (IQD) in the center of the dike, and an inclusion‐poor quartz diorite (QD) along the margins of the dike. New exposures of the intersection between the concentric Hess and radial Foy offset dikes provide an excellent opportunity to understand the relationship between the radial and concentric offset dikes and their internal phases. The goal was to constrain the timing of the dike emplacements relative to the impact and formation of the SIC. Results herein suggest that (1) the timing between the emplacement of the QD and IQD melts was geologically short, (2) the Hess and Foy dikes coexisted as melts at the same time and the intersection between them represents a mixture of the two, (3) the Foy dike has a slightly more evolved chemical composition than the Hess dike, and (4) the IQD melt from the Foy dike underwent some degree of chemical fractionation after its initial emplacement.     

Going Home: The Migration of Puerto Rican-Born Women from the United States to Puerto Rico*

This paper analyzes the migration of Puerto Rican-born women from the United States to Puerto Rico using longitudinal data. We hypothesize that sojourn length in the United States is a function of both structural (macro-level economic and cultural factors) and behavioral (micro-level life-cycle experiences and personal attributes) variables. We test these hypotheses by estimating a proportional hazards model. The parameter estimates of this model indicate that sojourn length in the United States, and thus the decision to return to Puerto Rico, is a function of wage trends and community characteristics on the mainland plus a number of individual attributes that include education, marriage, and childbirth.     

Volcano hazards program in the United States

Volcano monitoring and volcanic-hazards studies have received greatly increased attention in the United States in the past few years. Before 1980, the Volcanic Hazards Program was primarily focused on the active volcanoes of Kilauea and Mauna Loa, Hawaii, which have been monitored continuously since 1912 by the Hawaiian Volcano Observatory. After the reawakening and catastrophic eruption of Mount St. Helens in 1980, the program was substantially expanded as the government and general public became aware of the potential for eruptions and associated hazards within the conterminous United States. Integrated components of the expanded program include: volcanic-hazards assessment; volcano monitoring; fundamental research; and, in concert with federal, state, and local authorities, emergency-response planning.In 1980 the David A. Johnston Cascades Volcano Observatory was established in Vancouver, Washington, to systematically monitor the continuing activity of Mount St. Helens, and to acquire baseline data for monitoring the other, presently quiescent, but potentially dangerous Cascade volcanoes in the Pacific Northwest. Since June 1980, all of the eruptions of Mount St. Helens have been predicted successfully on the basis of seismic and geodetic monitoring.The largest volcanic eruptions, but the least probable statistically, that pose a threat to western conterminous United States are those from the large Pleistocene-Holocene volcanic systems, such as Long Valley caldera (California) and Yellowstone caldera (Wyoming), which are underlain by large magma chambers still potentially capable of producing catastrophic caldera-forming eruptions. In order to become better prepared for possible future hazards associated with such historically unpecedented events, detailed studies of these, and similar, large volcanic systems should be intensified to gain better insight into caldera-forming processes and to recognize, if possible, the precursors of caldera-forming eruptions.     

A study of objective methods for water mass analysis,applied to the Iceland Basin

Ocean Dynamics - Cluster analysis and principal component analysis have been applied to multi-parameter hydrographic data from the Iceland Basin. The results indicate that both objective methods...     

Vibrational nitrogen concentration in the ionosphere and its dependence on season and solar cycle

A fully time-dependent mathematical model, SUPIM, of the Earth’s plasmasphere is used in this investigation. The model solves coupled time-dependent equations of continuity, momentum and energy balance for the O⁺, H⁺, He⁺, N⁺₂, O⁺₂, NO⁺ ions and electrons; in the present study, the geomagnetic field is represented by an axial-centred dipole. Calculation of vibrationally excited nitrogen molecules, which has been incorporated into the model, is presented here. The enhanced model is then used to investigate the behaviour of vibrationally excited nitrogen molecules with F10.7 and solar EUV flux, during summer, winter and equinox conditions. The presence of vibrational nitrogen causes a reduction in the electron content. The diurnal peak in electron content increases linearly up to a certain value of F10.7, and above this value increases at a lesser rate, in agreement with previous observations and modelling work. The value of F10.7 at which this change in gradient occurs is reduced by the presence of vibrational nitrogen. Vibrational nitrogen is most effective at F-region altitudes during summer daytime conditions when a reduction in the electron density is seen. A lesser effect is seen at equinox, and in winter the effect is negligible. The summer reduction in electron density due to vibrational nitrogen therefore reinforces the seasonal anomaly.