Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Earth Observer”

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 955 records · Page 53Linked to original sources
USGS1988

Seismic design and engineering research at the U.S. Geological Survey

The Engineering Seismology Element of the USGS Earthquake Hazards Reduction Program is responsible for the coordination and operation of the National Strong Motion Network to collect, process, and disseminate earthquake strong-motion data; and, the development of improved methodologies to estimate and predict earthquake ground motion. Instrumental observations of strong ground shaking induced by damaging earthquakes and the corresponding response of man-made structures provide the basis for estimating the severity of shaking from future earthquakes, for earthquake-resistant design, and for understanding the physics of seismologic failure in the Earth's crust.

Earthquakes & Volcanoes (USGS)↗
USGS2016

Terrestrial nitrogen cycling in Earth system models revisited

Understanding the degree to which nitrogen (N) availability limits land carbon (C) uptake under global environmental change represents an unresolved challenge. First-generation ‘C-only’vegetation models, lacking explicit representations of N cycling,projected a substantial and increasing land C sink under rising atmospheric CO 2 concentrations. This prediction was questioned for not taking into account the potentially limiting effect of N availability, which is necessary for plant growth (Hungate et al.,2003). More recent global models include coupled C and N cycles in land ecosystems (C–N models) and are widely assumed to be more realistic. However, inclusion of more processes has not consistently improved their performance in capturing observed responses of the global C cycle (e.g. Wenzel et al., 2014). With the advent of a new generation of global models, including coupled C, N, and phosphorus (P) cycling, model complexity is sure to increase; but model reliability may not, unless greater attention is paid to the correspondence of model process representations ande mpirical evidence. It was in this context that the ‘Nitrogen Cycle Workshop’ at Dartington Hall, Devon, UK was held on 1–5 February 2016. Organized by I. Colin Prentice and Benjamin D. Stocker (Imperial College London, UK), the workshop was funded by the European Research Council,project ‘Earth system Model Bias Reduction and assessing Abrupt Climate change’ (EMBRACE). We gathered empirical ecologists and ecosystem modellers to identify key uncertainties in terrestrial C–N cycling, and to discuss processes that are missing or poorly represented in current models.

New Phytologist↗
USGS2020

The pervasive and multifaceted influence of biocrusts on water in the world’s drylands

The capture and use of water are critically important in drylands, which collectively constitute Earth's largest biome. Drylands will likely experience lower and more unreliable rainfall as climatic conditions change over the next century. Dryland soils support a rich community of microphytic organisms (biocrusts), which are critically important because they regulate the delivery and retention of water. Yet despite their hydrological significance, a global synthesis of their effects on hydrology is lacking. We synthesized 2,997 observations from 109 publications to explore how biocrusts affected five hydrological processes (times to ponding and runoff, early [sorptivity] and final [infiltration] stages of water flow into soil, and the rate or volume of runoff) and two hydrological outcomes (moisture storage, sediment production). We found that increasing biocrust cover reduced the time for water to pond on the surface (−40%) and commence runoff (−33%), and reduced infiltration (−34%) and sediment production (−68%). Greater biocrust cover had no significant effect on sorptivity or runoff rate/amount, but increased moisture storage (+14%). Infiltration declined most (−56%) at fine scales, and moisture storage was greatest (+36%) at large scales. Effects of biocrust type (cyanobacteria, lichen, moss, mixed), soil texture (sand, loam, clay), and climatic zone (arid, semiarid, dry subhumid) were nuanced. Our synthesis provides novel insights into the magnitude, processes, and contexts of biocrust effects in drylands. This information is critical to improve our capacity to manage dwindling dryland water supplies as Earth becomes hotter and drier.

Global Change Biology↗
USGS2020

Plate boundary localization, slip-rates and rupture segmentation of the Queen Charlotte Fault based on submarine tectonic geomorphology

Linking fault behavior over many earthquake cycles to individual earthquake behavior is a primary goal in tectonic geomorphology, particularly across an entire plate boundary. Here, we examine the 1150-km-long, right-lateral Queen Charlotte-Fairweather fault system using comprehensive multibeam bathymetry data acquired along the Queen Charlotte Fault (QCF) offshore southeastern Alaska and western British Columbia. Fine-scale analysis of tectonic geomorphology allowed us to identify and reconstruct 184 strike-slip piercing points over a 630 km stretch of the QCF. Age constraints from glacial recession and offshore sedimentation patterns yield a consistent slip-rate of ∼50–57 mm/yr since ∼17–12 ka, the fastest rate for a continent-ocean strike-slip fault on Earth. These slip-rates equal or exceed estimates of Pacific-North America (PA-NA) relative motion from global plate reconstructions, indicating that PA-NA motion is highly localized. The QCF cuts the seafloor along a narrow and unusually straight trace for its entire length and multiple fault traces are observed only at local step-overs. The geometry and behavior of the QCF over many earthquake cycles is simple and typical of mature faults with relatively homogeneous stress fields. Since the QCF is the primary PA-NA plate boundary, we used the trace of the QCF to define the small circle path for relative plate motion and computed the associated Euler pole. Predicted along-strike obliquity variations based on the new pole agree with observed tectonic geomorphology and suggest that previous global plate reconstructions overestimated the degree of oblique convergence along the QCF. We also find that subtle, long-wavelength (75–150 km) bends and discrete step-overs appear to define the endpoints of M>7 earthquakes, suggesting that obliquity and resultant fault geometry may control rupture segmentation and asperity development. Lastly, the agreement between predicted obliquity and tectonic geomorphology along the entire length of QCF compelled a reevaluation of regional tectonic models. In the north, the eastern Yakatat Terrane appears to be translating northwest with the Pacific plate, and slip transferred from the QCF to the Fairweather Fault results in ∼20 mm/yr of convergence along the southern St. Elias mountains. In the south, we predict a reduced rate of convergence along the QCF west of Haida Gwaii (∼5–6 mm/yr of shortening, on average) relative to previous studies. Our results support a model for transpression and strike-slip partitioning along the edge of a hot and weak Pacific Plate, leading to crustal thickening and growth of the Queen Charlotte Terrace to the west of Haida Gwaii.

Alaska, British Columbia↗
USGS2006

Study of volcano/ice interactions gains momentum

Observations of recent volcanic eruptions in Iceland and detailed studies of sub-glacially erupted deposits and the interaction of lava and pyroclastic flows with snow and ice have provided important new data that should lead to significant advances in the understanding of volcano/ice interaction on Earth and Mars. A conference on this subject, the first of its kind, recently brought together geologists, geophysicists, glaciologists, and planetary scientists studying various aspects of volcano-ice interaction.

Eos, Transactions, American Geophysical Union↗
USGS1981

Did ice streams carve martian outflow channels?

Outflow channels on Mars 1 are long sinuous linear depressions that occur mostly in the equatorial area (±30° lat.). They differ from small valley networks 2 by being larger and arising full born from chaotic terrains. Outflow channels resemble terrestrial stream beds, and their origin has generally been attributed to water 3–5 in catastrophic floods 6,7 or mudflows 8 . The catastrophic-flood hypothesis is derived primarily from the morphological similarities of martian outflow channels and features created by the catastrophic Spokane flood that formed the Washington scablands. These similarities have been documented extensively 3,6,7 , but differences of scale remain a major problem: martian channel features are on the average much larger than their proposed terrestrial analogues. We examine here the problem of channel origin from the perspective of erosional characteristics and the resultant landf orms created by former and present-day ice streams and glaciers on Earth. From morphologic comparisons, an ice-stream origin seems equally well suited to explain the occurrences and form of the outflow channels on Mars, and in contrast with the hydraulic hypothesis, ice streams and ice sheets produce terrestrial features of the same scale as those observed on Mars.

Nature↗
USGS1971

Aeromagnetic study of the midcontinent gravity high of central United States

A composite map of detailed aeromagnetic surveys over the midcontinent gravity high provides coverage of the 600-mi-long buried belt of mafic rocks of the Keweenawan Series from their outcrop localities in Minnesota and Wisconsin through Iowa and Nebraska. A map of the subsurface extent of the mafic rocks, based on the intricate magnetic patterns, shows that the rocks form a long, semicontinuous block, averaging 40 mi wide and consisting mainly of a sequence of layered flows. This sequence is probably fault-bounded and has been tilted up along the margins, where the linearity of the anomalies indicates steeper dips. The associated clastic rocks, indicated by a smoother magnetic pattern, occur in basins along both sides of the mafic belt and in grabens and a series of axial basins on the upper surface of the block. The well-defined outliers of flows marginal to the main block and the truncation of some of the outermost flow units along a diagonal boundary striking at an angle to them suggest that the present boundaries of the block are postdepositional structural features. The basins and the edges of the block appear to have controlled later, largely vertical movement in the overlying Paleozoic and younger sedimentary cover. Calculated models based on coincident magnetic and detailed gravity profiles along typical cross sections of the midcontinent gravity high show that the block of mafic rocks is steep-sided and as much as several miles thick. The free-air gravity anomaly, which consists of a large positive maximum flanked by minima, averages very close to zero, indicating that this major crustal feature is regionally compensated, although locally each of its components shows a large departure from equilibrium. Remanent magnetization is a primary factor in the interpretation of the magnetic data. Magnetic property studies of Keweenawan mafic rocks in the Lake Superior region show that remanent magnetization may be five times the magnetization induced by the present Earth's field and differs from it radically in direction. This magnetization was acquired before the flows were tilted into their present positions. A computed magnetic profile shows that a trough of flows with such a magnetization and inward-dipping limbs can account for the observed persistent lows along the western edge of the block, the relatively low magnetic values along the axis of the block, and the large positive anomaly along the eastern side of the block. Flows as much as 1 mi thick near the base of the sequence have a remanent magnetization with a nearly opposite polarity. This reverse polarity has been measured on both sides of Lake Superior and is probably also present farther south, particularly in Iowa where the outer units of the block in an area north of Des Moines give rise to a prominent magnetic low. The axis of this long belt of Keweenawan mafic rocks cuts discordantly through the prevailing east-west-trending fabric of the older Precambrian terrane from southern Kansas to Lake Superior. This belt has several major left-lateral offsets, one of which produces a complete hiatus in the vicinity of the 40th parallel where an east-west transcontinental rift or fracture zone has been proposed. The axial basins of clastic rocks are outlined by linear magnetic anomalies and show a concordant relation to the structure of the mafic flows. These basins are oriented at an angle to the main axis, suggesting that the entire feature originated as a major rift composed of a series of short, linear, en echelon segments with offsets similar to the transform faults characterizing the present mid-ocean rift system. This midcontinent rift may well have been part of a Keweenawan global rift system with initial offsets consisting of transform faults along pre-existing fractures, but apparently it never fully developed laterally into an ocean basin, and the upwelling mafic material was localized along a relatively narrow belt.

Geological Society of America Bulletin↗
USGS2012

Methods of practice and guidelines for using survey-grade global navigation satellite systems (GNSS) to establish vertical datum in the United States Geological Survey

Geodetic surveys have evolved through the years to the use of survey-grade (centimeter level) global positioning to perpetuate and post-process vertical datum. The U.S. Geological Survey (USGS) uses Global Navigation Satellite Systems (GNSS) technology to monitor natural hazards, ensure geospatial control for climate and land use change, and gather data necessary for investigative studies related to water, the environment, energy, and ecosystems. Vertical datum is fundamental to a variety of these integrated earth sciences. Essentially GNSS surveys provide a three-dimensional position x, y, and z as a function of the North American Datum of 1983 ellipsoid and the most current hybrid geoid model. A GNSS survey may be approached with post-processed positioning for static observations related to a single point or network, or involve real-time corrections to provide positioning "on-the-fly." Field equipment required to facilitate GNSS surveys range from a single receiver, with a power source for static positioning, to an additional receiver or network communicated by radio or cellular for real-time positioning. A real-time approach in its most common form may be described as a roving receiver augmented by a single-base station receiver, known as a single-base real-time (RT) survey. More efficient real-time methods involving a Real-Time Network (RTN) permit the use of only one roving receiver that is augmented to a network of fixed receivers commonly known as Continually Operating Reference Stations (CORS). A post-processed approach in its most common form involves static data collection at a single point. Data are most commonly post-processed through a universally accepted utility maintained by the National Geodetic Survey (NGS), known as the Online Position User Service (OPUS). More complex post-processed methods involve static observations among a network of additional receivers collecting static data at known benchmarks. Both classifications provide users flexibility regarding efficiency and quality of data collection. Quality assurance of survey-grade global positioning is often overlooked or not understood and perceived uncertainties can be misleading. GNSS users can benefit from a blueprint of data collection standards used to ensure consistency among USGS mission areas. A classification of GNSS survey qualities provide the user with the ability to choose from the highest quality survey used to establish objective points with low uncertainties, identified as a Level I, to a GNSS survey for general topographic control without quality assurance, identified as a Level IV. A Level I survey is strictly limited to post-processed methods, whereas Level II, Level III, and Level IV surveys integrate variations of a RT approach. Among these classifications, techniques involving blunder checks and redundancy are important, and planning that involves the assessment of the overall satellite configuration, as well as terrestrial and space weather, are necessary to ensure an efficient and quality campaign. Although quality indicators and uncertainties are identified in post-processed methods using CORS, the accuracy of a GNSS survey is most effectively expressed as a comparison to a local benchmark that has a high degree of confidence. Real-time and post-processed methods should incorporate these "trusted" benchmarks as a check during any campaign. Global positioning surveys are expected to change rapidly in the future. The expansion of continuously operating reference stations, combined with newly available satellite signals, and enhancements to the conterminous geoid, are all sufficient indicators for substantial growth in real-time positioning and quality thereof.

Techniques and Methods↗
USGS2012

GFDL's ESM2 global coupled climate-carbon Earth System Models. Part I: physical formulation and baseline simulation characteristics

We describe the physical climate formulation and simulation characteristics of two new global coupled carbon-climate Earth System Models, ESM2M and ESM2G. These models demonstrate similar climate fidelity as the Geophysical Fluid Dynamics Laboratory's previous CM2.1 climate model while incorporating explicit and consistent carbon dynamics. The two models differ exclusively in the physical ocean component; ESM2M uses Modular Ocean Model version 4.1 with vertical pressure layers while ESM2G uses Generalized Ocean Layer Dynamics with a bulk mixed layer and interior isopycnal layers. Differences in the ocean mean state include the thermocline depth being relatively deep in ESM2M and relatively shallow in ESM2G compared to observations. The crucial role of ocean dynamics on climate variability is highlighted in the El Niño-Southern Oscillation being overly strong in ESM2M and overly weak ESM2G relative to observations. Thus, while ESM2G might better represent climate changes relating to: total heat content variability given its lack of long term drift, gyre circulation and ventilation in the North Pacific, tropical Atlantic and Indian Oceans, and depth structure in the overturning and abyssal flows, ESM2M might better represent climate changes relating to: surface circulation given its superior surface temperature, salinity and height patterns, tropical Pacific circulation and variability, and Southern Ocean dynamics. Our overall assessment is that neither model is fundamentally superior to the other, and that both models achieve sufficient fidelity to allow meaningful climate and earth system modeling applications. This affords us the ability to assess the role of ocean configuration on earth system interactions in the context of two state-of-the-art coupled carbon-climate models.

Journal of Climate↗
USGS2014

Inundation of a barrier island (Chandeleur Islands, Louisiana, USA) during a hurricane: Observed water-level gradients and modeled seaward sand transport

Large geomorphic changes to barrier islands may occur during inundation, when storm surge exceeds island elevation. Inundation occurs episodically and under energetic conditions that make quantitative observations difficult. We measured water levels on both sides of a barrier island in the northern Chandeleur Islands during inundation by Hurricane Isaac. Wind patterns caused the water levels to slope from the bay side to the ocean side for much of the storm. Modeled geomorphic changes during the storm were very sensitive to the cross-island slopes imposed by water-level boundary conditions. Simulations with equal or landward sloping water levels produced the characteristic barrier island storm response of overwash deposits or displaced berms with smoother final topography. Simulations using the observed seaward sloping water levels produced cross-barrier channels and deposits of sand on the ocean side, consistent with poststorm observations. This sensitivity indicates that accurate water-level boundary conditions must be applied on both sides of a barrier to correctly represent the geomorphic response to inundation events. More broadly, the consequence of seaward transport is that it alters the relationship between storm intensity and volume of landward transport. Sand transported to the ocean side may move downdrift, or aid poststorm recovery by moving onto the beach face or closing recent breaches, but it does not contribute to island transgression or appear as an overwash deposit in the back-barrier stratigraphic record. The high vulnerability of the Chandeleur Islands allowed us to observe processes that are infrequent but may be important at other barrier islands.

Louisiana↗
USGS2017

Delayed seismicity rate changes controlled by static stress transfer

On 15 June 2010, a M w 5.7 earthquake occurred near Ocotillo, California, in the Yuha Desert. This event was the largest aftershock of the 4 April 2010 M w 7.2 El Mayor-Cucapah (EMC) earthquake in this region. The EMC mainshock and subsequent Ocotillo aftershock provide an opportunity to test the Coulomb failure hypothesis (CFS). We explore the spatiotemporal correlation between seismicity rate changes and regions of positive and negative CFS change imparted by the Ocotillo event. Based on simple CFS calculations we divide the Yuha Desert into three subregions, one triggering zone and two stress shadow zones. We find the nominal triggering zone displays immediate triggering, one stress shadowed region experiences immediate quiescence, and the other nominal stress shadow undergoes an immediate rate increase followed by a delayed shutdown. We quantitatively model the spatiotemporal variation of earthquake rates by combining calculations of CFS change with the rate-state earthquake rate formulation of Dieterich (1994), assuming that each subregion contains a mixture of nucleation sources that experienced a CFS change of differing signs. Our modeling reproduces the observations, including the observed delay in the stress shadow effect in the third region following the Ocotillo aftershock. The delayed shadow effect occurs because of intrinsic differences in the amplitude of the rate response to positive and negative stress changes and the time constants for return to background rates for the two populations. We find that rate-state models of time-dependent earthquake rates are in good agreement with the observed rates and thus explain the complex spatiotemporal patterns of seismicity.

California↗
USGS1981

Sudden death at the end of the Mesozoic

A paleoecological analysis of the fossil record before and after the Cretaceous/Tertiary boundary indicates that the widespread extinctions and biological stresses around the boundary are best explained in terms of a sudden, significant, but short temperature rise. L. Alvarez and co-authors, having found an enrichment in iridium at the same boundary, postulated that it was associated with the impact of an extraterrestrial body. If this body struck the ocean, the water injected into the atmosphere may have led to a transient increase in the global surface temperature. This temperature pulse may have been primarily responsible for the effects observed in the biosphere. The pattern of extinction of higher plant species suggests that splash down occurred in the northern Pacific-Bering Sea area.

Earth and Planetary Science Letters↗
USGS2009

Methane sources and production in the northern Cascadia margin gas hydrate system

The oceanographic and tectonic conditions of accretionary margins are well-suited for several potential processes governing methane generation, storage and release. To identify the relevant methane evolution pathways in the northern Cascadia accretionary margin, a four-site transect was drilled during Integrated Ocean Drilling Program Expedition 311. The δ13C values of methane range from a minimum value of − 82.2‰ on an uplifted ridge of accreted sediment near the deformation front (Site U1326, 1829 mbsl, meters below sea level) to a maximum value of − 39.5‰ at the most landward location within an area of steep canyons near the shelf edge (Site U1329, 946 mbsl). An interpretation based solely on methane isotope values might conclude the 13C-enrichment of methane indicates a transition from microbially- to thermogenically-sourced methane. However, the co-existing CO2 exhibits a similar trend of 13C-enrichment along the transect with values ranging from − 22.5‰ to +25.7‰. The magnitude of the carbon isotope separation between methane and CO2 (εc = 63.8 ± 5.8) is consistent with isotope fractionation during microbially mediated carbonate reduction. These results, in conjunction with a transect-wide gaseous hydrocarbon content composed of > 99.8% (by volume) methane and uniform δDCH4 values (− 172‰ ± 8) that are distinct from thermogenic methane at a seep located 60 km from the Expedition 311 transect, suggest microbial CO2 reduction is the predominant methane source at all investigated sites. The magnitude of the intra-site downhole 13C-enrichment of CO2 within the accreted ridge (Site U1326) and a slope basin nearest the deformation front (Site U1325, 2195 mbsl) is ~ 5‰. At the mid-slope site (Site U1327, 1304 mbsl) the downhole 13C-enrichment of the CO2 is ~ 25‰ and increases to ~ 40‰ at the near-shelf edge Site U1329. This isotope fractionation pattern is indicative of more extensive diagenetic alteration at sites with greater 13C-enrichment. The magnitude of the 13C-enrichment of CO2 correlates with decreasing sedimentation rates and a diminishing occurrence of stratigraphic gas hydrate. We suggest the decreasing sedimentation rates increase the exposure time of sedimentary organic matter to aerobic and anaerobic degradation, during burial, thereby reducing the availability of metabolizable organic matter available for methane production. This process is reflected in the occurrence and distribution of gas hydrate within the northern Cascadia margin accretionary prism. Our observations are relevant for evaluating methane production and the occurrence of stratigraphic gas hydrate within other convergent margins.

Earth and Planetary Science Letters↗
USGS2001

Earthquake stress drop and laboratory-inferred interseismic strength recovery

We determine the scaling relationships between earthquake stress drop and recurrence interval t r that are implied by laboratory-measured fault strength. We assume that repeating earthquakes can be simulated by stick-slip sliding using a spring and slider block model. Simulations with static/kinetic strength, time-dependent strength, and rate- and state-variable-dependent strength indicate that the relationship between loading velocity and recurrence interval can be adequately described by the power law V L ∝ t r n where n≈−1. Deviations from n=−1 arise from second order effects on strength, with n>−1 corresponding to apparent time-dependent strengthening and n<−1 corresponding to weakening. Simulations with rate and state-variable equations show that dynamic shear stress drop Δτ d scales with recurrence as d Δτ d / d ln t r ≤σ e ( b-a ), where σ e is the effective normal stress, μ=τ/σ e , and ( a-b )= d μ ss / d ln V is the steady-state slip rate dependence of strength. In addition, accounting for seismic energy radiation, we suggest that the static shear stress drop Δτ s scales as d Δτ s / d ln t r ≤σ e (1 +ζ)( b-a ), where ζ is the fractional overshoot. The variation of Δτ s with ln t r for earthquake stress drops is somewhat larger than implied by room temperature laboratory values of ζ and b-a . However, the uncertainty associated with the seismic data is large and the discrepancy between the seismic observations and the rate of strengthening predicted by room temperature experiments is less than an order of magnitude.

Journal of Geophysical Research B: Solid Earth↗
USGS1997

Stressing of the New Madrid seismic zone by a lower crust detachment fault

A new mechanical model for the cause of the New Madrid seismic zone in the central United States is analyzed. The model contains a subhorizontal detachment fault which is assumed to be near the domed top surface of locally thickened anomalous lower crust (“rift pillow”). Regional horizontal compression induces slip on the fault, and the slip creates a stress concentration in the upper crust above the rift pillow dome. In the coseismic stage of the model earthquake cycle, where the three largest magnitude 7–8 earthquakes in 1811–1812 are represented by a single model mainshock on a vertical northeast trending fault, the model mainshock has a moment equivalent to a magnitude 8 event. During the interseismic stage, corresponding to the present time, slip on the detachment fault exerts a right-lateral shear stress on the locked vertical fault whose failure produces the model mainshock. The sense of shear is generally consistent with the overall sense of slip of 1811–1812 and later earthquakes. Predicted rates of horizontal strain at the ground surface are about 10 −7 year −1 and are comparable to some observed rates. The model implies that rift pillow geometry is a significant influence on the maximum possible earthquake magnitude.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Drainage effects on the transient, near-surface hydrologic response of a steep hillslope to rainfall: Implications for slope stability, Edmonds, Washington, USA

Shallow landslides on steep (>25??) hillsides along Puget Sound have resulted in occasional loss of life and costly damage to property during intense or prolonged rainfall. As part of a larger project to assess landslide hazards in the Seattle area, the U.S. Geological Survey instrumented two coastal bluff sites in 2001 to observe the subsurface hydrologic response to rainfall. The instrumentation at one of these sites, near Edmonds, Washington, consists of two rain gauges, two water-content probes that measure volumetric water content at eight depths between 0.2 and 2.0 m, and two tensiometer nests that measure soil-water suction at six depths ranging from 0.2 to 1.5m. Measurements from these instruments are used to test one- and two-dimensional numerical models of infiltration and groundwater flow. Capillary-rise tests, performed in the laboratory on soil sample from the Edmonds site, are used to define the soil hydraulic properties for the wetting process. The field observations of water content and suction show an apparent effect of porosity variation with depth on the hydraulic response to rainfall. Using a range of physical properties consistent with our laboratory and field measurements, we perform sensitivity analyses to investigate the effects of variation in physical and hydraulic properties of the soil on rainfall infiltration, pore-pressure response, and, hence, slope stability. For a two-layer-system in which the hydraulic conductivity of the upper layer is at least 10 times greater than the conductivity of the lower layer, and the infiltration rate is greater than the conductivity of the lower layer, a perched water table forms above the layer boundary potentially destabilizing the upper layer of soil. Two-dimensional modeling results indicate that the addition of a simple trench drain to the same two-layer slope has differing effects on the hydraulic response depending on the initial pressure head conditions. For slope-parallel flow conditions, pressure head is significantly reduced near the drain; however, for transient, vertical infiltration in a partially saturated soil, conditions consistent with those observed during monitoring at the Edmonds site, the drain decreases the thickness of a perched water table by a small amount.

Natural Hazards and Earth System Sciences↗
USGS2003

Groundwater flow, heat transport, and water table position within volcanic edifices: Implications for volcanic processes in the Cascade Range

The position of the water table within a volcanic edifice has significant implications for volcano hazards, geothermal energy, and epithermal mineralization. We have modified the HYDROTHERM numerical simulator to allow for a free-surface (water table) upper boundary condition and a wide range of recharge rates, heat input rates, and thermodynamic conditions representative of continental volcano-hydrothermal systems. An extensive set of simulations was performed on a hypothetical stratovolcano system with unconfined groundwater flow. Simulation results suggest that the permeability structure of the volcanic edifice and underlying material is the dominant control on water table elevation and the distribution of pressures, temperatures, and fluid phases at depth. When permeabilities are isotropic, water table elevation decreases with increasing heat flux and increases with increasing recharge, but when permeabilities are anisotropic, these effects can be much less pronounced. Several conditions facilitate the ascent of a hydrothermal plume into a volcanic edifice: a sufficient source of heat and magmatic volatiles at depth, strong buoyancy forces, and a relatively weak topography-driven flow system. Further, the plume must be connected to a deep heat source through a pathway with a time-averaged effective permeability ???1 ?? 10-16 m2, which may be maintained by frequent seismicity. Topography-driven flow may be retarded by low permeability in the edifice and/or the lack of precipitation recharge; in the latter case, the water table may be relatively deep. Simulation results were compared with observations from the Quaternary stratovolcanoes along the Cascade Range of the western United States to infer hydrothermal processes within the edifices. Extensive ice caps on many Cascade Range stratovolcanoes may restrict recharge on the summits and uppermost flanks. Both the simulation results and limited observational data allow for the possibility that the water table beneath the stratovolcanoes is relatively deep.

Journal of Geophysical Research B: Solid Earth↗
USGS1993

Origin and significance of tourmaline-rich rocks in the Broken Hill district, Australia

Tourmaline-rich rocks are widespread minor lithologies within the Early Proterozoic Willyama Supergroup in the Broken Hill district, Australia. Tourmaline concentrations occur in strata-bound and local stratiform tourmalinites, clastic metasedimentary rocks, quartz-gahnite lode rocks, stratiform Pb-Zn-Ag sulfide ores, garnet quartzites, strata-bound scheelite deposits, quartz-tourmaline nodules, discordant quartz veins, and granitic pegmatites. Most of the tourmaline-rich rocks are within the Broken Hill Group that hosts the main Pb-Zn-Ag ores.At the Globe mine along the northeast end of the main lodes, tourmalinites are closely associated with Pb-Zn-Ag mineralization and in places are interbedded with Mn-rich garnet quartzites. Galena and other ore minerals occur locally in the cores of recrystallized tourmaline grains in these tourmalinites, indicating that tourmaline and sulfides were present together prior to deformation and high-grade metamorphism. Electron microprobe analyses of tourmalines intergrown with Fe sulfides at the Globe mine show Mg-rich compositions relative to tourmalines in sulfide-free assemblages from the same area, suggesting early (premetamorphic) introduction of boron and Mg enrichment of tourmaline by sulfide-silicate reactions during metamorphism.Combined field and geochemical data indicate that the district tourmalinites represent normal clastic sediments that were metasomatically altered by boron-rich hydrothermal fluids at or below the sediment-water interface. Whole-rock chemical analyses of 33 tourmaline-rich rocks show linear trends of data for major and trace elements that closely resemble the trends observed for unmineralized elastic metasedimentary rocks of the district. Average Fe/Al, Mg/Al, Na/Al, and Ti/Al molar ratios of the tourmaline-rich rocks and clastic metasediments are very similar; the average K/Al molar ratio of the tourmaline-rich rocks is significantly lower than that of the clastic metasediments, reflecting the loss of K during tourmalinite formation. Chondrite-normalized patterns of rare earth elements (REE) in the quartz-rich tourmalinites are generally similar to those of the clastic metasediments, except for minor depletions of light REE; local positive and negative Ce anomalies suggest tourmalinite formation in the presence of seawater or a seawater-derived pore fluid. The geochemical data imply relative immobility of Al, Ti, Cr, and heavy REE during hydrothermal alteration and later metamorphism. Boron isotope analyses of 52 tourmaline separates show a total range of delta 11 B values from -26.8 to -17.0 per mil. Fine-grained, euhedral, nonpoikilitic tourmalines from tourmalinites in the andalusite-muscovite zone in the northern part of the district (e.g., Black Prince mine) have delta 11 B values from -21 to -17 per mil, whereas coarse granoblastic and poikilo-blastic tourmalines from the sillimanite and two-pyroxene granulite zones in the southern part of the district (e.g., Globe mine) have delta 11 B values of-24 to -19 per mil. Tourmalines in strongly retrogressed tourmalinites have delta 11 B values from about -27 to -20 per mil. The observed variations in delta 11 B are consistent with prograde and retrograde metamorphic fractionation of boron isotopes, in which the fluid phase is preferentially enriched in the heavier isotope ( 11 B). Premetamorphic hydrothermal fluids that deposited the Black Prince tourmalinites had delta 11 B values of-8 to -5 per mil at 200 degrees to 300 degrees C, suggesting a boron source from nonmarine evaporite borates.Tourmalinites in the Broken Hill district apparently formed by the same submarine hydrothermal processes as the main Pb-Zn-Ag lodes and the siliceous ferromanganese protoliths of the garnet quartzites. In our model, the hydrothermal system(s) acquired abundant boron by leaching evaporitic borates within the Thackaringa Group, the stratigraphic sequence that underlies the Broken Hill Group and most of the tourmaline concentrations. We suggest that evaporites of the Thackaringa Group provided a source of readily extractable boron for formation of the tourmalinites and also the source of the fluoride, sulfur, and perhaps the carbonate in the main lodes; such evaporites may have been critical for increased metal chloride complexing and transport necessary for deposition of the high-grade Pb-Zn-Ag ores. The Broken Hill deposit may have formed contemporaneously with the Mount Isa and McArthur River Pb-Zn-Ag deposits in similar evaporite-bearing sequences during widespread Early Proterozoic continental rifting.

Economic Geology↗