Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research”

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 73 records · Page 4Linked to original sources

Change in atmospheric mineral aerosols in response to climate: Last glacial period, preindustrial, modern, and doubled carbon dioxide climates

Desert dust simulations generated by the National Center for Atmospheric Research's Community Climate System Model for the current climate are shown to be consistent with present day satellite and deposition data. The response of the dust cycle to last glacial maximum, preindustrial, modern, and doubled-carbon dioxide climates is analyzed. Only natural (non-land use related) dust sources are included in this simulation. Similar to some previous studies, dust production mainly responds to changes in the source areas from vegetation changes, not from winds or soil moisture changes alone. This model simulates a +92%, +33%, and -60% change in dust loading for the last glacial maximum, preindustrial, and doubled-carbon dioxide climate, respectively, when impacts of carbon dioxide fertilization on vegetation are included in the model. Terrestrial sediment records from the last glacial maximum compiled here indicate a large underestimate of deposition in continental regions, probably due to the lack of simulation of glaciogenic dust sources. In order to include the glaciogenic dust sources as a first approximation, we designate the location of these sources, and infer the size of the sources using an inversion method that best matches the available data. The inclusion of these inferred glaciogenic dust sources increases our dust flux in the last glacial maximum from 2.1 to 3.3 times current deposition. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗

Creep of water ices at planetary conditions: A compilation

Many constitutive laws for the flow of ice have been published since the advent of the Voyager explorations of the outer solar system. Conflicting data have occasionally come from different laboratories, and refinement of experimental techniques has led to the publication of laws that supersede earlier ones. In addition, there are unpublished data from ongoing research that also amend the constitutive laws. Here we compile the most current laboratory-derived flow laws for water ice phases I, II, III, V, and VI, and ice I mixtures with hard particulates. The rheology of interest is mainly that of steady state, and the conditions reviewed are the pressures and temperatures applicable to the surfaces and interiors of icy moons of the outer solar system. Advances in grain-size-dependent creep in ices I and II as well as in phase transformations and metastability under differential stress are also included in this compilation. At laboratory strain rates the several ice polymorphs are rheologically distinct in terms of their stress, temperature, and pressure dependencies but, with the exception of ice III, have fairly similar strengths. Hard particulates strengthen ice I significantly only at high particulate volume fractions. Ice III has the potential for significantly affecting mantle dynamics because it is much weaker than the other polymorphs and its region of stability, which may extend metastably well into what is nominally the ice II field, is located near likely geotherms of large icy moons.

Journal of Geophysical Research E: Planets↗

Bora event variability and the role of air-sea feedback

A two-way interacting high resolution numerical simulation of the Adriatic Sea using the Navy Coastal Ocean Model (NCOM) and Coupled Ocean/ Atmosphere Mesoscale Prediction System (COAMPS??) was conducted to improve forecast momentum and heat flux fields, and to evaluate surface flux field differences for two consecutive bora events during February 2003. (COAMPS?? is a registered trademark of the Naval Research Laboratory.) The strength, mean positions and extensions of the bora jets, and the atmospheric conditions driving them varied considerably between the two events. Bora 1 had 62% stronger heat flux and 51% larger momentum flux than bora 2. The latter displayed much greater diurnal variability characterized by inertial oscillations and the early morning strengthening of a west Adriatic barrier jet, beneath which a stronger west Adriatic ocean current developed. Elsewhere, surface ocean current differences between the two events were directly related to differences in wind stress curl generated by the position and strength of the individual bora jets. The mean heat flux bias was reduced by 72%, and heat flux RMSE reduced by 30% on average at four instrumented over-water sites in the two-way coupled simulation relative to the uncoupled control. Largest reductions in wind stress were found in the bora jets, while the biggest reductions in heat flux were found along the north and west coasts of the Adriatic. In bora 2, SST gradients impacted the wind stress curl along the north and west coasts, and in bora 1 wind stress curl was sensitive to the Istrian front position and strength. The two-way coupled simulation produced diminished surface current speeds of ???12% over the northern Adriatic during both bora compared with a one-way coupled simulation. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗

GCIP water and energy budget synthesis (WEBS)

As part of the World Climate Research Program's (WCRPs) Global Energy and Water-Cycle Experiment (GEWEX) Continental-scale International Project (GCIP), a preliminary water and energy budget synthesis (WEBS) was developed for the period 1996-1999 fromthe "best available" observations and models. Besides this summary paper, a companion CD-ROM with more extensive discussion, figures, tables, and raw data is available to the interested researcher from the GEWEX project office, the GAPP project office, or the first author. An updated online version of the CD-ROM is also available at http://ecpc.ucsd.edu/gcip/webs.htm/. Observations cannot adequately characterize or "close" budgets since too many fundamental processes are missing. Models that properly represent the many complicated atmospheric and near-surface interactions are also required. This preliminary synthesis therefore included a representative global general circulation model, regional climate model, and a macroscale hydrologic model as well as a global reanalysis and a regional analysis. By the qualitative agreement among the models and available observations, it did appear that we now qualitatively understand water and energy budgets of the Mississippi River Basin. However, there is still much quantitative uncertainty. In that regard, there did appear to be a clear advantage to using a regional analysis over a global analysis or a regional simulation over a global simulation to describe the Mississippi River Basin water and energy budgets. There also appeared to be some advantage to using a macroscale hydrologic model for at least the surface water budgets. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗

An eddy-resolving numerical model to study turbulent flow, sediment and bed evolution using detached eddy simulation in a lateral separation zone at the field-scale

Turbulence-resolving simulations elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field-scale, coupling a viscous flow and sediment transport solver using Detached Eddy Simulation techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative, non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars. tions elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field scale, coupling viscous flow and sediment transport solver using Detached Eddy Simulation (DES) techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in massive erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars.

Arizona↗

Submersible observations along the southern Juan de Fuca Ridge: 1984 Alvin program.

In September 1984, the research submersible Alvin provided direct observations of three major hydrothermal vent areas along the southernmost segment of the Juan de Fuca Ridge (JFR). The submersible operations focused on specific volcanologie, structural, and hydrothermal problems that had been identified during the preceding 4 years of photographic, dredging, acoustic imaging, and geophysical studies along a 12-km-long section of the ridge. A continuously maintained (from 1981 to the present) net of seafloor-anchored acoustic transponders allowed the observations from Alvin to be directly tied to all previous U.S. Geological Survey data sets and National Oceanic and Atmospheric Administration water column surveys from 1984 to the present. The three vent areas studied are the largest of at least six areas identified by previous deep-towed camera surveys that lie within a deep cleft, which marks the axis of symmetry of the JFR in this region. The cleft appears to be the locus of eruption for this segment of the JFR. The vent areas, at least in part, are localized near what appear to be previous volcanic eruptive centers marked by extensive lava lake collapse features adjacent to the cleft at these sites. Each hydrothermal area has several active discharge sites, and sulfide deposits occur as clusters (15–100 m 2 ) of small chimneys, individual large chimneys, or clusters of large branched chimneys. We review the dive program and present a brief synthesis of the geology of the vent sites together with sample and track line compilations.

Journal of Geophysical Research B: Solid Earth↗

Atmospheric chemistry of a 33-34 hour old volcanic cloud from Hekla Volcano (Iceland): Insights from direct sampling and the application of chemical box modeling

On 28 February 2000, a volcanic cloud from Hekla volcano, Iceland, was serendipitously sampled by a DC-8 research aircraft during the SAGE III Ozone Loss and Validation Experiment (SOLVE I). It was encountered at night at 10.4 km above sea level (in the lower stratosphere) and 33-34 hours after emission. The cloud is readily identified by abundant SO2 (???1 ppmv), HCl (???70 ppbv), HF (???60 ppbv), and particles (which may have included fine silicate ash). We compare observed and modeled cloud compositions to understand its chemical evolution. Abundances of sulfur and halogen species indicate some oxidation of sulfur gases but limited scavenging and removal of halides. Chemical modeling suggests that cloud concentrations of water vapor and nitric acid promoted polar stratospheric cloud (PSC) formation at 201-203 K, yielding ice, nitric acid trihydrate (NAT), sulfuric acid tetrahydrate (SAT), and liquid ternary solution H2SO4/H2O/HNO3 (STS) particles. We show that these volcanically induced PSCs, especially the ice and NAT particles, activated volcanogenic halogens in the cloud producing >2 ppbv ClOx. This would have destroyed ozone during an earlier period of daylight, consistent with the very low levels of ozone observed. This combination of volcanogenic PSCs and chlorine destroyed ozone at much faster rates than other PSCs that Arctic winter. Elevated levels of HNO3 and NOy in the cloud can be explained by atmospheric nitrogen fixation in the eruption column due to high temperatures and/or volcanic lightning. However, observed elevated levels of HOx remain unexplained given that the cloud was sampled at night. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗

Shrinking ponds in subarctic Alaska based on 1950-2002 remotely sensed images

Over the past 50 years, Alaska has experienced a warming climate with longer growing seasons, increased potential evapotranspiration, and permafrost warming. Research from the Seward Peninsula and Kenai Peninsula has demonstrated a substantial landscape-level trend in the reduction of surface water and number of closed-basin ponds. We investigated whether this drying trend occurred at nine other regions throughout Alaska. One study region was from the Arctic Coastal Plain where depp permafrost occurs continuously across the landscape. The other eight study regions were from the boreal forest regions where discontinuous permafrost occurs. Mean annual precipitation across the study regions ranged from 100 to over 700 min yr-1. We used remotely sensed imagery from the 1950s to 2002 to inventory over 10,000 closed-basin ponds from at least three periods from this time span. We found a reduction in the area and number of shallow, closed-basin ponds for all boreal regions. In contrast, the Arctic Coastal Plain region had negligible change in the area of closed-basin ponds. Since the 1950s, surface water area of closed-basin ponds included in this analysis decreased by 31 to 4 percent, and the total number of closed-basin ponds surveyed within each study region decreased from 54 to 5 percent. There was a significant increasing trend in annual mean temperature and potential evapotranspiration since the 1950s for all study regions. There was no significant trend in annual precipitation during the same period. The regional trend of shrinking ponds may be due to increased drainage as permafrost warms, or increased evapotranspiration during a warmer and extended growing season. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences↗

Topographically driven groundwater flow and the San Andreas heat flow paradox revisited

Evidence for a weak San Andreas Fault includes (1) borehole heat flow measurements that show no evidence for a frictionally generated heat flow anomaly and (2) the inferred orientation of ??1 nearly perpendicular to the fault trace. Interpretations of the stress orientation data remain controversial, at least in close proximity to the fault, leading some researchers to hypothesize that the San Andreas Fault is, in fact, strong and that its thermal signature may be removed or redistributed by topographically driven groundwater flow in areas of rugged topography, such as typify the San Andreas Fault system. To evaluate this scenario, we use a steady state, two-dimensional model of coupled heat and fluid flow within cross sections oriented perpendicular to the fault and to the primary regional topography. Our results show that existing heat flow data near Parkfield, California, do not readily discriminate between the expected thermal signature of a strong fault and that of a weak fault. In contrast, for a wide range of groundwater flow scenarios in the Mojave Desert, models that include frictional heat generation along a strong fault are inconsistent with existing heat flow data, suggesting that the San Andreas Fault at this location is indeed weak. In both areas, comparison of modeling results and heat flow data suggest that advective redistribution of heat is minimal. The robust results for the Mojave region demonstrate that topographically driven groundwater flow, at least in two dimensions, is inadequate to obscure the frictionally generated heat flow anomaly from a strong fault. However, our results do not preclude the possibility of transient advective heat transport associated with earthquakes.

Journal of Geophysical Research B: Solid Earth↗

Plastic faulting in ice

Plastic faulting is a brittle‐like failure phenomenon exhibited by water ice and several other rock types under confinement. It is suspected to be the mechanism of deep earthquakes and extreme cases of shear localization in shallow rocks. Unlike ordinary Coulombic failure, plastic faulting is characterized by a pressure‐independent failure strength and fault plane oriented 45° to maximum principal stress. To research the question of how the instability initiates, we conducted over 50 constant‐displacement‐rate experiments on polycrystalline ice (phases Ih and II) near the brittle‐to‐ductile (B‐D) transition, at confining pressures P = 0–300 MPa, applied strain rates = 5 × 10 −5 – 7 × 10 −3 s −1 , temperatures T = 105–233 K, and mean grain sizes d = 0.25–1.18 mm. We find that (1) the width of the B‐D transition in variable space is vanishingly narrow, to the point of appearing as a crossover, (2) a plastic fault plane, once formed, is not a zone of subsequent weakness, (3) distributed ice I→II phase transformation in small amounts (<1 vol%) shows no causal relationship to subsequent failure, and (4) plastic faulting also occurs in ice II. We hypothesize that the elusive nucleating “trigger” parallels that of metals and ceramics undergoing severe plastic deformation, wherein transient local structural rearrangement occurs, in turn causing material strength to drop to a level sufficiently low, in a volume sufficiently large, that adiabatic instability is nucleated. Our results do not require and often are inconsistent with phase transformation. Plastic faulting may therefore be available to all solids undergoing severe deformation, and its appearance in so few is simply the result of insufficiently extreme conditions.

Journal of Geophysical Research- Solid Earth↗

Responses of stream nitrate and dissolved organic carbon loadings to hydrological forcing and climate change in an upland forest of the northeast USA

[1] In coming decades, higher annual temperatures, increased growing season length, and increased dormant season precipitation are expected across the northeastern United States in response to anthropogenic forcing of global climate. We synthesized long-term stream hydrochemical data from the Sleepers River Research Watershed in Vermont, United States, to explore the relationship of catchment wetness to stream nitrate and DOC loadings. We modeled changes in growing season length and precipitation patterns to simulate future climate scenarios and to assess how stream nutrient loadings respond to climate change. Model results for the 2070&ndash;2099 time period suggest that stream nutrient loadings during both the dormant and growing seasons will respond to climate change. During a warmer climate, growing season stream fluxes (runoff +20%, nitrate +57%, and DOC +58%) increase as more precipitation (+28%) and quick flow (+39%) occur during a longer growing season (+43 days). During the dormant season, stream water and nutrient loadings decrease. Net annual stream runoff (+8%) and DOC loading (+9%) increases are commensurate with the magnitude of the average increase of net annual precipitation (+7%). Net annual stream water and DOC loadings are primarily affected by increased dormant season precipitation. In contrast, decreased annual loading of stream nitrate (&minus;2%) reflects a larger effect of growing season controls on stream nitrate and the effects of lengthened growing seasons in a warmer climate. Our findings suggest that leaching of nitrate and DOC from catchment soils will be affected by anthropogenic climate forcing, thereby affecting the timing and magnitude of annual stream loadings in the northeastern United States.

Vermont↗

Cataloging tectonic tremor energy radiation in the Cascadia subduction zone

For the past ∼12 years the Pacific Northwest Seismic Network has been automatically detecting and locating tectonic tremor across the Cascadia subduction zone, resulting in a catalog of more than 500,000 tremor epicenters to date, which has served as a valuable resource for tremor and slip research. This manuscript presents an updated methodology for routine tremor detection in Cascadia and a new catalog of over 180,000 tremor epicenters including amplitudes detected along the subduction zone margin from 2017 to 2021. The events are detected via cross-correlation of continuous vertical envelope data of 128 stations from northern California to northern Vancouver Island. The modified approach results in less scatter and a 55% increase in detected epicenters than previously observed, as well as a newly identified tremor source offset updip from the main tremor and slip region at the southern edge of the subduction zone. Radiated seismic energy in the 1.5–5 Hz band is used to assign epicenters an energy magnitude ( M eL ), which is calibrated to the M L of local earthquakes. Southern Cascadia is most active, but the highest tremor energy rates occur in northern Cascadia. Tremor in central Cascadia is systematically weaker and less frequent. Individual epicenter magnitudes range from ∼0.5–2 and spatiotemporally cluster into 1,060 swarms with cumulative M eL ranging from ∼0.8 to 3.7. The swarms reflect underlying slow slip events and occur with an earthquake-like energy distribution with a b value ∼1. Tremor epicenters, however, follow a tapered Gutenberg-Richter distribution with high b values, suggesting individual tremor bursts and their constituent low-frequency earthquakes are fault-dimension limited.

California, Oregon, Washington↗

Mars global digital dune database and initial science results

A new Mars Global Digital Dune Database (MGD 3 ) constructed using Thermal Emission Imaging System (THEMIS) infrared (IR) images provides a comprehensive and quantitative view of the geographic distribution of moderate‐ to large‐size dune fields (area >1 km 2 ) that will help researchers to understand global climatic and sedimentary processes that have shaped the surface of Mars. MGD 3 extends from 65°N to 65°S latitude and includes ∼550 dune fields, covering ∼70,000 km 2 , with an estimated total volume of ∼3,600 km 3 . This area, when combined with polar dune estimates, suggests moderate‐ to large‐size dune field coverage on Mars may total ∼800,000 km 2 , ∼6 times less than the total areal estimate of ∼5,000,000 km 2 for terrestrial dunes. Where availability and quality of THEMIS visible (VIS) or Mars Orbiter Camera narrow‐angle (MOC NA) images allow, we classify dunes and include dune slipface measurements, which are derived from gross dune morphology and represent the prevailing wind direction at the last time of significant dune modification. For dunes located within craters, the azimuth from crater centroid to dune field centroid (referred to as dune centroid azimuth) is calculated and can provide an accurate method for tracking dune migration within smooth‐floored craters. These indicators of wind direction are compared to output from a general circulation model (GCM). Dune centroid azimuth values generally correlate to regional wind patterns. Slipface orientations are less well correlated, suggesting that local topographic effects may play a larger role in dune orientation than regional winds.

Journal of Geophysical Research E: Planets↗

Comparison of atmospheric mercury speciation and deposition at nine sites across central and eastern North America

This study presents >5 cumulative years of tropospheric mercury (Hg) speciation measurements, over the period of 2003–2009, for eight sites in the central and eastern United States and one site in coastal Puerto Rico. The purpose of this research was to identify local and regional processes that impact Hg speciation and deposition (wet + dry) across a large swath of North America. Sites sampled were selected to represent both a wide range of mercury exposure and environmental conditions. Seasonal mean concentrations of elemental Hg (1.27 ± 0.31 to 2.94 ± 1.57 ng m−3; x ± s), reactive gaseous mercury (RGM; 1.5 ± 1.6 to 63.3 ± 529 pg m−3), and fine particulate Hg (1.2 ± 1.4 to 37.9 ± 492 pg m−3) were greatest at sites impacted by Hg point sources. Diel bin plots of Hgo and RGM suggest control by a variety of local/regional processes including impacts from Hg point sources and boundary layer/free tropospheric interactions as well as from larger‐scale processes affecting Hg speciation (i.e., input of the global Hg pool, RGM formed from oxidation of Hgo by photochemical compounds at coastal sites, and elemental Hg depletion during periods of dew formation). Comparison of wet Hg deposition (measured), RGM and fine particulate Hg dry deposition (calculated using a multiple resistance model), and anthropogenic point source emissions varied significantly between sites. Significant correlation between emission sources and dry deposition was observed but was highly dependant upon inclusion of data from two sites with exceptionally high deposition. Findings from this study highlight the importance of environmental setting on atmospheric Hg cycling and deposition rates.

Journal of Geophysical Research↗

Dating the Vostok ice core record by importing the Devils Hole chronology

The development of an accurate chronology for the Vostok record continues to be an open research question because these invaluable ice cores cannot be dated directly. Depth-to-age relationships have been developed using many different approaches, but published age estimates are inconsistent, even for major paleoclimatic events. We have developed a chronology for the Vostok deuterium paleotemperature record using a simple and objective algorithm to transfer ages of major paleoclimatic events from the radiometrically dated 500,000-year δ 18 O-paleotemperature record from Devils Hole, Nevada. The method is based only on a strong inference that major shifts in paleotemperature recorded at both locations occurred synchronously, consistent with an atmospheric teleconnection. The derived depth-to-age relationship conforms with the physics of ice compaction, and internally produces ages for climatic events 5.4 and 11.24 which are consistent with the externally assigned ages that the Vostok team needed to assume in order to derive their most recent chronology, GT4. Indeed, the resulting V-DH chronology is highly correlated with GT4 because of the unexpected correspondence even in the timing of second-order climatic events that were not constrained by the algorithm. Furthermore, the algorithm developed herein is not specific to this problem; rather, the procedure can be used whenever two paleoclimate records are proxies for the same physical phenomenon, and paleoclimatic conditions forcing the two records can be considered to have occurred contemporaneously. The ability of the algorithm to date the East Antarctic Dome Fuji core is also demonstrated.

Journal of Geophysical Research D: Atmospheres↗

Geological and thermal control of the hydrothermal system in northern Yellowstone Lake: Inferences from high resolution magnetic surveys

A multiscale magnetic survey of the northern basin of Yellowstone Lake was undertaken in 2016 as part of the Hydrothermal Dynamics of Yellowstone Lake Project (HD‐YLAKE)—a broad research effort to characterize the cause‐and‐effect relationships between geologic and environmental processes and hydrothermal activity on the lake floor. The magnetic survey includes lake surface, regional aeromagnetic, and near‐bottom autonomous underwater vehicle (AUV) data. The study reveals a strong contrast between the northeastern lake basin, characterized by a regional magnetic low punctuated by stronger local magnetic lows, many of which host hydrothermal vent activity, and the northwestern lake basin with higher‐amplitude magnetic anomalies and no obvious hydrothermal activity or punctuated magnetic lows. The boundary between these two regions is marked by a steep gradient in heat flow and magnetic values, likely reflecting a significant structure within the currently active ~20‐km‐long Eagle Bay‐Lake Hotel fault zone that may be related to the ~2.08‐Ma Huckleberry Ridge caldera rim. Modeling suggests that the broad northeastern magnetic low reflects both a shallower Curie isotherm and widespread hydrothermal activity that has demagnetized the rock. Along the western lake shoreline are sinuous‐shaped, high‐amplitude magnetic anomaly highs, interpreted as lava flow fronts of upper units of the West Thumb rhyolite. The AUV magnetic survey shows decreased magnetization at the periphery of the active Deep Hole hydrothermal vent. We postulate that lower magnetization in the outer zone results from enhanced hydrothermal alteration of rhyolite by hydrothermal condensates while the vapor‐dominated center of the vent is less altered.

Wyoming↗

Structure of the lower crust beneath the Carolina Trough, U.S. Atlantic continental margin

Data from three large-offset seismic profiles provide information on the crustal structure beneath the Carolina trough. The profiles, obtained by the U.S. Geological Survey, the Naval Oceanographic Research Development Agency, and the Scripps Institution of Oceanography in 1985, were oriented parallel to the trough and were located (1) seaward of the East Coast Magnetic Anomaly (ECMA), which is generally thought to represent the boundary between oceanic and continental crust; (2) along the axis of the trough between the ECMA and the hinge zone, which is thought to reflect the landward limit of highly stretched and altered transitional crust; and (3) along the Carolina platform landward of the basement hinge zone on crust thought to have been thinned only slightly during rifting. These data constrain the velocity structure of the lower crust and provide evidence for a thick lens of high-velocity (>7.1 km/s) lower crustal material that extends beneath the Carolina trough and the adjacent ocean basin. This lens reaches a maximum thickness of about 13 km beneath the deepest part of the trough, thins to about 5 km seaward of the ECMA, and is either very thin or absent landward of the hinge zone. It is interpreted to represent material that was underplated beneath and/or intruded into the crust during the late stage of continental rifting and that led to an anomalously thick plutonic layer during the early seafloor spreading phase. These data thus support the recent conclusions of White et al. (1987b) and Mutter et al. (1988) that the initiation of seafloor spreading is attended in many, if not most, cases by the generation of an anomalously large volume of melt.

Carolina Trough↗

Quantifying interdependencies in geyser eruptions at the Upper Geyser Basin, Yellowstone National Park

The Upper Geyser Basin at Yellowstone National Park (Wyoming, USA) harbors the greatest concentration of geysers worldwide. Research suggests that individual geysers are not isolated but rather are hydraulically connected in the subsurface with other geysers and thermal springs. To quantify such connections, we combined techniques from machine learning, causal inference, and dynamical systems to characterize the collective eruptive behavior of a set of 10 geysers over 18 months (April 2007 – September 2008) focusing on geyser-geyser interactions. Model predictions were up to 15 times more accurate when we sought to predict a geyser's eruption time series based on outflow channel temperatures from the network than based on its own time series alone, suggesting the existence of a complex interconnected subsurface groundwater system. On average, cone-type geysers had larger impacts on other geysers than did fountain-type geysers. Similarly, cone-type geysers were on average more insulated from other geysers. However, substantial unexplained variation remained after considering the cone versus fountain dichotomy. Distance between geysers also affected interactions: nearby geysers had stronger effects on focal geysers than did geysers located farther away. Collectively, results support the hypothesis of geyser interdependence at timescales of 5 min–10 days. Our analyses highlight the existence of quantifiable geyser-to-geyser interactions that can be resolved through pairwise and system-level analyses. These findings emphasize the subsurface interconnectedness of thermal features, provide information relevant to visitor experiences in Yellowstone National Park, and suggest strategies for exploring patterns of interdependence that may exist among other episodic geological phenomena.

Wyoming↗