Discrete‐Storm Water‐Table Fluctuation Method to Estimate Episodic Recharge

We have developed a method to identify and quantify recharge episodes, along with their associated infiltration‐related inputs, by a consistent, systematic procedure. Our algorithm partitions a time series of water levels into discrete recharge episodes and intervals of no episodic recharge. It correlates each recharge episode with a specific interval of rainfall, so storm characteristics such as intensity and duration can be associated with the amount of recharge that results. To be useful in humid climates, the algorithm evaluates the separability of events, so that those whose recharge cannot be associated with a single storm can be appropriately lumped together. Elements of this method that are subject to subjectivity in the application of hydrologic judgment are values of lag time, fluctuation tolerance, and master recession parameters. Because these are determined once for a given site, they do not contribute subjective influences affecting episode‐to‐episode comparisons. By centralizing the elements requiring scientific judgment, our method facilitates such comparisons by keeping the most subjective elements openly apparent, making it easy to maintain consistency. If applied to a period of data long enough to include recharge episodes with broadly diverse characteristics, the method has value for predicting how climatic alterations in the distribution of storm intensities and seasonal duration may affect recharge.     

Geophysical and Hydrologic Studies of Lake Seepage Variability

Variations in lake seepage were studied along a 130 m shoreline of Mirror Lake NH. Seepage was downward from the lake to groundwater; rates measured from 28 seepage meters varied from 0 to ?282 cm/d. Causes of this variation were investigated using electrical resistivity surveys and lakebed sediment characterization. Two‐dimensional (2D) resistivity surveys showed a transition in lakebed sediments from outwash to till that correlated with high‐ and low‐seepage zones, respectively. However, the 2D survey was not able to predict smaller scale variations within these facies. In the outwash, fast seepage was associated with permeability variations in a thin (2 cm) layer of sediments at the top of the lakebed. In the till, where seepage was slower than that in the outwash, a three‐dimensional resistivity survey mapped a point of high seepage associated with heterogeneity (lower resistivity and likely higher permeability). Points of focused flow across the sediment–water interface are difficult to detect and can transmit a large percentage of total exchange. Using a series of electrical resistivity geophysical methods in combination with hydrologic data to locate heterogeneities that affect seepage rates can help guide seepage meter placement. Improving our understanding of the causes and types of heterogeneity in lake seepage will provide better data for lake budgets and prediction of mass transfer of solutes or contaminants between lakes and groundwater.     

Controls of Wellbore Flow Regimes on Pump Effluent Composition

Where well water and formation water are compositionally different or heterogeneous, pump effluent composition will vary due to partial mixing and transport induced by pumping. Investigating influences of purging and sampling methodology on composition variability requires quantification of wellbore flow regimes and mixing. As a basis for this quantification, analytical models simulating Poiseuille flow were developed to calculate flow paths and travel times. Finite element modeling was used to incorporate influences of mixing. Parabolic velocity distributions within the screened interval accelerate with cumulative inflow approaching the pump intake while an annulus of inflowing formation water contracts uniformly to displace an axial cylinder of pre‐pumping well water as pumping proceeds. Increased dispersive mixing forms a more diffuse formation water annulus and the contribution of formation water to pump effluent increases more rapidly. Models incorporating viscous flow and diffusion scale mixing show that initially pump effluent is predominantly pre‐pumping well water and compositions vary most rapidly. After two screen volumes of pumping, 94% of pump effluent is inflowing formation water. Where the composition of formation water and pre‐pumping well water are likely to be similar, pump effluent compositions will not vary significantly and may be collected during early purging or with passive sampling. However, where these compositions are expected to be considerably different or heterogeneous, compositions would be most variable during early pumping, that is, when samples are collected during low‐flow sampling. Purging of two screen volumes would be required to stabilize the content and collect a sample consisting of 94% formation water.     

Modeling the effects of climate change projections on streamflow in the Nooksack River basin,Northwest Washington

The Nooksack River has its headwaters in the North Cascade Mountains and drains an approximately 2000 km² watershed in northwestern Washington State. The timing and magnitude of streamflow in a snowpack‐dominated drainage basin such as the Nooksack River basin are strongly influenced by temperature and precipitation. Projections of future climate made by general circulation models (GCMs) indicate increases in temperature and variable changes in precipitation for the Nooksack River basin. Understanding the response of the river to climate change is crucial for regional water resources planning because municipalities, tribes, and industry depend on the river for water use and for fish habitat. We combine three different climate scenarios downscaled from GCMs and the Distributed‐Hydrology‐Soil‐Vegetation Model to simulate future changes to timing and magnitude of streamflow in the higher elevations of the Nooksack River. Simulations of future streamflow and snowpack in the basin project a range of magnitudes, which reflects the variable meteorological changes indicated by the three GCM scenarios and the local natural variability employed in the modeling. Simulation results project increased winter flows, decreased summer flows, decreased snowpack, and a shift in timing of the spring melt peak and maximum snow water equivalent. These results are consistent with previous regional studies, but the magnitude of increased winter flows and total annual runoff is higher. Increases in temperature dominate snowpack declines and changes to spring and summer streamflow, whereas a combination of increases in temperature and precipitation control increased winter streamflow. Copyright © 2013 John Wiley & Sons, Ltd.     

Human impact on Holocene sediment dynamics in the Eastern Mediterranean – the example of the Roman harbour of Ephesus

During the past millennia, many erosion and accumulation processes have been modified by anthropogenic impact. This holds especially true for the environs of ancient settlements and their harbours along the Mediterranean coasts. Our multi‐proxy investigations in the Roman harbour and the harbour canal of Ephesus (western Turkey) reveals that humans have significantly triggered soil erosion during the last three millennia. Since the eighth century bc , and especially since the Hellenistic period, a high sedimentation rate indicates fast alluviation and delta progradation of the Küçük Menderes. Deforestation, agriculture (especially ploughing) and grazing (especially goats) were the main reasons for erosion of the river catchment area. One consequence was significant siltation of the Hellenistic/Roman harbour basin. This sediment trap archives the human impact, which was strongly enhanced from Hellenistic/Roman to Byzantine times (second/first centuries bc to the sixth/seventh centuries ad ), evidenced by high sedimentation rates, raised values of heavy metal contaminations [lead (Pb), copper (Cu)], the occurrence of fruit tree pollen and of intestinal parasites. From the middle to the end of the first millennium ad , the influence of Ephesus declined, which resulted in a decrease of human impact. Studies of several ancient settlements around the Mediterranean Sea tell a comparable story. They also confirm that during their most flourishing periods the human impact totally overprinted the climatic one. To detect the latter, geo‐bio‐archives of relatively pristine areas have to be investigated in detail. Copyright © 2016 John Wiley & Sons, Ltd.     

Polar mesosphere summer echoes with ESRAD,Kiruna, Sweden: Variations and trends over 1997–2008

Measurements of Polar Mesosphere Summer Echoes (PMSE) with the MST radar ESRAD, in northern Sweden, from 1997 to 2008, were used to study diurnal, day-to-day and year-to-year variations of PMSE. The PMSE occurrence rate and volume reflectivity on a daily scale show predominantly semidiurnal variations with small interannual variability in the shape of the diurnal curves. Day-to-day and interannual variations of PMSE are found to correlate with geomagnetic activity while they do not correlate with mesopause temperature or solar activity. No statistically significant trends in PMSE occurrence rate and length of PMSE season were detected over the observation interval.     

A Simplified Kinetic Element Formulation for the Rotation of a Perturbed Mass-Asymmetric Rigid Body

Euler's equations, describing the rotation of an arbitrarily torqued mass asymmetric rigid body, are scaled using linear transformations that lead to a simplified set of first order ordinary differential equations without the explicit appearance of the principal moments of inertia. These scaled differential equations provide trivial access to an analytical solution and two constants of integration for the case of torque-free motion. Two additional representations for the third constant of integration are chosen to complete two new kinetic element sets that describe an osculating solution using the variation of parameters. The elements' physical representations are amplitudes and either angular displacement or initial time constant in the torque-free solution. These new kinetic elements lead to a considerably simplified variation of parameters solution to Euler's equations. The resulting variational equations are quite compact. To investigate error propagation behaviour of these new variational formulations in computer simulations, they are compared to the unmodified equations without kinematic coupling but under the influence of simulated gravity-gradient torques.     

Shock compression of deuterium to 900 Kbar

An experimental program to compress liquid deuterium with a two-stage gun is described. A xenon standard was used to check the adequacy of the EOS of the Dural container in the 7000 (1.7 Mbar) to 3000° K (0.0 Mbar) range. The agreement is reasonable. Further work on this standard is in progress to improve the accuracy. The target design parameters were checked with some sound speed measurements which confirm the liquid-like behavior of Dural at the temperatures, pressures, and strain rates of these experiments. The uncertainties of the results are analyzed and tabulated with the data. Comparison between the molecular and the metallic equations of state indicates a small volume change on metallization due to the high compressibility of the molecular phase vis-a-vis the metallic phase.     

Frequency-dependent Shear-wave Q Models for the Crust of China and Surrounding Regions

— Starting with fundamental-mode Rayleigh-wave attenuation coefficient values (_R) predicted by previously determined frequency-independent models of shear-wave Q (Q), we have obtained frequency-dependent Q models that explain measured values of _R as well as of Lg coda Q and its frequency dependence at 1 Hz (Q_o and , respectively) for China and some adjacent regions. The process combines trial-and-error selection of a model for the depth distribution of the frequency dependence parameter () for Q with a formal inversion for the depth distribution of Q at 1 Hz. Fifteen of the derived models have depth distributions of that are constant, or nearly constant, between the surface and a depth of 30 km. distributions that vary with depth in the upper 30 km are necessary to explain the remaining seven models. values for the depth-independent models vary between 0.4 and 0.7 everywhere except in the western portion of the Tibetan Plateau where they range between about 0.1 and 0.3 for three paths. These low values lie in a region where Q_Lg and crustal Q are very low and suggest that they should also be low for high-frequency propagation. The models in which varies with depth all show a decrease in that value ranging between 0.55 and 0.8 in the upper 15 km of the crust and (with two exceptions where =0.0) between 0.3 and 0.55 in the depth range 15–30 km. The distribution of values between 0.6 and 0.8 (the higher part of the range) in the upper crust indicates that high-frequency waves will propagate most efficiently, relative to low-frequency waves, in a band that includes, and strikes north-northeastward from the path between event 212/97 and KMI to the path between event 180/95 and station HIA in the north.Acknowledgement. We thank Lianli Cong for providing his code for plotting crustal Q models and Robert Herrmann for writing the mode summation code for computing Lg synthetics used in this study. Our work benefited from helpful discussions with Jack Xie at Lamont-Doherty Earth Observatory of Columbia University. This research was sponsored by the Defense Threat Reduction Agency Contract No. DTRA-01-00-C-0213.