SEARCH · Search USGS

Results for “Journal of Geophysics Research”

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

SourceResearch connectionFrom yearThrough year

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

At least 865 records · Page 48Linked to original sources

USGS2012

Dissolved organic carbon and chromophoric dissolved organic matter properties of rivers in the USA

Dissolved organic carbon (DOC) concentration and chromophoric dissolved organic matter (CDOM) parameters were measured over a range of discharge in 30 U.S. rivers, covering a diverse assortment of fluvial ecosystems in terms of watershed size and landscape drained. Relationships between CDOM absorption at a range of wavelengths (a 254 , a 350 , a 440 ) and DOC in the 30 watersheds were found to correlate strongly and positively for the majority of U.S. rivers. However, four rivers (Colorado, Colombia, Rio Grande and St. Lawrence) exhibited statistically weak relationships between CDOM absorption and DOC. These four rivers are atypical, as they either drain from the Great Lakes or experience significant impoundment of water within their watersheds, and they exhibited values for dissolved organic matter (DOM) parameters indicative of autochthonous or anthropogenic sources or photochemically degraded allochthonous DOM and thus a decoupling between CDOM and DOC. CDOM quality parameters in the 30 rivers were found to be strongly correlated to DOM compositional metrics derived via XAD fractionation, highlighting the potential for examining DOM biochemical quality from CDOM measurements. This study establishes the ability to derive DOC concentration from CDOM absorption for the majority of U.S. rivers, describes characteristics of riverine systems where such an approach is not valid, and emphasizes the possibility of examining DOM composition and thus biogeochemical function via CDOM parameters. Therefore, the usefulness of CDOM measurements, both laboratory-based analyses and in situ instrumentation, for improving spatial and temporal resolution of DOC fluxes and DOM dynamics in future studies is considerable in a range of biogeochemical studies.

Journal of Geophysical Research G: Biogeosciences↗

USGS2012

Uncovering the nonadiabatic response of geosynchronous electrons to geomagnetic disturbance

We describe an energy spectrum method for scaling electron integral flux, which is measured at a constant energy, to phase space density at a constant value of the first adiabatic invariant which removes much of the variation due to reversible adiabatic effects. Applying this method to nearly a solar cycle (1995 - 2006) of geosynchronous electron integral flux (E>2.0MeV) from the GOES satellites, we see that much of the diurnal variation in electron phase space density at constant energy can be removed by the transformation to phase space density at constant μ (4000 MeV/G). This allows us a clearer picture of underlying non-adiabatic electron population changes due to geomagnetic activity. Using scaled phase space density, we calculate the percentage of geomagnetic storms resulting in an increase, decrease or no change in geosynchronous electrons as 38%, 7%, and 55%, respectively. We also show examples of changes in the electron population that may be different than the unscaled fluxes alone suggest. These examples include sudden electron enhancements during storms which appear during the peak of negative Dst for μ-scaled phase space density, contrary to the slow increase seen during the recovery phase for unscaled phase space density for the same event.

Journal of Geophysical Research↗

USGS2012

Problem of the Love‐Gannon relation between the asymmetric disturbance field and Dst

Love and Gannon (2009) discovered that statistically, over a fifty year period the difference in the dawn and dusk disturbance‐field H component at low latitudes (hourly averaged) is linearly proportional to Dst. If the difference is designated by δ DD in units of nT/R E , then the Love‐Gannon (L‐G) relation is δ DD = − 0.2 Dst. At any time departures from the relation can be large. Nonetheless, the relation is evident for all values of Dst and persists throughout magnetic storms, both the main phase and the recovery phase. The Love‐Gannon discovery presents a problem to current understanding of the relation between the causes of δ DD and Dst because the dawn dusk asymmetry in the disturbance field is presumably governed by a long‐established magnetosphere‐ionosphere coupling theory which predicts a characteristic time scale (the shielding time) of less than an hour whereas the characteristic time scale for Dst (the ring current decay time) is more like ten hours. Thus, without forcing both time scales toward each other to the limits of their ranges, a linear proportionality between δ DD and Dst cannot be derived from the current understanding of the causes of the asymmetry and the ring current. This conclusion is the paper's main contribution. In addition, we attempt to get around the conflict of time scales by looking at other possibilities for generating δ DD that depend directly on the ring current. The most promising of these is the possibility that the ring current decay mechanism creates a quasi‐permanent, local‐time modification of the ring current compared to what it would be in the absence of the decay mechanism and that this modification causes a field‐aligned current that closes through the ionosphere and generates the asymmetry δ DD . This idea has the virtue of coupling the asymmetry directly to the ring current and of accounting for the persistence of the L‐G proportionality through the recovery phase of magnetic storms.

Journal of Geophysical Research A: Space Physics↗

USGS2012

Juan de Fuca slab geometry and its relation to Wadati-Benioff zone seismicity

A new model of the subducted Juan de Fuca plate beneath western North America allows first-order correlations between the occurrence of Wadati-Benioff zone earthquakes and slab geometry, temperature, and hydration state. The geo-referenced 3D model, constructed from weighted control points, integrates depth information from earthquake locations and regional seismic velocity studies. We use the model to separate earthquakes that occur in the Cascadia forearc from those that occur within the underlying Juan de Fuca plate and thereby reveal previously obscured details regarding the spatial distribution of earthquakes. Seismicity within the slab is most prevalent where the slab is warped beneath northwestern California and western Washington suggesting that slab flexure, in addition to expected metamorphic dehydration processes, promotes earthquake occurrence within the subducted oceanic plate. Earthquake patterns beneath western Vancouver Island are consistent with slab dehydration processes. Conversely, the lack of slab earthquakes beneath western Oregon is consistent with an anhydrous slab. Double-differenced relocated seismicity resolves a double seismic zone within the slab beneath northwestern California that strongly constrains the location of the plate interface and delineates a cluster of seismicity 10 km above the surface that includes the 1992 M7.1 Mendocino earthquake. We infer that this earthquake ruptured a surface within the Cascadia accretionary margin above the Juan de Fuca plate. We further speculate that this earthquake is associated with a detached fragment of former Farallon plate. Other subsurface tectonic elements within the forearc may have the potential to generate similar damaging earthquakes.

British Columbia;California;Oregon;Washington↗

USGS2012

Laboratory triggering of stick-slip events by oscillatory loading in the presence of pore fluid with implications for physics of tectonic tremor

The physical mechanism by which the low-frequency earthquakes (LFEs) that make up portions of tectonic (also called non-volcanic) tremor are created is poorly understood. In many areas of the world, tectonic tremor and LFEs appear to be strongly tidally modulated, whereas ordinary earthquakes are not. Anomalous seismic wave speeds, interpreted as high pore fluid pressure, have been observed in regions that generate tremor. Here we build upon previous laboratory studies that investigated the response of stick-slip on artificial faults to oscillatory, tide-like loading. These previous experiments were carried out using room-dry samples of Westerly granite, at one effective stress. Here we augment these results with new experiments on Westerly granite, with the addition of varying effective stress using pore fluid at two pressures. We find that raising pore pressure, thereby lowering effective stress can significantly increase the degree of correlation of stick-slip to oscillatory loading. We also find other pore fluid effects that become important at higher frequencies, when the period of oscillation is comparable to the diffusion time of pore fluid into the fault. These results help constrain the conditions at depth that give rise to tidally modulated LFEs, providing confirmation of the effective pressure law for triggering and insights into why tremor is tidally modulated while earthquakes are at best only weakly modulated.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Heat flow in vapor dominated areas of the Yellowstone Plateau volcanic field: implications for the thermal budget of the Yellowstone Caldera

Characterizing the vigor of magmatic activity in Yellowstone requires knowledge of the mechanisms and rates of heat transport between magma and the ground surface. We present results from a heat flow study in two vapor dominated, acid-sulfate thermal areas in the Yellowstone Caldera, the 0.11 km 2 Obsidian Pool Thermal Area (OPTA) and the 0.25 km 2 Solfatara Plateau Thermal Area (SPTA). Conductive heat flux through a low permeability layer capping large vapor reservoirs is calculated from soil temperature measurements at >600 locations and from laboratory measurements of soil properties. The conductive heat output is 3.6 ± 0.4 MW and 7.5 ± 0.4 MW from the OPTA and the SPTA, respectively. The advective heat output from soils is 1.3 ± 0.3 MW and 1.2 ± 0.3 MW from the OPTA and the SPTA, respectively and the heat output from thermal pools in the OPTA is 6.8 ± 1.4 MW. These estimates result in a total heat output of 11.8 ± 1.4 MW and 8.8 ± 0.4 MW from OPTA and SPTA, respectively. Focused zones of high heat flux in both thermal areas are roughly aligned with regional faults suggesting that faults in both areas serve as conduits for the rising acid vapor. Extrapolation of the average heat flux from the OPTA (103 ± 2 W·m −2 ) and SPTA (35 ± 3 W·m −2 ) to the ~35 km 2 of vapor dominated areas in Yellowstone yields 3.6 and 1.2 GW, respectively, which is less than the total heat output transported by steam from the Yellowstone Caldera as estimated by the chloride inventory method (4.0 to 8.0 GW).

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Externally triggered renewed bubble nucleation in basaltic magma: the 12 October 2008 eruption at Halema‘uma‘u Overlook vent, Kīlauea, Hawai‘i, USA

From October 2008 until present, dozens of small impulsive explosive eruptions occurred from the Overlook vent on the southeast side of Halema‘uma‘u Crater, at Kīlauea volcano, USA. These eruptions were triggered by rockfalls from the walls of the volcanic vent and conduit onto the top of the lava column. Here we use microtextural observations and data from clasts erupted during the well-characterized 12 October 2008 explosive eruption at Halema‘uma‘u to extend existing models of eruption triggering. We present a potential mechanism for this eruption by combining microtextural observations with existing geophysical and visual data sets. We measure the size and number density of bubbles preserved in juvenile ejecta using 2D images and X-ray microtomography. Our data suggest that accumulations of large bubbles with diameters of >50 μ m to at least millimeters existed at shallow levels within the conduit prior to the 12 October 2008 explosion. Furthermore, a high number density of small bubbles <50 μ m is measured in the clasts, implying very rapid nucleation of bubbles. Visual observations, combined with preexisting geophysical data, suggest that the impact of rockfalls onto the magma free surface induces pressure changes over short timescales that (1) nucleated new additional bubbles in the shallow conduit leading to high number densities of small bubbles and (2) expanded the preexisting bubbles driving upward acceleration. The trigger of eruption and bubble nucleation is thus external to the degassing system.

Hawaii↗

USGS2012

Hydrate morphology: Physical properties of sands with patchy hydrate saturation

The physical properties of gas hydrate-bearing sediments depend on the volume fraction and spatial distribution of the hydrate phase. The host sediment grain size and the state of effective stress determine the hydrate morphology in sediments; this information can be used to significantly constrain estimates of the physical properties of hydrate-bearing sediments, including the coarse-grained sands subjected to high effective stress that are of interest as potential energy resources. Reported data and physical analyses suggest hydrate-bearing sands contain a heterogeneous, patchy hydrate distribution, whereby zones with 100% pore-space hydrate saturation are embedded in hydrate-free sand. Accounting for patchy rather than homogeneous hydrate distribution yields more tightly constrained estimates of physical properties in hydrate-bearing sands and captures observed physical-property dependencies on hydrate saturation. For example, numerical modeling results of sands with patchy saturation agree with experimental observation, showing a transition in stiffness starting near the series bound at low hydrate saturations but moving toward the parallel bound at high hydrate saturations. The hydrate-patch size itself impacts the physical properties of hydrate-bearing sediments; for example, at constant hydrate saturation, we find that conductivity (electrical, hydraulic and thermal) increases as the number of hydrate-saturated patches increases. This increase reflects the larger number of conductive flow paths that exist in specimens with many small hydrate-saturated patches in comparison to specimens in which a few large hydrate saturated patches can block flow over a significant cross-section of the specimen.

Journal of Geophysical Research B: Solid Earth↗

USGS2013

Clustering of GPS velocities in the Mojave Block, southeastern California

We find subdivisions within the Mojave Block using cluster analysis to identify groupings in the velocities observed at GPS stations there. The clusters are represented on a fault map by symbols located at the positions of the GPS stations, each symbol representing the cluster to which the velocity of that GPS station belongs. Fault systems that separate the clusters are readily identified on such a map. The most significant representation as judged by the gap test involves 4 clusters within the Mojave Block. The fault systems bounding the clusters from east to west are 1) the faults defining the eastern boundary of the Northeast Mojave Domain extended southward to connect to the Hector Mine rupture, 2) the Calico-Paradise fault system, 3) the Landers-Blackwater fault system, and 4) the Helendale-Lockhart fault system. This division of the Mojave Block is very similar to that proposed by Meade and Hager. However, no cluster boundary coincides with the Garlock Fault, the northern boundary of the Mojave Block. Rather, the clusters appear to continue without interruption from the Mojave Block north into the southern Walker Lane Belt, similar to the continuity across the Garlock Fault of the shear zone along the Blackwater-Little Lake fault system observed by Peltzer et al. Mapped traces of individual faults in the Mojave Block terminate within the block and do not continue across the Garlock Fault [Dokka and Travis, ].

California↗

USGS2012

Coupled atmosphere-ocean-wave simulations of a storm event over the Gulf of Lion and Balearic Sea

The coastal areas of the North-Western Mediterranean Sea are one of the most challenging places for ocean forecasting. This region is exposed to severe storms events that are of short duration. During these events, significant air-sea interactions, strong winds and large sea-state can have catastrophic consequences in the coastal areas. To investigate these air-sea interactions and the oceanic response to such events, we implemented the Coupled Ocean-Atmosphere-Wave-Sediment Transport Modeling System simulating a severe storm in the Mediterranean Sea that occurred in May 2010. During this event, wind speed reached up to 25 m.s-1 inducing significant sea surface cooling (up to 2°C) over the Gulf of Lion (GoL) and along the storm track, and generating surface waves with a significant height of 6 m. It is shown that the event, associated with a cyclogenesis between the Balearic Islands and the GoL, is relatively well reproduced by the coupled system. A surface heat budget analysis showed that ocean vertical mixing was a major contributor to the cooling tendency along the storm track and in the GoL where turbulent heat fluxes also played an important role. Sensitivity experiments on the ocean-atmosphere coupling suggested that the coupled system is sensitive to the momentum flux parameterization as well as air-sea and air-wave coupling. Comparisons with available atmospheric and oceanic observations showed that the use of the fully coupled system provides the most skillful simulation, illustrating the benefit of using a fully coupled ocean-atmosphere-wave model for the assessment of these storm events.

Journal of Geophysical Research↗

USGS2012

Currents, drag, and sediment transport induced by a tsunami

We report observations of water surface elevation, currents, and suspended sediment concentration (SSC) from a 10-m deep site on the inner shelf in northern Monterey Bay during the arrival of the 2010 Chile tsunami. Velocity profiles were measured from 3.5 m above the bed (mab) to the surface at 2 min intervals, and from 0.1 to 0.7 mab at 1 Hz. SSC was determined from the acoustic backscatter of the near-bed profiler. The initial tsunami waves were directed cross shore and had a period of approximately 16 min. Maximum wave height was 1.1 m, and maximum current speed was 0.36 m/s. During the strongest onrush, near-bed velocities were clearly influenced by friction and a logarithmic boundary layer developed, extending more than 0.3 mab. We estimated friction velocity and bed shear stress from the logarithmic profiles. The logarithmic structure indicates that the flow can be characterized as quasi-steady at these times. At other phases of the tsunami waves, the magnitude of the acceleration term was significant in the near-bed momentum equation, indicating unsteady flow. The maximum tsunami-induced bed shear stress (0.4 N/m 2 ) exceeded the critical shear stress for the medium-grained sand on the seafloor. Cross-shore sediment flux was enhanced by the tsunami. Oscillations of water surface elevation and currents continued for several days. The oscillations were dominated by resonant frequencies, the most energetic of which was the fundamental longitudinal frequency of Monterey Bay. The maximum current speed (hourly-timescale) in 18 months of observations occurred four hours after the tsunami arrived.

California↗

USGS2012

Bathymetric controls on sediment transport in the Hudson River estuary: Lateral asymmetry and frontal trapping

Analyses of field observations and numerical model results have identified that sediment transport in the Hudson River estuary is laterally segregated between channel and shoals, features frontal trapping at multiple locations along the estuary, and varies significantly over the spring-neap tidal cycle. Lateral gradients in depth, and therefore baroclinic pressure gradient and stratification, control the lateral distribution of sediment transport. Within the saline estuary, sediment fluxes are strongly landward in the channel and seaward on the shoals. At multiple locations, bottom salinity fronts form at bathymetric transitions in width or depth. Sediment convergences near the fronts create local maxima in suspended-sediment concentration and deposition, providing a general mechanism for creation of secondary estuarine turbidity maxima at bathymetric transitions. The lateral bathymetry also affects the spring-neap cycle of sediment suspension and deposition. In regions with broad, shallow shoals, the shoals are erosional and the channel is depositional during neap tides, with the opposite pattern during spring tides. Narrower, deeper shoals are depositional during neaps and erosional during springs. In each case, the lateral transfer is from regions of higher to lower bed stress, and depends on the elevation of the pycnocline relative to the bed. Collectively, the results indicate that lateral and along-channel gradients in bathymetry and thus stratification, bed stress, and sediment flux lead to an unsteady, heterogeneous distribution of sediment transport and trapping along the estuary rather than trapping solely at a turbidity maximum at the limit of the salinity intrusion.

Journal of Geophysical Research↗

USGS2012

Along-shelf current variability on the Catalan inner-shelf (NW Mediterranean)

We examine the circulation over the inner shelf of the Catalan Sea using observations of currents obtained from three ADCPs within the inner-shelf (24 and 50 m depth) during March-April 2011. The along-shelf current fluctuations during that period are mainly controlled by the local wind stress on short time scales and by remote pressure gradients on synoptic time scales. Different forcing mechanisms are involved in the along-shelf momentum balance. During storm conditions, wind stress, sea level gradients and the non-linear terms dominate the balance. During weak wind conditions, the momentum balance is controlled by the pressure gradient, while during periods of moderate wind in the presence of considerable stratification, the balance is established between the Coriolis and wind stress terms. Vertical variations of velocity are affected by the strong observed density gradient. The increased vertical shear is accompanied by the development of stratified conditions due to local heating when the wind is not able to counteract (and destroy) stratification. The occasional influence of the Besòs river plume is observed in time scales of hours to days in a limited area in front of Barcelona. The area affected by the plume depends on the vertical extend of the fresher layer, the fast river discharge peak, and the relaxation of cross-shore velocities after northeast storm events. This contribution provides a first interpretation of the inner-shelf dynamics in the Catalan Sea.

Journal of Geophysical Research↗

USGS2012

Modeling thermal dynamics of active layer soils and near-surface permafrost using a fully coupled water and heat transport model

Thawing and freezing processes are key components in permafrost dynamics, and these processes play an important role in regulating the hydrological and carbon cycles in the northern high latitudes. In the present study, we apply a well-developed soil thermal model that fully couples heat and water transport, to simulate the thawing and freezing processes at daily time steps across multiple sites that vary with vegetation cover, disturbance history, and climate. The model performance was evaluated by comparing modeled and measured soil temperatures at different depths. We use the model to explore the influence of climate, fire disturbance, and topography (north- and south-facing slopes) on soil thermal dynamics. Modeled soil temperatures agree well with measured values for both boreal forest and tundra ecosystems at the site level. Combustion of organic-soil horizons during wildfire alters the surface energy balance and increases the downward heat flux through the soil profile, resulting in the warming and thawing of near-surface permafrost. A projection of 21st century permafrost dynamics indicates that as the climate warms, active layer thickness will likely increase to more than 3 meters in the boreal forest site and deeper than one meter in the tundra site. Results from this coupled heat-water modeling approach represent faster thaw rates than previously simulated in other studies. We conclude that the discussed soil thermal model is able to well simulate the permafrost dynamics and could be used as a tool to analyze the influence of climate change and wildfire disturbance on permafrost thawing.

Journal of Geophysical Research D: Atmospheres↗

USGS2012

A robust method to forecast volcanic ash clouds

Ash clouds emanating from volcanic eruption columns often form trails of ash extending thousands of kilometers through the Earth's atmosphere, disrupting air traffic and posing a significant hazard to air travel. To mitigate such hazards, the community charged with reducing flight risk must accurately assess risk of ash ingestion for any flight path and provide robust forecasts of volcanic ash dispersal. In response to this need, a number of different transport models have been developed for this purpose and applied to recent eruptions, providing a means to assess uncertainty in forecasts. Here we provide a framework for optimal forecasts and their uncertainties given any model and any observational data. This involves random sampling of the probability distributions of input (source) parameters to a transport model and iteratively running the model with different inputs, each time assessing the predictions that the model makes about ash dispersal by direct comparison with satellite data. The results of these comparisons are embodied in a likelihood function whose maximum corresponds to the minimum misfit between model output and observations. Bayes theorem is then used to determine a normalized posterior probability distribution and from that a forecast of future uncertainty in ash dispersal. The nature of ash clouds in heterogeneous wind fields creates a strong maximum likelihood estimate in which most of the probability is localized to narrow ranges of model source parameters. This property is used here to accelerate probability assessment, producing a method to rapidly generate a prediction of future ash concentrations and their distribution based upon assimilation of satellite data as well as model and data uncertainties. Applying this method to the recent eruption of Eyjafjallajökull in Iceland, we show that the 3 and 6 h forecasts of ash cloud location probability encompassed the location of observed satellite-determined ash cloud loads, providing an efficient means to assess all of the hazards associated with these ash clouds.

Journal of Geophysical Research D: Atmospheres↗

USGS2012

Improving the accuracy of S0 2 column densities and emission rates obtained from upward-looking UV-spectroscopic measurements of volcanic plumes by taking realistic radiative transfer into account

Sulfur dioxide (SO 2 ) is monitored using ultraviolet (UV) absorption spectroscopy at numerous volcanoes around the world due to its importance as a measure of volcanic activity and a tracer for other gaseous species. Recent studies have shown that failure to take realistic radiative transfer into account during the spectral retrieval of the collected data often leads to large errors in the calculated emission rates. Here, the framework for a new evaluation method which couples a radiative transfer model to the spectral retrieval is described. In it, absorption spectra are simulated, and atmospheric parameters are iteratively updated in the model until a best match to the measurement data is achieved. The evaluation algorithm is applied to two example Differential Optical Absorption Spectroscopy (DOAS) measurements conducted at Kilauea volcano (Hawaii). The resulting emission rates were 20 and 90% higher than those obtained with a conventional DOAS retrieval performed between 305 and 315 nm, respectively, depending on the different SO 2 and aerosol loads present in the volcanic plume. The internal consistency of the method was validated by measuring and modeling SO2 absorption features in a separate wavelength region around 375 nm and comparing the results. Although additional information about the measurement geometry and atmospheric conditions is needed in addition to the acquired spectral data, this method for the first time provides a means of taking realistic three-dimensional radiative transfer into account when analyzing UV-spectral absorption measurements of volcanic SO 2 plumes.

Journal of Geophysical Research D: Atmospheres↗

USGS2012

Spatial distribution of damage around faults in the Joe Lott Tuff Member of the Mount Belknap Volcanics, Utah: A mechanical analog for faulting in pyroclastic deposits on Mars

Volcanic ash is thought to comprise a large fraction of the Martian equatorial layered deposits and much new insight into the process of faulting and related fluid flow in these deposits can be gained through the study of analogous terrestrial tuffs. This study identifies a set of fault-related processes that are pertinent to understanding the evolution of fault systems in fine-grained, poorly indurated volcanic ash by investigating exposures of faults in the Miocene-aged Joe Lott Tuff Member of the Mount Belknap Volcanics, Utah. The porosity and granularity of the host rock are found to control the style of localized strain that occurs prior to and contemporaneous with faulting. Deformation bands occur in tuff that was porous and granular at the time of deformation, while fractures formed where the tuff lost its porous and granular nature due to silicic alteration. Non-localized deformation of the host rock is also prominent and occurs through compaction of void space, including crushing of pumice clasts. Significant off-fault damage of the host rock, resembling fault pulverization, is recognized adjacent to one analog fault and may reflect the strain rate dependence of the resulting fault zone architecture. These findings provide important new guidelines for future structural analyses and numerical modeling of faulting and subsurface fluid flow through volcanic ash deposits on Mars.

Utah↗

USGS2012

Soil genesis on the island of Bermuda in the Quaternary: the importance of African dust transport and deposition

The origin of terra rossa, red or reddish-brown, clay-rich soils overlying high-purity carbonate substrates, has intrigued geologists and pedologists for decades. Terra rossa soils can form from accumulation of insoluble residues during dissolution of the host limestones, addition of volcanic ash, or addition of externally derived, long-range-transported (LRT) aeolian particles. We studied soils and paleosols on high-purity, carbonate aeolianites of Quaternary age on Bermuda, where terra rossa origins have been debated for more than a century. Potential soil parent materials on this island include sand-sized fragments of local volcanic bedrock, the LRT, fine-grained (<20 μm) component of distal loess from the lower Mississippi River Valley, and LRT dust from Africa. These parent materials can be characterized geochemically using trace elements that are immobile in the soil-forming environment. Results indicate that local volcanic bedrock on Bermuda has Sc-Th-La, Cr-Ta-Nd, and Eu/Eu*, La N /Yb N , Gd N /Yb N that can be distinguished from African dust and lower Mississippi River valley loess. Bermuda soils have Sc-Th-La, Cr-Ta-Nd, and Eu/Eu*, La N /Yb N , Gd N /Yb N that indicate derivation from a combination of LRT dust from Africa and local volcanic bedrock. Our results indicate that soils on islands in a very broad latitudinal belt of the western Atlantic margin have been influenced by African LRT dust inputs over much of the past –500 ka.

Journal of Geophysical Research↗