Variable near-surface deformation along the Commerce segment of the Commerce geophysical lineament, southeast Missouri to southern Illinois, USA

Recent studies have demonstrated a plausible link between surface and near-surface tectonic features and the vertical projection of the Commerce geophysical lineament (CGL). The CGL is a 5- to 10-km-wide zone of basement magnetic and gravity anomalies traceable for more than 600 km, extending from Arkansas through southeast Missouri and southern Illinois and into Indiana. Twelve kilometers of high-resolution seismic reflection data, collected at four sites along a 175-km segment of the CGL projection, are interpreted to show varying amounts of deformation involving Tertiary and some Quaternary sediments. Some of the locally anomalous geomorphic features in the northern Mississippi embayment region (i.e., paleoliquefaction features, anomalous directional changes in stream channels, and areas of linear bluff escarpments) overlying the CGL can be correlated with specific faults and/or narrow zones of deformed (faulted and folded) strata that are imaged on high-resolution seismic reflection data. There is an observable change in near-surface deformation style and complexity progressing from the southwest to the northeast along the trace of the CGL. The seismic reflection data collaborate mapping evidence which suggests that this region has undergone a complex history of deformation, some of which is documented to be as young as Quaternary, during multiple episodes of reactivation under varying stress fields. This work, along with that of other studies presented in this volume, points to the existence of at least one major crustal feature outside the currently defined zone of seismic activity (New Madrid Seismic Zone) that should be considered as a significant potential source zone for seismogenic activity within the midcontinent region of the United States.     

Vulnerability of headwater catchment resources to incidences of <Superscript>210</Superscript>Pb excess and <Superscript>137</Superscript>Cs radionuclide fallout

Recent identification of elevated excess ²¹⁰Pb (≤302.6 mBq L⁻¹) and ¹³⁷Cs (≤111.3 mBq L⁻¹) activity in drinking water wells up to 20 m depth indicates some transport of airborne radionuclide fallout beyond soils in the Shaker Village catchment, Maine. Estimated airborne mass loading ²¹⁰Pb_ex fluxes of about 0.9 mBq m⁻³, canvass this headwater catchment and may be sufficient to pose risks to unprotected shallow wells. Inventories of ²¹⁰Pb_ex and ¹³⁷Cs in pond sediments indicate maximum median activities of 943 mBq g⁻¹ and 40.0 mBq g⁻¹, respectively. Calculated ²¹⁰Pb_ex fluxes in the catchment soils range from 0.62–0.78 Bq cm⁻² year⁻¹ and yield a mean residence time of near 140 years. Measured ¹³⁷Cs activity up to 51.1 mBq g⁻¹ occurs in sediments at least to 5 m depth. Assumed particle transport in groundwater with apparent ⁸⁵Kr ages less than 5 years BP (2005) may explain the correlation between these particle-reactive radionuclides and elevated activity in some drinking water wells.     

Geotechnical properties of lower trench inner-slope sediments

Geotechnical property analyses on sediments recovered by DSDP drilling within select convergent margins were integrated to determine the variability in these properties along the trench inner-slope and to delineate the effects of convergence on these properties. Consolidation states range from very underconsolidated to highly overconsolidated, with preconsolidation pressures exceeding 46,000 kPa in Quaternary sediments. Underconsolidated sediments are attributed to:

1. (1) high sedimentation rates.

2. (2) low sediment permeability relative to the length of the drainage path; (3) laterally applied stresses; and (4) induced pore water pressures resulting from the subduction of pelagic sediments with high water content. Factors contributing to the state of overconsolidation include: (1) tectonically induced overpressures; (2) removal of overburden by mass movement processes.

3. (3) low sediment accumulation rates.

4. (4) age.

Vertical gradients of index properties also vary greatly with the maximum gradient associated with overconsolidated sediments. Values of porosity generally exceed 30% at all margins studied, suggesting this porosity represents the minimum attainable solely by the effects of convergence.Geotechnical property results and site-specific parameters suggest two end members should be considered in the general geotechnical property model for convergent margin sediments. Clastic dominated margins form the basis for the initial geotechnical property model. These margins are characterized by thick trench sediment sequences which are folded and faulted and become progressively more deformed upslope. Convergence rates are low, and no well-developed faults or hrst and graben structures are evident on the downbending oceanic plate. Highly overconsolidated sediments and index property values which change rapidly with depth result from the slow progressive deformation of trench and lower slope sediments. Examples of clastic-dominated margins include the Nankai Trough, Aleutian Trench, and Washington continental margin.Pelagic-dominated margins have significantly different geotechnical properties and site-specific parameters. These margins are characterized by thin clastic trench sediment sequences overlying pelagic sediments on a rapidly converging oceanic plate. Well-developed faults and/or hrst and graben structures are evident on the downbending oceanic plate. Sediment deformation along the trench and lower slope appears limited. These sediments are underconsolidated and index property values change gradually but often irregularly with depth. Japan Trench and Middle America Trench-Guatemala margins are considered pelagic-dominated. Middle America Trench-Mexico appears to represent an intermediate case having characteristics of both margin types.     

A landscape analysis of the candelaria watershed in mexico: Insights into paleoclimates affecting upland horticulture in the southern yucatan peninsula semi-karst

Biocultural systems adapt to global climate change through its regional manifestations. Subsistence customs are the cultural interface between regional climate and culture at large. Swidden horticulturalists in Campeche, Mexico report that dry April followed by early onset of the wet season enhances the productivity of upland tropical gardens, or milpas. To relate regional seasonality of moisture to global climate, growing season discharge for Mexico's Candelaria River from 1958 to 1990 was analyzed relative to global average temperature. Analysis of covariance revealed a statistically significant relationship (p < 0.001). Further analysis showed that hot global climate eliminates the dry season, which lowers milpa productivity by preventing burning of the slash. Cold global climate delays the wet season and planting, also at cost of productivity. Intermediate global temperature fosters optimal wet—dry season combinations. Productivity of milpas is therefore directly related to global climate through the intervening mechanism of seasonality of moisture. A regression model reflecting these findings is used to retrodict paleohydrology for the last 3000 years. The pattern of ascendancy and decline of ancient southern Maya lowland urban centers is reviewed in the perspective of changing hydrological conditions. The model indicates that fluorescence occurred with optimal balance of wet and dry season duration and catastrophes unfolded during extended wet or dry periods. We suggest that the southern Maya lowlands have had a precipitous record of urban development and collapse in part because of complex interactions of global climate and upland horticulture. The most productive conditions for milpa issue from an inherently unstable overlay of global climate on a relatively narrow band of partially developed karst (semikarst) geological formations.     

Multiple P–T–d–t paths reveal the evolution of the final Nuna assembly in northeast Australia

The final assembly of the Mesoproterozoic supercontinent Nuna was marked by the collision of Laurentia and Australia at 1.60 Ga, which is recorded in the Georgetown Inlier of NE Australia. Here, we decipher the metamorphic evolution of this final Nuna collisional event using petrostructural analysis, major and trace element compositions of key minerals, thermodynamic modelling, and multi-method geochronology. The Georgetown Inlier is characterised by deformed and metamorphosed 1.70–1.62 Ga sedimentary and mafic rocks, which were intruded by c. 1.56 Ga old S-type granites. Garnet Lu–Hf and monazite U–Pb isotopic analyses distinguish two major metamorphic events (M1 at c. 1.60 Ga and M2 at c. 1.55 Ga), which allows at least two composite fabrics to be identified at the regional scale—c. 1.60 Ga S1 (consisting in fabrics S1a and S1b) and c. 1.55 Ga S2 (including fabrics S2a and S2b). Also, three tectono-metamorphic domains are distinguished: (a) the western domain, with S1 defined by low-P (LP) greenschist facies assemblages; (b) the central domain, where S1 fabric is preserved as medium-P (MP) amphibolite facies relicts, and locally as inclusion trails in garnet wrapped by the regionally dominant low-P amphibolite facies S2 fabric; and (c) the eastern domain dominated by upper amphibolite to granulite facies S2 foliation. In the central domain, 1.60 Ga MP–medium-T (MT) metamorphism (M1) developed within the staurolite–garnet stability field, with conditions ranging from 530–550°C at 6–7 kbar (garnet cores) to 620–650°C at 8–9 kbar (garnet rims), and it is associated with S1 fabric. The onset of 1.55 Ga LP–high-T (HT) metamorphism (M2) is marked by replacement of staurolite by andalusite (M2a/D2a), which was subsequently pseudomorphed by sillimanite (M2b/D2b) where granite and migmatite are abundant. P–T conditions ranged from 600 to 680°C and 4–6 kbar for the M2b sillimanite stage. 1.60 Ga garnet relicts within the S2 foliation highlight the progressive obliteration of the S1 fabric by regional S2 in the central zone during peak M2 metamorphism. In the eastern migmatitic complex, partial melting of paragneiss and amphibolite occurred syn- to post- S2, at 730–770°C and 6–8 kbar, and at 750–790°C and 6 kbar, respectively. The pressure–temperature–deformation–time paths reconstructed for the Georgetown Inlier suggest a c. 1.60 Ga M1/D1 event recorded under greenschist facies conditions in the western domain and under medium-P and medium-T conditions in the central domain. This event was followed by the regional 1.56–1.54 Ga low-P and high-T phase (M2/D2), extensively recorded in the central and eastern domains. Decompression between these two metamorphic events is ascribed to an episode of exhumation. The two-stage evolution supports the previous hypothesis that the Georgetown Inlier preserves continental collisional and subsequent thermal perturbation associated with granite emplacement.     

Solute and geothermal flux monitoring using electrical conductivity in the Madison,Firehole, and Gibbon Rivers,Yellowstone National Park

The thermal output from the Yellowstone magma chamber can be estimated from the Cl flux in the major rivers in Yellowstone National Park; and by utilizing continuous discharge and electrical conductivity measurements the Cl flux can be calculated. The relationship between electrical conductivity and concentrations of Cl and other geothermal solutes (Na, SO₄, F, HCO₃, SiO₂, K, Li, B, and As) was quantified at monitoring sites along the Madison, Gibbon, and Firehole Rivers, which receive discharge from some of the largest and most active geothermal areas in Yellowstone. Except for some trace elements, most solutes behave conservatively and the ratios between geothermal solute concentrations are constant in the Madison, Gibbon, and Firehole Rivers. Hence, dissolved concentrations of Cl, Na, SO₄, F, HCO₃, SiO₂, K, Li, Ca, B and As correlate well with conductivity (R² > 0.9 for most solutes) and most exhibit linear trends. The 2011 flux for Cl, SO₄, F and HCO₃ determined using automated conductivity sensors and discharge data from nearby USGS gaging stations is in good agreement with those of previous years (1983–1994 and 1997–2008) at each of the monitoring sites. Continuous conductivity monitoring provides a cost- and labor-effective alternative to existing protocols whereby flux is estimated through manual collection of numerous water samples and subsequent chemical analysis. Electrical conductivity data also yield insights into a variety of topics of research interest at Yellowstone and elsewhere: (1) Geyser eruptions are easily identified and the solute flux quantified with conductivity data. (2) Short-term heavy rain events can produce conductivity anomalies due to dissolution of efflorescent salts that are temporarily trapped in and around geyser basins during low-flow periods. During a major rain event in October 2010, 180,000 kg of additional solute was measured in the Madison River. (3) The output of thermal water from the Gibbon River appears to have increased by about 0.2%/a in recent years, while the output of thermal water for the Firehole River shows a decrease of about 10% from 1983 to 2011. Confirmation of these trends will require continuing Cl flux monitoring over the coming decades.     

A combinatorial optimization scheme for parameter structure identification in ground water modeling

This research develops a methodology for parameter structure identification in ground water modeling. For a given set of observations, parameter structure identification seeks to identify the parameter dimension, its corresponding parameter pattern and values. Voronoi tessellation is used to parameterize the unknown distributed parameter into a number of zones. Accordingly, the parameter structure identification problem is equivalent to finding the number and locations as well as the values of the basis points associated with the Voronoi tessellation. A genetic algorithm (GA) is allied with a grid search method and a quasi-Newton algorithm to solve the inverse problem. GA is first used to search for the near-optimal parameter pattern and values. Next, a grid search method and a quasi-Newton algorithm iteratively improve the GA's estimates. Sensitivities of state variables to parameters are calculated by the sensitivity-equation method. MODFLOW and MT3DMS are employed to solve the coupled flow and transport model as well as the derived sensitivity equations. The optimal parameter dimension is determined using criteria based on parameter uncertainty and parameter structure discrimination. Numerical experiments are conducted to demonstrate the proposed methodology, in which the true transmissivity field is characterized by either a continuous distribution or a distribution that can be characterized by zones. We conclude that the optimized transmissivity zones capture the trend and distribution of the true transmissivity field.     

Helium isotope geochemistry of some volcanic rocks from Saint Helena

³He/⁴He ratios have been measured for olivine and clinopyroxene phenocrysts in 7–15 m.y. old basaltic lavas from the island of St. Helena. Magmatic helium was effectively resolved from post-eruptive radiogenic helium by employing various extraction techniques, includingin vacuo crushing, and stepwise heating or fusion of the powders following crushing. The inherited³He/⁴He ratio at St. Helena is 4.3–5.9 R_A. Helium isotope disequilibrium is present within the phenocrysts, with lower³He/⁴He upon heating and fusion of the powders following crushing, due to radiogenic ingrowth or to -particle implantation from the surrounding(U + Th)-rich lavas.

A single crushing analysis for clinopyroxene in a basalt from Tubuaii gave³He/⁴He= 7.1 R_A.³He/⁴He ratios at St. Helena and Tubuaii (HIMU hotspots characterized by radiogenic Pb isotope signatures) are similar to³He/⁴He ratios previously measured at Tristan da Cunha and Gough Island (EM hotspots characterized by low²⁰⁶Pb/²⁰⁴Pb). Overall, the HeSrPb isotope systematics at these islands are consistent with a mantle origin as contiguous, heterogeneous materials, such as recycled crust and/or lithosphere.³He/⁴He ratios at HIMU hotspots are similar to mantle xenoliths which display nearly the entire range of Pb isotope compositions found at ocean islands, and are only slightly less than values found in mid-ocean ridge basalts (7–9 R_A). This suggests that the recycled materials were injected into the mantle within the last 10⁹ yrs.     

Impact of historic land‐use change on sediment delivery to a Chesapeake Bay subestuarine delta

Historic land use in the Chesapeake Bay drainage basin induced large fluxes of fluvial sediment to subestuarine tributaries. Stratigraphic and palaeoecologic analyses of deltaic deposits may be used to infer changes on the landscape, but are not sufficient to quantify past sediment supply. When viewed as an inverse boundary‐value problem, reconstruction of the sediment supply function may be achieved by combining deltaic sedimentation chronologies with an equation governing delta progradation. We propose that the diffusion equation is appropriate for simulating delta progradation and obtaining the sediment supply function provided a suitable diffusion constant (D) can be determined. Three new methods for estimating D are presented for the case of estuarine deltas. When the inverse boundary‐value technique was applied to Otter Point Creek, a tidal freshwater delta at the head of Bush River in upper Chesapeake Bay, D values ranged from 3763 to 6199 m² a^?1. Delta growth simulations showed a 1740–1760 initial pulse, a 1760–1780 erosive/redistributive interval, a 1780–1920 growth period, and a 1920‐present erosive/redistributive era. Coupling of simulated delta elevations with an empirical plant habitat predictive equation allowed for comparison of predicted versus actual relative habitat areas. Also, the model yielded reconstructed watershed erosion rates and stream suspended sediment concentrations that could be useful for development of water quality regulations. Copyright © 2001 John Wiley & Sons, Ltd.