USGS Science⌕ Search

SEARCH · USGS Science

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 523 records · Page 29Linked to original sources

Postseismic relaxation associated with transient creep rheology

[1] Perfettini and Avouac (2004) postulated that both the aftershock rate (assumed proportional to the local stressing rate) and the postseismic relaxation are driven by the loading imposed by postseismic slip on the brittle creep fault zone (BCFZ), the downdip extension of the fault zone below the coseismic rupture. I explore the consequences of that hypothesis for a long, strike-slip fault in the case where the BCFZ rheology is compatible with ordinary transient creep (creep strain proportional to log e (1 + t / τ 2 )). Because the important relaxation occurs near the bottom of the coseismic rupture, I calculate the postearthquake response with a model in which the BCFZ is represented by a viscoelastic half-space below the coseismic rupture. I find that both the predicted postseismic relaxation and the cumulative number of aftershocks can be approximated by the same temporal dependence N MO ( t ) = a MO (1−(1 + t / τ ) 1− p )/( p − 1), where t is the time after the earthquake and a MO , τ , and p are the constants chosen to fit either data set. Notice that d N MO ( t )/d t = ( a MO / τ )/(1 + t / τ ) p is the modified Omori law used to describe the rate of aftershock occurrence. Thus, the modified Omori law can be understood as a consequence of the Perfettini–Avouac hypothesis (aftershocks driven by slip on the BCFZ) and a BCFZ rheology compatible with ordinary transient creep. Moreover, the temporal dependence N MO ( t ) has been shown to fit postseismic surface deformation following at least 9 earthquakes. I also show that the conventional, one-dimensional, spring-block model of a BFCZ with a rheology compatible with ordinary transient creep leads to the same temporal dependence ( N MO ( t )).

Journal of Geophysical Research B: Solid Earth↗

Hydrothermal fluid flow and deformation in large calderas: Inferences from numerical simulations

[1] Inflation and deflation of large calderas is traditionally interpreted as being induced by volume change of a discrete source embedded in an elastic or viscoelastic half-space, though it has also been suggested that hydrothermal fluids may play a role. To test the latter hypothesis, we carry out numerical simulations of hydrothermal fluid flow and poroelastic deformation in calderas by coupling two numerical codes: (1) TOUGH2 [Pruess et al., 1999], which simulates flow in porous or fractured media, and (2) BIOT2 [Hsieh, 1996], which simulates fluid flow and deformation in a linearly elastic porous medium. In the simulations, high-temperature water (350°C) is injected at variable rates into a cylinder (radius 50 km, height 3–5 km). A sensitivity analysis indicates that small differences in the values of permeability and its anisotropy, the depth and rate of hydrothermal injection, and the values of the shear modulus may lead to significant variations in the magnitude, rate, and geometry of ground surface displacement, or uplift. Some of the simulated uplift rates are similar to observed uplift rates in large calderas, suggesting that the injection of aqueous fluids into the shallow crust may explain some of the deformation observed in calderas.

Journal of Geophysical Research B: Solid Earth↗

Three-dimensional geophysical mapping of rock alteration and water content at Mount Adams, Washington: Implications for lahar hazards

[1] Hydrothermally altered rocks, particularly if water saturated, can weaken stratovolcanoes, thereby increasing the potential for catastrophic sector collapses that can lead to far-traveled, destructive debris flows. Evaluating the hazards associated with such alteration is difficult because alteration has been mapped on few active volcanoes and the distribution and intensity of subsurface alteration are largely unknown on any active volcano. At Mount Adams, some Holocene debris flows contain abundant hydrothermal minerals derived from collapse of the altered edifice. Intense hydrothermal alteration significantly reduces the resistivity and magnetization of volcanic rock, and therefore hydrothermally altered rocks can be identified with helicopter electromagnetic and magnetic measurements. Electromagnetic and magnetic data, combined with geological mapping and rock property measurements, indicate the presence of appreciable thicknesses of hydrothermally altered rock in the central core of Mount Adams north of the summit. We identify steep cliffs at the western edge of this zone as the likely source for future large debris flows. In addition, the electromagnetic data identified water in the brecciated core of the upper 100–200 m of the volcano. Water helps alter the rocks, reduces the effective stress, thereby increasing the potential for slope failure, and acts, with entrained melting ice, as a lubricant to transform debris avalanches into lahars. Therefore knowing the distribution of water is also important for hazard assessments. Our results demonstrate that high-resolution geophysical and geological observations can yield unprecedented views of the three-dimensional distribution of altered rock and shallow pore water aiding evaluation of the debris avalanche hazard.

Journal of Geophysical Research B: Solid Earth↗

Integrating laboratory creep compaction data with numerical fault models: A Bayesian framework

[1] We developed a robust Bayesian inversion scheme to plan and analyze laboratory creep compaction experiments. We chose a simple creep law that features the main parameters of interest when trying to identify rate-controlling mechanisms from experimental data. By integrating the chosen creep law or an approximation thereof, one can use all the data, either simultaneously or in overlapping subsets, thus making more complete use of the experiment data and propagating statistical variations in the data through to the final rate constants. Despite the nonlinearity of the problem, with this technique one can retrieve accurate estimates of both the stress exponent and the activation energy, even when the porosity time series data are noisy. Whereas adding observation points and/or experiments reduces the uncertainty on all parameters, enlarging the range of temperature or effective stress significantly reduces the covariance between stress exponent and activation energy. We apply this methodology to hydrothermal creep compaction data on quartz to obtain a quantitative, semiempirical law for fault zone compaction in the interseismic period. Incorporating this law into a simple direct rupture model, we find marginal distributions of the time to failure that are robust with respect to errors in the initial fault zone porosity.

Journal of Geophysical Research B: Solid Earth↗

Constraining fault constitutive behavior with slip and stress heterogeneity

We study how enforcing self-consistency in the statistical properties of the preshear and postshear stress on a fault can be used to constrain fault constitutive behavior beyond that required to produce a desired spatial and temporal evolution of slip in a single event. We explore features of rupture dynamics that (1) lead to slip heterogeneity in earthquake ruptures and (2) maintain these conditions following rupture, so that the stress field is compatible with the generation of aftershocks and facilitates heterogeneous slip in subsequent events. Our three-dimensional fmite element simulations of magnitude 7 events on a vertical, planar strike-slip fault show that the conditions that lead to slip heterogeneity remain in place after large events when the dynamic stress drop (initial shear stress) and breakdown work (fracture energy) are spatially heterogeneous. In these models the breakdown work is on the order of MJ/m2, which is comparable to the radiated energy. These conditions producing slip heterogeneity also tend to produce narrower slip pulses independent of a slip rate dependence in the fault constitutive model. An alternative mechanism for generating these confined slip pulses appears to be fault constitutive models that have a stronger rate dependence, which also makes them difficult to implement in numerical models. We hypothesize that self-consistent ruptures could also be produced by very narrow slip pulses propagating in a self-sustaining heterogeneous stress field with breakdown work comparable to fracture energy estimates of kJ/M2. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

A global search inversion for earthquake kinematic rupture history: Application to the 2000 western Tottori, Japan earthquake

[1] We present a two-stage nonlinear technique to invert strong motions records and geodetic data to retrieve the rupture history of an earthquake on a finite fault. To account for the actual rupture complexity, the fault parameters are spatially variable peak slip velocity, slip direction, rupture time and risetime. The unknown parameters are given at the nodes of the subfaults, whereas the parameters within a subfault are allowed to vary through a bilinear interpolation of the nodal values. The forward modeling is performed with a discrete wave number technique, whose Green's functions include the complete response of the vertically varying Earth structure. During the first stage, an algorithm based on the heat-bath simulated annealing generates an ensemble of models that efficiently sample the good data-fitting regions of parameter space. In the second stage (appraisal), the algorithm performs a statistical analysis of the model ensemble and computes a weighted mean model and its standard deviation. This technique, rather than simply looking at the best model, extracts the most stable features of the earthquake rupture that are consistent with the data and gives an estimate of the variability of each model parameter. We present some synthetic tests to show the effectiveness of the method and its robustness to uncertainty of the adopted crustal model. Finally, we apply this inverse technique to the well recorded 2000 western Tottori, Japan, earthquake ( Mw 6.6); we confirm that the rupture process is characterized by large slip (3-4 m) at very shallow depths but, differently from previous studies, we imaged a new slip patch (2-2.5 m) located deeper, between 14 and 18 km depth.

Tottori↗

Reconstructed historical land cover and biophysical parameters for studies of land-atmosphere interactions within the eastern United States

Over the past 350 years, the eastern half of the United States experienced extensive land cover changes. These began with land clearing in the 1600s, continued with widespread deforestation, wetland drainage, and intensive land use by 1920, and then evolved to the present-day landscape of forest regrowth, intensive agriculture, urban expansion, and landscape fragmentation. Such changes alter biophysical properties that are key determinants of land-atmosphere interactions (water, energy, and carbon exchanges). To understand the potential implications of these land use transformations, we developed and analyzed 20-km land cover and biophysical parameter data sets for the eastern United States at 1650, 1850, 1920, and 1992 time slices. Our approach combined potential vegetation, county-level census data, soils data, resource statistics, a Landsat-derived land cover classification, and published historical information on land cover and land use. We reconstructed land use intensity maps for each time slice and characterized the land cover condition. We combined these land use data with a mutually consistent set of biophysical parameter classes, to characterize the historical diversity and distribution of land surface properties. Time series maps of land surface albedo, leaf area index, a deciduousness index, canopy height, surface roughness, and potential saturated soils in 1650, 1850, 1920, and 1992 illustrate the profound effects of land use change on biophysical properties of the land surface. Although much of the eastern forest has returned, the average biophysical parameters for recent landscapes remain markedly different from those of earlier periods. Understanding the consequences of these historical changes will require land-atmosphere interactions modeling experiments.

Journal of Geophysical Research D: Atmospheres↗

Nature and origin of the hematite-bearing plains of Terra Meridiani based on analyses of orbital and Mars Exploration rover data sets

The ∼5 km of traverses and observations completed by the Opportunity rover from Endurance crater to the Fruitbasket outcrop show that the Meridiani plains consist of sulfate‐rich sedimentary rocks that are largely covered by poorly‐sorted basaltic aeolian sands and a lag of granule‐sized hematitic concretions. Orbital reflectance spectra obtained by Mars Express OMEGA over this region are dominated by pyroxene, plagioclase feldspar, crystalline hematite (i.e., concretions), and nano‐phase iron oxide dust signatures, consistent with Pancam and Mini‐TES observations. Mössbauer Spectrometer observations indicate more olivine than observed with the other instruments, consistent with preferential optical obscuration of olivine features in mixtures with pyroxene and dust. Orbital data covering bright plains located several kilometers to the south of the landing site expose a smaller areal abundance of hematite, more dust, and a larger areal extent of outcrop compared to plains proximal to the landing site. Low‐albedo, low‐thermal‐inertia, windswept plains located several hundred kilometers to the south of the landing site are predicted from OMEGA data to have more hematite and fine‐grained olivine grains exposed as compared to the landing site. Low calcium pyroxene dominates spectral signatures from the cratered highlands to the south of Opportunity. A regional‐scale model is presented for the formation of the plains explored by Opportunity, based on a rising ground water table late in the Noachian Era that trapped and altered local materials and aeolian basaltic sands. Cessation of this aqueous process led to dominance of aeolian processes and formation of the current configuration of the plains.

Journal of Geophysical Research E: Planets↗

Active dust devils in Gusev crater, Mars: Observations from the Mars Exploration Rover Spirit

A full dust devil "season" was observed from Spirit from 10 March 2005 (sol 421, first active dust devil observed) to 12 December 2005 (sol 691, last dust devil seen); this corresponds to the period Ls 173.2?? to 339.5??, or the southern spring and summer on Mars. Thermal Emission Spectrometer data suggest a correlation between high surface temperatures and a positive thermal gradient with active dust devils in Gusev and that Spirit landed in the waning stages of a dust devil season as temperatures decreased. 533 active dust devils were observed, enabling new characterizations; they ranged in diameter from 2 to 276 m, with most in the range of 10-20 m in diameter, and occurred from about 0930 to 1630 hours local true solar time (with the maximum forming around 1300 hours) and a peak occurrence in southern late spring (Ls ??? 250??). Horizontal speeds of the dust devils ranged from <1 to 21 m/s, while vertical wind speeds within the dust devils ranged from 0.2 to 8.8 m/s. These data, when combined with estimates of the dust content within the dust devils, yield dust fluxes of 3.95 ?? 10-9 to 4.59-4 kg/m2/s. Analysis of the dust devil frequency distribution over the inferred dust devil zone within Gusev crater yields ???50 active dust devils/km2/sol, suggesting a dust loading into the atmosphere of ???19 kg/km2/sol. This value is less than one tenth the estimates by Cantor et al. (2001) for regional dust storms on Mars. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

Erosion rates at the Mars Exploration Rover landing sites and long-term climate change on Mars

Erosion rates derived from the Gusev cratered plains and the erosion of weak sulfates by saltating sand at Meridiani Planum are so slow that they argue that the present dry and desiccating environment has persisted since the Early Hesperian. In contrast, sedimentary rocks at Meridiani formed in the presence of groundwater and occasional surface water, and many Columbia Hills rocks at Gusev underwent aqueous alteration during the Late Noachian, approximately coeval with a wide variety of geomorphic indicators that indicate a wetter and likely warmer environment. Two-toned rocks, elevated ventifacts, and perched and undercut rocks indicate localized deflation of the Gusev plains and deposition of an equivalent amount of sediment into craters to form hollows, suggesting average erosion rates of ???0.03 nm/yr. Erosion of Hesperian craters, modification of Late Amazonian craters, and the concentration of hematite concretions in the soils of Meridiani yield slightly higher average erosion rates of 1-10 nm/yr in the Amazonian. These erosion rates are 2-5 orders of magnitude lower than the slowest continental denudation rates on Earth, indicating that liquid water was not an active erosional agent. Erosion rates for Meridiani just before deposition of the sulfate-rich sediments and other eroded Noachian areas are comparable with slow denudation rates on Earth that are dominated by liquid water. Available data suggest the climate change at the landing sites from wet and likely warm to dry and desiccating occurred sometime between the Late Noachian and the beginning of the Late Hesperian (3.7-3.5 Ga). Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

Evidence for montmorillonite or its compositional equivalent in Columbia Hills, Mars

During its exploration of the Columbia Hills, the Mars Exploration Rover "Spirit" encountered several similar samples that are distinctly different from Martian meteorites and known Gusev crater soils, rocks, and sediments. Occurring in a variety of contexts and locations, these "Independence class" samples are rough-textured, iron-poor (equivalent FeO ??? 4 wt%), have high Al/Si ratios, and often contain unexpectedly high concentrations of one or more minor or trace elements (including Cr, Ni, Cu, Sr, and Y). Apart from accessory minerals, the major component common to these samples has a compositional profile of major and minor elements which is similar to the smectite montmorillonite, implicating this mineral, or its compositional equivalent. Infrared thermal emission spectra do not indicate the presence of crystalline smectite. One of these samples was found spatially associated with a ferric sulfate-enriched soil horizon, possibly indicating a genetic relationship between these disparate types of materials. Compared to the nearby Wishstone and Watchtower class rocks, major aqueous alteration involving mineral dissolution and mobilization with consequent depletions of certain elements is implied for this setting and may be undetectable by remote sensing from orbit because of the small scale of the occurrences and obscuration by mantling with soil and dust.

Journal of Geophysical Research E: Planets↗

Spectrophotometric properties of materials observed by Pancam on the Mars Exploration Rovers: 2. Opportunity

The Panoramic Camera (Pancam) on the Mars Exploration Rover Opportunity acquired visible/near‐infrared multispectral observations of soils and rocks under varying viewing and illumination geometries that were modeled using radiative transfer theory to improve interpretations of the microphysical and surface scattering nature of materials in Meridiani Planum. Nearly 25,000 individual measurements were collected of rock and soil units identified by their color and morphologic properties over a wide range of phase angles (0–150°) at Eagle crater, in the surrounding plains, in Endurance crater, and in the plains between Endurance and Erebus craters through Sol 492. Corrections for diffuse skylight incorporated sky models based on observations of atmospheric opacity throughout the mission. Disparity maps created from Pancam stereo images allowed inclusion of local facet orientation estimates. Outcrop rocks overall exhibited the highest single scattering albedos (≤0.9 at 753 nm), and most spherule‐rich soils exhibited the lowest (≤0.6 at 753 nm). Macroscopic roughness among outcrop rocks varied but was typically larger than spherule‐rich soils. Data sets with sufficient phase angle coverage (resulting in well‐constrained Hapke parameters) suggested that models using single‐term and two‐term Henyey‐Greenstein phase functions exhibit a dominantly broad backscattering trend for most undisturbed spherule‐rich soils. Rover tracks and other compressed soils exhibited forward scattering, while outcrop rocks were intermediate in their scattering behaviors. Some phase functions exhibited wavelength‐dependent trends that may result from variations in thin deposits of airfall dust that occurred during the mission.

Journal of Geophysical Research E: Planets↗

Overview of the Opportunity Mars Exploration Rover Mission to Meridiani Planum: Eagle Crater to Purgatory Ripple

The Mars Exploration Rover Opportunity touched down at Meridiani Planum in January 2004 and since then has been conducting observations with the Athena science payload. The rover has traversed more than 5 km, carrying out the first outcrop‐scale investigation of sedimentary rocks on Mars. The rocks of Meridiani Planum are sandstones formed by eolian and aqueous reworking of sand grains that are composed of mixed fine‐grained siliciclastics and sulfates. The siliciclastic fraction was produced by chemical alteration of a precursor basalt. The sulfates are dominantly Mg‐sulfates and also include Ca‐sulfates and jarosite. The stratigraphic section observed to date is dominated by eolian bedforms, with subaqueous current ripples exposed near the top of the section. After deposition, interaction with groundwater produced a range of diagenetic features, notably the hematite‐rich concretions known as “blueberries.” The bedrock at Meridiani is highly friable and has undergone substantial erosion by wind‐transported basaltic sand. This sand, along with concretions and concretion fragments eroded from the rock, makes up a soil cover that thinly and discontinuously buries the bedrock. The soil surface exhibits both ancient and active wind ripples that record past and present wind directions. Loose rocks on the soil surface are rare and include both impact ejecta and meteorites. While Opportunity's results show that liquid water was once present at Meridiani Planum below and occasionally at the surface, the environmental conditions recorded were dominantly arid, acidic, and oxidizing and would have posed some significant challenges to the origin of life.

Journal of Geophysical Research E: Planets↗

Visible and near-infrared multispectral analysis of rocks at Meridiani Planum, Mars, by the Mars Exploration Rover Opportunity

Multispectral measurements in the visible and near infrared of rocks at Meridiani Planum by the Mars Exploration Rover Opportunity's Pancam are described. The Pancam multispectral data show that the outcrops of the Burns formation consist of two main spectral units which in stretched 673, 535, 432 nm color composites appear buff- and purple-colored. These units are referred to as the HFS and LFS spectral units based on higher and lower values of 482 to 535 nm slope. Spectral characteristics are consistent with the LFS outcrop consisting of less oxidized, and the HFS outcrop consisting of more oxidized, iron-bearing minerals. The LFS surfaces are not as common and appear, primarily, at the distal ends of outcrop layers and on steep, more massive surfaces, locations that are subject to greater eolian erosion. Consequently, the HFS surfaces are interpreted as a weathering rind. Further inherent spectral differences between layer's and between different outcrop map units, both untouched and patches abraded by the rover's Rock Abrasion Tool, are also described. Comparisons of the spectral parameters of the Meridiani outcrop with a set of laboratory reflectance measurements of Fe3+-bearing minerals show that the field of outcrop measurements plots near the fields of hematite, ferrihydrite, poorly crystalline goethite, and schwertmannite. Rind and fracture fill materials, observed intermittently at outcrop exposures, are intermediate in their spectral character between both the HFS and LFS spectral classes and other, less oxidized, surface materials (basaltic sands, spherules, and cobbles). Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

Dust deposition on the Mars Exploration Rover Panoramic Camera (Pancam) calibration targets

The Panoramic Camera (Pancam) on the Mars Exploration Rover mission has acquired in excess of 20,000 images of the Pancam calibration targets on the rovers. Analysis of this data set allows estimates of the rate of deposition and removal of aeolian dust on both rovers. During the first 150-170 sols there was gradual dust accumulation on the rovers but no evidence for dust removal. After that time there is ample evidence for both dust removal and dust deposition on both rover decks. We analyze data from early in both rover missions using a diffusive reflectance mixing model. Assuming a dust settling rate proportional to the atmospheric optical depth, we derive spectra of optically thick layers of airfall dust that are consistent with spectra from dusty regions on the Martian surface. Airfall dust reflectance at the Opportunity site appears greater than at the Spirit site, consistent with other observations. We estimate the optical depth of dust deposited on the Spirit calibration target by sol 150 to be 0.44 ?? 0.13. For Opportunity the value was 0.39 ?? 0.12. Assuming 80% pore space, we estimate that the dust layer grew at a rate of one grain diameter per ???100 sols on the Spirit calibration target. On Opportunity the rate was one grain diameter per ???125 sols. These numbers are consistent with dust deposition rates observed by Mars Pathfinder taking into account the lower atmospheric dust optical depth during the Mars Pathfinder mission. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

Coordinated analyses of orbital and spirit rover data to characterize surface materials on the cratered plains of Gusev Crater, Mars

Comparison of the Mars Exploration Rover Spirit's Pancam (0.4 to 1.0 ??m) and Mars Express Observatoire pour la Mineralogie l'Eau, les Glaces et l'Activite?? (OMEGA) (0.4 to 2.5 ??m) spectral reflectance data over Spirit's traverses shows that Gusev cratered plains are dominated by nanophase ferric-oxide-rich dust covering weakly altered basaltic sands. This interpretation is also consistent with both observations from OMEGA data covering plains beyond the traverse region and interpretations of data from the other payload instruments on the Spirit Rover. OMEGA observations of relatively low albedo regions where dust has presumably been stripped by dust devils show negative spectral reflectance slopes from 1.5 to 2.5 ??m and moderately masked spectral features which are indicative of olivine or pyroxene. High-albedo regions north and south of the Spirit landing site have flat spectral reflectance slopes and few spectral features, although all spectra have a nanophase ferric-oxide absorption edge between 0.4 and 0.75 ??m. Comparison of THEMIS-derived thermal inertia values with OMEGA-derived spectral parameters shows that although the dust cover can be optically thick (0.4 to 2.5 ??m wavelength region) in some areas, it is not thick enough (???1 cm) to mask the thermal inertia of the underlying substrate for areas included in this study. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

Sorted bed forms as self-organized patterns: 2. complex forcing scenarios

We employ a numerical model to study the development of sorted bed forms under a variety of hydrodynamic and sedimentary conditions. Results indicate that increased variability in wave height decreases the growth rate of the features and can potentially give rise to complicated, a priori unpredictable, behavior. This happens because the system responds to a change in wave characteristics by attempting to self-organize into a patterned seabed of different geometry and spacing. The new wavelength might not have enough time to emerge before a new change in wave characteristics occurs, leading to less regular seabed configurations. The new seabed configuration is also highly dependent on the preexisting morphology, which further limits the possibility of predicting future behavior. For the same reasons, variability in the mean current magnitude and direction slows down the growth of features and causes patterns to develop that differ from classical sorted bed forms. Spatial variability in grain size distribution and different types of net sediment aggradation/degradation can also result in the development of sorted bed forms characterized by a less regular shape. Numerical simulations qualitatively agree with observed geometry (spacing and height) of sorted bed forms. Also in agreement with observations is that at shallower depths, sorted bed forms are more likely to be affected by changes in the forcing conditions, which might also explain why, in shallow waters, sorted bed forms are described as ephemeral features. Finally, simulations indicate that the different sorted bed form shapes and patterns observed in the field might not necessarily be related to diverse physical mechanisms. Instead, variations in sorted bed form characteristics may result from variations in local hydrodynamic and/or sedimentary conditions.

Journal of Geophysical Research F: Earth Surface↗

Modeling englacial radar attenuation at Siple Dome, West Antarctica, using ice chemistry and temperature data

The radar reflectivity of an ice-sheet bed is a primary measurement for discriminating between thawed and frozen beds. Uncertainty in englacial radar attenuation and its spatial variation introduces corresponding uncertainty in estimates of basal reflectivity. Radar attenuation is proportional to ice conductivity, which depends on the concentrations of acid and sea-salt chloride and the temperature of the ice. We synthesize published conductivity measurements to specify an ice-conductivity model and find that some of the dielectric properties of ice at radar frequencies are not yet well constrained. Using depth profiles of ice-core chemistry and borehole temperature and an average of the experimental values for the dielectric properties, we calculate an attenuation rate profile for Siple Dome, West Antarctica. The depth-averaged modeled attenuation rate at Siple Dome (20.0 ?? 5.7 dB km-1) is somewhat lower than the value derived from radar profiles (25.3 ?? 1.1 dB km-1). Pending more experimental data on the dielectric properties of ice, we can match the modeled and radar-derived attenuation rates by an adjustment to the value for the pure ice conductivity that is within the range of reported values. Alternatively, using the pure ice dielectric properties derived from the most extensive single data set, the modeled depth-averaged attenuation rate is 24.0 ?? 2.2 dB km-1. This work shows how to calculate englacial radar attenuation using ice chemistry and temperature data and establishes a basis for mapping spatial variations in radar attenuation across an ice sheet. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗