USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Physics of the Earth”

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 109 records · Page 6Linked to original sources
USGS2022

Multi-decadal simulation of marsh topography evolution under sea level rise and episodic sediment loads

Coastal marsh within Mediterranean climate zones is exposed to episodic watershed runoff and sediment loads that occur during storm events. Simulating future marsh accretion under sea level rise calls for attention to: (a) physical processes acting over the time scale of storm events and (b) biophysical processes acting over time scales longer than storm events. Using the upper Newport Bay in Southern California as a case study, we examine the influence of event-scale processes on simulated change in marsh topography by comparing: (a) a biophysical model that integrates with an annual time step and neglects event-scale processes (BP-Annual), (b) a physical model that resolves event-scale processes but neglects biophysical interactions (P-Event), and (c) a biophysical model that resolves event-scale physical processes and biophysical processes at annual and longer time scales (BP-Event). A calibrated BP-Event model shows that large (>20-year return period) episodic storm events are major drivers of marsh accretion, depositing up to 30 cm of sediment in one event. Greater deposition is predicted near fluvial sources and tidal channels and less on marshes further from fluvial sources and tidal channels. In contrast, the BP-Annual model poorly resolves spatial structure in marsh accretion as a consequence of neglecting event-scale processes. Furthermore, the P-Event model significantly overestimates marsh accretion as a consequence of neglecting marsh surface compaction driven by annual scale biophysical processes. Differences between BP-Event and BP-Annual models translate up to 20 cm per century in marsh surface elevation.

California↗
USGS2022

Detection of aseismic slip and poroelastic reservoir deformation at the North Brawley Geothermal Field from 2009 to 2019

The North Brawley Geothermal Field, located within the Brawley Seismic Zone of Southern California, presents a case study for understanding seismic hazards linked to fluid injection and geothermal energy extraction. An earthquake swarm near the geothermal field in 2012 included two earthquakes with magnitudes greater than 5 and was potentially preceded by a years-long aseismic slip transient. To better understand ground deformation around the geothermal field, including its evolution with time and its physical mechanisms, we analyze deformation before, during, and after the swarm using ground- and satellite-based geodetic techniques between 2009 and 2019. We integrate observations from GNSS, Sentinel-1, TerraSAR-X, UAVSAR, and leveling surveys into a single deformation history. Modeling of this new collection of observations at the North Brawley Geothermal Field provides evidence for 80% more pre-swarm aseismic slip than previously recognized from 2009 to 2012. During the 2012 Brawley swarm, our geodetic slip inversions closely match the results of seismic waveform inversions from the swarm events. After the 2012 swarm, surface deformation is dominated by poroelastic deformation of a shallow fluid reservoir at <1 km depth rather than fault slip. The deformation history and seismicity catalogs at North Brawley suggest a cessation of fault-related slip during the ∼7 years after the 2012 earthquake swarm.

California↗
USGS2018

Consequences of abrading bed load on vertical and lateral bedrock erosion in a curved experimental channel

In this study, we conducted multiple physical experiments to estimate the efficacy and spatial pattern of erosion by abrading sediment moving through a simple U‐shaped channel bend with erodible bed and banks. The experiments showed that in the bend, lateral abrasion followed a monotonically increasing linear relationship with sediment feed rate. However, vertical incision had a more complex relation with the sediment feed rate, with an initial increase in abrasion as the feed rate increased followed by a decrease in abrasion of the bed as cover effects became dominant at higher feed rates. Bank erosion was large in places where the width and the lateral slope of the point bar were relatively large. On the other hand, in places where the width of the point bar was smaller, the bedrock bed was eroded primarily along the boundary of the point bar, resulting in a bedrock bench near the outer bank.

Shimanto River↗
USGS2002

A bibliography of terrain modeling (geomorphometry), the quantitative representation of topography: Supplement 4.0

Terrain modeling, the practice of ground-surface quantification, is an amalgam of Earth science, mathematics, engineering, and computer science. The discipline is known variously as geomorphometry (or simply morphometry), terrain analysis, and quantitative geomorphology. It continues to grow through myriad applications to hydrology, geohazards mapping, tectonics, sea-floor and planetary exploration, and other fields. Dating nominally to the co-founders of academic geography, Alexander von Humboldt (1808, 1817) and Carl Ritter (1826, 1828), the field was revolutionized late in the 20th Century by the computer manipulation of spatial arrays of terrain heights, or digital elevation models (DEMs), which can quantify and portray ground-surface form over large areas (Maune, 2001). Morphometric procedures are implemented routinely by commercial geographic information systems (GIS) as well as specialized software (Harvey and Eash, 1996; Köthe and others, 1996; ESRI, 1997; Drzewiecki et al., 1999; Dikau and Saurer, 1999; Djokic and Maidment, 2000; Wilson and Gallant, 2000; Breuer, 2001; Guth, 2001; Eastman, 2002). The new Earth Surface edition of the Journal of Geophysical Research, specializing in surficial processes, is the latest of many publication venues for terrain modeling. This is the fourth update of a bibliography and introduction to terrain modeling (Pike, 1993, 1995, 1996, 1999) designed to collect the diverse, scattered literature on surface measurement as a resource for the research community. The use of DEMs in science and technology continues to accelerate and diversify (Pike, 2000a). New work appears so frequently that a sampling must suffice to represent the vast literature. This report adds 1636 entries to the 4374 in the four earlier publications1. Forty-eight additional entries correct dead Internet links and other errors found in the prior listings. Chronicling the history of terrain modeling, many entries in this report predate the 1999 supplement. Coverage is representative from about 1800 through early–mid 2002. Papers increasingly are published exclusively or in duplicate on the Internet's World Wide Web; the dates given here for Web addresses (URLs) that lack a print publication indicate a Web site's last update or my last access of it. The bibliography is arranged alphabetically and thus is not readily summarized. This introduction cites about 500 entries, a third of them grouped under 24 morphometric topics, as a guide to the listing's contents. Continuing the practice of previous bibliographies in the series to provide more information on a few applications (see summary of past topics in Pike, 2000a), this report elaborates further on topographic data, putative new parameters, tectonic geomorphology/neo-orometry, biogeography, ice-cap morphometry, results from the Mars Global DEM, landslide-hazard mapping, terrain modeling as physics, Hack's law, and broad-scale computer visualization. The literature of some of these subjects is large, and none of the summaries is intended to more than introduce the topic and comment on some of the current contributions of terrain modeling. Closing the essay is a discussion of pre-1900 papers that trace the evolution of ridge-line and watercourse quantification by descriptive geometry, as well as comments on some new books and an on-line bulletin board.

Open-File Report↗
USGS2012

Tektites in Cretaceous-Tertiary boundary rocks on Haiti and their bearing on the Alvarez impact extinction hypothesis

Relic tektites are associated with a Pt-group metal abundance anomaly and shocked minerals in a thin marl bed that marks the K-T boundary on Haiti. The presence of these three impact-produced materials at the precise K-T boundary enormously strengthens the Alvarez impact extinction hypothesis. The tektites occur in smectite spherules that have external shapes typical of tektites. Their chemical and physical properties are broadly similar to those of other tektite groups, except that the Haitian tektites have lower Si and higher Fe, Ca, and Na. On average, they contain more Sc, V, Cu, Zn, Ga, Sr, Sn, and Ba and less Cr, Ni, Co, B, Mn, and Hf than Other tektite groups. Amounts of rare earth elements (REE) in the tektites indicate that their progenitor materials were not melted mafic or ultramafic oceanic crust; rather they were sedimentary deposits having a bulk composition of andesite. Rare tektites contain unusually high amounts of CaO (∼20%) and S (0.4%), and these data suggest that some target materials consisted of CaSO 4 . Anhydrite beds occur in the subsurface at two candidate impact sites (Chicxulub and Manson). Sm-Nd isotopic data for the tektites indicate that the melted precursor sediments were most likely deposited less than −400 m.y. between the Late Cretaceous and Silurian periods. Major chemical changes accompanied the diagenetic change of glass to smectite. The Haitian tektites are the first datable impact products in K-T boundary rocks, and 40 Ar- 39 Ar ages of the glass show that the K-T boundary and impact event are coeval at 64.5±0.1 Ma.

Journal of Geophysical Research↗
USGS2011

Detecting aseismic strain transients from seismicity data

Aseismic deformation transients such as fluid flow, magma migration, and slow slip can trigger changes in seismicity rate. We present a method that can detect these seismicity rate variations and utilize these anomalies to constrain the underlying variations in stressing rate. Because ordinary aftershock sequences often obscure changes in the background seismicity caused by aseismic processes, we combine the stochastic Epidemic Type Aftershock Sequence model that describes aftershock sequences well and the physically based rate‐ and state‐dependent friction seismicity model into a single seismicity rate model that models both aftershock activity and changes in background seismicity rate. We implement this model into a data assimilation algorithm that inverts seismicity catalogs to estimate space‐time variations in stressing rate. We evaluate the method using a synthetic catalog, and then apply it to a catalog of M ≥ 1.5 events that occurred in the Salton Trough from 1990 to 2009. We validate our stressing rate estimates by comparing them to estimates from a geodetically derived slip model for a large creep event on the Obsidian Buttes fault. The results demonstrate that our approach can identify large aseismic deformation transients in a multidecade long earthquake catalog and roughly constrain the absolute magnitude of the stressing rate transients. Our method can therefore provide a way to detect aseismic transients in regions where geodetic resolution in space or time is poor.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

Volcanic plume height measured by seismic waves based on a mechanical model

In August 2008 an unmonitored, largely unstudied Aleutian volcano, Kasatochi, erupted catastrophically. Here we use seismic data to infer the height of large eruptive columns such as those of Kasatochi based on a combination of existing fluid and solid mechanical models. In so doing, we propose a connection between a common, observable, short-period seismic wave amplitude to the physics of an eruptive column. To construct a combined model, we estimate the mass ejection rate of material from the vent on the basis of the plume height, assuming that the height is controlled by thermal buoyancy for a continuous plume. Using the estimated mass ejection rate, we then derive the equivalent vertical force on the Earth through a momentum balance. Finally, we calculate the far-field surface waves resulting from the vertical force. The model performs well for recent eruptions of Kasatochi and Augustine volcanoes if v , the velocity of material exiting the vent, is 120-230 m s -1 . The consistency between the seismically inferred and measured plume heights indicates that in these cases the far-field ~1 s seismic energy radiated by fluctuating flow in the volcanic jet during the eruption is a useful indicator of overall mass ejection rates. Thus, use of the model holds promise for characterizing eruptions and evaluating ash hazards to aircraft in real time on the basis of far-field short-period seismic data. This study emphasizes the need for better measurements of eruptive plume heights and a more detailed understanding of the full spectrum of seismic energy radiated coeruptively.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Compositional trends in aeolian dust along a transect across the southwestern United States

Aeolian dust strongly influences ecology and landscape geochemistry over large areas that span several desert ecosystems of the southwestern United States. This study evaluates spatial and temporal variations and trends of the physical and chemical properties of dust in the southwestern United States by examining dust deposited in natural depressions on high isolated surfaces along a transect from the Mojave Desert to the central Colorado Plateau. Aeolian dust is recognized in these depressions on the basis of textural, chemical, isotopic, and mineralogical characteristics and comparisons of those characteristics to the underlying bedrock units. Spatial and temporal trends suggest that although local dust sources are important to the accumulated material in these depressions, Mojave Desert dust sources may also contribute. Depth trends in the depressions suggest that Mojave sources may have contributed more dust to the Colorado Plateau recently than in the past. These interpretations point to the important roles of far-traveled aeolian dust for landscape geochemistry and imply future changes to soil geochemistry under changing conditions in far-distant dust source areas. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2008

Distributed energy balance modeling of South Cascade Glacier, Washington and assessment of model uncertainty

We have developed a physically based, distributed surface energy balance model to simulate glacier mass balance under meteorological and climatological forcing. Here we apply the model to estimate summer ablation on South Cascade Glacier, Washington, for the 2004 and 2005 mass balance seasons. To arrive at optimal mass balance simulations, we investigate and quantify model uncertainty associated with selecting from a range of physical parameter values that are not commonly measured in glaciological mass balance field studies. We optimize the performance of the model by varying values for atmospheric transmissivity, the albedo of surrounding topography, precipitation-elevation lapse rate, surface roughness for turbulent exchange of momentum, and snow albedo aging coefficient. Of these the snow aging parameter and precipitation lapse rates have the greatest influence on the modeled ablation. We examined model sensitivity to varying parameters by performing an additional 103 realizations with parameters randomly chosen over a ??5% range centered about the optimum values. The best fit suite of model parameters yielded a net balance of -1.69??0.38 m water equivalent (WE) for the 2004 water year and -2.10??0.30 m WE up to 11 September 2005. The 2004 result is within 3% of the measured value. These simulations account for 91% and 93% of the variance in measured ablation for the respective years. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2008

Constitutive relationships and physical basis of fault strength due to flash heating

We develop a model of fault strength loss resulting from phase change at asperity contacts due to flash heating that considers a distribution of contact sizes and nonsteady state evolution of fault strength with displacement. Laboratory faulting experiments conducted at high sliding velocities, which show dramatic strength reduction below the threshold for bulk melting, are well fit by the model. The predicted slip speed for the onset of weakening is in the range of 0.05 to 2 m/s, qualitatively consistent with the limited published observations. For this model, earthquake stress drops and effective shear fracture energy should be linearly pressure-dependent, whereas the onset speed may be pressure-independent or weakly pressure-dependent. On the basis of the theory, flash weakening is expected to produce large dynamic stress drops, small effective shear fracture energy, and undershoot. Estimates of the threshold slip speed, stress drop, and fracture energy are uncertain due to poor knowledge of the average ontact dimension, shear zone thickness and gouge particle size at seismogenic depths. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Seismicity and deformation induced by magma accumulation at three basaltic volcanoes

We analyzed the evolution of volcano‐tectonic (VT) seismicity and deformation at three basaltic volcanoes (Kilauea, Mauna Loa, Piton de la Fournaise) during phases of magma accumulation. We observed that the VT earthquake activity displays an accelerating evolution at the three studied volcanoes during the time of magma accumulation. At the same times, deformation rates recorded at the summit of Kilauea and Mauna Loa volcanoes were not accelerating but rather tend to decay. To interpret these observations, we propose a physical model describing the evolution of pressure produced by the accumulation of magma into a reservoir. This variation of pressure is then used to force a simple model of damage, where damage episodes are equivalent to earthquakes. This model leads to an exponential increase of the VT activity and to an exponential decay of the deformation rate during accumulation phases. Seismicity and deformation data are well fitted by such an exponential model. The time constant, deduced from the exponential increase of the seismicity, is in agreement with the time constant predicted by the model of magma accumulation. This VT activity can thus be a direct indication of the accumulation of magma at depth, and therefore can be seen as a long‐term precursory phenomenon, at least for the three studied basaltic volcanoes. Unfortunately, it does not allow the prediction of the onset of future eruptions, as no diverging point (i.e., critical time) is present in the model.

Journal of Geophysical Research B: Solid Earth↗
USGS2010

A universal approximation to grain size from images of non-cohesive sediment

The two-dimensional spectral decomposition of an image of sediment provides a direct statistical estimate, grid-by-number style, of the mean of all intermediate axes of all single particles within the image. We develop and test this new method which, unlike existing techniques, requires neither image processing algorithms for detection and measurement of individual grains, nor calibration. The only information required of the operator is the spatial resolution of the image. The method is tested with images of bed sediment from nine different sedimentary environments (five beaches, three rivers, and one continental shelf), across the range 0.1 mm to 150 mm, taken in air and underwater. Each population was photographed using a different camera and lighting conditions. We term it a “universal approximation” because it has produced accurate estimates for all populations we have tested it with, without calibration. We use three approaches (theory, computational experiments, and physical experiments) to both understand and explore the sensitivities and limits of this new method. Based on 443 samples, the root-mean-squared (RMS) error between size estimates from the new method and known mean grain size (obtained from point counts on the image) was found to be ±≈16%, with a 95% probability of estimates within ±31% of the true mean grain size (measured in a linear scale). The RMS error reduces to ≈11%, with a 95% probability of estimates within ±20% of the true mean grain size if point counts from a few images are used to correct bias for a specific population of sediment images. It thus appears it is transferable between sedimentary populations with different grain size, but factors such as particle shape and packing may introduce bias which may need to be calibrated for. For the first time, an attempt has been made to mathematically relate the spatial distribution of pixel intensity within the image of sediment to the grain size.

Journal of Geophysical Research F: Earth Surface↗
USGS2020

SPEAR: The next generation GFDL modeling system for seasonal to multidecadal prediction and projection

We document the development and simulation characteristics of the next generation modeling system for seasonal to decadal prediction and projection at the Geophysical Fluid Dynamics Laboratory (GFDL). SPEAR ( S eamless System for P rediction and EA rth System R esearch) is built from component models recently developed at GFDL—the AM4 atmosphere model, MOM6 ocean code, LM4 land model, and SIS2 sea ice model. The SPEAR models are specifically designed with attributes needed for a prediction model for seasonal to decadal time scales, including the ability to run large ensembles of simulations with available computational resources. For computational speed SPEAR uses a coarse ocean resolution of approximately 1.0° (with tropical refinement). SPEAR can use differing atmospheric horizontal resolutions ranging from 1° to 0.25°. The higher atmospheric resolution facilitates improved simulation of regional climate and extremes. SPEAR is built from the same components as the GFDL CM4 and ESM4 models but with design choices geared toward seasonal to multidecadal physical climate prediction and projection. We document simulation characteristics for the time mean climate, aspects of internal variability, and the response to both idealized and realistic radiative forcing change. We describe in greater detail one focus of the model development process that was motivated by the importance of the Southern Ocean to the global climate system. We present sensitivity tests that document the influence of the Antarctic surface heat budget on Southern Ocean ventilation and deep global ocean circulation. These findings were also useful in the development processes for the GFDL CM4 and ESM4 models.

Journal of Advances in Modeling Earth Systems↗
USGS2025

Advancing current understanding of Martian impact-generated hydrothermal systems through novel coupled modeling: Insights from Gale, Jezero, and other craters

Impact-generated hydrothermal systems are considered potentially habitable environments on Mars, Earth, and other planetary bodies for microbial life. However, there is an ongoing debate regarding what geological features on Mars provide definitive evidence for such systems. Although earlier studies have modeled hydrothermal processes in Martian craters, they often lacked integration with shock physics hydrocodes to constrain initial impact conditions. The importance of this two-code coupling was demonstrated by successfully replicating alteration signatures in the Earth's Haughton impact structure. In this study, we use a similar two-code approach, combining the iSALE hydrocode with the HYDROTHERM hydrothermal model to simulate the full evolution of impact-generated hydrothermal systems. We apply this method to craters the size of Jezero (∼50 km) and Gale (∼154 km) in diameter. Although Jezero's interior is largely buried, our results align with hypothesized hydrothermal vents and alteration minerals near central uplifts in similarly sized exposed craters, such as Toro and Auki. Furthermore, our models correspond to alteration patterns observed by the Curiosity in the lower layers of Mount Sharp, which may represent remnants of impact-driven hydrothermal activity. A key finding is that these systems may persist much longer than previously estimated. Our simulations suggest that a Jezero-sized system could remain habitable for thermophiles for approximately 720,000 years, whereas a Gale-sized system could persist for nearly 2 million years. Additionally, simulations under unsaturated crustal conditions reveal that air-dominated near-surface layers can suppress vertical fluid flow, enabling deep subsurface alteration without producing detectable mineral signatures at the surface.

Journal of Geophysical Research - Planets↗
USGS2018

Can the earth be dated from decay of its magnetic field?.

Thomas G. Barnes, Emeritus Professor of Physics at the University of Texas, El Paso, and a creationist, argues that the geomagnetic field was created by unknown processes at the time of creation of the earth and has been decaying irreversibly and exponentially, with a half-life of about 1400 years, since then. He calculates that the field would have been impossibly large in 8000 B.C. and concludes that the earth, therefore, is less than 10,000 years old. He denies that the earth's field has reversed polarity in the past, and attempts to refute the dynamo theory of the origin of the field. Barnes' arguments and conclusions are commonly cited in creationist literature as definitive proof that the earth is very young. Barnes is wrong, and has ignored or misrepresented much of the data on the earth's magnetic field. Barnes calculations are based on observatory measurements showing that the earth's dipole field strength has decreased since 1835. But these same measurements also show a corresponding increase in the strength of the nondipole field so that the total field energy external to the core has remained about constant. Paleomagnetic evidence shows that the earth's field has existed for more than three billion years, and that the dipole field both fluctuates in strength and irregularly reverses polarity. The earth's field is generated by a fluid dynamo in the iron-nickel core; all of the elements necessary for a dynamo, including fluid motion and energy sources, are present. Contrary to Barnes' assertions, there are no theoretical reasons why a dynamo cannot exist in the earth. There are no properties of the magnetid field that can be used to place an upper limit on the earth's age.

Journal of Geological Education↗
USGS1969

Oxygen isotope fractionation in divalent metal carbonates

Equilibrium fractionation factors for the distribution of 18 O between alkaline‐earth carbonates and water have been measured over the temperature range 0–500°C. The fractionation factors α"> α α can be represented by the equations CaCO3–H2O, 1000lnα = 2.78(106T−2)− 3.39,"> CaCO 3 – H 2 O , 1000 ln α = 2.78(10 6 T −2 )− 3.39, CaCO3–H2O, 1000lnα = 2.78(106T−2)− 3.39, SrCO3–H2O, 1000lnα = 2.69(106T−2)− 3.74,"> SrCO 3 – H 2 O , 1000 ln α = 2.69(10 6 T −2 )− 3.74, SrCO3–H2O, 1000lnα = 2.69(106T−2)− 3.74, BaCO3–H2O, 1000lnα = 2.57(106T−2)− 4.73."> BaCO 3 – H 2 O , 1000 ln α = 2.57(10 6 T −2 )− 4.73. BaCO3–H2O, 1000lnα = 2.57(106T−2)− 4.73. Measurements on MnCO 3 , CdCO 3 , and PbCO 3 were made at isolated temperatures. A statistical‐mechanical calculation of the isotopic partition function ratios gives reasonably good agreement with experiment. Both cationic size and mass are important in isotopic fractionation, the former predominantly in its effect on the internal vibrations of the anion, the latter in its effect on the lattice vibrations.

Journal of Chemical Physics↗
USGS2001

Digital terrain modelling and industrial surface metrology - Converging crafts

Quantitative characterisation of surface form, increasingly from digital 3-D height data, is cross-disciplinary and can be applied at any scale. Thus, separation of industrial-surface metrology from its Earth-science counterpart, (digital) terrain modelling, is artificial. Their growing convergence presents an opportunity to develop in surface morphometry a unified approach to surface representation. This paper introduces terrain modelling and compares it with metrology, noting their differences and similarities. Examples of potential redundancy among parameters illustrate one of the many issues common to both disciplines. ?? 2001 Elsevier Science Ltd. All rights reserved.

International Journal of Machine Tools and Manufac↗
USGS2013

Bayesian inversion of data from effusive volcanic eruptions using physics-based models: Application to Mount St. Helens 2004--2008

Physics-based models of volcanic eruptions can directly link magmatic processes with diverse, time-varying geophysical observations, and when used in an inverse procedure make it possible to bring all available information to bear on estimating properties of the volcanic system. We develop a technique for inverting geodetic, extrusive flux, and other types of data using a physics-based model of an effusive silicic volcanic eruption to estimate the geometry, pressure, depth, and volatile content of a magma chamber, and properties of the conduit linking the chamber to the surface. A Bayesian inverse formulation makes it possible to easily incorporate independent information into the inversion, such as petrologic estimates of melt water content, and yields probabilistic estimates for model parameters and other properties of the volcano. Probability distributions are sampled using a Markov-Chain Monte Carlo algorithm. We apply the technique using GPS and extrusion data from the 2004–2008 eruption of Mount St. Helens. In contrast to more traditional inversions such as those involving geodetic data alone in combination with kinematic forward models, this technique is able to provide constraint on properties of the magma, including its volatile content, and on the absolute volume and pressure of the magma chamber. Results suggest a large chamber of >40 km 3 with a centroid depth of 11–18 km and a dissolved water content at the top of the chamber of 2.6–4.9 wt%.

Washington↗