SEARCH · USGS Science

Results for “Geophysical Journal”

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

SourceResearch connectionFrom yearThrough year

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

At least 973 records · Page 54Linked to original sources

USGS2005

Generation and verification of theoretical rating curves in the Whitewater River basin, Kansas

[1] A new method for generating stage-discharge relations (rating curves) for geomorphically stable channels is presented and applied to two streams in the Whitewater River basin, Kansas. The approach converts measurements of stage into discharge using a fluid mechanically based model. The model does not use empirical roughness coefficients, such as Manning coefficients, but rather determines channel roughness from field measurements of the (1) channel geometry, (2) the physical roughness of the bed, banks, and floodplain, and (3) the vegetation density on the banks and floodplain. These measurements are used to calculate explicitly the drag on the small-scale topographic features on the boundary, the drag on the stems and branches of woody vegetation, and the friction on the bed, banks, and floodplain. The theoretical rating curves produced by the model for two study reaches, which are near U.S. Geological Survey (USGS) streamflow-gauging stations, are in good agreement with direct measurements of discharge made by the USGS. Our method has the potential of providing accurate estimates of stream flows less expensively than conventional gauging methods. In addition, the method can be used to obtain more accurate discharge estimates than conventional indirect methods for determining discharge, which are based on estimates of Manning's coefficient.

Journal of Geophysical Research F: Earth Surface↗

USGS2006

Effects of tributary debris on the longitudinal profile of the Colorado River in Grand Canyon

The Colorado River in Grand Canyon has long been known as a "rapids-and-pools" river, with the rapids owing their existence primarily to tributary debris flows. The debris flows deposit subaerial debris fans that constrict the channel laterally and, when they enter the river, raise the bed elevation. The rapids are short-wavelength (???0.1 to ???1 km), small-amplitude (??????5 m) convexities in the river's longitudinal profile, arising from the shallow gradient in the upstream pool and the steep gradient through the rapid itself. Analysis of the entire longitudinal profile through Grand Canyon reveals two long-wavelength (???100 km), large-amplitude (15-30 m) river profile convexities: the eastern canyon convexity between river mile (RM) 30 and RM 80 and the western canyon convexity between RM 150 and RM 250. Convexities of intermediate scale are also identified in the longitudinal profile. These longer-wavelength, larger-amplitude convexities have strong spatial correlations with high rates of debris flow occurrence, high densities of Holocene debris fans, the largest debris fans along the river, and alluvial thicknesses of 10 m or more. River profile convexities are unstable and require an active and powerful geologic process to maintain them, in this case the abundant, frequent, and voluminous Holocene debris flow activity in Grand Canyon. At all wavelengths the most likely cause for these river profile convexities is Holocene aggradation of the riverbed beneath them, driven by the coarse particles of tributary debris flows. Large enough debris flows will slow river flow for kilometers upstream, causing it to drop much of its suspended load. Integrated over time and all of the tributary point source contributions, this process will build short-wavelength convexities into long-wavelength convexities. For most if not all of the Holocene the Colorado River has been dissipating most of its energy in the rapids and expending the remainder in transporting fine sediment through Grand Canyon, with little or no regional incision of bedrock.

Journal of Geophysical Research F: Earth Surface↗

USGS2005

Regulation of landslide motion by dilatancy and pore pressure feedback

A new mathematical model clarifies how diverse styles and rates of landslide motion can result from regulation of Coulomb friction by dilation or contraction of water‐saturated basal shear zones. Normalization of the model equations shows that feedback due to coupling between landslide motion, shear zone volume change, and pore pressure change depends on a single dimensionless parameter α, which, in turn, depends on the dilatancy angle ψ and the intrinsic timescales for pore pressure generation and dissipation. If shear zone soil contracts during slope failure, then α < 0, and positive pore pressure feedback and runaway acceleration are inevitable. If the shear zone dilates, then α > 0, and negative feedback permits slow, steady landslide motion to occur while positive pore pressure is supplied by rain infiltration. Steady state slip velocities v 0 obey v 0 = −( K /ψ) p * e , where K is the hydraulic conductivity and p * e is the normalized (dimensionless) negative pore pressure generated by dilation. If rain infiltration and attendant pore pressure growth continue unabated, however, their influence ultimately overwhelms the stabilizing influence of negative p * e . Then, unbounded landslide acceleration occurs, accentuated by an instability that develops if ψ diminishes as landslide motion proceeds. Nonetheless, numerical solutions of the model equations show that slow, nearly steady motion of a clay‐rich landslide may persist for many months as a result of negative pore pressure feedback that regulates basal Coulomb friction. Similarly stabilized motion is less likely to occur in sand‐rich landslides that are characterized by weaker negative feedback.

Journal of Geophysical Research F: Earth Surface↗

USGS2005

Transient rheology of the upper mantle beneath central Alaska inferred from the crustal velocity field following the 2002 Denali earthquake

The M7.9 2002 Denali earthquake, Alaska, is one of the largest strike-slip earthquakes ever recorded. The postseismic GPS velocity field around the 300-km-long rupture is characterized by very rapid horizontal velocity up to ???300 mm/yr for the first 0.1 years and slower but still elevated horizontal velocity up to ???100 mm/yr for the succeeding 1.5 years. I find that the spatial and temporal pattern of the displacement field may be explained by a transient mantle rheology. Representing the regional upper mantle as a Burghers body, I infer steady state and transient viscosities of ??1 = 2.8 ?? 1018 Pa s and ??2 = 1.0 ?? 1017 Pa s, respectively, corresponding to material relaxation times of 1.3 and 0.05 years. The lower crustal viscosity is poorly constrained by the considered horizontal velocity field, and the quoted mantle viscosities assume a steady state lower crust viscosity that is 7??1. Systematic bias in predicted versus observed velocity vectors with respect to a fixed North America during the first 3-6 months following the earthquake is reduced when all velocity vectors are referred to a fixed site. This suggests that the post-Denali GPS time series for the first 1.63 years are shaped by a combination of a common mode noise source during the first 3-6 months plus viscoelastic relaxation controlled by a transient mantle rheology.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

The U-tube: A novel system for acquiring borehole fluid samples from a deep geologic CO2 sequestration experiment

A novel system has been deployed to obtain geochemical samples of water and gas, at in situ pressure, during a geologic CO 2 sequestration experiment conducted in the Frio brine aquifer in Liberty County, Texas. Project goals required high-frequency recovery of representative and uncontaminated aliquots of a rapidly changing two-phase fluid (supercritical CO 2 and brine) fluid from 1.5 km depth. The data sets collected, using both the liquid and gas portions of the downhole samples, provide insights into the coupled hydrogeochemical issues affecting CO 2 sequestration in brine-filled formations. While the basic premise underlying the U-tube sampler is not new, the system is unique because careful consideration was given to the processing of the recovered two-phase fluids. In particular, strain gauges mounted beneath the high-pressure surface sample cylinders measured the ratio of recovered brine to supercritical CO 2 . A quadrupole mass spectrometer provided real-time gas analysis for perfluorocarbon and noble gas tracers that were injected along with the CO 2 . The U-tube successfully acquired frequent samples, facilitating accurate delineation of the arrival of the CO 2 plume, and on-site analysis revealed rapid changes in geochemical conditions.

Journal of Geophysical Research Solid Earth↗

USGS2006

Seismic characteristics of central Brazil crust and upper mantle: A deep seismic refraction study

[1] A two‐dimensional model of the Brazilian central crust and upper mantle was obtained from the traveltime interpretation of deep seismic refraction data from the Porangatu and Cavalcante lines, each approximately 300 km long. When the lines were deployed, they overlapped by 50 km, forming an E‐W transect approximately 530 km long across the Tocantins Province and western São Francisco Craton. The Tocantins Province formed during the Neoproterozoic when the São Francisco, the Paranapanema, and the Amazon cratons collided, following the subduction of the former Goiás ocean basin. Average crustal V P and V P /V S ratios, Moho topography, and lateral discontinuities within crustal layers suggest that the crust beneath central Brazil can be associated with major geological domains recognized at the surface. The Moho is an irregular interface, between 36 and 44 km deep, that shows evidences of first‐order tectonic structures. The 8.05 and 8.23 km s −1 P wave velocities identify the upper mantle beneath the Porangatu and Cavalcante lines, respectively. The observed seismic features allow for the identification of (1) the crust has largely felsic composition in the studied region, (2) the absence of the mafic‐ultramafic root beneath the Goiás magmatic arc, and (3) block tectonics in the foreland fold‐and‐thrust belt of the northern Brasília Belt during the Neoproterozoic. Seismic data also suggested that the Bouguer gravimetric discontinuities are mainly compensated by differences in mass distribution within the lithospheric mantle. Finally, the Goiás‐Tocantins seismic belt can be interpreted as a natural seismic alignment related to the Neoproterozoic mantle domain.

Journal of Geophysical Research B: Solid Earth↗

USGS2006

V P and V S structure of the Yellowstone hot spot from teleseismic tomography: Evidence for an upper mantle plume

The movement of the lithosphere over a stationary mantle magmatic source, often thought to be a mantle plume, explains key features of the 16 Ma Yellowstone–Snake River Plain volcanic system. However, the seismic signature of a Yellowstone plume has remained elusive because of the lack of adequate data. We employ new teleseismic P and S wave traveltime data to develop tomographic images of the Yellowstone hot spot upper mantle. The teleseismic data were recorded with two temporary seismograph arrays deployed in a 500 km by 600 km area centered on Yellowstone. Additional data from nearby regional seismic networks were incorporated into the data set. The V P and V S models reveal a strong low-velocity anomaly from ∼50 to 200 km directly beneath the Yellowstone caldera and eastern Snake River Plain, as has been imaged in previous studies. Peak anomalies are −2.3% for V P and −5.5% for V S . A weaker, anomaly with a velocity perturbation of up to −1.0% V P and −2.5% V S continues to at least 400 km depth. This anomaly dips 30° from vertical, west-northwest to a location beneath the northern Rocky Mountains. We interpret the low-velocity body as a plume of upwelling hot, and possibly wet rock, from the mantle transition zone that promotes small-scale convection in the upper ∼200 km of the mantle and long-lived volcanism. A high-velocity anomaly, 1.2% V P and 1.9% V S , is located at ∼100 to 250 km depth southeast of Yellowstone and may represent a downwelling of colder, denser mantle material.

Idaho, Montana, Nevada, Oregon, Utah, Wyoming↗

USGS2006

Distribution, 14 C chronology, and paleomagnetism of latest Pleistocene and Holocene lava flows at Haleakala volcano, Island of Maui, Hawai'i: A revision of lava flow hazard zones

New mapping and 60 new radiocarbon ages define the age and distribution of latest Pleistocene and Holocene (past 13,000 years) lava flows at Haleakalā volcano, Island of Maui. Paleomagnetic directions were determined for 118 sites, of which 89 are in lava flows younger than 13,000 years. The paleomagnetic data, in conjunction with a reference paleosecular variation (PSV) curve for the Hawaiian Islands, are combined with our knowledge of age limitations based on stratigraphic control to refine age estimates for some of the undated lava flows. The resulting volumetric rate calculations indicate that within analytical error, the extrusion rate has remained nearly constant during the past 13,000 years, in the range 0.05–0.15 km 3 /kyr, only about half the long-term rate required to produce the postshield strata emplaced in the past &sim;1 Myr. Haleakalā's eruptive frequency is similar to that of Hualālai volcano on the Island of Hawai‘i, but its lava flows cover substantially less area per unit time. The reduced rates of lava coverage indicate a lower volcanic hazard than in similar zones at Hualālai.

Hawaii↗

USGS2007

Generation of fine hydromagmatic ash by growth and disintegration of glassy rinds

The deposits of mafic hydromagmatic eruptions are more fine grained and variable in vesicularity than dry magmatic deposits. Blocky, equant shapes of many hydromagmatic clasts also contrast with droplet, thread, and bubble wall morphology of dry magmatic fragments. Small (<∼180 μ m), blocky hydromagmatic pyroclasts have traditionally been interpreted to result from discrete vapor explosions, although such explosions tend to occur only under certain conditions. This paper considers a process of hydromagmatic ash formation that involves repeated growth and disintegration of glassy rinds on pyroclast surfaces as they deform within turbulent flows. This process, termed “turbulent shedding”, may occur during the expansion phase of vapor explosions or during turbulent but nonexplosive mixing of magma with water, steam, or water sprays. The occurrence of turbulent shedding and the resulting fragment sizes depend on the timescale for rind growth and the timescale between disturbances that remove or disintegrate glassy rinds. Turbulent shedding is directly observable in some small littoral jets at Kilauea. Calculations suggest that, in the presence of liquid water or water sprays, glassy rinds having a thickness of microns to millimeters should form in milliseconds to seconds. This is similar to the timescale between turbulent velocity fluctuations that can shred lava globules and remove such rinds. The fraction of a deposit consisting of fine ash should increase with the duration of this process: Large‐scale Surtseyan jets generate hundreds or thousands of shedding events; bubble bursts or tephra jets at Kilauea's coast may produce only a few.

Journal of Geophysical Research B: Solid Earth↗

USGS2006

Inference of postseismic deformation mechanisms of the 1923 Kanto earthquake

Coseismic slip associated with the M7.9, 1923 Kanto earthquake is fairly well understood, involving slip of up to 8 m along the Philippine Sea-Honshu interplate boundary under Sagami Bay and its onland extension. Postseismic deformation after the 1923 earthquake, however, is relatively poorly understood. We revisit the available deformation data in order to constrain possible mechanisms of postseismic deformation and to examine the consequences for associated stress changes in the surrounding crust. Data from two leveling lines and one tide gage station over the first 7-8 years postseismic period are of much greater amplitude than the corresponding expected interseismic deformation during the same period, making these data suitable for isolating the signal from postseismic deformation. We consider both viscoelastic models of asthenosphere relaxation and afterslip models. A distributed coseismic slip model presented by Pollitz et al. (2005), combined with prescribed parameters of a viscoelastic Earth model, yields predicted postseismic deformation that agrees with observed deformation on mainland Honshu from Tokyo to the Izu peninsula. Elsewhere (southern Miura peninsula; Boso peninsula), the considered viscoelastic models fail to predict observed deformation, and a model of ???1 in shallow afterslip in the offshore region south of the Boso peninsula, with equivalent moment magnitude Mw = 7.0, adequately accounts for the observed deformation. Using the distributed coseismic slip model, layered viscoelastic structure, and a model of interseismic strain accumulation, we evaluate the post-1923 stress evolution, including both the coseismic and accumulated postseismic stress changes and those stresses contributed by interseismic loading. We find that if account is made for the varying tectonic regime in the region, the occurrence of both immediate (first month) post-1923 crustal aftershocks as well as recent regional crustal seismicity is consistent with the predicted stress pattern. This suggests that the influence of the 1923 earthquake on regional seismicity is fairly predictable and has persisted for at least seven decades following the earthquake.

Journal of Geophysical Research B: Solid Earth↗

USGS2006

Tectonic stressing in California modeled from GPS observations

What happens in the crust as a result of geodetically observed secular motions? In this paper we find out by distorting a finite element model of California using GPS-derived displacements. A complex model was constructed using spatially varying crustal thickness, geothermal gradient, topography, and creeping faults. GPS velocity observations were interpolated and extrapolated across the model and boundary condition areas, and the model was loaded according to 5-year displacements. Results map highest differential stressing rates in a 200-km-wide band along the Pacific-North American plate boundary, coinciding with regions of greatest seismic energy release. Away from the plate boundary, GPS-derived crustal strain reduces modeled differential stress in some places, suggesting that some crustal motions are related to topographic collapse. Calculated stressing rates can be resolved onto fault planes: useful for addressing fault interactions and necessary for calculating earthquake advances or delays. As an example, I examine seismic quiescence on the Garlock fault despite a calculated minimum 0.1-0.4 MPa static stress increase from the 1857 M???7.8 Fort Tejon earthquake. Results from finite element modeling show very low to negative secular Coulomb stress growth on the Garlock fault, suggesting that the stress state may have been too low for large earthquake triggering. Thus the Garlock fault may only be stressed by San Andreas fault slip, a loading pattern that could explain its erratic rupture history.

Journal of Geophysical Research B: Solid Earth↗

USGS2006

Aftershock decay, productivity, and stress rates in Hawaii: Indicators of temperature and stress from magma sources

We examined dozens of aftershock sequences in Hawaii in terms of Gutenberg-Richter and modified Omori law parameters. We studied p, the rate of aftershock decay; A p , the aftershock productivity, defined as the observed divided by the expected number of aftershocks; and c, the time delay when aftershock rates begin to fall. We found that for earthquakes shallower than 20 km, p values >1.2 are near active magma centers. We associate this high decay rate with higher temperatures and faster stress relaxation near magma reservoirs. Deep earthquakes near Kilauea's inferred magma transport path show a range of p values, suggesting the absence of a large, deep magma reservoir. Aftershock productivity is >4.0 for flank earthquakes known to be triggered by intrusions but is normal (0.25 to 4.0) for isolated main shocks. We infer that continuing, post-main shock stress from the intrusion adds to the main shock's stress step and causes higher A p . High A p in other zones suggests less obvious intrusions and pulsing magma pressure near Kilauea's feeding conduit. We calculate stress rates and stress rate changes from pre-main shock and aftershock rates. Stress rate increased after many intrusions but decreased after large M7–8 earthquakes. Stress rates are highest in the seismically active volcano flanks and lowest in areas far from volcanic centers. We found sequences triggered by intrusions tend to have high A p , high (>0.10 day) c values, a stress rate increase, and sometimes a peak in aftershock rate hours after the main shock. We interpret these values as indicating continuing intrusive stress after the main shock.

Hawai'i↗

USGS2006

Dislocation pileup as a representation of strain accumulation on a strike-slip fault

The conventional model of strain accumulation on a vertical transform fault is a discrete screw dislocation in an elastic half-space with the Burgers vector of the dislocation increasing at the rate of relative plate motion. It would be more realistic to replace that discrete dislocation by a dislocation distribution, presumably a pileup in which the individual dislocations are in equilibrium. The length of the pileup depends upon the applied stress and the amount of slip that has occurred at depth. I argue here that the dislocation pileup (the transition on the fault from no slip to slip at the full plate rate) occupies a substantial portion of the lithosphere thickness. A discrete dislocation at an adjustable depth can reproduce the surface deformation profile predicted by a pileup so closely that it will be difficult to distinguish between the two models. The locking depth (dislocation depth) of that discrete dislocation approximation is substantially (???30%) larger than that (depth to top of the pileup) in the pileup model. Thus, in inverting surface deformation data using the discrete dislocation model, the locking depth in the model should not be interpreted as the true locking depth. Although dislocation pileup models should provide a good explanation of the surface deformation near the fault trace, that explanation may not be adequate at greater distances from the fault trace because approximating the expected horizontally distributed deformation at subcrustal depths by uniform slip concentrated on the fault is not justified.

Journal of Geophysical Research B: Solid Earth↗

USGS2006

A new 1649-1884 catalog of destructive earthquakes near Tokyo and implications for the long-term seismic process

In order to assess the long-term character of seismicity near Tokyo, we construct an intensity-based catalog of damaging earthquakes that struck the greater Tokyo area between 1649 and 1884. Models for 15 historical earthquakes are developed using calibrated intensity attenuation relations that quantitatively convey uncertainties in event location and magnitude, as well as their covariance. The historical catalog is most likely complete for earthquakes M ??? 6.7; the largest earthquake in the catalog is the 1703 M ??? 8.2 Genroku event. Seismicity rates from 80 years of instrumental records, which include the 1923 M = 7.9 Kanto shock, as well as interevent times estimated from the past ???7000 years of paleoseismic data, are combined with the historical catalog to define a frequency-magnitude distribution for 4.5 ??? M ??? 8.2, which is well described by a truncated Gutenberg-Richter relation with a b value of 0.96 and a maximum magnitude of 8.4. Large uncertainties associated with the intensity-based catalog are propagated by a Monte Carlo simulation to estimations of the scalar moment rate. The resulting best estimate of moment rate during 1649-2003 is 1.35 ?? 1026 dyn cm yr-1 with considerable uncertainty at the 1??, level: (-0.11, + 0.20) ?? 1026 dyn cm yr-1. Comparison with geodetic models of the interseismic deformation indicates that the geodetic moment accumulation and likely moment release rate are roughly balanced over the catalog period. This balance suggests that the extended catalog is representative of long-term seismic processes near Tokyo and so can be used to assess earthquake probabilities. The resulting Poisson (or time-averaged) 30-year probability for M ??? 7.9 earthquakes is 7-11%.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

Modeling and measuring the nocturnal drainage flow in a high-elevation, subalpine forest with complex terrain

The nocturnal drainage flow of air causes significant uncertainty in ecosystem CO 2 , H 2 O, and energy budgets determined with the eddy covariance measurement approach. In this study, we examined the magnitude, nature, and dynamics of the nocturnal drainage flow in a subalpine forest ecosystem with complex terrain. We used an experimental approach involving four towers, each with vertical profiling of wind speed to measure the magnitude of drainage flows and dynamics in their occurrence. We developed an analytical drainage flow model, constrained with measurements of canopy structure and SF 6 diffusion, to help us interpret the tower profile results. Model predictions were in good agreement with observed profiles of wind speed, leaf area density, and wind drag coefficient. Using theory, we showed that this one‐dimensional model is reduced to the widely used exponential wind profile model under conditions where vertical leaf area density and drag coefficient are uniformly distributed. We used the model for stability analysis, which predicted the presence of a very stable layer near the height of maximum leaf area density. This stable layer acts as a flow impediment, minimizing vertical dispersion between the subcanopy air space and the atmosphere above the canopy. The prediction is consistent with the results of SF 6 diffusion observations that showed minimal vertical dispersion of nighttime, subcanopy drainage flows. The stable within‐canopy air layer coincided with the height of maximum wake‐to‐shear production ratio. We concluded that nighttime drainage flows are restricted to a relatively shallow layer of air beneath the canopy, with little vertical mixing across a relatively long horizontal fetch. Insight into the horizontal and vertical structure of the drainage flow is crucial for understanding the magnitude and dynamics of the mean advective CO 2 flux that becomes significant during stable nighttime conditions and are typically missed during measurement of the turbulent CO 2 flux. The model and interpretation provided in this study should lead to research strategies for the measurement of these advective fluxes and their inclusion in the overall mass balance for CO 2 at this site with complex terrain.

Journal of Geophysical Research D: Atmospheres↗

USGS2006

Globally synchronous ice core volcanic tracers and abrupt cooling during the last glacial period

We perform a Monte Carlo pattern recognition analysis of the coincidence between three regional volcanic histories from ice coring of Greenland and Antarctica over the period 2 to 45 ka, using SO4 anomalies in Greenland and East Antarctica determined by continuous core chemistry, together with West Antarctic volcanic ash layers determined by remote optical borehole logging and core assays. We find that the Antarctic record of volcanism correlates with Glacial abrupt climate change at a 95% to >99.8% (???3??) significance level and that volcanic depositions at the three locations match at levels exceeding 3??, likely indicating that many common horizons represent single eruptive events which dispersed material world wide. These globally coincident volcanics were associated with abrupt cooling, often simultaneous with onsets or sudden intensifications of millennial cold periods. The striking agreement between sites implies that the consistency of current timescales obtained by isotopic and glaciological dating methods is better than estimated. Copyright 2006 by the American Geogphysical Union.

Journal of Geophysical Research D: Atmospheres↗

USGS2006

Atmospheric mercury emissions from substrates and fumaroles associated with three hydrothermal systems in the western United States

This paper quantifies atmospheric mercury (Hg) emissions from substrates and fumaroles associated with three hydrothermal systems: Lassen Volcanic Center, California (LVC); Yellowstone Caldera, Wyoming (YC); and Dixie Valley, Nevada (DV). Substrate Hg fluxes were measured using field chamber methods at thermal and nonthermal sites. The highest Hg fluxes (up to 541 ng m −2 h −1 ) were measured at thermal active areas. Fluxes from altered and unaltered nonthermal sites were <10 ng m −2 h −1 and were comparable to those measured in natural low-Hg background regions for YC and DV, and at LVC they were slightly elevated. Similarly, reactive gaseous mercury concentrations were higher in thermal active areas. Using a Geographic Information System framework for scaling, estimated area-average Hg emissions from substrates were 0.9–3.8 ng m −2 h −1 at LVC, 0.8–2.8 ng m −2 h −1 at YC, and 0.4–0.5 ng m −2 h −1 at DV. At LVC, nonthermal areas were the primary substrate source of atmospheric Hg (>98%). At YC, substrate Hg emissions were dominated (50 to 90%) by acidically altered thermal areas. Substrate emissions at DV were low and primarily from nonthermal areas (66% to 75%). Fumarole emissions at LVC (91–146 kg yr −1 ) and YC (0.18–1.6 kg yr −1 for Mud Volcano) were estimated by applying Hg:H 2 O and Hg:CO 2 ratios in hydrothermal gas samples to H 2 O and CO 2 emissions. Applying total area-average emissions from substrates and thermal features at LVC, YC, and DV to similar systems across the conterminous United States, yearly atmospheric Hg emissions from active hydrothermal systems are projected to be 1.3–2.1 Mg.

Journal of Geophysical Research D: Atmospheres↗

USGS2006

Characterization and petrologic interpretation of olivine-rich basalts at Gusev Crater, Mars

Rocks on the floor of Gusev crater are basalts of uniform composition and mineralogy. Olivine, the only mineral to have been identified or inferred from data by all instruments on the Spirit rover, is especially abundant in these rocks. These picritic basalts are similar in many respects to certain Martian meteorites (olivine-phyric shergottites). The olivine megacrysts in both have intermediate compositions, with modal abundances ranging up to 20-30%. Associated minerals in both include low-calcium and high-calcium pyroxenes, plagioclase of intermediate composition, iron-titanium-chromium oxides, and phosphate. These rocks also share minor element trends, reflected in their nickel-magnesium and chromium-magnesium ratios. Gusev basalts and shergottites appear to have formed from primitive magmas produced by melting an undepleted mantle at depth and erupted without significant fractionation. However, apparent differences between Gusev rocks and shergottites in their ages, plagioclase abundances, and volatile contents preclude direct correlation. Orbital determinations of global olivine distribution and compositions by thermal emission spectroscopy suggest that olivine-rich rocks may be widespread. Because weathering under acidic conditions preferentially attacks olivine and disguises such rocks beneath alteration rinds, picritic basalts formed from primitive magmas may even be a common component of the Martian crust formed during ancient and recent times.

Journal of Geophysical Research E: Planets↗