USGS Science⌕ Search

SEARCH · USGS Science

Results for “Mechanics of Materials”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 523 records · Page 29Linked to original sources

Thermokarst and thaw-related landscape dynamics -- an annotated bibliography with an emphasis on potential effects on habitat and wildlife

Permafrost has warmed throughout much of the Northern Hemisphere since the 1980s, with colder permafrost sites warming more rapidly (Romanovsky and others, 2010; Smith and others, 2010). Warming of the near-surface permafrost may lead to widespread terrain instability in ice-rich permafrost in the Arctic and the Subarctic, and may result in thermokarst development and other thaw-related landscape features (Jorgenson and others, 2006; Gooseff and others, 2009). Thermokarst and other thaw-related landscape features result from varying modes and scales of permafrost thaw, subsidence, and removal of material. An increase in active-layer depth, water accumulation on the soil surface, permafrost degradation and associated retreat of the permafrost table, and changes to lake shores and coastal bluffs act and interact to create thermokarst and other thaw-related landscape features (Shur and Osterkamp, 2007). There is increasing interest in the spatial and temporal dynamics of thermokarst and other thaw-related features from diverse disciplines including landscape ecology, hydrology, engineering, and biogeochemistry. Therefore, there is a need to synthesize and disseminate knowledge on the current state of near-surface permafrost terrain. The term "thermokarst" originated in the Russian literature, and its scientific use has varied substantially over time (Shur and Osterkamp, 2007). The modern definition of thermokarst refers to the process by which characteristic landforms result from the thawing of ice-rich permafrost or the melting of massive ice (van Everdingen, 1998), or, more specifically, the thawing of ice-rich permafrost and (or) melting of massive ice that result in consolidation and deformation of the soil surface and formation of specific forms of relief (Shur, 1988). Jorgenson (2013) identifies 23 distinct thermokarst and other thaw-related features in the Arctic, Subarctic, and Antarctic based primarily on differences in terrain condition, ground-ice volume, and heat and mass transfer processes. Typical Arctic thermokarst landforms include thermokarst lakes, collapsed pingos, sinkholes, and pits. Thermokarst is differentiated from thermal erosion, which refers to the erosion of the land surface by thermal and mechanical processes (Mackay, 1970; van Everdingen, 1998). Typical thermal erosional features include thermo-erosional gullies. Thermal abrasion is further differentiated from thermokarst and thermal erosion by association with the reworking of ocean, river, and lake bluffs (Are, 1988). Typical thermo-abrasion features include erosional niches at the base of bluffs. Thermal denudation is another distinct term that refers to the effect of incoming solar energy on the thaw of frozen slopes and permafrost bodies that subsequently become transported downhill by gravity (Shur and Osterkamp, 2007). Active layer detachment slides and thaw slumps are typical thermal denudation features. Shur and Osterkamp (2007) noted that these various transport processes may occur together with thermokarst or in instances that would not be considered thermokarst. This compilation of references regarding thermokarst and other thaw-related features is focused on the Arctic and the Subarctic. References were drawn from North America as well as Siberia. English-language literature mostly was targeted, with 167 references annotated in version 1.0; however, an additional 28 Russian-language references were taken from Shur and Osterkamp (2007) and are provided at the end of this document. This compilation may be missing key references and inevitably will become outdated soon after publication. We hope that this document, version 1.0, will serve as the foundation for a comprehensive compilation of thermokarst and permafrost-terrain stability references, and that it will be updated continually over the coming years.

Open-File Report↗

Magmatic interactions as recorded in plagioclase phenocrysts of Chaos Crags, Lassen Volcanic Center, California

The silicic lava domes of Chaos Crags in Lassen Volcanic National Park contain a suite of variably quenched, hybrid basaltic andesite magmatic inclusions. The inclusions represent thorough mixing between rhyodacite and basalt recharge liquids accompanied by some mechanical disaggregation of the inclusions resulting in crystals mixing into the rhyodacite host preserved by quenching on dome emplacement. 87Sr/86Sr ratios (~0.7037-0.7038) of the inclusions are distinctly lower than those of the host rhyodacite (~0.704-0.7041), which are used to fingerprint the origin of mineral components and to monitor the mixing and mingling process. Chemical, isotopic, and textural characteristics indicate that the inclusions are hybrid magmas formed from the mixing and undercooling of recharge basaltic magma with rhyodacitic magma. All the host magma phenocrysts (biotite, plagioclase, hornblende and quartz crystals) also occur in the inclusions, where they are rimmed by reaction products. Compositional and strontium isotopic data from cores of unresorbed plagioclase crystals in the host rhyodacite, partially resorbed plagioclase crystals enclosed within basaltic andesite inclusions, and partially resorbed plagioclase crystals in the rhyodacitic host are all similar. Rim 87Sr/86Sr ratios of the partially resorbed plagioclase crystals in both inclusions and host are lower and close to those of the whole-rock hybrid basaltic andesite values. This observation indicates that some crystals originally crystallized in the silicic host, were partially resorbed and subsequently overgrown in the hybrid basaltic andesite magma, and then some of these partially resorbed plagioclase crystals were recycled back into the host rhyodacite. Textural evidence, in the form of sieve zones and major dissolution boundaries of the resorbed plagioclase crystals, indicates immersion of crystals into a hotter, more calcic magma. The occurrence of partially resorbed plagioclase together with plagioclase microlites and olivine crystals reflects disaggregation of inclusions and mingling of this material into the silicic host. These processes are commonplace in some orogenic magma systems and may be elucidated by isotopic microsampling and analysis of the plagioclases crystallizing from them.

Journal of Petrology↗

The future of landslides’ past—A framework for assessing consecutive landsliding systems

Landslides often happen where they have already occurred in the past. The potential of landslides to reduce or enhance conditions for further landsliding has long been recognized and has often been reported, but the mechanisms and spatial and temporal scales of these processes have previously received little specific attention. Despite a preponderance of qualitative and anecdotal evidence, there has been limited analysis. As a result, there is little consensus on the meaning of terms such as landslide repetition, recurrence, and reactivation. This source of confusion is evident when such terms are also used to describe systems where landsliding is prevalent but unrelated to landslide history. Recent findings, partly based on a rare multi-temporal landslide inventory for an area in Italy, show that the impacts of earlier landslides affect a substantial fraction of landslides, that landslides following earlier landslides differ from those that do not, and that accounting for the effect of previous landslides can improve susceptibility assessments. These findings await confirmation in other landslide prone landscapes but show that consecutive landsliding deserves more attention, which requires consistent terminology. No such terminology is presently available, and we therefore propose it in this manuscript. We use the term 'uncorrelated landsliding' to describe situations where landslides are common, but where a correlation with environmental variables such as terrain steepness is not implied. We propose 'correlated landsliding' to describe situations where landslides are common and correlations with environmental variables exist, and 'path-dependent landsliding' to describe situations where causal relations exist between consecutive landslides, for instance when landslides occur at the scarp of previous landslides. These are situations where past landslides impact future landslides. Within the path-dependent category, we distinguish three subcategories based on the spatial distance between earlier and later landslides: 'reactivation' or 'continuation' if essentially the same material recommences or continues to slide, 'local activation' if an earlier slide causes changes in a local hillslope that cause a later slide, and 'remote activation' if an earlier slide causes changes elsewhere in the landscape that cause a later landslide. We use this proposed set of terms to outline some prominent knowledge gaps and potential research questions.

Landslides↗

Theoretical basis of the borehole deepening method of absolute stress measurement

Knowledge of the initial state of stress in rocks provides a key to the solution of many problems in rock mechanics. The initial state of stress is part of the basic data required for rational design of structures in rock, since its redistribution when engineering activities are conducted is a primary load on the rock portion of the engineering system. Knowledge of initial and subsequent changes in the state of stress is also relevant in the study of active faults and crustal strain. The purposes of stress measurements are not generally well served by methods which disclose only one of the stress components at a point. Recently, interest has centered on overcoring methods which though expensive and difficult have the important advantage of yielding data on three or more stress components. In a previous paper, the authors introduced a new method--termed borehole deepening--based on measurement of the deformations that occur close to the bottom of a borehole when it is deepened by further drilling without overcoring . The proposed method is at a disadvantage when compared with overcoring procedures in that the quantities to be measured are somewhat smaller and the measuring environment is more hostile. However, the borehole deepening method presents three distinct advantages in comparison with overcoring procedures: (1) Using the same borehole, both the initial and the subsequent changes in the state of stress can be obtained by borehole deepening. Stress changes with time are very significant for stability monitoring as well as in predicting earthquake activity in regions close to known active faults and volcanoes. (2) The necessity of completely unloading a portion of the rock mass where deformation measurements are being carried out in the overcoring methods is ill-suited to nonlinearly elastic materials such as rocks. The borehole deepening method, on the other hand, is based on a change of stress from one level to another, rather than complete strain relief and is therefore less affected by nonlinear behavior. (3) Numerous stress measurements can be carried out economically in a single borehole. Thus, the regional stress distribution can be more readily obtained. Such information is considerably more significant than a single stress measurement at a point; considering the variability of rock masses and the magnitudes and direction of the stresses acting on them from one region to another and within the regions themselves, the reliability of a single measurement when extrapolated for the entire region is not very reasonable. This advantage of more numerous measurements makes it possible to predetermine the initial state of stress with confidence before actual access and construction. It also makes it possible to consider studying regional tectonic systems. Integrated regional measurements of in-situ state of stress would provide the necessary information for plotting the variation of crustal stresses; this would be of great interest in evaluating the several theories of orogenesis and in delineating tectonically active regions. This chapter is concerned with the theoretical basis for the borehole deepening method--the "software" for converting deformation measurements into stress values.

Conference Paper↗

Sediment characterization and dynamics in Lake Pontchartrain, Louisiana

Lake Pontchartrain in southeastern Louisiana is the largest of several shallow estuaries that together cover over 15,000 km 2 . Wetlands, forests, and large urban areas surround the lake. Primary transport mechanisms of sediments to Lake Pontchartrain include urban runoff, major diversions of the Mississippi River, discharge from streams along the north and west shores, and tidal circulation. Sediments deposited in Lake Pontchartrain are subjected to resuspension and mixing by natural and human activities. Bioturbation and water turbulence throughout the lake are the major mixing agents, and mechanical shell dredging has reworked much of the lake bottom over the last century. Sediment characterization through direct sampling and geophysical surveys indicates that these processes continually rework the top meter of sediment. The lake receives discharge from roadways and industrial and agricultural sources. Contaminants from these sources accumulate in the lake sediments and are an important contributor to the degradation of the estuary. Decline in populations of various benthic organisms, such as shrimp and clams, has been documented in the lake. To characterize the health of this important estuary, the U.S. Geological Survey (USGS) conducted a comprehensive evaluation of the geology, geomorphology, coastal processes, and environmental condition of the Pontchartrain Basin from 1994 to 1997. This report presents an assessment of sediment distribution and quality using a multidisciplinary approach to characterize the influence of various physical and chemical parameters: nearsurface stratigraphy, major trace metal concentrations (Cu, Pb, Zn, and Ni), and short-lived radionuclides ( 210 Pb, 7 Be, and 137 Cs). The results are compared with water-circulation patterns to determine high-resolution sedimentation patterns in the lake. The data show a significant increase in trace metals in the top 1 m of lake sediments. Above this horizon, pollen analysis indicates a correlation with land clearing in the area, a proxy for increasing human development of the surrounding landscape and an increase in surface run-off. The data also show that the top meter of sediment undergoes frequent resuspension during high-energy circulation events and via circulation gyres in the lake. This regular turnover does not allow stratification of recently deposited sediments, restricting the sequestration of contaminated material that enters the lake.

Louisiana↗

Hydrogeology and ground-water quality at a land reclamation site, Neshaminy State Park, Pennsylvania

At Neshaminy State park, the most important aquifer is the informally named 'Trenton gravel' of Pleistocene age, which consists of poorly sorted sand and gravel. This is underlain by less permeable crystalline rock that limits the downward movement of water. Up to 5 feet of Holocene (or perhaps Pleistocene) alluvium consisting of clay and silt was deposited above the Trenton gravel, but much of the surficial material is dredge spoil, mostly sand and silt from the Delaware River. Normally, ground water discharges to the Delaware River and Neshaminy Creek. However, because the Trenton gravel is confined by the alluvium and is hydraulically connected to the Delaware River throughout its entire thickness due to channel dredging, hydraulic-head variations in the river are readily transmitted to wells. Thus, flow is toward the Delaware River and Neshaminy Creek at low tide, but the river and and creek recharge the aquifer as the tide rises. In September 1981 and October 1983, infertile dredge spoil from the Delaware River at Neshaminy State Park was mixed with composted sewage sludge and seeded to grasses to produce a soil product that has not adversely affected ground-water quality as of March 1985. Analyses for pH, specific conductance, nutrients, and heavy metals in ground water at the beginning of the study indicate that the chemistry of the ground water in the dredge spoil at the site differs from that in the Trenton gravel. Analyses of ground-water samples collected after the first two sludge applications (120 tons per acre and 450 tons per acre), indicate that no significant change occurred in the chemistry of the samples from the Trenton gravel, whereas organic nitrogen increased temporarily in ground water from the dredge spoil 6 months after the larger of the two sludge applications, but quickly returned to background levels. The lack of chemical change with time in the ground water implies either that little of the more than 100 inches of precipitation that fell from April 1983 through March 1985 reached the water table or, more likely, that a mechanism exists beneath the soil- factory site that retards or prevents the downard migration of contaminants.

Pennsylvania↗

Retardations in fault creep rates before local moderate earthquakes along the San Andreas fault system, central California

Records of shallow aseismic slip (fault creep) obtained along parts of the San Andreas and Calaveras faults in central California demonstrate that significant changes in creep rates often have been associated with local moderate earthquakes. An immediate postearthquake increase followed by gradual, long-term decay back to a previous background rate is generally the most obvious earthquake effect on fault creep. This phenomenon, identified as aseismic afterslip, usually is characterized by above-average creep rates for several months to a few years. In several cases, minor step-like movements, called coseismic slip events, have occurred at or near the times of mainshocks. One extreme case of coseismic slip, recorded at Cienega Winery on the San Andreas fault 17.5 km southeast of San Juan Bautista, consisted of 11 mm of sudden displacement coincident with earthquakes of ML=5.3 and ML=5.2 that occurred 2.5 minutes apart on 9 April 1961. At least one of these shocks originated on the main fault beneath the winery. Creep activity subsequently stopped at the winery for 19 months, then gradually returned to a nearly steady rate slightly below the previous long-term average. The phenomena mentioned above can be explained in terms of simple models consisting of relatively weak material along shallow reaches of the fault responding to changes in load imposed by sudden slip within the underlying seismogenic zone. In addition to coseismic slip and afterslip phenomena, however, pre-earthquake retardations in creep rates also have been observed. Onsets of significant, persistent decreases in creep rates have occurred at several sites 12 months or more before the times of moderate earthquakes. A 44-month retardation before the 1979 ML=5.9 Coyote Lake earthquake on the Calaveras fault was recorded at the Shore Road creepmeter site 10 km northwest of Hollister. Creep retardation on the San Andreas fault near San Juan Bautista has been evident in records from one creepmeter site for the past 5 years. Retardations with durations of 21 and 19 months also occurred at Shore Road before the 1974 and 1984 earthquakes of ML=5.2 and ML=6.2, respectively. Although creep retardation remains poorly understood, several possible explanations have been discussed previously. (1) Certain onsets of apparent creep retardation may be explained as abrupt terminations of afterslip generated from previous moderate-mainshock sequences. (2) Retardations may be related to significant decreases in the rate of seismic and/or aseismic slip occurring within or beneath the underlying seismogenic zone. Such decreases may be caused by changes in local conditions related to growth of asperities, strain hardening, or dilatancy, or perhaps by passage of stress-waves or other fluctuations in driving stresses. (3) Finally, creep rates may be lowered (or increased) by stresses imposed on the fault by seismic or aseismic slip on neighboring faults. In addition to causing creep-rate increases or retardations, such fault interactions occasionally may trigger earthquakes. Regardless of the actual mechanisms involved and the current lack of understanding of creep retardation, it appears that shallow fault creep is sensitive to local and regional effects that promote or accompany intermediate-term preparation stages leading to moderate earthquakes. A strategy for more complete monitoring of fault creep, wherever it is known to occur, therefore should be assigned a higher priority in our continuing efforts to test various hypotheses concerning the mechanical relations between seismic and aseismic slip. ?? 1988 Birkha??user Verlag.

Pure and Applied Geophysics PAGEOPH↗

Rockslide-debris avalanche of May 18, 1980, Mount St. Helens Volcano, Washington

This report provides a detailed picture of the rockslide-debris avalanche of the May 18, 1980, eruption of Mount St. Helens volcano. It provides a characterization of the deposit, a reinterpretation of the details of the first minutes of the eruption of May 18, and insight into the transport mechanism of the mass movement. Details of the rockslide event, as revealed by eyewitness photographs, are correlated with features of the deposit. The photographs show three slide blocks in the rockslide movement. Slide block I was triggered by a magnitude 5.1 earthquake at 8:32 a.m. Pacific Daylight Time (P.D.T.). An exploding cryptodome burst through slide block II to produce the 'blast surge.' Slide block III consisted of many discrete failures that were carried out in continuing pyroclastic currents generated from the exploding cryptodome. The cryptodome continued to depressurize after slide block III, producing a blast deposit that rests on top of the debris-avalanche deposit. The hummocky 2.5 cubic kilometer debris-avalanche deposit consists of block facies (pieces of the pre-eruption Mount St. Helens transported relatively intact) and matrix facies (a mixture of rocks from the old mountain and cryptodome dacite). Block facies is divided into five lithologic units. Matrix facies was derived from the explosively generated current of slide block III as well as from disaggregation and mixing of debris-avalanche blocks. The mean density of the old cone was measured to be abut 20 percent greater than the mean density of the avalanche deposit. Density in the deposit does not decrease with distance which suggests that debris-avalanche blocks were dilated at the mountain, rather than during transport. Various grain-size parameters that show that clast size converges about a mean with distance suggest mixing during transport. The debris-avalanche flow can be considered a grain flow, where particles -- either debris-avalanche blocks or the clasts within the blocks -- collided and created dispersive stress normal to the movement of material. The dispersive stress preserved the dilation of the material and allowed it to flow.

Washington↗

Volcanic ash activates the NLRP3 inflammasome in murine and human macrophages

Volcanic ash is a heterogeneous mineral dust that is typically composed of a mixture of amorphous (glass) and crystalline (mineral) fragments. It commonly contains an abundance of the crystalline silica (SiO 2 ) polymorph cristobalite. Inhalation of crystalline silica can induce inflammation by stimulating the NLRP3 inflammasome, a cytosolic receptor complex that plays a critical role in driving inflammatory immune responses. Ingested material results in the assembly of NLRP3, ASC, and caspase-1 with subsequent secretion of the interleukin-1 family cytokine IL-1β. Previous toxicology work suggests that cristobalite-bearing volcanic ash is minimally reactive, calling into question the reactivity of volcanically derived crystalline silica, in general. In this study, we target the NLRP3 inflammasome as a crystalline silica responsive element to clarify volcanic cristobalite reactivity. We expose immortalized bone marrow-derived macrophages of genetically engineered mice and primary human peripheral blood mononuclear cells (PBMCs) to ash from the Soufrière Hills volcano as well as representative, pure-phase samples of its primary componentry (volcanic glass, feldspar, cristobalite) and measure NLRP3 inflammasome activation. We demonstrate that respirable Soufrière Hills volcanic ash induces the activation of caspase-1 with subsequent release of mature IL-1β in a NLRP3 inflammasome-dependent manner. Macrophages deficient in NLRP3 inflammasome components are incapable of secreting IL-1β in response to volcanic ash ingestion. Cellular uptake induces lysosomal destabilization involving cysteine proteases. Furthermore, the response involves activation of mitochondrial stress pathways leading to the generation of reactive oxygen species. Considering ash componentry, cristobalite is the most reactive pure-phase with other components inducing only low-level IL-1β secretion. Inflammasome activation mediated by inhaled ash and its potential relevance in chronic pulmonary disease was further evidenced in PBMC using the NLRP3 small-molecule inhibitor CP-456,773 (CRID3, MCC950). Our data indicate the functional activation of the NLRP3 inflammasome by volcanic ash in murine and human macrophages in vitro . Cristobalite is identified as the apparent driver, thereby contesting previous assertions that chemical and structural imperfections may be sufficient to abrogate the reactivity of volcanically derived cristobalite. This is a novel mechanism for the stimulation of a pro-inflammatory response by volcanic particulate and provides new insight regarding chronic exposure to environmentally occurring particles.

Frontiers in Immunology↗

Timescales of carbon turnover in soils with mixed crystalline mineralogies

Organic matter–mineral associations stabilize much of the carbon (C) stored globally in soils. Metastable short-range-order (SRO) minerals such as allophane and ferrihydrite provide one mechanism for long-term stabilization of organic matter in young soil. However, in soils with few SRO minerals and a predominance of crystalline aluminosilicate or Fe (and Al) oxyhydroxide, C turnover should be governed by chemisorption with those minerals. Here, we correlate mineral composition from soils containing small amounts of SRO minerals with mean turnover time (TT) of C estimated from radiocarbon ( 14 C) in bulk soil, free light fraction and mineral-associated organic matter. We varied the mineral amount and composition by sampling ancient soils formed on different lithologies in arid to subhumid climates in Kruger National Park (KNP), South Africa. Mineral contents in bulk soils were assessed using chemical extractions to quantify Fe oxyhydroxides and SRO minerals. Because of our interest in the role of silicate clay mineralogy, particularly smectite (2 : 1) and kaolinite (1 : 1), we separately quantified the mineralogy of the clay-sized fraction using X-ray diffraction (XRD) and measured 14 C on the same fraction. Density separation demonstrated that mineral associated C accounted for 40–70 % of bulk soil organic C in A and B1 horizons for granite, nephelinite and arid-zone gabbro soils, and > 80 % in other soils. Organic matter strongly associated with the isolated clay-sized fraction represented only 9–47 % of the bulk soil C. The mean TT of C strongly associated with the clay-sized fraction increased with the amount of smectite (2 : 1 clays); in samples with > 40 % smectite it averaged 1020 ± 460 years. The C not strongly associated with clay-sized minerals, including a combination of low-density C, the C associated with minerals of sizes between 2 µm and 2 cm (including Fe oxyhydroxides as coatings), and C removed from clay-sized material by 2 % hydrogen peroxide had TTs averaging 190 ± 190 years in surface horizons. Summed over the bulk soil profile, we found that smectite content correlated with the mean TT of bulk soil C across varied lithologies. The SRO mineral content in KNP soils was generally very low, except for the soils developed on gabbros under more humid climate that also had very high Fe and C contents with a surprisingly short, mean C TTs. In younger landscapes, SRO minerals are metastable and sequester C for long timescales. We hypothesize that in the KNP, SRO minerals represent a transient stage of mineral evolution and therefore lock up C for a shorter time. Overall, we found crystalline Fe-oxyhydroxides (determined as the difference between Fe in dithionate citrate and oxalate extractions) to be the strongest predictor for soil C content, while the mean TT of soil C was best predicted from the amount of smectite, which was also related to more easily measured bulk properties such as cation exchange capacity or pH. Combined with previous research on C turnover times in 2 : 1 vs. 1 : 1 clays, our results hold promise for predicting C inventory and persistence based on intrinsic timescales of specific carbon–mineral interactions.

SOIL↗

Use of an intact core and stable-metal isotopes to examine leaching characteristics of a fluvial tailings deposit

The upper Arkansas River south of Leadville, Colorado, USA, contains deposits of fluvial tailings from historical mining operations in the Leadville area. These deposits are possible non-point sources of acid and metal contamination to surface- and ground-water systems. We used stable-metal isotopes to help ascertain metal retention and release mechanisms that influence metal transport from the deposits to shallow ground-water systems. To accomplish this, we excavated an intact core from a small fluvial tailings deposit and performed laboratory column experiments to examine the amount of metals leaching through the core. Deionised water was continuously applied to the top of the core, and effluent was collected from the bottom of the core for a period of about two months. Three saturation conditions were used, including unsaturated, partially saturated, and fully saturated. Reducing conditions developed upon partial and complete saturation of the core. During the unsaturated leaching phase, core effluents were acidic (pH 2.8 - 3.5) and contained elevated metal concentrations. During a portion of the unsaturated leaching phase, stable-metal isotope spikes and conservative tracers were simultaneously applied to the top of the core. The isotopes provide a means to delineate metal behaviour within the core. During the unsaturated leaching phase, the isotope spikes did not exhibit breakthrough with simultaneously added conservative tracers. However, some of the stable-metal isotope spikes did break through when reducing conditions occurred in the core (during the partial- and completesaturation conditions). In this paper, we use Cd as an example of the utility of stable-metal isotopes in geochemical studies. In the case of Cd, after the core was partially saturated, the 111 Cd spike was released as evidenced by a change in the Cd isotope ratios in the effluent. This release continued during the fully saturated leaching phase, however, the total Cd concentration did not increase. These results suggest that the 111 Cd spike was retained inside the core during the unsaturated leaching phase, and only partially released as reducing conditions developed. Results from this core-leaching experiment indicate there is a large reservoir of water-soluble material within the fluvial tailings deposit, which yields elevated metal concentrations and high acidity, and which may degrade adjacent ground- and surface-water quality. Use of stable metal isotopes in this study facilitated the determination of different metal-retention processes, metal-release processes, and metal sources in the fluvial tailings deposit in response to changing geochemical conditions.

Conference Paper↗

Physical data of soil profiles formed on late Quaternary marine terraces near Santa Cruz, California

The marine terraces in and around Santa Cruz, California, represent a set of well-preserved terraces formed as a product of geology, sea level, and climate. A marine terrace begins as a wave cut platform. Eustatic sea level changes, seacliff erosion, and tectonic uplift work together to generate marine terraces. "When a wave-cut platform is raised (due to tectonic activity) above sea level and cliffed by wave action it becomes a marine terrace" (Bradley, 1957, p. 424). During glacial periods, eustatic sea level is estimated to have dropped by 150 meters (Fairbanks, 1989). Cliff retreat measured from aerial photographs between 1930 and 1980 vary from 0.0 to 0.2 m yr–1 (Best and Griggs, 1991). Estimates of uplift rates along the Santa Cruz coastline vary from 0.10 to 0.48 m kyr–1 (Bradley and Griggs, 1976; Weber and others, 1999). Uplift mechanisms include coseismic uplift associated both with a reverse component of slip on the steeply SW dipping Loma Prieta fault in the restraining bend of the San Andreas Fault and a small component of reverse slip on the steeply SE dipping San Gregorio fault (Anderson and Menking 1994). Previous work studying physical properties on these terraces include Pinney and others (in press) and Aniku (1986) and Bowman and Estrada (1980). Sedimentary deposits of the marine terraces are a mixture of terrestrial and marine sediments but generally consist of a sheet of marine deposits overlying the old platform and a wedge of nonmarine deposits banked against the old sea cliff (Bradley, 1957). Bedrock underlying the terraces in the Santa Cruz area is generally either Santa Margarita Sandstone or Santa Cruz Mudstone. The Santa Margarita Sandstone represents an upper Miocene, transgressive, tidally dominated marine-shelf deposit with crossbedded sets of sand and gravel and horizontally stratified and bioturbated invertebrate-fossils beds (Phillips, 1990). The siliceous Santa Cruz Mudstone, of late Miocene age, conformably overlies the Santa Margarita Sandstone. The Santa Cruz Mudstone is a thin to medium-bedded siliceous mudstone with nonsiliceous mudstone and siltstone and minor amounts of sandstone. The siliceous nature implies organic deposition in a quiescent, deep-water environment. Bedrock is mantled by 1–4 meters of medium to coarse-grained regressive beach sediment and fluvial deposits from the Ben Lomond Mountains. Terrace age increases with elevation above sea level, and weathering of primary minerals increases with age. The suite of soils formed on the terraces is referred to as a soil chronosequence. Soil chronosequences, important tools in characterizing natural weathering rates, are defined as a group of soils that differ in age and therefore in duration of weathering but have similar climatic conditions, vegetation, geomorphic position, and parent material (Jenny, 1941; Birkland, 1999). Soils are frequently useful indicators of geomorphic age (Muhs, 1982; Switzer and others, 1988) and are a function of pedogenic and/or eolian processes. Some aspects of soil development can be episodic but when viewed on large time scales can be perceived as continuous (Switzer and others, 1988). The age of the soil may be constrained by the age of the deposit, since soil formation generally commences when deposition has ceased (Birkland, 1999). Dating of the terraces provides an unprecedented opportunity to study weathering and soil-formation rates (Perg and others, 2001; Hanks and others, 1984; Bradley and Griggs, 1976; Bradley and Addicott, 1968; Bradley, 1956). Ages of the terraces recently dated by cosmogenic radionuclide are, starting with the youngest, 65, 92, 137, 139, and 226 k.y. (Perg and others, 2001). However, these ages are much younger than recent radiometric dates on mollusk shells (Muhs, U.S. Geological Survey, personal communication, 2002; Bradley and Addicott, 1968). For this study, soils were sampled on five terraces. Terrace one in the Lighthouse Field along Westcliff in Santa Cruz was the last site selected, and this report contains minimal data on this terrace. Sites on the second, third, and fourth terraces are located in Wilder Ranch, Santa Cruz, California. Site five is on private property north of Wilder Ranch. Careful consideration was taken in selecting field sites, choosing locations in a topographically flat area to avoid effects of erosion, and trying to keep parent material similar. This report contains physical properties of the soil profiles on four of the five marine terraces near Santa Cruz, California, excluding the youngest terrace in all tables except 6 and 7. Data includes field descriptions, bulk density, grain size analyses, weight percent magnetic fraction, and the soil development index. Soil properties are important when trying to understand the chemistry of a given profile or when comparing profiles. Grain size constrains the movement of water in a profile, thus controlling movement of chemicals and weathering rates. Bulk density is a useful property to calculate chemical inventory. Quantifying the magnetic fraction aids in understanding the Fe inventory for these soils. The soil development index is a semi-quantitative way to define the degree of development of a soil profile. This is a useful way to compare development of profiles for this chronosequence or compare the Santa Cruz terraces to a suit of other terraces or another chronosequence.

California↗

A note on the effect of fault gouge composition on the stability of frictional sliding

The frictional properties of fault gouge have been studied at confining pressures to 6 kbars. If the gouge is composed of strong materials such as crushed granite or quartz sand, the frictional strength is high, and violent stick-slip occurs at confining pressures above approximately 1.5 kbars. If the gouge is composed of minerals such as illite, kaolinite, chlorite, or antigorite, which have weak bonding forces between the structural layers, the frictional strength is slightly lower, but violent stick-slip still occurs under high confining pressure. The expanding clays, montmorillonite and vermiculite, which have free water between their structural layers, slide stably at confining pressures as high as 6.25 kbars and exhibit low friction. A similar stable behavior with lowered strength is observed in water-saturated quartz sand when the water is confined within the fault zone during deformation. The results of this series of experiments support water being the stabilizing influence when it is either (1) trapped within or between rocks of low permeability and can provide a high pore pressure when the rocks are deformed, or (2) loosely bonded in a mineral structure, as in the hydrated clays, where it can produce a pseudo-pore pressure when the clay is compressed. In both these cases, the effective stress can be reduced and the deformation stabilized. ?? 1977.

International Journal of Rock Mechanics and Mining↗

Poroelastic rebound along the Landers 1992 earthquake surface rupture

Maps of surface displacement following the 1992 Landers, California, earthquake, generated by interferometric processing of ERS-1 synthetic aperture radar (SAR) images, reveal effects of various postseismic deformation processes along the 1992 surface rupture. The large-scale pattern of the postseismic displacement field includes large lobes, mostly visible on the west side of the fault, comparable in shape with the lobes observed in the coseismic displacement field. This pattern and the steep displacement gradient observed near the Emerson-Camp Rock fault cannot be simply explained by afterslip on deep sections of the 1992 rupture. Models show that horizontal slip occurring on a buried dislocation in a Poisson's material produces a characteristic quadripole pattern in the surface displacement field with several centimeters of vertical motion at distances of 10–20 km from the fault, yet this pattern is not observed in the postseismic interferograms. As previously proposed to explain local strain in the fault step overs [ Peltzer et al ., 1996b], we argue that poroelastic rebound caused by pore fluid flow may also occur over greater distances from the fault, compensating the vertical ground shift produced by fault afterslip. Such a rebound is explained by the gradual change of the crustal rocks' Poisson's ratio value from undrained (coseismic) to drained (postseismic) conditions as pore pressure gradients produced by the earthquake dissipate. Using the Poisson's ratio values of 0.27 and 0.31 for the drained and undrained crustal rocks, respectively, elastic dislocation models show that the combined contributions of afterslip on deep sections of the fault and poroelastic rebound can account for the range change observed in the SAR data and the horizontal displacement measured at Global Positioning System (GPS) sites along a 60-km-long transect across the Emerson fault [ Savage and Svarc , 1997]. Using a detailed surface slip distribution on the Homestead Valley, Kickapoo, and Johnson Valley faults, we modeled the poroelastic rebound in the Homestead Valley pull apart. A Poisson's ratio value of 0.35 for the undrained gouge rocks in the fault zone is required to account for the observed surface uplift in the 3.5 years following the earthquake. This large value implies a seismic velocity ratio V p /V s of 2.1, consistent with the observed low V s values of fault zone guided waves at shallow depth [ Li et al ., 1997]. The SAR data also reveal postseismic creep along shallow patches of the Eureka Peak and Burnt Mountain faults with a characteristic decay time of 0.8 years. Coseismic, dilatant hardening (locking process) followed by post-seismic, pore pressure controlled fault creep provide a plausible mechanism to account for the decay time of the observed slip rate along this section of the fault.

Journal of Geophysical Research B: Solid Earth↗

Application of a hydrodynamic and sediment transport model for guidance of response efforts related to the Deepwater Horizon oil spill in the Northern Gulf of Mexico along the coast of Alabama and Florida

U.S. Geological Survey (USGS) scientists have provided a model-based assessment of transport and deposition of residual Deepwater Horizon oil along the shoreline within the northern Gulf of Mexico in the form of mixtures of sand and weathered oil, known as surface residual balls (SRBs). The results of this USGS research, in combination with results from other components of the overall study, will inform operational decisionmaking. The results will provide guidance for response activities and data collection needs during future oil spills. In May 2012 the U.S. Coast Guard, acting as the Deepwater Horizon Federal on-scene coordinator, chartered an operational science advisory team to provide a science-based review of data collected and to conduct additional directed studies and sampling. The goal was to characterize typical shoreline profiles and morphology in the northern Gulf of Mexico to identify likely sources of residual oil and to evaluate mechanisms whereby reoiling phenomena may be occurring (for example, burial and exhumation and alongshore transport). A steering committee cochaired by British Petroleum Corporation (BP) and the National Oceanic and Atmospheric Administration (NOAA) is overseeing the project and includes State on-scene coordinators from four States (Alabama, Florida, Louisiana, and Mississippi), trustees of the U.S. Department of the Interior (DOI), and representatives from the U.S. Coast Guard. This report presents the results of hydrodynamic and sediment transport models and developed techniques for analyzing potential SRB movement and burial and exhumation along the coastline of Alabama and Florida. Results from these modeling efforts are being used to explain the complexity of reoiling in the nearshore environment and to broaden consideration of the different scenarios and difficulties that are being faced in identifying and removing residual oil. For instance, modeling results suggest that larger SRBs are not, under the most commonly observed low-energy wave conditions, likely to move very far alongshore. This finding suggests that SRBs from one source location may not (outside of storm conditions) be redistributed to other up or down coast locations. This information can guide operational response decisions. In addition, because SRBs are less mobile compared with sand, they are likely to become buried and unburied under normal sand transport processes thereby lengthening the time SRBs may take to move onshore. The rate of onshore movement was not specifically addressed by this study, yet the results resolve the cross-shore domain and cross-shore variations in alongshore transport that are relevant to achieving the primary objectives. Furthermore, during infrequent events (for example, winter storms and severe meteorological events such as Hurricane Isaac of August 2012), energy is shown to be sufficient to move a greater range of SRB sizes and potentially expose and break up submerged oil mats. When SRBs do move alongshore, the models indicate that there are regions that are more conducive to accumulation of SRB material than others. Accumulation can occur where there are reversals and decelerations in alongshore currents and where forces created by shear stress drops below critical thresholds to maintain or initiate SRB movement. In addition, flow and SRB mobility patterns around inlets indicate patterns in hydrodynamic forces that influence redistribution of SRBs and the surface oil that mixed with sediment to form oil mats in the first place.

Alabama;Florida↗

Topography of closed depressions, scarps, and grabens in the North Tharsis region of Mars: implications for shallow crustal discontinuities and graben formation

Using Viking Orbiter images, detailed photoclinometric profiles were obtained across 10 irregular depressions, 32 fretted fractures, 49 troughs and pits, 124 solitary scarps, and 370 simple grabens in the north Tharsis region of Mars. These data allow inferences to be made on the shallow crustal structure of this region. The frequency modes of measured scarp heights correspond with previous general thickness estimates of the heavily cratered and ridged plains units. The depths of the flat-floored irregular depressions (55-175 m), fretted fractures (85-890 m), and troughs and pits (60-1620 m) are also similar to scarp heights (thicknesses) of the geologic units in which these depressions occur, which suggests that the depths of these flat-floored features were controlled by erosional base levels created by lithologic contacts. Although the features have a similar age, both their depths and their observed local structural control increase in the order listed above, which suggests that the more advanced stages of associated fracturing facilitated the development of these depressions by increasing permeability. If a ground-ice zone is a factor in development of these features, as has been suggested, our observation that the depths of these features decrease with increasing latitude suggests that either the thickness of the ground-ice zone does not increase poleward or the depths of the depressions were controlled by the top of the ground-ice zone whose depth may decrease with latitude. Deeper discontinuities are inferred from fault-intersection depths of 370 simple grabens (assuming 60° dipping faults that initiate at a mechanical discontinuity) in Tempe Terra and Alba Patera and from the depths of the large, flat-floored troughs in Tempe Terra. The frequency distributions of these fault-intersection and large trough depths show a concentration at 1.0-1.6 km depth, similar to data obtained for Syria, Sinai, and Lunae Plana. The consistency of these depth data over such a large region of western Mars suggests that a discontinuity or a process that transcends local and regional geology is responsible for the formation of these features. If this discontinuity is represented by the base of the cryosphere, its uniform depth over 55° of latitude suggests that the cryosphere did not thicken poleward. Alternatively, the concentration of depths at 1.0-1.6 km may represent the upper level of noneruptive dike ascent (lateral dike propagation) of Mars, which is controlled by gravity and atmospheric pressure and magma and country-rock characteristics, and was probably controlled, in part, by ground ice. Fault-intersection depths in the north Tharsis region locally extend down to a depth of 5-7 km. The depth data between 2 and 3 km are attributed to the discontinuity at the interface of megaregolith and basement or to the upper limit of noneruptive dike ascent of magma with a high volatile content. Intersection depths greater than 3 km, which were found at Alba Patera, may be due to the megaregolith-basement discontinuity, which was buried and depressed by volcanic loading, or to the upper level of noneruptive dike ascent of magma with a low volatile content. The near absence of narrow simple grabens with fault-initiation depths less than 0.6-1.0 km in this study area, as well as in most of western Mars, suggests that this depth represents the minimum depth that normal faults can initiate; at shallower depths tension cracks or joints would form instead. This hypothesis is supported by the application of the Griffith failure criterion to this minimum depth of normal fault initiation, which suggests that shallow crustal materials have a tensile strength of 2-4 MPa throughout most of western Mars, in close agreement with previous estimates of tensile strength of martian basaltic rock.

Icarus↗

River channel patterns: Braided, meandering, and straight

Channel pattern is used to describe the plan view of a reach of river as seen from an airplane, and includes meandering, braiding, or relatively straight channels. Natural channels characteristically exhibit alternating pools or deep reaches and riffles or shallow reaches, regardless of the type of pattern. The length of the pool or distance between riffles in a straight channel equals the straight line distance between successive points of inflection in the wave pattern of a meandering river of the same width. The points of inflection are also shallow points and correspond to riffles in the straight channel. This distance, which is half the wavelength of the meander, varies approximately as a linear function of channel width. In the data we analysed the meander wavelength, or twice the distance between successive riffles, is from 7 to 12 times the channel width. It is concluded that the mechanics which may lead to meandering operate in straight channels. River braiding is characterized by channel division around alluvial islands. The growth of an island begins as the deposition of a central bar which results from sorting and deposition of the coarser fractions of the load which locally cannot be transported. The bar grows downstream and in height by continued deposition on its surface, forcing the water into the flanking channels, which, to carry the flow, deepen and cut laterally into the original banks. Such deepening locally lowers the water surface and the central bar emerges as an island which becomes stabilized by vegetation. Braiding was observed in a small river in a laboratory. Measurements of the adjustments of velocity, depth, width, and slope associated with island development lead to the conclusion that braiding is one of the many patterns which can maintain quasi-equilibrium among discharge, load, and transporting ability. Braiding does not necessarily indicate an excess of total load. Channel cross section and pattern are ultimately controlled by the discharge and load provided by the drainage basin. It is important, therefore, to develop a picture of how the several variables involved in channel shape interact to result in observed channel characteristics. Such a rationale is summarized as follows: Channel width appears to be primarily a function of near-bankfull discharge, in conjunction with the inherent resistance of bed and bank to scour. Excessive width increases the shear on the bed at the expense of that on the bank and the reverse is true for very narrow widths. Because at high stages width adjustment can take place rapidly and with the evacuation or deposition of relatively small volumes of debris, achievement of a ,relatively stable width at high flow is a primary adjustment to which the further interadjustments between depth, velocity, slope, and roughness tend to accommodate. Channel roughness, to the extent that it is determined by particle size, is an independent factor related to the drainage basin rather than to the channel. Roughness in streams carrying fine material, however, is also a function of the dunes or other characteristics of bed configuration. Where roughness is independently determined as well as discharge and load, these studies indicate that a particular slope is associated with the roughness. At the width determined by the discharge, velocity and depth must be adjusted to satisfy quasi-equilibrium in accord with the particular slope. But if roughness also is variable, depending on the transitory configuration of the bed, then a number of combinations of velocity, depth, and slope will satisfy equilibrium. An increase in load at constant discharge, width, and caliber of load tends to be associated with an increasing slope if the roughness (dune or bed configuration) changes with the load. In the laboratory river an increase of load at constant discharge, width, and caliber resulted in progressive aggradation of long reaches of channel at constant slope. The adjustments of several variables tending toward the establishment of quasi-equilibrium in river channels lead to the different channel patterns observed in nature. For example, the data indicate that at a given discharge, meanders occur at smaller values of slope than do’ braids. Further, at the same slope braided channels are associated with higher bankfull discharges than are meanders. An additional example is provided by the division of discharge around islands in braided rivers which produces numerous small channels. The changes in slope, roughness, and channel shape which accompany this division are in accord with quasi-equilibrium adjustments observed in the comparison of large and small rivers.

Professional Paper↗

Approaches to modeling weathered regolith

Sustainable soils are a requirement for maintaining human civilizations ( Carter and Dale 1974 ; Lal 1989 ). However, as the “most complicated biomaterial on the planet” ( Young and Crawford 2004 ), soils represent one of the most difficult systems to understand and model with respect to chemical, physical, and biological coupling over time (Fig. 1 ). Despite the complexity of these interactions, certain patterns in soil properties and development are universally observed and have been used in soil science as a means for classification. Elemental, mineralogical, or isotopic concentrations in soils plotted versus depth beneath the land surface comprise such patterns. Soil depth profiles are often reported for solid soil materials, and, less frequently, for solutes in soil pore waters. These profiles cross a large range in spatial scales that traditionally have been studied by different disciplines. For example, shallow, biologically active horizons are commonly defined as the soil zone in agronomic studies whereas the mobile layer of the regolith is referred to as soil in geomorphological studies. In contrast, many geochemical studies target chemical weathering to tens or even hundreds of meters in depth, sometimes extending the definition of “soils” to include the entire regolith down to parent bedrock or alluvium. Soil profiles also exhibit a large range in temporal scales ( Amundson 2004 ; Brantley 2008b ). Solid-state profiles document chemical and mineralogical changes integrated over the time scales of evolution of regolith from protolith. This “geologic time” can vary from tens to hundreds of years for weathered material developed on moraines deposited by active glaciers ( Anderson et al. 1997 ), to millions or possibly hundreds of millions of years of regolith evolution as documented in laterites and bauxites on stable cratons ( Nahon 1986 ). In contrast, solute profiles reflect much shorter time scales corresponding to the residence time of the soil water which commonly ranges from days to decades ( Stonestrom et al. 1998 ). Factors impacting soil minerals can therefore be related to geologically old processes while those impacting pore waters are related to contemporary processes. We first discuss a geochemical frame work for modeling soil profiles, including a simple scheme that depends on the extent of enrichment or depletion. Such profiles are comprised of reaction fronts affected by chemical, hydrologic, geologic and biologic processes that control soil evolution. We then present a hierarchy of models that have been used to interpret both solid state and solute compositions in regolith. The more simple approaches to model depletion in soils, using analytical models, are first described. The most elementary of these is a linear model that calculates rate constants from the slopes of either solid or solute weathering gradients: these rate constants represent lumped parameters that describe weathering in terms of an integrated reaction rate. Two other analytical models are then presented that have been used to fit solid state elemental profiles with exponential and sigmoidal functions. All of these analytical approaches are derived for models of soils as containing a limited number of components, phases, and species. At a more complex level, numerical models are then presented to elucidate how kinetic and transport parameters as well as chemical, hydrologic, and physical soil data can be incorporated. We consider two forms of these models, first relatively simple spreadsheet calculators and then more sophisticated multi-component, multi-phase reactive-transport numerical codes. Our treatment of reactive transport modeling is relatively cursory, in recognition of the treatment in the chapter by Steefel and Maher (2009 , this volume). Because these models incorporate more phases, components, and species than the other approaches and explicitly model the more fundamental reaction mechanisms involved, they generally have a greater need for parameterization. In our conclusion section, we discuss how this hierarchy of approaches can yield generalizations about soils that are often complementary.

Reviews in Mineralogy and Geochemistry↗