Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Geophysical Research - Oceans”

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 487 records · Page 27Linked to original sources
USGS2006

Gas slug ascent through changes in conduit diameter: Laboratory insights into a volcano-seismic source process in low-viscosity magmas

Seismic signals generated during the flow and degassing of low-viscosity magmas include long-period (LP) and very-long-period (VLP) events, whose sources are often attributed to dynamic fluid processes within the conduit. We present the results of laboratory experiments designed to investigate whether the passage of a gas slug through regions of changing conduit diameter could act as a suitable source mechanism. A vertical, liquid-filled glass tube featuring a concentric diameter change was used to provide canonical insights into potentially deep or shallow seismic sources. As gas slugs ascend the tube, we observe systematic pressure changes varying with slug size, liquid depth, tube diameter, and liquid viscosity. Gas slugs undergoing an abrupt flow pattern change upon entering a section of significantly increased tube diameter induce a transient pressure decrease in and above the flare and an associated pressure increase below it, which stimulates acoustic and inertial resonant oscillations. When the liquid flow is not dominantly controlled by viscosity, net vertical forces on the apparatus are also detected. The net force is a function of the magnitude of the pressure transients generated and the tube geometry, which dictates where, and hence when, the traveling pressure pulses can couple into the tube. In contrast to interpretations of related volcano-seismic data, where a single downward force is assumed to result from an upward acceleration of the center of mass in the conduit, our experiments suggest that significant downward forces can result from the rapid deceleration of relatively small volumes of downward-moving liquid. Copyright 2006 by the American Geophysical Union.

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

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 ∼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

Uplift and subsidence associated with the great Aceh-Andaman earthquake of 2004

Rupture of the Sunda megathrust on 26 December 2004 produced broad regions of uplift and subsidence. We define the pivot line separating these regions as a first step in defining the lateral extent and the downdip limit of rupture during that great Mw ??? 9.2 earthquake. In the region of the Andaman and Nicobar islands we rely exclusively on the interpretation of satellite imagery and a tidal model. At the southern limit of the great rupture we rely principally on field measurements of emerged coral microatolls. Uplift extends from the middle of Simeulue Island, Sumatra, at ??? 2.5??N, to Preparis Island, Myanmar (Burma), at ??? 14.9??N. Thus the rupture is ??? 1600 km long. The distance from the pivot line to the trench varies appreciably. The northern and western Andaman Islands rose, whereas the southern and eastern portion of the islands subsided. The Nicobar Islands and the west coast of Aceh province, Sumatra, subsided. Tilt at the southern end of the rupture is steep; the distance from 1.5 m of uplift to the pivot line is just 60 km. Our method of using satellite imagery to recognize changes in elevation relative to sea surface height and of using a tidal model to place quantitative bounds on coseismic uplift or subsidence is a novel approach that can be adapted to other forms of remote sensing and can be applied to other subduction zones in tropical regions. Copyright 2006 by the American Geophysical Union.

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↗
USGS2006

Analysis of gravity and topography in the GLIMPSE study region: Isostatic compensation and uplift of the Sojourn and Hotu Matua Ridge systems

The Gravity Lieations Intraplate Melting Petrologic and Seismic Expedition (GLIMPSE) Experiment investigated the formation of a series of non-hot spot, intraplate volcanic ridges in the South Pacific and their relationship to cross-grain gravity lineaments detected by satellite altimetry. Using shipboard gravity measurements and a simple model of surface loading of a thin elastic plate, we estimate effective elastic thicknesses ranging from ???2 km beneath the Sojourn Ridge to a maximum of 10 km beneath the Southern Cross Seamount. These elastic thicknesses are lower than predicted for the 3-9 Ma seafloor on which the volcanoes lie, perhaps due to reheating and thinning of the plate during emplacement. Anomalously low apparent densities estimated for the Matua and Southern Cross seamounts 2050 and 2250 kg m-3, respectively, probably are artifacts caused by the assumption of only surface loading, ignoring the presence of subsurface loading in the form of underplated crust and/or low-density mantle. Using satellite free-air gravity and shipboard bathymetry, we calculate the age-detrended, residual mantle Bouguer anomaly (rMBA). The rMBA corrects the free-air anomaly for the direct effects of topography, including the thickening of the crust beneath the seamounts and volcanic ridges due to surface loading of the volcanic edifices. There are broad, negative rMBA anomalies along the Sojourn and Brown ridges and the Hotu Matua seamount chain that extend nearly to the East Pacific Rise. These negative rMBA anomalies connect to negative free-air anomalies in the western part of the study area that have been recognized previously as the beginnings of the cross-grain gravity lineaments. Subtracting the topographic effects of surface loading by the ridges and seamounts from the observed topography reveals that the ridges are built on broad bands of anomalously elevated seafloor. This swell topography and the negative rMBA anomalies contradict the predictions of lithospheric cracking models for the origin of gravity lineaments and associated volcanic ridges, favoring models with a dynamic mantle component such as small-scale convection or channelized asthenospheric return flow. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Imaging the transition from Aleutian subduction to Yakutat collision in central Alaska, with local earthquakes and active source data

In southern and central Alaska the subduction and active volcanism of the Aleutian subduction zone give way to a broad plate boundary zone with mountain building and strike-slip faulting, where the Yakutat terrane joins the subducting Pacific plate. The interplay of these tectonic elements can be best understood by considering the entire region in three dimensions. We image three-dimensional seismic velocity using abundant local earthquakes, supplemented by active source data. Crustal low-velocity correlates with basins. The Denali fault zone is a dominant feature with a change in crustal thickness across the fault. A relatively high-velocity subducted slab and a low-velocity mantle wedge are observed, and high Vp/Vs beneath the active volcanic systems, which indicates focusing of partial melt. North of Cook Inlet, the subducted Yakutat slab is characterized by a thick low-velocity, high-Vp/Vs, crust. High-velocity material above the Yakutat slab may represent a residual older slab, which inhibits vertical flow of Yakutat subduction fluids. Alternate lateral flow allows Yakutat subduction fluids to contribute to Cook Inlet volcanism and the Wrangell volcanic field. The apparent northeast edge of the subducted Yakutat slab is southwest of the Wrangell volcanics, which have adakitic composition consistent with melting of this Yakutat slab edge. In the mantle, the Yakutat slab is subducting with the Pacific plate, while at shallower depths the Yakutat slab overthrusts the shallow Pacific plate along the Transition fault. This region of crustal doubling within the shallow slab is associated with extremely strong plate coupling and the primary asperity of the Mw 9.2 great 1964 earthquake. Copyright 2006 by the American Geophysical Union.

Alaska↗
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

Cloud-to-ground lightning and surface rainfall in warm-season Florida thunderstorms

Relationships between cloud-to-ground (CG) lightning and surface rainfall have been examined in nine isolated, warm-season thunderstorms on the east coast of central Florida. CG flashes and the associated rain volumes were measured as a function of time in storm-centered reference frames that followed each storm over a network of rain gauges. Values of the storm-average rain volume per CG flash ranged from 0.70 ?? 104 to 6.4 ?? 104 m3/CG flash, with a mean (and standard deviation) of 2.6 ?? 104 ?? 2.1 ?? 104 m3/CG flash. Values of the rain volume concurrent with CG flashes ranged from 0.11 ?? 104 to 4.9 ?? 104 m3/CG flash with a mean of 2.1 ?? 104 ?? 2.0 ?? 104 m3/CG flash. The lag-time between the peak CG flash rate and the peak rainfall rate (using 5 min bins), and the results of a lag correlation analysis, show that surface rainfall tends to follow the lightning (positive lag) by up to 20 min in six storms. In one storm the rainfall preceded the lightning by 5 min, and two storms had nonsignificant lags. Values of the lagged rain volume concurrent with CG flashes ranged from 0.43 ?? 104 to 4.9 ?? 104 m3/CG flash, and the mean was 1.9 ?? 104 ?? 1.7 ?? 104 m3/CG flash. For the five storms that produced 12 or more flashes and had significant lags, a plot of the optimum lag time versus the total number of CG flashes shows a linear trend (R2 = 0.56). The number of storms is limited, but the lag results do indicate that large storms tend to have longer lags. A linear fit to the lagged rain volume vs. the number of concurrent CG flashes has a slope of 1.9 ?? 104 m3/CG flash (R2 = 0.83). We conclude that warm-season Florida thunderstorms produce a roughly constant rain volume per CG flash and that CG lightning can be used to estimate the location and intensity of convective rainfall in that weather regime. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗
USGS2006

Crater gradation in Gusev crater and Meridiani Planum, Mars

The Mars Exploration Rovers investigated numerous craters in Gusev crater and Meridiani Planum during the first ???400 sols of their missions. Craters vary in size and preservation state but are mostly due to secondary impacts at Gusev and primary impacts at Meridiani. Craters at both locations are modified primarily by eolian erosion and infilling and lack evidence for modification by aqueous processes. Effects of gradation on crater form are dependent on size, local lithology, slopes, and availability of mobile sediments. At Gusev, impacts into basaltic rubble create shallow craters and ejecta composed of resistant rocks. Ejecta initially experience eolian stripping, which becomes weathering-limited as lags develop on ejecta surfaces and sediments are trapped within craters. Subsequent eolian gradation depends on the slow production of fines by weathering and impacts and is accompanied by minor mass wasting. At Meridiani the sulfate-rich bedrock is more susceptible to eolian erosion, and exposed crater rims, walls, and ejecta are eroded, while lower interiors and low-relief surfaces are increasingly infilled and buried by mostly basaltic sediments. Eolian processes outpace early mass wasting, often produce meters of erosion, and mantle some surfaces. Some small craters were likely completely eroded/buried. Craters >100 m in diameter on the Hesperian-aged floor of Gusev are generally more pristine than on the Amazonian-aged Meridiani plains. This conclusion contradicts interpretations from orbital views, which do not readily distinguish crater gradation state at Meridiani and reveal apparently subdued crater forms at Gusev that may suggest more gradation than has occurred. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2005

Flow strength of highly hydrated Mg- and Na-sulfate hydrate salts, pure and in mixtures with water ice, with application to Europa

We selected two Europan-ice-shell candidate highly hydrated sulfate salts for a laboratory survey of ductile flow properties: MgSO4 ?? 7H2O (epsomite) and Na2SO4 ?? 10H2O (mirabilite), called MS7 and NS10, respectively. Polycrystalline samples in pure form and in mixtures with water ice I were tested using our cryogenic high-pressure creep apparatus at temperatures 232 ??? T ??? 294 K, confining pressures P = 50 and 100 MPa, and strain rates 4 ?? 10-8 ??? ???dot;e ??? 7 ?? 10-5 s-1. Grain size of NS10 samples was > 100 ??m. The flow strength ?? of pure MS7 was over 100 times that of polycrystalline ice I at comparable conditions; that of pure NS10 over 20 times that of ice. In terms of the creep law ???dot;e = A??n e-Q/RT, where R is the gas constant, we determine parameter values of A = 1012.1 MPa-ns-1, n = 5.4, and Q = 128 kJ/mol for pure NS10. Composites of ice I and NS10 of volume fraction ?? NS10 have flow strength ??c = [??NS10??NS10J + (1 - ?? NS10)??iceIJ]1/J where J ??? -0.5, making the effect on the flow of ice with low volume fractions of NS10 much like that of virtually undeformable hard rock inclusions. Being much stronger and denser than ice, massive sulfate inclusions in the warmer, ductile layer of the Europan ice shell are less likely to be entrained in convective ice flow and more likely to be drawn to the base of the ice shell by gravitational forces and eventually expelled. With only smaller, dispersed sulfate inclusions, at probable sulfate ?? < 0.2, the shell may be treated rheologically as pure, polycrystalline ice, with boundary conditions perhaps influenced by the high density and low thermal conductivity of the hydrated salts. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
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↗