Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Interfaces”

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 901 records · Page 50Linked to original sources
USGS2024

Urbanization and water management control stream water quality along a mountain to plains transition

Urbanization can have substantial effects on water quality due to altered hydrology and introduction of constituents to water bodies. In arid and semi-arid environments, streams are further stressed by dewatering as a result of diversions. We conducted a high-resolution synoptic survey of two streams in Colorado, USA that transition abruptly from granitic/metamorphic forested mountains to sedimentary urbanized plains and are both highly managed for water supply, yet differ in degree of urbanization. A synoptic mass balance approach developed for mine drainage applications was adapted to elucidate effects of urbanization, geology, and diversions on stream chemistry during baseflow conditions. Urbanization was a more important driver of stream concentrations than geology. The urban area was a strong source of bromide, calcium, chloride, and manganese, while lanthanum and dissolved organic carbon were primarily sourced from the mountains. A majority of streamflow was removed by diversions near the mountains/plains interface. Groundwater accounted for 31% of the subsequent flow increase to the urbanized stream, and delivered at least 33% of chloride loading. Constituents that were primarily urban-derived (bromide, calcium, chloride, and manganese) were 2–3 times higher in the urban region due to diversions; without diversions, stream water quality would have largely retained characteristics of forested streams through the urban reach. This study provides insights into processes that affect water quality in highly managed streams of the semi-arid western USA.

Colorado↗
USGS2024

Transdisciplinary research supports the sustainability of barrier island systems threatened by climate change

The management of developed barrier islands is often piece-meal and reactionary despite the complex, dynamic nature of these systems, and sustainable practices will become increasingly difficult due to heightened pressures of climate change. Adaptation actions, including nature-based solutions, need to be thoroughly evaluated prior to implementation to understand system-wide impacts and avoid maladaptation. Anarde et al. (2024a), ( https://doi.org/10.1029/2023ef003672 ), Anarde et al. (2024b), ( https://doi.org/10.1029/2023ef004200 ) is the latest important contribution in a growing body of transdisciplinary research that more robustly evaluates the complex physical process-and-response relationship of barrier systems via sophisticated numerical modeling approaches that also interface with socioeconomic models to inform coastal management actions in response to mitigating coastal risk. This new research indicates the importance of coordinated system-scale barrier island management, as strategies to reduce coastal hazard risk in one location will directly affect adjacent communities. Further, this work demonstrates that reducing barrier management interventions may actually promote barrier recovery and sustainability in the face of sea level rise. In addition, recent advances in the analysis and application of remotely sensed data from satellites and oblique aerial photography provide scientists an unprecedented opportunity to track coastal evolution over a wide range of spatial and temporal scales at minimal cost. As sea level rise and changing storm patterns challenge the sustainable management of barrier island systems, integrating these advanced, transdisciplinary tools will enable scientists and coastal practitioners to more thoroughly evaluate coastal adaptation options, efficiently invest limited resources to mitigate coastal hazard risk for communities, support healthy ecosystems, and reduce system-wide impacts.

Earth's Future↗
USGS2025

Physicochemical properties and bioreactivity of sub-10 µm geogenic particles: Comparison of volcanic ash and desert dust

Exposure to ambient particulate matter (PM) with an aerodynamic diameter of <10 μm (PM 10 ) is a well-established health hazard. There is increasing evidence that geogenic (Earth-derived) particles can induce adverse biological effects upon inhalation, though there is high variability in particle bioreactivity that is associated with particle source and physicochemical properties. In this study, we investigated physicochemical properties and biological reactivity of volcanic ash from the April 2021 eruption of La Soufrière volcano, St. Vincent, and two desert dust samples: a standardized test dust from Arizona and an aeolian Gobi Desert dust sampled in China. We determined particle size, morphology, mineralogy, surface texture and chemistry in sub-10 μm material to investigate associations between particle physicochemical properties and observed bioreactivity. We assessed cellular responses (cytotoxic and pro-inflammatory effects) to acute particle exposures (24 hr) in monocultures at the air-liquid interface using two types of cells of the human airways: BEAS-2B bronchial epithelial cells and A549 alveolar type II epithelial cells. In acellular assays, we also assessed particle oxidative potential and the presence of microorganisms. The results showed that volcanic ash and desert dust exhibit intrinsically different particle morphology, surface textures and chemistry, and variable mineralogical content. We found that Gobi Desert dust is more bioreactive than freshly erupted volcanic ash and Arizona test dust, which is possibly linked to the presence of microorganisms (bacteria) and/or nanoscale elongated silicate minerals (potentially clay such as illite or vermiculite) on particle surfaces.

GeoHealth↗
USGS2025

Dynamic feedbacks between river meandering and landsliding in northwestern Washington glacial terraces

Landsliding in river valleys poses unique risks for cascading hazards and can damage infrastructure and cause fatalities. In postglacial valleys, many landslides are posited to occur in relation to lateral river erosion, but the dynamics of fluvial-hillslope interactions are not well understood. Here, we investigate a section of the Nooksack River in western Washington State where the channel is flanked by landslide-prone glacial terraces similar to those that failed in the 2014 State Route 530 “Oso” landslide. We map 216 landslides through time across 17 aerial imagery data sets (1933–2022) and analyze them in relation to river meandering and curvature. We observe dynamic feedbacks between lateral river meandering and valley-adjacent landsliding. Terrace lateral retreat rates of up to 25 m/year owing to combined fluvial erosion and slope failure occur on pinned, outer meander bends immediately downstream from peaks in river curvature (>0.0075 1/m); these locations are predisposed to both shallow and deep-seated landslides. Deep-seated landslides extending 17%–32% of the active valley width into the floodplain can displace the river away from the floodplain margin and change the channel planform. River-displacing landslides relocate meanders up- or downstream, thereby conditioning the location of subsequent landslides. This conceptual model of coupled landslide-driven meander displacement and valley-adjacent landsliding is exemplified across western Washington river systems. The distance between up- and downstream valley-adjacent landsliding scales with valley width, meander wavelength, and terrace height. Our results can advance our understanding of the river-hillslope interface in landscape evolution and can be used to inform hazard management in river corridors.

Washington↗
USGS2024

Integrating depth measurements from gaging stations with image archives for spectrally based remote sensing of river bathymetry

Remote sensing can be an effective tool for mapping river bathymetry, but the need for direct measurements to calibrate image-derived depth estimates impedes broader application of this approach. One way to circumvent the need for field campaigns dedicated to calibration is to capitalize upon existing data. In this study, we introduce a framework for Bathymetric Mapping using Gage Records and Image Databases (BaMGRID). This workflow involves retrieving depth measurements made during gaging station site visits, downloading archived multispectral images, and then combining these two data sets to establish a relationship between depth and reflectance. We developed a processing chain that involves using application programming interfaces to obtain both depth measurements made during site visits and images centered on the gage and then linking depth to reflectance via an optimal band ratio analysis (OBRA) algorithm modified for small sample sizes. Applying this workflow to selected gages within two river basins indicated that depth retrieval from multispectral satellite images could be highly accurate, but with variable results from one image to the next at a given site. High resolution aerial photography was less conducive to bathymetric mapping in one of the basin considered. Of the four predictors of depth retrieval performance we evaluated (mean and standard deviation of depth, width, and an index of water clarity), only width was consistently significantly correlated with OBRA R 2 ( p < 0.026). Currently, BaMGRID is best-suited for site-by-site analysis to support practical applications at the reach scale; continuous, basin-wide mapping of river bathymetry will require additional research.

Willamette River basin, Delaware River basin↗
USGS2026

Hydrological control on soil redox condition and carbon loss of coastal wetland under sea-kevel rise

Coastal wetlands are critical carbon sinks with their biogeochemical and ecological functioning shaped by dynamic hydrological conditions that are increasingly influenced by climate change. A key unresolved question is how hydrologic flow, vegetation response, and rising sea levels interact to regulate soil redox condition and carbon loss in coastal wetlands. Using a field-tested hydrological–biogeochemical–ecological modeling framework, we reveal how the interplay between terrestrial groundwater discharge and tidal fluctuations generates complex groundwater flow patterns at the terrestrial–aquatic interface, and how these patterns modulate soil redox conditions, in turn influencing soil organic matter decomposition and carbon loss. Notably, rising sea levels suppress soil CO 2 emissions while reducing lateral dissolved carbon losses, thereby enhancing litter carbon sequestration under anoxic conditions. As vegetation responds to sea-level rise and carbon inputs diminish, litter carbon subsequently declines. These findings underscore a critical hydrological control on carbon cycling, advancing our understanding of coastal ecosystem resilience in a warming world.

Louisiana↗
USGS2026

From collision to subduction: Thermal-kinematic inversions constrain plate boundary structure and dip-slip activity in southwestern New Zealand

Constraints on three-dimensional fault system geometry and kinematics are required to understand both long-term tectonics and contemporary seismic hazard. In many cases, these parameters must be inferred from a combination of different data or analyses. Here, we invert low-temperature thermochronological data to quantitatively parameterize both the crustal-scale geometry and late Cenozoic dip-slip activity of the plate boundary system in southwest New Zealand. We find that the modern transpressive plate boundary system has accommodated the bulk of Australian-Pacific Plate convergence since establishment in the middle Miocene-early Pliocene. Shortening has occurred via (a) dextral-reverse slip on the moderate angle, listric central Alpine Fault décollement, and (b) reverse-slip on either the sharply curved, anti-listric northern Puysegur subduction interface or low-angle offshore thrust fault system that splays from the high-angle southern Alpine Fault at depth. We resolve late Cenozoic shortening rates on these structures that are indistinguishable, and interpret these to reflect the transfer of convergent deformation outboard of the predominantly strike-slip southern Alpine Fault. Our work helps elucidate the transition from strain localized oblique continental collision to strain partitioned subduction. We also provide an estimate of the across-strike width of the central-southern Alpine Fault section boundary at seismogenic depths. This is the first quantitative constraint on the geometry of a structural feature inferred to have arrested several major prehistoric earthquakes. This study furthers understanding of fault structure and dip-slip kinematics along the greatest sources of seismic hazard in and offshore of southwest New Zealand.

South Island↗
USGS2026

Constraining the permeability and outer-rise hydration at the Central America margin

Subduction zones are important drivers of Earth's geochemical evolution; the presence of water within the subduction system has substantial implications for geohazards, including the style of seismicity and the volume of magma generated at arc volcanoes. We utilize estimates of pore fluid distributions in the outer-rise region of the Central American margin imaged using controlled-source electro-magnetic surveys and novel two-dimensional forward models simulating poro-elastic deformation of slab bending to constrain the background permeability structure of the Cocos Plate and within outer-rise fault damage zones. The model results suggest that a relatively low background permeability structure is required to match the observed porosity at depths greater than 3 km (less than ~10 −21 m 2 . Similarly, increasing the permeability within the damage zones of outer-rise faults by more than a factor of 10 also over-predicts observed porosity distributions at depths greater than 3 km. However, the simulations consistently underpredict pore fluid distributions between depths of 1–3 km, which may reflect nonlinear elastic or anelastic deformation processes not included in the simulations or inherent uncertainty in the conversion of resistivity measurements to porosity distributions. The best-fitting depth-permeability relationship is relatively low, which has important implications for the timing of pore fluid release down-dip of the trench, as it restricts shallow fluid release, likely leading to higher pore fluid pressures, more widespread hydraulic fracturing, and deeper release of pore fluid along the plate interface.

Journal of Geophysical Research, Solid Earth↗
USGS2026

Corundum discovered by SuperCam and the Perseverance rover at Jezero crater, Mars

Mars is primarily composed of mafic mineral assemblages and their alteration products, but small, scattered rocks strewn across the landscape offer clues to greater petrological diversity. While traversing the Jezero crater rim, the Perseverance rover encountered several plagioclase-rich light-toned float rocks. SuperCam identified the distinctive signature of corundum (α-Al 2 O 3 ) in these rocks using time-resolved luminescence spectroscopy. Two strong peaks (692.7 and 694.1 nm) with millisecond lifetimes, and additional supporting lines, are consistent with Cr 3+ substitution for Al 3+ in corundum. Corundum forms in Al-rich, Si-depleted environments through magmatic or metamorphic processes. Given the rocks' small size, association with plagioclase, and location on the crater rim, we interpret the most plausible formation scenario to be impact induced metamorphism at the interface of a felsic and a mafic/ultramafic member with the likely action of fluids at some stage, although other possibilities are not excluded.

Geophysical Research Letters↗
USGS2012

Tide gage measurements of uplift along the south coast of Alaska

Annual mean sea levels along the south coast of Alaska are used to measure uplift along the Alaska‐Aleutian subduction zone. Oceanographic effects are removed from the observed annual mean sea levels by subtracting a correction that is proportional to the sea level fluctuations observed in southeast Alaska. That correction is effective in reducing fluctuations in the observed, annual mean sea level as far west as the tip of Alaska peninsula. Additional corrections to remove the eustatic rise in sea level and the apparent fall in sea level due to postglacial isostatic rebound of the land are introduced. This corrected sea level record should provide a measure of tectonic subsidence. In the area affected by the 1964 Alaska earthquake, postseismic uplift occurs where coseismic subsidence was observed, and postseismic subsidence occurs where coseismic uplift was observed. The immediate postseismic response is damped out within the first decade, and the subsequent uplift rates appear to be steady over the 1974–1989 interval. However, some of those rates seem to be too high to be sustained over the ∼1000 year earthquake recurrence interval appropriate to this area if the interseismic deformation is only to recover the coseismic displacement. Thus a long‐term ( ∼100 years) relaxation in uplift rates is postulated. The immediate (time constant ∼5 years) postseismic relaxation is attributed to postseismic slip on the plate interface directly downdip from the coseismic rupture. The long‐term (time constant ∼100 years) relaxation is attributed to flow in the asthenosphere.

Alaska↗
USGS2012

Basin and range crustal and upper mantle structure, northwest to central Nevada

We present an interpretation of the crustal and uppermost mantle structure of the Basin and Range of northwestern Nevada based on seismic refraction/wide‐angle reflection, near‐vertical reflection, and gravity data. In comparison to most previous estimates, we find that the crust is somewhat thicker (32–36 km versus 22–30 km), and the uppermost mantle velocity is somewhat higher (8.0 km/s versus 7.3–7.9 km/s). Along our transects, the crust is thinnest (32 km) in the Carson Sink‐Buena Vista Valley region and increases by 2–4 km to the west and east, respectively. There is considerable complexity throughout the crust where velocities range from of 2.5 km/s at the surface to 7.4 km/s in the lowermost crust. Variations in velocity and structure of the upper crustal layers reveal apparent basement velocity depressions (areas of lower velocities extending up to 10 km in depth) that underlie some surface ranges as well as the basins. The middle crust rises from about 20 km beneath central Nevada to within 12 km of the surface beneath the area of thinnest crust and is characterized by a modest (∼0.1 km/s) change in velocity and low‐velocity gradients. These midcrustal layers mark the onset of high crustal reflectivity and the apparent limiting depth to which Basin and Range faults can be traced in near‐vertical reflection profiles, suggesting that these midcrustal layers represent the transition between the brittle and ductile zones of the crust. The lower crust is more structurally complex, with layers thickening and thinning in a systematic manner with the upper crustal layers; generally, where there are velocity depressions in the upper crust, the lower crust is thickest and shallowest. The geometry of these lower crustal layers (derived from refraction modeling) coincides with changes in the crustal reflectivity, determined from the Consortium of Continental Reflection Profiling reflection data. The lower crustal layer is unusually high in velocity (7.4 km/s) and is likely the layer identified as mantle in some previous studies. We do not identify the 7.4 km/s layer as mantle because (1) there is an underlying layer with a velocity (8.0 km/s) that is more consistent with the worldwide average velocity for the upper mantle, and (2) the 7.4 km/s layer does not correspond to the “reflection” Moho. Gravity modeling and comparison to existing seismic models show a general consensus in many aspects with respect to crustal structure. This new model forms the basis for speculation on some of the processes associated with rifting of the Basin and Range Province. One such process, lithospheric magmatism, is inferred from the strong attenuation of transmitted seismic waves, which occurs at the same interface at which high‐amplitude, bright spot reflections originate. Unlike previous models, the overall structure and velocity of the crust and uppermost mantle of our new model are similar to other regions worldwide which have undergone high degrees of extension.

Nevada↗
USGS2010

Paleomagnetism and rotation constraints for the middle Miocene southwestern Nevada volcanic field

Middle Miocene rocks of the southwestern Nevada volcanic field (SWNVF) lie across the projection of the Walker Lane belt within the Basin and Range province and thus provide an interesting opportunity to test for late Cenozoic vertical-axis rotation. Paleomagnetic data from individual ash flow sheets document no significant relative vertical-axis rotation among localities within central SWNVF, an area of relatively low stratal tilts and widely spaced faults. A time-averaged mean paleomagnetic direction (D = 351.4°, I = 52.7°, α 95 = 4.5°) calculated from data from numerous separate rock units suggests that the central SWNVF underwent minimal counterclockwise vertical-axis rotation (R = −7.1° ± 6.6°) with respect to the North American craton. No clockwise vertical-axis rotation is found to support projection of dextral faults of the Walker Lane beneath the central SWNVF. Clockwise rotation of variable magnitude is common at numerous sites from southern and western margins of the field. These clockwise rotations probably reflect dextral shear strain developed at the interface between the little extended central SWNVF block and more strongly extended areas to the south and southwest of the field. Negligible rotation of 11.45-Ma to 13.25-Ma tuffs relative to the central SWNVF was found at the southeast margin of the field where 90° clockwise rotation at the northwest termination of the Las Vegas Valley shear zone had been postulated. Any clockwise rotation in this area must predate 13.25 Ma, and thus dextral shear within this part of the Walker Lane belt was not synchronous or connected across the southern margin of the field. Small counterclockwise vertical-axis rotation relative to the craton, as found for the central SWNVF block, might be a regional feature in the western Great Basin.

Tectonics↗
USGS1996

Cyclic stressing and seismicity at strongly coupled subduction zones

We use the finite element method to analyze stress variations in and near a strongly coupled subduction zone during an earthquake cycle. Deformation is assumed to be uniform along strike (plane strain on a cross section normal to the trench axis), and periodic earthquake slip is imposed consistent with the long-term rate of plate convergence and degree of coupling. Simulations of stress and displacement rate fields represent periodic fluctuations in time superimposed on an average field. The oceanic plate, descending slab, and continental lithosphere are assumed here to respond elastically to these fluctuations, and the remaining mantle under and between plates is assumed to respond as Maxwell viscoelastic. In the first part of the analysis we find that computed stress fluctuations in space and time are generally consistent with observed earthquake mechanism variations with time since a great thrust event. In particular, trench-normal extensional earthquakes tend to occur early in the earthquake cycle toward the outer rise but occur more abundantly late in the cycle in the subducting slab downdip of the main thrust zone. Compressional earthquakes, when they occur at all, have the opposite pattern. Our results suggest also that the actual timing of extensional outer rise events is controlled by the rheology of the shallow aseismic portion of the thrust interface. The second part of the analysis shows the effects of mantle relaxation on the rate of ground surface deformation during the earthquake cycle. Models without relaxation predict a strong overall compressional strain rate in the continental plate above the main thrust zone, with the strain rate constant between mainshocks. However with significant relaxation present, a localized region of unusually low compressional, or even slightly extensional, strain rate develops along the surface of the continental plate above and somewhat inland from the downdip edge of the locked main thrust zone. The low strain rate starts in the middle or late part of the cycle, depending on position. This result suggests that the negligible or small contraction measured on the Shumagin Islands, Alaska, during 1980 to 1991, may not invalidate an interpretation of that region as being a moderately coupled subduction zone. In contrast, mantle relaxation causes only modest temporal nonuniformity of uplift rates in the overriding plate and of extensional stress rates in the subducting plate, even when the Maxwell time is an order of magnitude less than the recurrence interval.

Journal of Geophysical Research B: Solid Earth↗
USGS2010

State-dependent anisotrophy: Comparison of quasi-analytical solutions with stochastic results for steady gravity drainage

Anisotropy in large-scale unsaturated hydraulic conductivity of layered soils changes with the moisture state. Here, state-dependent anisotropy is computed under conditions of large-scale gravity drainage. Soils represented by Gardner's exponential function are perfectly stratified, periodic, and inclined. Analytical integration of Darcy’s law across each layer results in a system of nonlinear equations that is solved iteratively for capillary suction at layer interfaces and for the Darcy flux normal to layering. Computed fluxes and suction profiles are used to determine both upscaled hydraulic conductivity in the principal directions and the corresponding “state-dependent” anisotropy ratio as functions of the mean suction. Three groups of layered soils are analyzed and compared with independent predictions from the stochastic results of Yeh et al. (1985b). The small-perturbation approach predicts appropriate behaviors for anisotropy under nonarid conditions. However, the stochastic results are limited to moderate values of mean suction; this limitation is linked to a Taylor series approximation in terms of a group of statistical and geometric parameters. Two alternative forms of the Taylor series provide upper and lower bounds for the state-dependent anisotropy of relatively dry soils.

Water Resources Research↗
USGS2010

Hydrologic and geochemical effects on oxygen uptake in bottom sediments of an effluent-dominated river

More than 95% of the water in the South Platte River downstream from the largest wastewater treatment plant serving the metropolitan Denver, Colorado, area consists of treated effluent during some periods of low flow. Fluctuations in effluent-discharge rates caused daily changes in river stage that promoted exchange of water between the river and bottom sediments. Groundwater discharge measurements indicated fluxes of water across the sediment-water interface as high as 18 m 3 s −1 km −1 . Laboratory experiments indicated that downward movement of surface water through bottom sediments at velocities comparable to those measured in the field (median rate ≈0.005 cm s −1 ) substantially increased dissolved oxygen uptake rates in bottom sediments (maximum rate 212 ± 10 μmol O 2 L −1 h −1 ) compared with rates obtained when no vertical advective flux was generated (maximum rate 25 ± 8.8 μmol O 2 L −1 h −1 ). Additions of dissolved ammonium to surface waters generally increased dissolved oxygen uptake rates relative to rates measured in experiments without ammonium. However, the magnitude of the advective flux through bottom sediments had a greater effect on dissolved oxygen uptake rates than did the availability of ammonium. Results from this study indicated that efforts to improve dissolved oxygen dynamics in effluent-dominated rivers might include stabilizing daily fluctuations in river stage.

Colorado↗
USGS1996

Phytoplankton bloom dynamics in coastal ecosystems: A review with some general lessons from sustained investigation of San Francisco Bay, California

Phytoplankton blooms are prominent features of biological variability in shallow coastal ecosystems such as estuaries, lagoons, bays, and tidal rivers. Long-term observation and research in San Francisco Bay illustrates some patterns of phytoplankton spatial and temporal variability and the underlying mechanisms of this variability. Blooms are events of rapid production and accumulation of phytoplankton biomass that are usually responses to changing physical forcings originating in the coastal ocean (e.g., tides), the atmosphere (wind), or on the land surface (precipitation and river runoff). These physical forcings have different timescales of variability, so algal blooms can be short-term episodic events, recurrent seasonal phenomena, or rare events associated with exceptional climatic or hydrologic conditions. The biogeochemical role of phytoplankton primary production is to transform and incorporate reactive inorganic elements into organic forms, and these transformations are rapid and lead to measurable geochemical change during blooms. Examples include the depletion of inorganic nutrients (N, P, Si), supersaturation of oxygen and removal of carbon dioxide, shifts in the isotopic composition of reactive elements (C, N), production of climatically active trace gases (methyl bromide, dimethylsulfide), changes in the chemical form and toxicity of trace metals (As, Cd, Ni, Zn), changes in the biochemical composition and reactivity of the suspended particulate matter, and synthesis of organic matter required for the reproduction and growth of heterotrophs, including bacteria, zooplankton, and benthic consumer animals. Some classes of phytoplankton play special roles in the cycling of elements or synthesis of specific organic molecules, but we have only rudimentary understanding of the forces that select for and promote blooms of these species. Mounting evidence suggests that the natural cycles of bloom variability are being altered on a global scale by human activities including the input of toxic contaminants and nutrients, manipulation of river flows, and translocation of species. This hypothesis will be a key component of our effort to understand global change at the land-sea interface. Pursuit of this hypothesis will require creative approaches for distinguishing natural and anthropogenic sources of phytoplankton population variability, as well as recognition that the modes of human disturbance of coastal bloom cycles operate interactively and cannot be studied as isolated processes.

California↗
USGS1997

Multiphase groundwater flow near cooling plutons

We investigate groundwater flow near cooling plutons with a computer program that can model multiphase flow, temperatures up to 1200°C, thermal pressurization, and temperature‐dependent rock properties. A series of experiments examines the effects of host‐rock permeability, size and depth of pluton emplacement, single versus multiple intrusions, the influence of a caprock, and the impact of topographically driven groundwater flow. We also reproduce and evaluate some of the pioneering numerical experiments on flow around plutons. Host‐rock permeability is the principal factor influencing fluid circulation and heat transfer in hydrothermal systems. The hottest and most steam‐rich systems develop where permeability is of the order of 10 −15 m 2 . Temperatures and life spans of systems decrease with increasing permeability. Conduction‐dominated systems, in which permeabilities are ≤10 −16 m 2 , persist longer but exhibit relatively modest increases in near‐surface temperatures relative to ambient conditions. Pluton size, emplacement depth, and initial thermal conditions have less influence on hydrothermal circulation patterns but affect the extent of boiling and duration of hydrothermal systems. Topographically driven groundwater flow can significantly alter hydrothermal circulation; however, a low‐permeability caprock effectively decouples the topographically and density‐driven systems and stabilizes the mixing interface between them thereby defining a likely ore‐forming environment.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Nonlinear refraction and reflection travel time tomography

We develop a rapid nonlinear travel time tomography method that simultaneously inverts refraction and reflection travel times on a regular velocity grid. For travel time and ray path calculations, we apply a wave front method employing graph theory. The first-arrival refraction travel times are calculated on the basis of cell velocities, and the later refraction and reflection travel times are computed using both cell velocities and given interfaces. We solve a regularized nonlinear inverse problem. A Laplacian operator is applied to regularize the model parameters (cell slownesses and reflector geometry) so that the inverse problem is valid for a continuum. The travel times are also regularized such that we invert travel time curves rather than travel time points. A conjugate gradient method is applied to minimize the nonlinear objective function. After obtaining a solution, we perform nonlinear Monte Carlo inversions for uncertainty analysis and compute the posterior model covariance. In numerical experiments, we demonstrate that combining the first arrival refraction travel times with later reflection travel times can better reconstruct the velocity field as well as the reflector geometry. This combination is particularly important for modeling crustal structures where large velocity variations occur in the upper crust. We apply this approach to model the crustal structure of the California Borderland using ocean bottom seismometer and land data collected during the Los Angeles Region Seismic Experiment along two marine survey lines. Details of our image include a high-velocity zone under the Catalina Ridge, but a smooth gradient zone between. Catalina Ridge and San Clemente Ridge. The Moho depth is about 22 km with lateral variations. Copyright 1998 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗