Search USGS⌕ Search

SEARCH · Search USGS

Results for “Earth Surface Dynamics”

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 181 records · Page 10Linked to original sources
USGS2013

Runoff-generated debris flows: observations and modeling of surge initiation, magnitude, and frequency

Runoff during intense rainstorms plays a major role in generating debris flows in many alpine areas and burned steeplands. Yet compared to debris flow initiation from shallow landslides, the mechanics by which runoff generates a debris flow are less understood. To better understand debris flow initiation by surface water runoff, we monitored flow stage and rainfall associated with debris flows in the headwaters of two small catchments: a bedrock-dominated alpine basin in central Colorado (0.06 km 2 ) and a recently burned area in southern California (0.01 km 2 ). We also obtained video footage of debris flow initiation and flow dynamics from three cameras at the Colorado site. Stage observations at both sites display distinct patterns in debris flow surge characteristics relative to rainfall intensity (I). We observe small, quasiperiodic surges at low I; large, quasiperiodic surges at intermediate I; and a single large surge followed by small-amplitude fluctuations about a more steady high flow at high I. Video observations of surge formation lead us to the hypothesis that these flow patterns are controlled by upstream variations in channel slope, in which low-gradient sections act as “sediment capacitors,” temporarily storing incoming bed load transported by water flow and periodically releasing the accumulated sediment as a debris flow surge. To explore this hypothesis, we develop a simple one-dimensional morphodynamic model of a sediment capacitor that consists of a system of coupled equations for water flow, bed load transport, slope stability, and mass flow. This model reproduces the essential patterns in surge magnitude and frequency with rainfall intensity observed at the two field sites and provides a new framework for predicting the runoff threshold for debris flow initiation in a burned or alpine setting.

California;Colorado↗
USGS2018

Using stereo satellite imagery to account for ablation, entrainment, and compaction in volume calculations for rock avalanches on Glaciers: Application to the 2016 Lamplugh Rock Avalanche in Glacier Bay National Park, Alaska

The use of preevent and postevent digital elevation models (DEMs) to estimate the volume of rock avalanches on glaciers is complicated by ablation of ice before and after the rock avalanche, scour of material during rock avalanche emplacement, and postevent ablation and compaction of the rock avalanche deposit. We present a model to account for these processes in volume estimates of rock avalanches on glaciers. We applied our model by calculating the volume of the 28 June 2016 Lamplugh rock avalanche in Glacier Bay National Park, Alaska. We derived preevent and postevent 2‐m resolution DEMs from WorldView satellite stereo imagery. Using data from DEM differencing, we reconstructed the rock avalanche and adjacent surfaces at the time of occurrence by accounting for elevation changes due to ablation and scour of the ice surface, and postevent deposit changes. We accounted for uncertainties in our DEMs through precise coregistration and an assessment of relative elevation accuracy in bedrock control areas. The rock avalanche initially displaced 51.7 ± 1.5 Mm3 of intact rock and then scoured and entrained 13.2 ± 2.2 Mm3 of snow and ice during emplacement. We calculated the total deposit volume to be 69.9 ± 7.9 Mm3. Volume estimates that did not account for topographic changes due to ablation, scour, and compaction underestimated the deposit volume by 31.0–46.8 Mm3. Our model provides an improved framework for estimating uncertainties affecting rock avalanche volume measurements in glacial environments. These improvements can contribute to advances in the understanding of rock avalanche hazards and dynamics.

Alaska↗
USGS2015

Crustal stress and structure at Kīlauea Volcano inferred from seismic anisotropy

Seismic anisotropy, measured through shear wave splitting (SWS) analysis, can be indicative of the state of stress in Earth's crust. Changes in SWS at Kīlauea Volcano, Hawai‘i, associated with the onset of summit eruptive activity in 2008 hint at the potential of the technique for tracking volcanic activity. To use SWS observations as a monitoring tool, however, it is important to understand the cause of seismic anisotropy at the volcano throughout the eruptive cycle. To address this need, we analyzed SWS results from across Kīlauea in combination with macroscopic surface structures (mapped fractures, faults, and fissures) and stress orientations inferred from fault plane solutions. Seismic anisotropy seems to be due to pervasive aligned structures in most regions of the volcano. The upper East and Southwest Rift Zones, however, show a bimodality in stress and SWS, suggesting a stress discontinuity with depth, perhaps related to magma conduits that trend obliquely to the dominant structure. Other areas in and around Kīlauea Caldera display principal stresses of similar magnitudes, indicating that small stress perturbations can rotate the maximum horizontal compressive stress direction by up to 90 ° . In these locations, static structures generally control SWS, but dynamic conditions due to magmatic activity can override the structural control. Monitoring of SWS may therefore provide important signs of impending volcanism.

Hawaii↗
USGS2026

The path to FAIR research models: Lessons learned

Numerical modeling of Earth surface processes emerged as an important scientific tool in the late 1960s to mid-1970s, driven by the development of finite element methods in computer science. These advancements, initially applied in civil engineering, enabled scientists to simulate complex geological phenomena. At that time, models were often only described in publications, access was limited to researchers with direct connections to the developers, and the code was rarely documented for reuse, limiting their application beyond the original research context. The FAIR principles (Findability, Accessibility, Interoperability, and Reusability) as applied to data began to take shape in the 21st century with the rise of open science, digital repositories, and standardized data sharing frameworks. In the late 2010s, grassroots movements began to apply some of the FAIRness goals to numerical models. Subsequently, more formalized FAIR model principles were developed that addressed the specific needs of the scientific modeling community, resulting in the formulation of the FAIR principles for research software (FAIR4RS). In this study, we examine the development and implementation of strategies by two geoscience research infrastructures – the CSDMS (Community Surface Dynamics Modeling System) Model Repository and the U.S. Geological Survey Model Catalog – to enhance the FAIRness of models guided by FAIR4RS. Some of the development and implementation efforts described predate the formalization of FAIR and FAIR4RS principles, making this an ongoing and adaptive process. We evaluate the temporal progression towards increased FAIR4RS alignment across three phases of research infrastructure development: prototype, refinement, and growth and iteration. Although certain principles were more straightforward to implement early in prototypes of the catalog infrastructures, others required broader community collaboration during refinement, and some continue to pose practical challenges in the growth and iteration phase. By tracing these dynamics, our aim is to provide insights that can guide other modeling initiatives in effectively adopting FAIR4RS principles within their communities.

Geoscientific Model Development↗
USGS2014

A geologic approach to field methods in fluvial geomorphology

A geologic approach to field methods in fluvial geomorphology is useful for understanding causes and consequences of past, present, and possible future perturbations in river behavior and floodplain dynamics. Field methods include characterizing river planform and morphology changes and floodplain sedimentary sequences over long periods of time along a longitudinal river continuum. Techniques include topographic and bathymetric surveying of fluvial landforms in valley bottoms and describing floodplain sedimentary sequences through coring, trenching, and examining pits and exposures. Historical sediment budgets that include floodplain sedimentary records can characterize past and present sources and sinks of sediment along a longitudinal river continuum. Describing paleochannels and floodplain vertical accretion deposits, estimating long-term sedimentation rates, and constructing historical sediment budgets can assist in management of aquatic resources, habitat, sedimentation, and flooding issues.

Developments in Earth Surface Processes↗
USGS2010

Inelastic off-fault response and three-dimensional dynamics of earthquake rupture on a strike-slip fault

Large dynamic stress off the fault incurs an inelastic response and energy loss, which contributes to the fracture energy, limiting the rupture and slip velocity. Using an explicit finite element method, we model three-dimensional dynamic ruptures on a vertical strike-slip fault in a homogeneous half-space. The material is subjected to a pressure-dependent Drucker-Prager yield criterion. Initial stresses in the medium increase linearly with depth. Our simulations show that the inelastic response is confined narrowly to the fault at depth. There the inelastic strain is induced by large dynamic stresses associated with the rupture front that overcome the effect of the high confining pressure. The inelastic zone increases in size as it nears the surface. For material with low cohesion (~5 MPa) the inelastic zone broadens dramatically near the surface, forming a "flowerlike" structure. The near-surface inelastic strain occurs in both the extensional and the compressional regimes of the fault, induced by seismic waves ahead of the rupture front under a low confining pressure. When cohesion is large (~10 MPa), the inelastic strain is significantly reduced near the surface and confined mostly to depth. Cohesion, however, affects the inelastic zone at depth less significantly. The induced shear microcracks show diverse orientations near the surface, owing to the low confining pressure, but exhibit mostly horizontal slip at depth. The inferred rupture-induced anisotropy at depth has the fast wave direction along the direction of the maximum compressive stress.

Journal of Geophysical Research B: Solid Earth↗
USGS2021

Temporal seismic velocity variations: Recovery following from the 2019 Mw 7.1 Ridgecrest, California earthquake

We investigated seismic velocity changes ( dv / v ) associated with the 2019 Ridgecrest earthquake sequence with high‐frequency autocorrelations of ambient seismic noise data. Daily autocorrelation functions were computed for the entirety of 2019 and the first quarter of 2020 for broadband stations within the region, including the temporary broadband stations installed during the aftershock deployment. Travel time shifts in the daily autocorrelation functions, relative to the mean autocorrelation waveform, were computed to produce dv / v time series, which are sensitive to the evolving material properties of the shallow crust surrounding the Ridgecrest fault zone (RFZ). A short‐term velocity drop follows the M w 7.1 earthquake at stations in the vicinity of the rupture surface, while those greater than 50 km away showed no such drop. The maximum, absolute changes in seismic velocity are proportional to the logarithm of distance from the fault rupture and to the peak dynamic strain experienced during the earthquake. Near the areas of the highest coseismic slip within the RFZ, seismic velocities recovered over 3 months. However, in the vicinity of the nearby Garlock fault, where triggered slip manifested, and north of the RFZ, seismic velocities recovered within a month. We interpret the seismic velocity changes and their recovery to be largely due to changes in the physical properties of the shallow crust, such as fault zone damage recovery caused by the earthquake rupture process and in response to the large dynamic stresses of passing seismic waves from the mainshock.

California↗
USGS2021

Miocene phosphatization of rocks from the summit of Rio Grande Rise, Southwest Atlantic Ocean

Marine phosphorites are an important part of the oceanic phosphorus cycle and are related to the effects of long-term global climate changes. We use petrography, mineralogy, rare earth elements contents, and 87 Sr/ 86 Sr-determined carbonate fluorapatite (CFA) and calcite ages to investigate the paragenesis and history of phosphatization of carbonate sediments, limestones, ferromanganese crusts, and ironstones from the summit of Rio Grande Rise (RGR), Southwest Atlantic Ocean. Phosphatization of all the rock types occurred throughout the Miocene from 20.2 to 6.8 million years ago (Ma), and occasionally during the Quaternary, mainly through the cementation of carbonate sediments by cryptocrystalline CFA, likely involving the dissolution of the smaller size fraction of foraminifera-nannofossil ooze. Porosity/permeability and abundance of fine calcite material were important factors determining the intensity of phosphatization of the various rock types. Phosphatization was initiated during a transition to a more dynamic circulation system in the South Atlantic Ocean, which remobilized phosphorus from deeper waters and increased primary productivity that culminated with the middle-Miocene Climatic Optimum between ∼17 and 14.8 Ma. The relatively shallow-water depth of RGR summit during the Miocene provided proximity to the oxygen minimum zone, a reservoir for reactive phosphorus, especially during periods of enhanced phosphorus cycling spurred by surface primary productivity. The cessation of phosphatization at RGR resulted from a rapidly cooling and dry climate that characterized the Miocene-Pliocene transition. Our results support previous observations that periods of broadly intensified ocean circulation and local hydrodynamic changes were the key paleoceanographic links to phosphorite formation.

Paleoceanography and Paleoclimatology↗
USGS2013

Large wood budget and transport dynamics on a large river using radio telemetry

Despite the abundance of large wood (LW) river studies there is still a lack of understanding of LW transport dynamics on large low gradient rivers. This study used 290 radio frequency identification tagged (RFID) LW and 54 metal (aluminum) tagged LW, to quantify the percent of in-channel LW that moves per year and what variables play a role in LW transport dynamics. Aluminum tags were installed and monitored on LW in-transit during the rising limb of a flood, the mean distance traveled by those pieces during the week was 13.3 river kilometers (km) with a maximum distance of 72 km. RFID tagged LW moved a mean of 11.9 km/yr with a maximum observed at 101.1 km/yr. Approximately 41% of LW low on the bank moves per year. The high rate of transport and distance traveled is likely due to the lack of interaction between LW floating in the channel and the channel boundaries, caused primarily by the width of the channel relative to length of the LW. Approximately 80% of the RFID tags moved past a fixed reader during the highest 20% of river stage per year. LW transport and logjam dynamics are complicated at high flows as pieces form temporary jams that continually expand and contract. Unlike most other studies, key members that create a logjam were defined more by stability than jam size or channel/hydrologic conditions. Finally, using an existing geomorphic database for the river, and data from this study, we were able to develop a comprehensive LW budget showing that 5% of the in-channel LW population turns over each year (input from mass wasting and fluvial erosion equals burial, decomposition, and export out of system) and another 16% of the population moving within the system.

North Carolina;Virginia↗
USGS2005

Teleseismic body waves from dynamically rupturing shallow thrust faults: Are they opaque for surface-reflected phases?

We investigate whether a shallow-dipping thrust fault is prone to waveslip interactions via surface-reflected waves affecting the dynamic slip. If so, can these interactions create faults that are opaque to radiated energy? Furthermore, in this case of a shallow-dipping thrust fault, can incorrectly assuming a transparent fault while using dislocation theory lead to underestimates of seismic moment? Slip time histories are generated in three-dimensional dynamic rupture simulations while allowing for varying degrees of wave-slip interaction controlled by fault-friction models. Based on the slip time histories, P and SH seismograms are calculated for stations at teleseismic distances. The overburdening pressure caused by gravity eliminates mode I opening except at the tip of the fault near the surface; hence, mode I opening has no effect on the teleseismic signal. Normalizing by a Haskell-like traditional kinematic rupture, we find teleseismic peak-to-peak displacement amplitudes are approximately 1.0 for both P and SH waves, except for the unrealistic case of zero sliding friction. Zero sliding friction has peak-to-peak amplitudes of 1.6 for P and 2.0 for SH waves; the fault slip oscillates about its equilibrium value, resulting in a large nonzero (0.08 Hz) spectral peak not seen in other ruptures. These results indicate wave-slip interactions associated with surface-reflected phases in real earthquakes should have little to no effect on teleseismic motions. Thus, Haskell-like kinematic dislocation theory (transparent fault conditions) can be safety used to simulate teleseismic waveforms in the Earth.

Bulletin of the Seismological Society of America↗
USGS2019

Seasonal, tidal, and geomorphic controls on sediment export to Amazon River tidal floodplains

Mainstem–floodplain material exchange in the tidal freshwater reach of ma jor rivers may lead to significant sequestration of riverine sediment, but this zone remains understudied compared to adjacent fluvial and marine environ ments. This knowledge gap prompts investigation of floodplain-incising tidalchannels found along the banks of tidal rivers and their role in facilitating water and suspended-sediment fluxes between mainstem and floodplain. To evaluate this role, and how it evolves along the tidal river and with time, we measured water level, flow velocity, temperature, and suspended-sediment concentration (SSC) in four tidal channels along the tidal Amazon River, Brazil. Eleven deployments were made during low, rising, high, and falling seasonal Amazon discharge. Generally, channels export high-SSC water from the mainstem to the tidal floodplain on flood tides and transfer low-SSC water back to the mainstem on ebbs. Along the length of the tidal river, the interaction between tidal and seasonal water-level variations and channel–floodplain morphology is a primary control on tidal-channel sediment dynamics. Close to the river mouth, where tides are large, this interaction produces transient flow features and current induced sediment resuspension, but the importance of these processes decreases with distance upstream. Although the magnitude of the exchange of water and sediment between mainstem and floodplain via tidal channels is a small percentage of the total mainstem discharge in this large tidal-river system, tidal channels are important conduits for material flux between these two environments. This flux is critical to resisting floodplain submergence during times of rising sea level.

Amazon River↗
USGS1998

From polar wander to dynamic planet: A tribute to Keith Runcorn

The evolution of Keith Runcorn's ideas from the static, elastic Earth of Jeffreys to a dynamic, convecting planet are presented based on discussions as his colleague over 25 years from 1963 to his retirement in 1988. Keith reached the concept of a dynamic planet by way of polar wander and continental drift using palaeomagnetism with great zeal. It took some time for Keith to convince himself of the reality of continental drift. Once convinced, he became an evangelist converting others and enthusiastically pursuing possible mechanisms for explaining it, homing in on mantle convection. To establish the nature and history of mantle convection his interests ranged from the world rift system and satellite gravity anomalies to the radiometric age peaks. A great step forward occurred when we realised that a region of uprising convection was not necessary under all the rifts as had been commonly advocated in the 1960s and that the mid-Atlantic, African and Indian Ocean rifts could be equally well explained by one large region of upwelling mantle convection. It was also realised that the plate convergence zones (island arcs, trenches and deep focus earthquakes) were much better correlated with the satellite gravity anomalies and it was much easier to locate the possible regions of downwelling mantle convection. Now, seismic tomography helps to establish the nature of mantle convection and it appears that relations among the Earth's surface features, the geoid anomalies and peturbations of mantle seismic velocities are near to being established. In the next few years a far better and accurate picture of the geometry of mantle convection so enthusiastically advocated by Keith Runcorn is likely to be seen.

Physics and Chemistry of the Earth↗
USGS2025

Modeling byproduct and coproduct mine production and mineral substitution using multidimensional supply curves: Application to the Cu-Co-Ni system and beyond

Rapid demand growth is expected for many metals used in the energy transition. Many of these metals are byproducts of other commodities. Byproduct production’s price response is tied to host mineral economics, complicating its supply dynamics. Moreover, many of these metals are used in applications where the material properties desired are difficult to substitute; effectively, limiting how quickly demand can adapt to changes in commodity price. Previous work has demonstrated the interconnectivity of jointly produced mineral commodities from the supply side, where the copper–cobalt–nickel system was used and demand was assumed independent across commodities. Studies to understand byproduct-coproduct market interconnectivity on the demand side are limited, while studies on the interconnectivity of supply and demand simultaneously are even more so. We propose a modification to the multicommodity supply curve method to enable inter-commodity effects on demand simultaneous with supply. In batteries, high cobaltprices may push consumers to transition to high-nickel chemistries, causing the nickel demand surface to decrease with nickel price but increase with cobaltprice, creating a two-dimensional demand surface. Below cross-price elasticities of 0.05, inter-commodity effects were found to be negligible, potentially permitting exclusion of these effects for many commodities. This additional demand curve complexity introduces potential computation challenges alongside the capacity to model many interrelated commodity systems such as rare earth elements, ferroalloys, country-oriented subsidies or restrictions, and bifurcated sustainable metals markets. By presenting the work done on multicommodity supply surfaces to date and potential new directions, this work aims to catalyze the next round of innovative approaches to modeling jointly produced commodities.

Conference Paper↗
USGS2015

Crustal deformation in the New Madrid seismic zone and the role of postseismic processes

Global Navigation Satellite System data across the New Madrid seismic zone (NMSZ) in the central United States over the period from 2000 through 2014 are analyzed and modeled with several deformation mechanisms including the following: (1) creep on subsurface dislocations, (2) postseismic frictional afterslip and viscoelastic relaxation from the 1811–1812 and 1450 earthquakes in the NMSZ, and (3) regional strain. In agreement with previous studies, a dislocation creeping at about 4 mm/yr between 12 and 20 km depth along the downdip extension of the Reelfoot fault reproduces the observations well. We find that a dynamic model of postseismic frictional afterslip from the 1450 and February 1812 Reelfoot fault events can explain this creep. Kinematic and dynamic models involving the Cottonwood Grove fault provide minimal predictive power. This is likely due to the smaller size of the December 1811 event on the Cottonwood Grove fault and a distribution of stations better suited to constrain localized strain across the Reelfoot fault. Regional compressive strain across the NMSZ is found to be less than 3 × 10 −9 /yr. If much of the present-day surface deformation results from afterslip, it is likely that many of the earthquakes we see today in the NMSZ are aftershocks from the 1811–1812 New Madrid earthquakes. Despite this conclusion, our results are consistent with observations and models of intraplate earthquake clustering. Given this and the recent paleoseismic history of the region, we suggest that seismic hazard is likely to remain significant.

Journal of Geophysical Research B: Solid Earth↗
USGS2023

A global perspective on bacterial diversity in the terrestrial deep subsurface

While recent efforts to catalogue Earth’s microbial diversity have focused upon surface and marine habitats, 12–20 % of Earth’s biomass is suggested to exist in the terrestrial deep subsurface, compared to ~1.8 % in the deep subseafloor. Metagenomic studies of the terrestrial deep subsurface have yielded a trove of divergent and functionally important microbiomes from a range of localities. However, a wider perspective of microbial diversity and its relationship to environmental conditions within the terrestrial deep subsurface is still required. Our meta-analysis reveals that terrestrial deep subsurface microbiota are dominated by Betaproteobacteria, Gammaproteobacteria and Firmicutes , probably as a function of the diverse metabolic strategies of these taxa. Evidence was also found for a common small consortium of prevalent Betaproteobacteria and Gammaproteobacteria operational taxonomic units across the localities. This implies a core terrestrial deep subsurface community, irrespective of aquifer lithology, depth and other variables, that may play an important role in colonizing and sustaining microbial habitats in the deep terrestrial subsurface. An in silico contamination-aware approach to analysing this dataset underscores the importance of downstream methods for assuring that robust conclusions can be reached from deep subsurface-derived sequencing data. Understanding the global panorama of microbial diversity and ecological dynamics in the deep terrestrial subsurface provides a first step towards understanding the role of microbes in global subsurface element and nutrient cycling.

Microbiology↗
USGS2007

Dynamic triggering

A growing body of evidence demonstrates that dynamic stresses propagating as seismic waves from large earthquakes are capable of triggering additional earthquakes ranging from aftershocks in the near-field (within one or two source dimensions of the mainshock epicenter) to remotely triggered earthquakes at distances exceeding 10 000 km. Most of the triggered earthquakes are small (generally M ≤ 3) except within the near field, where dynamic stresses may trigger slip on subadjacent fault segments leading to complex rupture comprising of several large earthquakes of comparable magnitude. Crustal surface waves with periods of 15–30 s and peak dynamic stresses greater than ∼0.01 MPa seem to be most efficient in triggering remote seismicity. Current models for dynamic triggering fall under two broad groups: one appealing to Coulomb failure with various friction laws, and the other appealing to the activation of crustal fluids either hydrous or magmatic. No single model appears capable of accounting for the wide variation observed in the nature of triggered activity. Spatial sampling of dynamic triggering on a global scale is still woefully inadequate because of the limited distribution of adequate seismic networks. From the limited data currently available, it appears that extensional stress regimes hosting geothermal and volcanic activity are more susceptible to remote dynamic triggering than compressional stress regimes, although remote triggering is not limited to extensional regimes. Instances of remote triggering in the few areas with continuous, high-resolution deformation instrumentation (all volcanic or geothermal areas) include distinctive deformation transients, suggesting that the locally triggered seismicity in these areas may be a secondary response to a more fundamental aseismic process that likely involves some form of fluid transport or phase change. Recent evidence for triggering by solid Earth tides and ocean loading in convergent plate margins provides a low-frequency, low-amplitude reference point for the spectrum of stresses capable of dynamic triggering. Remaining challenges include establishing better sampling of the distribution of triggered seismicity and better constraints on physical models for the triggering process.

Book chapter↗
USGS2014

Ball-and-socket tectonic rotation during the 2013 M w 7.7 Balochistan earthquake

The September 2013 M w 7.7 Balochistan earthquake ruptured a ∼200-km-long segment of the curved Hoshab fault in southern Pakistan with 10±0.2 m of peak sinistral and ∼1.7±0.8 m of dip slip. This rupture is unusual because the fault dips 60±15° towards the focus of a small circle centered in northwest Pakistan, and, despite a 30° increase in obliquity along strike, the ratios of strike and dip slip remain relatively uniform. Surface displacements and geodetic and teleseismic source inversions quantify a bilateral rupture that propagated rapidly at shallow depths from a transtensional jog near the northern end of the rupture. Static friction prior to rupture was unusually weak (μ<0.05), and friction may have approached zero during dynamic rupture. Here we show that the inward-dipping Hoshab fault defines the northern rim of a structural unit in southeast Makran that rotates – akin to a 2-D ball-and-socket joint – counter-clockwise in response to India's penetration into the Eurasian plate. This rotation accounts for complexity in the Chaman fault system and, in principle, reduces seismic potential near Karachi; nonetheless, these findings highlight deficiencies in strong ground motion equations and tectonic models that invoke Anderson–Byerlee faulting predictions.

Balochistan↗
USGS2010

Two major Cenozoic episodes of phosphogenesis recorded in equatorial Pacific seamount deposits

Seamount phosphorites have been recognized since the 1950s, but this is the first study to provide an in depth exploration of the origin and history of these widespread deposits. Representative samples from equatorial Pacific Cretaceous seamounts were analyzed for chemical, mineralogical, and stable isotope compositions. The phosphorites occur in a wide variety of forms, but most commonly carbonate fluorapatite (CFA) replaced middle Eocene and older carbonate sediment in a deep water environment (>1000 m). Element ratios distinguish seamount phosphorites from continental margin, plateau, and insular phosphorites. Uranium and thorium contents are low and total rare earth element (REE) contents are generally high. REE ratios and shale-normalized patterns demonstrate that the REEs and host CFA were derived from seawater. Strontium isotopic compositions compared with inferred Cenozoic seawater curves define two major episodes of Cenozoic phosphatization: Late Eocene/early Oligocene (39–34 Ma) and late Oligocene/early Miocene (27–21 Ma); three minor events are also indicated. The major episodes occurred at times of climate transition, the first from a nonglacial to glacial earth and the second from a predominantly glacial to warm earth. The paleoceanographic conditions that existed at those times initiated and sustained development of phosphorite by accumulation of dissolved phosphorus in the deep sea during relatively stable climatic conditions when oceanic circulation was sluggish. Fluctuations in climate, sealevel, and upwelling that accompanied the climate transitions may have driven cycles of enrichment and depletion of the deep-sea phosphorus reservoir. As temperature gradients in the oceans increased, Antarctic glaciation expanded and oceanic circulation and upwelling intensified. Expansion and intensification of the oxygen minimum zone may have increased the capacity for midwater storage of phosphorus supplied by dynamic upwelling around seamounts; however, the bottom waters never became anoxic during the phosphogenic episodes. Fluctuations in the CCD and lysocline, CO 2 fluxes, and changes in bottom water circulation and temperatures may have bathed the seamount carbonates in more corrosive waters which, coupled with increased supplies of dissolved phosphorus, promoted replacement processes. The late Eocene/early Oligocene phosphogenic episode recorded in seamount deposits is not matched by large phosphorite deposits in the geologic record, whereas the late Oligocene/early Miocene episode and middle Miocene event are matched by large deposits distributed globally. The seamount phosphorites are exposed at the surface of the seamounts and have been for most of the Neogene and Oligocene. The phosphorites do not show signs of etching that would indicate substantial undersaturation of seawater phosphate with respect to CFA. Mass balance calculations indicate that about 5.4–19 × 10 12 g of P 2 O 5 are locked up in equatorial Pacific seamount phosphorites. That amount is equivalent to about 2-7 years of the present annual input from rivers.

Paleoceanography and Paleoclimatology↗