Headwater stream length dynamics across four physiographic provinces of the Appalachian Highlands

Understanding patterns of expansion, contraction, and disconnection of headwater stream length in diverse settings is invaluable for the effective management of water resources as well as for informing research in the hydrology, ecology, and biogeochemistry of temporary streams. More accurate mapping of the stream network and quantitative measures of flow duration in the vast headwater regions facilitate implementation of water quality regulation and other policies to protect waterways. We determined the length and connectivity of the wet stream and geomorphic channel network in 3 forested catchments (<75 ha) in each of 4 physiographic provinces of the Appalachian Highlands: the New England, Appalachian Plateau, Valley and Ridge, and Blue Ridge. We mapped wet stream length 7 times at each catchment to characterize flow conditions between exceedance probabilities of <5% and >90% of the mean daily discharge. Stream network dynamics reflected geologic controls at both regional and local scales. Wet stream length was most variable at two Valley and Ridge catchments on a shale scarp slope and changed the least in the Blue Ridge. The density and source area of flow origins differed between the crystalline and sedimentary physiographic provinces, as the Appalachian Plateau and Valley and Ridge had fewer origins with much larger contributing areas than New England and the Blue Ridge. However, the length and surface connectivity of the wet stream depended on local lithology, geologic structure, and the distribution of surficial deposits such as boulders, glacially derived material, and colluival debris or sediment valley fills. Several proxies indicate the magnitude of stream length dynamics, including bankfull channel width, network connectivity, the base flow index, and the ratio of geomorphic channel to wet stream length. Consideration of geologic characteristics at multiple spatial scales is imperative for future investigations of flow intermittency in headwaters.     

Imaging Faults and Shear Zones Using Receiver Functions

The geometry of faults at seismogenic depths and their continuation into the ductile zone is of interest for a number of applications ranging from earthquake hazard to modes of lithospheric deformation. Teleseismic passive source imaging of faults and shear zones can be useful particularly where faults are not outlined by local seismicity. Passive seismic signatures of faults may arise from abrupt changes in lithology or foliation orientation in the upper crust, and from mylonitic shear zones at greater depths. Faults and shear zones with less than near-vertical dip lend themselves to detection with teleseismic mode-converted waves (receiver functions) provided that they have either a contrast in isotropic shear velocity (V _s), or a contrast in orientation or strength of anisotropic compressional velocity (V _p). We introduce a detection method for faults and shear zones based on receiver functions. We use synthetic seismograms to demonstrate common features of dipping isotropic interfaces and contrasts in dipping foliation that allows determination of their strike and depth without making further assumptions about the model. We proceed with two applications. We first image a Laramide thrust fault in the western U.S. (the Wind River thrust fault) as a steeply dipping isotropic velocity contrast in the middle crust near the surface trace of the fault; further downdip and across the range, where basin geometry suggests the fault may sole into a subhorizontal shear zone, we identify a candidate shear zone signal from midcrustal depths. The second application is the use of microstructural data from exhumed ductile shear zones in Scotland and in the western Canadian Shield to predict the character of seismic signatures of present-day deep crustal shear zones. Realistic anisotropy in observed shear fabrics generates a signal in receiver functions that is comparable in amplitude to first-order features like the Moho. Observables that can be robustly constrained without significant tradeoffs are foliation strike and the depth of the foliation contrast. We find that an anisotropy of only a few percent in the shear zone is sufficient to generate a strong signal, but that the shear zone width is required to be >2 km for typical frequencies used in receiver function analysis to avoid destructive interference due to the signals from the boundaries of the shear zone.     

Thermal precursors in satellite images of the 1999 eruption of Shishaldin Volcano

Shishaldin Volcano, Unimak Island Alaska, began showing signs of thermal unrest in satellite images on 9 February 1999. A thermal anomaly and small steam plume were detected at the summit of the volcano in short-wave thermal infrared AVHRR (advanced very high resolution radiometer) satellite data. This was followed by over 2 months of changes in the observed thermal character of the volcano. Initially, the thermal anomaly was only visible when the satellite passed nearly directly over the volcano, suggesting a hot source deep in the central crater obscured from more oblique satellite passes. The "zenith angle" needed to see the anomaly increased with time, presumably as the thermal source rose within the conduit. Based on this change, an ascent rate of ca. 14 m per day for the thermal source was estimated, until it reached the summit on around 21 March. It is thought that Strombolian activity began around this time. The precursory activity culminated in a sub-Plinian eruption on 19 April, ejecting ash to over 45,000 ft. (13,700 m). The thermal energy output through the precursory period was calculated based on geometric constraints unique to Shishaldin. These calculations show fluctuations that can be tied to changes in the eruptive character inferred from seismic records and later geologic studies. The remote location of this volcano made satellite images a necessary observation tool for this eruption. To date, this is the longest thermal precursory activity preceding a sub-Plinian eruption recorded by satellite images in the region. This type of thermal monitoring of remote volcanoes is central in the efforts of the Alaska Volcano Observatory to provide timely warnings of volcanic eruption, and mitigate their associated hazards to air-traffic and local residents.     

Salt ramps: Wind-induced depositional features on Tunisian playas

Semicircular and crescent-shaped accumulations of salt crystals developed on salt crusts were measured on three Tunisian playas in September 1990. These features have been termed salt ramps. Their morphology and chemistry suggest that they are formed in the late stages of shallow ephemeral lake desiccation in playa basins. They form by salt precipitation from shallow brine lakes that are blown across salt-encrusted playa surfaces by the wind. Moreover, they appear to be short-lived features and their degradation is related to the flooding of playas with less saline water, and possibly rainfall and deflation.     

Forecasting Natural Aquifer Discharge Using a Numerical Model and Convolution

If the nature of groundwater sources and sinks can be determined or predicted, the data can be used to forecast natural aquifer discharge. We present a procedure to forecast the relative contribution of individual aquifer sources and sinks to natural aquifer discharge. Using these individual aquifer recharge components, along with observed aquifer heads for each January, we generate a 1‐year, monthly spring discharge forecast for the upcoming year with an existing numerical model and convolution. The results indicate that a forecast of natural aquifer discharge can be developed using only the dominant aquifer recharge sources combined with the effects of aquifer heads (initial conditions) at the time the forecast is generated. We also estimate how our forecast will perform in the future using a jackknife procedure, which indicates that the future performance of the forecast is good (Nash‐Sutcliffe efficiency of 0.81). We develop a forecast and demonstrate important features of the procedure by presenting an application to the Eastern Snake Plain Aquifer in southern Idaho.     

Dynamic response analysis of solitary flexible rocking bodies: modeling and behavior under pulse‐like ground excitation

Results obtained for rigid structures suggest that rocking can be used as seismic response modification strategy. However, actual structures are not rigid: structural elements where rocking is expected to occur are often slender and flexible. Modeling of the rocking motion and impact of flexible bodies is a challenging task. A non‐linear elastic viscously damped zero‐length spring rocking model, directly usable in conventional finite element software, is presented in this paper. The flexible rocking body is modeled using a conventional beam‐column element with distributed masses. This model is verified by comparing its pulse excitation response to the corresponding analytical solution and validated by overturning analysis of rocking blocks subjected to a recorded ground motion excitation. The rigid rocking block model provides a good approximation of the seismic response of solitary flexible columns designed to uplift when excited by pulse‐like ground motions. Guidance for development of rocking column models in ordinary finite element software is provided. Copyright © 2014 John Wiley & Sons, Ltd.     

Assessing characteristics of Mediterranean explosive cyclones for different data resolution

A comparison of two objective climatologies of explosive cyclones in the Mediterranean region is performed. The results are derived from two different mean sea-level pressure reanalysis data resolutions, but from the same assimilation model, in order to quantify the pure impact of higher resolution on the identification and characteristics of explosive cyclones, when the assimilation model is the same. The explosive cyclones were identified with the aid of the Melbourne University automatic cyclone finding and tracking scheme over a 40-year period, using the 6-hourly analyses of ERA-40 mean sea-level pressure (MSLP) on: (a) 2.5?×?2.5 and (b) 1?×?1 latitude–longitude grid. The comparison of the two datasets revealed the significant role of the increase in spatial resolution of MSLP data on the identification and tracking process, and the number of the explosive cyclones in the high-resolution dataset is almost four times greater than the respective one in the lower resolution dataset. However, the comparison of explosive cyclone characteristics, including spatial and temporal variations of explosive deepening, revealed differences in the geographical distribution of the location of the maximum explosive deepening and average explosive cyclone Laplacian of the central pressure. These differences are due to the identification in the higher resolution set of smaller scale and secondary explosives along the strongly baroclinic northern Mediterranean boundaries, south of the Alps and the Pyrenees. Explosive deepening appears a bias to the daytime period from 12 to 18 Coordinated Universal Time (UTC) for both datasets, which is more prominent in the LR dataset. Statistically significant difference of pressure tendency between the two datasets appear for the daytime period from 06 to 12 UTC, accounting for better representation of orographic forcing in the HR dataset.     

Geostatistical Borehole Image‐Based Mapping of Karst‐Carbonate Aquifer Pores

Quantification of the character and spatial distribution of porosity in carbonate aquifers is important as input into computer models used in the calculation of intrinsic permeability and for next‐generation, high‐resolution groundwater flow simulations. Digital, optical, borehole‐wall image data from three closely spaced boreholes in the karst‐carbonate Biscayne aquifer in southeastern Florida are used in geostatistical experiments to assess the capabilities of various methods to create realistic two‐dimensional models of vuggy megaporosity and matrix‐porosity distribution in the limestone that composes the aquifer. When the borehole image data alone were used as the model training image, multiple‐point geostatistics failed to detect the known spatial autocorrelation of vuggy megaporosity and matrix porosity among the three boreholes, which were only 10 m apart. Variogram analysis and subsequent Gaussian simulation produced results that showed a realistic conceptualization of horizontal continuity of strata dominated by vuggy megaporosity and matrix porosity among the three boreholes.     

Climate-related variations in crustal trace elements in Dome C (East Antarctica) ice during the past 672 kyr

Cr, Fe, Rb, Ba and U were determined by inductively coupled plasma sector field mass spectrometry (ICP-SFMS) in various sections of the 3,270 m deep ice core recently drilled at Dome C on the high East Antarctic plateau as part of the EPICA program. The sections were dated from 263 kyr bp (depth of 2,368 m) to 672 kyr bp (depth of 3,062 m). When combined with the data previously obtained by Gabrielli and co-workers for the upper 2,193 m of the core, it gives a detailed record for these elements during a 672-kyr period from the Holocene back to Marine Isotopic Stage (MIS) 16.2. Concentrations and fallout fluxes of all elements are found to be highly variable with low values during the successive interglacial periods and much higher values during the coldest periods of the last eight climatic cycles. Crustal enrichment factors indicates that rock and soil dust is the dominant source for Fe, Rb, Ba and U whatever the period and for Cr during the glacial maxima. The relationship between Cr, Fe, Rb, Ba and U concentrations and the deuterium content of the ice appears to be similar before and after the Mid-Brunhes Event (MBE, around 430 kyr bp). Mean concentration values observed during the successive interglacials from the Holocene to MIS 15.5 appear to vary from one interglacial to another at least for part of the elements. Concentrations observed during the successive glacial maxima suggest a decreasing trend from the most recent glacial maxima (MIS 2.2 and 4.2) to the oldest glacial maxima such as MIS 14.2, 14.4 and 16.2, which could be linked with changes in the size distribution of dust particles transported from mid-latitude areas to the East Antarctic ice cap.