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.     

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.     

Alternating field characteristics of pseudo-single-domain (2–14 μm) and multidomain magnetite

We report normalized AF demagnetization curves of anhysteretic remanences (ARM's) produced by 1-, 10- and 40-Oe steady fields and of saturation isothermal remanence (IRM_s) in a suite of dispersed, unannealed magnetite powders with median sizes of 2, 4, 6, 10 and 14 μm (pseudo-single-domain or PSD size range) and 100 μm (multidomain or MD size). Interpreted in the light of the domain structure test first proposed by Lowrie and Fuller [12], the relative stability trend of curves for the 2 μm sample is of single-domain (SD) type, the 1-Oe ARM being most resistant to demagnetization followed by the 10-Oe and 40-Oe ARM's and IRM_s. For the 100-μm sample, the trend is exactly reversed and is of MD-type. In the 4–14 μm samples, hitherto undescribed transitional trends between SD-type and MD-type occur. At 6 μm, 1-Oe, 10-Oe and 40-Oe ARM's preserve an SD-type trend but for all AF's > 75 Oe, IRM_s is more resistant than any of these remanences. At 10 μm, this trend is unmistakable, and only at 14 μm do the 1-Oe, 10-Oe and 40-Oe ARM curves merge. We conclude (1) that the Lowrie-Fuller test distinguishes between small MD grains enhanced by PSD remanence and large MD grains lacking PSD remanence, rather than between SD and MD structures per se, and (2) that in the PSD transition region from 6 to 14 μm in magnetite, IRM_s changes over to MD-type relative stability around 6 μm, whereas 10-Oe and 40-Oe ARM's achieve an MD-type trend around 14 μm, in accord with the predicted field dependence of the PSD threshold size.Our theoretical interpretation assumes that the intrinsic (internal field) coercive force spectra of weak-field and strong-field remanences are identical but that the observed (external field) spectrum is shifted to lower fields as a result of the internal demagnetizing field — NJ_r of the remanence J_r. The effect is slight for weak-field J_r's but substantial for IRM_s. Since all coercivities, high as well as low, are shifted, the result of the Lowrie-Fuller test is determined simply by the shape of the intrinsic coercivity spectrum or the corresponding AF demagnetization curve. Depending on the model of self-demagnetization used, either subexponential or sublinear AF decay curves of weak-field remanence will automatically lead to an MD-type trend, whereas by either model the decay curves that characterize SD and PSD remanences (decaying slowly initially and then more rapidly) will always produce and SD-type trend.     

Formation and pathways of North Atlantic Deep Water in a coupled ice–ocean model of the Arctic–North Atlantic Oceans

We investigate the formation process and pathways of deep water masses in a coupled ice–ocean model of the Arctic and North Atlantic Oceans. The intent is to determine the relative roles of these water masses from the different source regions (Arctic Ocean, Nordic Seas, and Subpolar Atlantic) in the meridional overturning circulation. The model exhibits significant decadal variability in the deep western boundary current and the overturning circulation. We use detailed diagnostics to understand the process of water mass formation in the model and the resulting effects on the North Atlantic overturning circulation. Particular emphasis is given to the multiple sources of North Atlantic Deep Water, the dominant deep water masses of the world ocean. The correct balance of Labrador Sea, Greenland Sea and Norwegian Sea sources is difficult to achieve in climate models, owing to small-scale sinking and convection processes. The global overturning circulation is described as a function of potential temperature and salinity, which more clearly signifies dynamical processes and clarifies resolution problems inherent to the high latitude oceans. We find that fluxes of deep water masses through various passages in the model are higher than observed estimates. Despite the excessive volume flux, the Nordic Seas overflow waters are diluted by strong mixing and enter the Labrador Sea at a lighter density. Through strong subpolar convection, these waters along with other North Atlantic water masses are converted into the densest waters [similar density to Antarctic Bottom Water (AABW)] in the North Atlantic. We describe the diminished role of salinity in the Labrador Sea, where a shortage of buoyant surface water (or excess of high salinity water) leads to overly strong convection. The result is that the Atlantic overturning circulation in the model is very sensitive to the surface heat flux in the Labrador Sea and hence is correlated with the North Atlantic Oscillation. As strong subpolar convection is found in other models, we discuss broader implications.     

Lightning optical pulse statistics from storm overflights during the Altus Cumulus Electrification Study

The Altus Cumulus Electrification Study (ACES) was conducted during the month of August 2002 in an area near Key West, Florida. One of the goals of this uninhabited aerial vehicle (UAV) study was to collect time resolved optical pulse data from thunderstorms. During the month long campaign, we acquired 5256 lightning generated optical pulses. Most of these observations were made while quite close to the top of the storms. We divided our data into two amplitude groups based on prior NASA U2 aircraft optical data and the pulse characteristics. The group of strong pulses with radiance greater than 2.1 mW m^{− 2} sr^{− 1} had mean and median 10–10% optical pulse widths of 770 and 740 μs, respectively, 50–50% pulse widths of 399 and 355 μs, respectively, and 10–90% risetimes of 292 and 260 μs. These values are very similar to the previous U2 based optical results The other group of pulses consisting of slightly more than a quarter of the total pulses observed had radiances less than the minimum values detected in the U2 study. These weak pulses were narrower than the strong pulses with 50–50% mean and median values of 199 and 160 μs, respectively. Only 12% of the flashes observed contained only weak pulses. The source of these weak pulses is unknown, but we suspect that some are artifacts of the close proximity of the aircraft to cloud top.