Search USGSSearch

SEARCH · Search USGS

Results for “Geophysical Journal”

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 289 records · Page 16Linked to original sources

Earthquake-by-earthquake fold growth above the Puente Hills blind thrust fault, Los Angeles, California: Implications for fold kinematics and seismic hazard

Boreholes and high-resolution seismic reflection data collected across the forelimb growth triangle above the central segment of the Puente Hills thrust fault (PHT) beneath Los Angeles, California, provide a detailed record of incremental fold growth during large earthquakes on this major blind thrust fault. These data document fold growth within a discrete kink band that narrows upward from ???460 m at the base of the Quaternary section (200-250 m depth) to <150 m at 2.5 m depth, with most growth during the most recent folding event occurring within a zone only ???60 m wide. These observations, coupled with evidence from petroleum industry seismic reflection data, demonstrate that most (>82% at 250 m depth) folding and uplift occur within discrete kink bands, thereby enabling us to develop a paleoseismic history of the underlying blind thrust fault. The borehole data reveal that the youngest part of the growth triangle in the uppermost 20 m comprises three stratigraphically discrete growth intervals marked by southward thickening sedimentary strata that are separated by intervals in which sediments do not change thickness across the site. We interpret the intervals of growth as occurring after the formation of now-buried paleofold scarps during three large PHT earthquakes in the past 8 kyr. The intervening intervals of no growth record periods of structural quiescence and deposition at the regional, near-horizontal stream gradient at the study site. Minimum uplift in each of the scarp-forming events, which occurred at 0.2-2.2 ka (event Y), 3.0-6.3 ka (event X), and 6.6-8.1 ka (event W), ranged from ???1.1 to ???1.6 m, indicating minimum thrust displacements of ???2.5 to 4.5 m. Such large displacements are consistent with the occurrence of large-magnitude earthquakes (Mw > 7). Cumulative, minimum uplift in the past three events was 3.3 to 4.7 m, suggesting cumulative thrust displacement of ???7 to 10.5 m. These values yield a minimum Holocene slip rate for the PHT of ???0.9 to 1.6 mm/yr. The borehole and seismic reflection data demonstrate that dip within the kink band is acquired incrementally, such that older strata that have been deformed by more earthquakes dip more steeply than younger strata. Specifically, strata dip 0.4?? at 4 m depth, 0.7?? at 20 m depth, 8?? at 90 m, 16?? at 110 m, and 17?? at 200 m. Moreover, structural restorations of the borehole data show that the locus of active folding (the anticlinal active axial surface) does not extend to the surface in exactly the same location from earthquake to earthquake. Rather, that the axial surfaces migrate from earthquake to earthquake, reflecting a component of fold growth by kink band migration. The incremental acquisition of bed dip in the growth triangle may reflect some combination of fold growth by limb rotation in addition to kink band migration, possibly through a component of trishear or shear fault bend folding. Alternatively, the component of limb rotation may result from curved hinge fault bend folding, and/or the mechanical response of loosely consolidated granular sediments in the shallow subsurface to folding at depth. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth

Mars Exploration Rover Navigation Camera in-flight calibration

The Navigation Camera (Navcam) instruments on the Mars Exploration Rover (MER) spacecraft provide support for both tactical operations as well as scientific observations where color information is not necessary: large-scale morphology, atmospheric monitoring including cloud observations and dust devil movies, and context imaging for both the thermal emission spectrometer and the in situ instruments on the Instrument Deployment Device. The Navcams are a panchromatic stereoscopic imaging system built using identical charge-coupled device (CCD) detectors and nearly identical electronics boards as the other cameras on the MER spacecraft. Previous calibration efforts were primarily focused on providing a detailed geometric calibration in line with the principal function of the Navcams, to provide data for the MER navigation team. This paper provides a detailed description of a new Navcam calibration pipeline developed to provide an absolute radiometric calibration that we estimate to have an absolute accuracy of 10% and a relative precision of 2.5%. Our calibration pipeline includes steps to model and remove the bias offset, the dark current charge that accumulates in both the active and readout regions of the CCD, and the shutter smear. It also corrects pixel-to-pixel responsivity variations using flat-field images, and converts from raw instrument-corrected digital number values per second to units of radiance (W m-2 nm-1 sr-1), or to radiance factor (I/F). We also describe here the initial results of two applications where radiance-calibrated Navcam data provide unique information for surface photometric and atmospheric aerosol studies. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research E: Planets

Climatic effects of 30 years of landscape change over the Greater Phoenix, Arizona, region: 1. Surface energy budget changes

This paper is part 1 of a two-part study that evaluates the climatic effects of recent landscape change for one of the nation's most rapidly expanding metropolitan complexes, the Greater Phoenix, Arizona, region. The region's landscape evolution over an approximate 30-year period since the early 1970s is documented on the basis of analyses of Landsat images and land use/land cover (LULC) data sets derived from aerial photography (1973) and Landsat (1992 and 2001). High-resolution, Regional Atmospheric Modeling System (RAMS), simulations (2-km grid spacing) are used in conjunction with consistently defined land cover data sets and associated biophysical parameters for the circa 1973, circa 1992, and circa 2001 time periods to quantify the impacts of intensive land use changes on the July surface temperatures and the surface radiation and energy budgets for the Greater Phoenix region. The main findings are as follows: since the early 1970s the region's landscape has been altered by a significant increase in urban/suburban land area, primarily at the expense of decreasing plots of irrigated agriculture and secondarily by the conversion of seminatural shrubland. Mean regional temperatures for the circa 2001 landscape were 0.12??C warmer than the circa 1973 landscape, with maximum temperature differences, located over regions of greatest urbanization, in excess of 1??C. The significant reduction in irrigated agriculture, for the circa 2001 relative to the circa 1973 landscape, resulted in dew point temperature decreases in excess of 1??C. The effect of distinct land use conversion themes (e.g., conversion from irrigated agriculture to urban land) was also examined to evaluate how the most important conversion themes have each contributed to the region's changing climate. The two urbanization themes studied (from an initial landscape of irrigated agriculture and seminatural shrubland) have the greatest positive effect on near-surface temperature, increasing maximum daily temperatures by 1??C. Overall, sensible heat flux differences between the circa 2001 and circa 1973 landscapes result in a 1 W m -2 increase in domain-wide sensible heating, and a similar order of magnitude decrease in latent heating, highlighting the importance of surface repartitioning in establishing near-surface temperature trends. In part 2 of this study, we address the role of the surface budget changes on the mesoscale dynamics/thermodynamics, in context of the large-scale environment. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres

Transient and steady state creep response of ice I and magnesium sulfate hydrate eutectic aggregates

Using uniaxial compression creep experiments, we characterized the transient and steady state deformation behaviors of eutectic aggregates of system ice I and MgSO4 11H2O (MS11; meridianiite), which has significance because of its likely presence on moons of the outer solar system. Synthetic samples of eutectic liquid bulk composition, which produce eutectic colonies containing 0.35-0.50 volume fraction MS11, were tested as functions of colony size and lamellar spacing, temperature (230-250 K), and confining pressure (0.1 and 50 MPa) to strains ???0.2. Up to a differential stress of 6 MPa, the ice I-MS11 aggregates display an order of magnitude higher effective viscosity and higher stress sensitivity than do aggregates of pure polycrystalline ice at the same conditions. The creep data and associated microstructural observations demonstrate, however, that the aggregates are additionally more brittle than pure ice, approaching rate-independent plasticity that includes rupture of the hydrate phase at 6-8 MPa, depending on the scale of the microstructure. Microstructures of deformed samples reveal forms of semibrittle flow in which the hydrate phase fractures while the ice phase deforms plastically. Semibrittle flow in the icy shell of a planetary body would truncate the lithospheric strength envelope and thereby decrease the depth to the brittle-ductile transition by 55% and reduce the failure limit for compressional surface features from 10 to ???6 MPa. A constitutive equation that includes eutectic colony boundary sliding and intracolony flow is used to describe the steady state rheology of the eutectic aggregates. Copyright ?? 2011 by the American Geophysical Union.

Journal of Geophysical Research E: Planets

Flow of variably fluidized granular masses across three-dimensional terrain I. Coulomb mixture theory

Rock avalanches, debris flows, and related phenomena consist of grain-fluid mixtures that move across three-dimensional terrain. In all these phenomena the same basic forces, govern motion, but differing mixture compositions, initial conditions, and boundary conditions yield varied dynamics and deposits. To predict motion of diverse grain-fluid masses from initiation to deposition, we develop a depth-averaged, threedimensional mathematical model that accounts explicitly for solid- and fluid-phase forces and interactions. Model input consists of initial conditions, path topography, basal and internal friction angles of solid grains, viscosity of pore fluid, mixture density, and a mixture diffusivity that controls pore pressure dissipation. Because these properties are constrained by independent measurements, the model requires little or no calibration and yields readily testable predictions. In the limit of vanishing Coulomb friction due to persistent high fluid pressure the model equations describe motion of viscous floods, and in the limit of vanishing fluid stress they describe one-phase granular avalanches. Analysis of intermediate phenomena such as debris flows and pyroclastic flows requires use of the full mixture equations, which can simulate interaction of high-friction surge fronts with more-fluid debris that follows. Special numerical methods (described in the companion paper) are necessary to solve the full equations, but exact analytical solutions of simplified equations provide critical insight. An analytical solution for translational motion of a Coulomb mixture accelerating from rest and descending a uniform slope demonstrates that steady flow can occur only asymptotically. A solution for the asymptotic limit of steady flow in a rectangular channel explains why shear may be concentrated in narrow marginal bands that border a plug of translating debris. Solutions for static equilibrium of source areas describe conditions of incipient slope instability, and other static solutions show that nonuniform distributions of pore fluid pressure produce bluntly tapered vertical profiles at the margins of deposits. Simplified equations and solutions may apply in additional situations identified by a scaling analysis. Assessment of dimensionless scaling parameters also reveals that miniature laboratory experiments poorly simulate the dynamics of full-scale flows in which fluid effects are significant. Therefore large geophysical flows can exhibit dynamics not evident at laboratory scales.

Journal of Geophysical Research B: Solid Earth

Slicing up the San Francisco Bay Area: Block kinematics and fault slip rates from GPS-derived surface velocities

Observations of surface deformation allow us to determine the kinematics of faults in the San Francisco Bay Area. We present the Bay Area velocity unification (BA??VU??, "bay view"), a compilation of over 200 horizontal surface velocities computed from campaign-style and continuous Global Positioning System (GPS) observations from 1993 to 2003. We interpret this interseismic velocity field using a three-dimensional block model to determine the relative contributions of block motion, elastic strain accumulation, and shallow aseismic creep. The total relative motion between the Pacific plate and the rigid Sierra Nevada/Great Valley (SNGV) microplate is 37.9 ?? 0.6 mm yr-1 directed toward N30.4??W ?? 0.8?? at San Francisco (??2??). Fault slip rates from our preferred model are typically within the error bounds of geologic estimates but provide a better fit to geodetic data (notable right-lateral slip rates in mm yr-1: San Gregorio fault, 2.4 ?? 1.0; West Napa fault, 4.0 ?? 3.0; zone of faulting along the eastern margin of the Coast Range, 5.4 ?? 1.0; and Mount Diablo thrust, 3.9 ?? 1.0 of reverse slip and 4.0 ?? 0.2 of right-lateral strike slip). Slip on the northern Calaveras is partitioned between both the West Napa and Concord/ Green Valley fault systems. The total convergence across the Bay Area is negligible. Poles of rotation for Bay Area blocks progress systematically from the North America-Pacific to North America-SNGV poles. The resulting present-day relative motion cannot explain the strike of most Bay Area faults, but fault strike does loosely correlate with inferred plate motions at the time each fault initiated. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth

Debris-bed friction of hard-bedded glaciers

[1] Field measurements of debris-bed friction on a smooth rock tablet at the bed of Engabreen, a hard-bedded, temperate glacier in northern Norway, indicated that basal ice containing 10% debris by volume exerted local shear traction of up to 500 kPa. The corresponding bulk friction coefficient between the dirty basal ice and the tablet was between 0.05 and 0.08. A model of friction in which nonrotating spherical rock particles are held in frictional contact with the bed by bed-normal ice flow can account for these measurements if the power law exponent for ice flowing past large clasts is 1. A small exponent (n < 2) is likely because stresses in ice are small and flow is transient. Numerical calculations of the bed-normal drag force on a sphere in contact with a flat bed using n = 1 show that this force can reach values several hundred times that on a sphere isolated from the bed, thus drastically increasing frictional resistance. Various estimates of basal friction are obtained from this model. For example, the shear traction at the bed of a glacier sliding at 20 m a-1 with a geothermally induced melt rate of 0.006 m a-1 and an effective pressure of 300 kPa can exceed 100 kPa. Debris-bed friction can therefore be a major component of sliding resistance, contradicting the common assumption that debris-bed friction is negligible. Copyright 2005 by the American Geophysical Union.

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

Diffuse flow hydrothermal manganese mineralization along the active Mariana and southern Izu-Bonin arc system, western Pacific

Abundant ferromanganese oxides were collected along 1200 km of the active Izu-Bonin-Mariana arc system. Chemical compositions and mineralogy show that samples were collected from two deposit types: Fe-Mn crusts of mixed hydrogenetic/hydrothermal origin and hydrothermal Mn oxide deposits; this paper addresses only the second type. Mn oxides cement volcaniclastic and biogenic sandstone and breccia layers (Mn sandstone) and form discrete dense stratabound layers along bedding planes and within beds (stratabound Mn). The Mn oxide was deposited within coarse-grained sediments from diffuse flow systems where precipitation occurred below the seafloor. Deposits were exposed at the seabed by faulting, mass wasting, and erosion. Scanning electron microscopy and microprobe analyses indicate the presence of both amorphous and crystalline 10 ?? and 7 ?? manganate minerals, the fundamental chemical difference being high water contents in the amorphous Mn oxides. Alternation of amorphous and crystalline laminae occurs in many samples, which likely resulted from initial rapid precipitation of amorphous Mn oxides from waxing pulses of hydrothermal fluids followed by precipitation of slow forming crystallites during waning stages. The chemical composition is characteristic of a hydrothermal origin including strong fractionation between Fe (mean 0.9 wt %) and Mn (mean 48 wt %) for the stratabound Mn, generally low trace metal contents, and very low rare earth element and platinum group element contents. However, Mo, Cd, Zn, Cu, Ni, and Co occur in high concentrations in some samples and may be good indicator elements for proximity to the heat source or to massive sulfide deposits. For the Mn sandstones, Fe (mean-8.4%) and Mn (12.4%) are not significantly fractionated because of high Fe contents in the volcaniclastic material. However, the proportion of hydrothermal Fe (nondetrital Fe) to total Fe is remarkably constant (49-58%) for all the sample groups, regardless of the degree of Mn mineralization. Factor analyses indicate various mixtures of two dominant components: hydrothermal Mn oxide for the stratabound Mn and detrital aluminosilicate for the Mn-cemented sandstone; and two minor components, hydrothermal Fe oxyhydroxide and biocarbonate/biosilica. Our conceptual model shows that Mn mineralization was produced by hydrothermal convection cells within arc volcanoes and sedimentary prisms that occur along, the flanks and within calderas. The main source of hydrothermal fluid was seawater that penetrated through fractures, faults, and permeable volcanic edifices. The fluids were heated by magma, enriched in metals by leaching of basement rocks and sediments, and mixed with magmatic fluids and gases. Dikes and sills may have been another source of heat that drove small-scale circulation within sedimentary prisms. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth

Evaluating the effects of historical land cover change on summertime weather and climate in New Jersey: Land cover and surface energy budget changes

The 19th-century agrarian landscape of New Jersey (NJ) and the surrounding region has been extensively transformed to the present-day land cover by urbanization, reforestation, and localized areas of deforestation. This study used a mesoscale atmospheric numerical model to investigate the sensitivity of the warm season climate of NJ to these land cover changes. Reconstructed 1880s-era and present-day land cover data sets were used as surface boundary conditions for a set of simulations performed with the Regional Atmospheric Modeling System (RAMS). Three-member ensembles with historical and present-day land cover were compared to examine the sensitivity of surface air and dew point temperatures, rainfall, and the individual components of the surface energy budget to these land cover changes. Mean temperatures for the present-day landscape were 0.3-0.6??C warmer than for the historical landscape over a considerable portion of NJ and the surrounding region, with daily maximum temperatures at least 1.0??C warmer over some of the highly urbanized locations. Reforested regions, however, were slightly cooler. Dew point temperatures decreased by 0.3-0.6??C, suggesting drier, less humid near-surface air for the present-day landscape. Surface warming was generally associated with repartitioning of net radiation from latent to sensible heat flux, and conversely for cooling. While urbanization was accompanied by strong surface albedo decreases and increases in net shortwave radiation, reforestation and potential changes in forest composition have generally increased albedos and also enhanced landscape heterogeneity. The increased deciduousness of forests may have further reduced net downward longwave radiation. Copyright 2008 by the American Geophysical Union.

New Jersey

Equilibrium shoreline response of a high wave energy beach

Four years of beach elevation surveys at Ocean Beach, San Francisco, California, are used to extend an existing equilibrium shoreline change model, previously calibrated with fine sand and moderate energy waves, to medium sand and higher-energy waves. The shoreline, characterized as the cross-shore location of the mean high water contour, varied seasonally by between 30 and 60 m, depending on the alongshore location. The equilibrium shoreline change model relates the rate of horizontal shoreline displacement to the hourly wave energy E and the wave energy disequilibrium, the difference between E and the equilibrium wave energy that would cause no change in the present shoreline location. Values for the model shoreline response coefficients are tuned to fit the observations in 500 m alongshore segments and averaged over segments where the model has good skill and the estimated effects of neglected alongshore sediment transport are relatively small. Using these representative response coefficients for 0.3 mm sand from Ocean Beach and driving the model with much lower-energy winter waves observed at San Onofre Beach (also 0.3 mm sand) in southern California, qualitatively reproduces the small seasonal shoreline fluctuations at San Onofre. This consistency suggests that the shoreline model response coefficients depend on grain size and may be constant, and thus transportable, between sites with similar grain size and different wave climates. The calibrated model response coefficients predict that for equal fluctuations in wave energy, changes in shoreline location on a medium-grained (0.3 mm) beach are much smaller than on a previously studied fine-grained (0.2 mm) beach. Copyright ?? 2011 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans

InSAR observations of aseismic slip associated with an earthquake swarm in the Columbia River flood basalts

In 2009 a swarm of small shallow earthquakes occurred within the basalt flows of the Columbia River Basalt Group (CRBG). The swarm occurred within a dense seismic network in the U.S. Department of Energys Hanford Site. Data from the seismic network along with interferometric synthetic aperture radar (InSAR) data from the European Space Agencys (ESA) ENVISAT satellite provide insight into the nature of the swarm. By modeling the InSAR deformation data we constructed a model that consists of a shallow thrust fault and a near horizontal fault. We suggest that the near horizontal lying fault is a bedding-plane fault located between basalt flows. The geodetic moment of the modeled fault system is about eight times the cumulative seismic moment of the swarm. Precise location estimates of the swarm earthquakes indicate that the area of highest slip on the thrust fault, ???70mm of slip less than ???0.5km depth, was not located within the swarm cluster. Most of the slip on the faults appears to have progressed aseismically and we suggest that interbed sediments play a central role in the slip process. Copyright 2011 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth

Bora event variability and the role of air-sea feedback

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

Journal of Geophysical Research C: Oceans

Quantity-activity relationship of denitrifying bacteria and environmental scaling in streams of a forested watershed

The spatial variability of subreach denitrification rates in streams was evaluated with respect to controlling environmental conditions, molecular examination of denitrifying bacteria, and dimensional analysis. Denitrification activities ranged from 0 and 800 ng-N gsed-1 d-1 with large variations observed within short distances (<50 m) along stream reaches. A log-normal probability distribution described the range in denitrification activities and was used to define low (16% of the probability distributibn), medium (68%), and high (16%) denitrification potential groups. Denitrifying bacteria were quantified using a competitive polymerase chain reaction (cPCR) technique that amplified the nirK gene that encodes for nitrite reductase. Results showed a range of nirK quantities from 103 to 107 gene-copy-number gsed.-1 A nonparametric statistical test showed no significant difference in nirK quantifies among stream reaches, but revealed that samples with a high denitrification potential had significantly higher nirK quantities. Denitrification activity was positively correlated with nirK quantities with scatter in the data that can be attributed to varying environmental conditions along stream reaches. Dimensional analysis was used to evaluate denitrification activities according to environmental variables that describe fluid-flow properties, nitrate and organic material quantities, and dissolved oxygen flux. Buckingham's pi theorem was used to generate dimensionless groupings and field data were used to determine scaling parameters. The resulting expressions between dimensionless NO3- flux and dimensionless groupings of environmental variables showed consistent scaling, which indicates that the subreach variability in denitrification rates can be predicted by the controlling physical, chemical, and microbiological conditions. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences

Biogeochemical legacy of prescribed fire in a giant sequoia - Mixed conifer forest: A 16-year record of watershed balances

The effects of prescription burning on watershed balances of major ions in mixed conifer forest were examined in a 16-year paired catchment study in Sequoia National Park, California. The objective was to determine whether fire-related changes in watershed balances persist as long as estimated low-end natural fire-return intervals (???10 years), and whether cumulative net export caused by fire could deplete nutrient stocks between successive fires. Inputs (wet + dry deposition) and outputs (stream export) of N, S, Cl-, HCO3-, Ca2+, Mg2+, Na+, K+, H+, and SiO2 were measured for 7 years preceding, and 9 years following, a prescribed burn of one of the catchments. After fire, runoff coefficients increased by 7% (in dry years) to 35% (in wet years). Inorganic N was elevated in stream water for 3 years after fire. Increased export of water, SO42-, Cl-,SiO2, and base cations continued through the end of the study. Pools and processes attributed to fire led to the cumulative loss, per hectare, of 1.2 kg N, 16 kg S, 25 kg Cl-, 130 kg Ca2+, 19 kg Mg2+, 71 kg Na+, 29 kg K+ and 192 kg Si, above that predicted by prefire regression equations relating export in the paired catchments. This additional export equaled <1% of the N, up to one-third of the Ca and Mg, and up to three-fourths of the K, contained in the forest floor prior to combustion. Changes in watershed balances indicated that low-end natural fire-return intervals may prevent complete reaccumulation of several elements between fires. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences

Variational analysis of drifter positions and model outputs for the reconstruction of surface currents in the central Adriatic during fall 2002

In this paper we present an application of a variational method for the reconstruction of the velocity field in a coastal flow in the central Adriatic Sea, using in situ data from surface drifters and outputs from the ROMS circulation model. The variational approach, previously developed and tested for mesoscale open ocean flows, has been improved and adapted to account for inhomogeneities on boundary current dynamics over complex bathymetry and coastline and for weak Lagrangian persistence in coastal flows. The velocity reconstruction is performed using nine drifter trajectories over 45 d, and a hierarchy of indirect tests is introduced to evaluate the results as the real ocean state is not known. For internal consistency and impact of the analysis, three diagnostics characterizing the particle prediction and transport, in terms of residence times in various zones and export rates from the boundary current toward the interior, show that the reconstruction is quite effective. A qualitative comparison with sea color data from the MODIS satellite images shows that the reconstruction significantly improves the description of the boundary current with respect to the ROMS model first guess, capturing its main features and its exchanges with the interior when sampled by the drifters. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans

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

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

Journal of Geophysical Research F: Earth Surface

Sediment and nutrient delivery from thermokarst features in the foothills of the North Slope, Alaska: Potential impacts on headwater stream ecosystems

Permafrost is a defining characteristic of the Arctic environment. However, climate warming is thawing permafrost in many areas leading to failures in soil structure called thermokarst. An extensive survey of a 600 km2 area in and around the Toolik Lake Natural Research Area (TLNRA) revealed at least 34 thermokarst features, two thirds of which were new since ???1980 when a high resolution aerial survey of the area was done. Most of these thermokarst features were associated with headwater streams or lakes. We have measured significantly increased sediment and nutrient loading from thermokarst features to streams in two well-studied locations near the TLNRA. One small thermokarst gully that formed in 2003 on the Toolik River in a 0.9 km2 subcatchment delivered more sediment to the river than is normally delivered in 18 years from 132 km2 in the adjacent upper Kuparuk River basin (a long-term monitoring reference site). Ammonium, nitrate, and phosphate concentrations downstream from a thermokarst feature on Imnavait Creek increased significantly compared to upstream reference concentrations and the increased concentrations persisted over the period of sampling (1999-2005). The downstream concentrations were similar to those we have used in a long-term experimental manipulation of the Kuparuk River and that have significantly altered the structure and function of that river. A subsampling of other thermokarst features from the extensive regional survey showed that concentrations of ammonium, nitrate, and phosphate were always higher downstream of the thermokarst features. Our previous research has shown that even minor increases in nutrient loading stimulate primary and secondary production. However, increased sediment loading could interfere with benthic communities and change the responses to increased nutrient delivery. Although the terrestrial area impacted by thermokarsts is limited, the aquatic habitat altered by these failures can be extensive. If warming in the Arctic foothills accelerates thermokarst formation, there may be substantial and wide-spread impacts on arctic stream ecosystems that are currently poorly understood. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences