USGS Science⌕ Search

SEARCH · USGS Science

Results for “The Earth Observer”

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 1,765 records · Page 98Linked to original sources

Interannual variation of rare earth element abundances in corals from northern coast of the South China Sea and its relation with sea-level change and human activities

Here we present interannual rare earth element (REE) records spanning the last two decades of the 20th century in two living Porites corals, collected from Longwan Bay, close to the estuarine zones off Wanquan River of Hainan Island and Hong Kong off the Pearl River Delta of Guangdong Province in the northern South China Sea. The results show that both coral REE contents (0.5–40 ng g −1 in Longwan Bay and 2–250 ng g −1 in Hong Kong for La–Lu) are characterized with a declining trend, which are significantly negative correlated with regional sea-level rise (9.4 mm a −1 from 1981 to 1996 in Longwan Bay, 13.7 mm a −1 from 1991 to 2001 in Hong Kong). The REE features are proposed to be resulted from seawater intrusion into the estuaries in response to contemporary sea-level rise. However, the tendency for the coral Er/Nd time series at Hong Kong site is absent and there is no significant relation between Er/Nd and total REEs as found for the coral at Longwan Bay site. The observations are likely attributed to changes of the water discharge and sediment load of Pearl River, which have been significantly affected by intense human activities, such as the construction of dams/reservoirs and riverbed sediment mining, in past decades. The riverine sediment load/discharge ratio of the Pearl River decreased sharply with a rate of 0.02 kg m −3 a −1 , which could make significant contribution to the declining trend of coral REE. We propose that coastal corals in Longwan Bay and similar unexplored sites with little influences of river discharge and anthropogenic disruption are ideal candidates to investigate the influence of sea-level change on seawater/coral REE.

South China Sea↗

Trace element geochemistry of altered volcanic ash layers (tonsteins) in late Permian coal-bearing formations of eastern Yunnan and western Guizhou Provinces, China

Trace element compositions were determined (by instrumental neutron activation analysis; INAA) in 30 samples of synsedimentary volcanic ash-derived tonsteins and detrital claystones from coal seams within the late Permian coal-bearing formation of eastern Yunnan and western Guizhou Provinces, China. The characteristics of trace-element geochemistry in the tonsteins can be distinguished from those of detrital claystones because of the former's unique volcanic-ash origin. The detrital claystones are characterized by their relatively high content of V, Ti, Sc, Cr, Co and Ni, relatively low content of Th and U, Th/U ratio, and small negative Eu anomaly (Eu/Eu* 0.63-0.93). Overall, these trace element characteristics are consistent with a mafic source similar to the composition of basalt rocks in the erosional region on the western edge of the study area. In contrast, the tonsteins are low in V, Ti, Sc, Cr, Co and Ni contents and have a high Th/U ratio with a distinct negative Eu anomaly (Eu/Eu* normally in the range of 0.2-0.4), consistent with a silicic magmatic source. Within the group of tonsteins, those from the lower section (P2.1) of the coal-bearing formation are relatively high in Ti, Zr, Hf, Nb, Ta and rare earth elements (REE), as compared to those from the middle and upper sections (P2.2+3). In trace-element discrimination diagrams (scatter plots) of Hf-Ta, Ti-Ta, Ti-V, Hf-Sc, Lu-Hf and Lu-Th, tonsteins from the P2.1 horizon always fall in isolated distribution areas, separate from the tonsteins of the P2.2+3 horizon. These results suggest that the source materials of tonsteins from the two separate horizons were probably derived from volcanic ash falls of two distinctly different natures. Based on a comparison of the concentrations and assemblages of trace elements between various magmatic rocks, the source materials of tonsteins from P2.1 horizon were mostly composed of calc-alkalic, silica-poor volcanic ash (similar to rhyodacitic magma), whereas those from P2.+3 were apparently more siliceous and K-rich (rhyolitic magma). Thus, tonsteins from the two different horizons are characterized by unique geochemical properties, which remain constant over a wide lateral extent. Integration of trace-elemental compositions with mineralogical and textural observations makes possible the establishment of tonstein stratigraphy, thus, facilitating more precise and reliable coal-seam correlations. (C) 2000 Elsevier Science B.V. All rights reserved.Trace element compositions were determined in 30 samples of synsedimentary volcanic ash-derived tonsteins and detrital claystones from coal seams within the late Permian coal-bearing formation of eastern Yunnan and western Guizhou Provinces, China. The characteristics of trace-element geochemistry in the tonsteins can be distinguished from those of detrital claystones because of the former's unique volcanic-ash origin. The detrital claystones are characterized by their relatively high content of V, Ti, Sc, CR, Co and Ni, relatively low content of Th and U, Th/U ratio, and small negative Eu anomaly.

International Journal of Coal Geology↗

Determination and prediction of rare earth element eeochemical associations in acid mine drainage treatment wastes

Acid mine drainage (AMD) has been proposed by various researchers as a novel source of rare earth elements (REE), a group of elements that include critical metals for clean energy and modern technologies. REE tend to be sequestered in the Fe-Al-Mn-rich solids produced during the treatment of AMD. These solids are typically managed as waste, but could be a low-cost, readily available REE source. Here, results from field sampling, solids characterization, and geochemical modeling are presented to identify the mechanism(s) of REE attenuation and determine the minerals/solid phases in AMD solids that are enriched in REE. This study reveals that solids produced from low-pH AMD that was passively treated by limestone contain elevated concentrations of REE with Al, Fe, and/or Mn. AMD solid characterization via sequential extraction demonstrated that Al and Mn oxides were more abundant than Fe oxides and that the REEs are mainly associated with Al/Mn phases. Additionally, sequential extractions demonstrate that for the AMD solids evaluated, acidic and/or reducing extractions are required to mobilize the REE. Finally, the “CausticTitrationREYs.exe” geochemical equilibrium model demonstrated in this study indicates that the observed dissolved REE attenuation can be explained via surface complexation on Fe, Al, and Mn oxides/hydroxides and not by REE compound precipitation. The model accurately predicts the pH dependent removal of dissolved REE and that Al and Mn oxides/hydroxides are largely responsible for dissolved REE removal for the systems evaluated. The modeling results are consistent with the characterization results that show that Al and Mn hydroxides are important hosting phases of REEs in AMD treatment systems. The results presented here can be used to identify conditions favorable for accumulation of REE-enriched AMD solids and possible chemical treatment(s) to mobilize REE. The geochemical model can be applied to active and/or passive AMD treatment systems to predict REE attenuation with Fe, Al, and Mn during treatment and what phases may be enriched in REE. This information can be used to engineer AMD systems to produce specific phases enriched in REE. The recovery of REE from AMD solids is an opportunity to transform the environmental and economical challenge of polluted mine drainage into an asset.

Conference Paper↗

Impact of precipitation and increasing temperatures on drought trends in eastern Africa

In eastern Africa droughts can cause crop failure and lead to food insecurity. With increasing temperatures, there is an a priori assumption that droughts are becoming more severe. However, the link between droughts and climate change is not sufficiently understood. Here we investigate trends in long-term agricultural drought and the influence of increasing temperatures and precipitation deficits. Using a combination of models and observational datasets, we studied trends, spanning the period from 1900 (to approximate pre-industrial conditions) to 2018, for six regions in eastern Africa in four drought-related annually averaged variables: soil moisture, precipitation, temperature, and evaporative demand ( E 0 ). In standardized soil moisture data, we found no discernible trends. The strongest influence on soil moisture variability was from precipitation, especially in the drier or water-limited study regions; temperature and E 0 did not demonstrate strong relations to soil moisture. However, the error margins on precipitation trend estimates are large and no clear trend is evident, whereas significant positive trends were observed in local temperatures. The trends in E 0 are predominantly positive, but we do not find strong relations between E 0 and soil moisture trends. Nevertheless, the E 0 trend results can still be of interest for irrigation purposes because it is E 0 that determines the maximum evaporation rate. We conclude that until now the impact of increasing local temperatures on agricultural drought in eastern Africa is limited and we recommend that any soil moisture analysis be supplemented by an analysis of precipitation deficit.

Earth Systems Science Dynamics↗

Nearly frictionless faulting by unclamping in long-term interaction models

In defiance of direct rock-friction observations, some transform faults appear to slide with little resistance. In this paper finite element models are used to show how strain energy is minimized by interacting faults that can cause long-term reduction in fault-normal stresses (unclamping). A model fault contained within a sheared elastic medium concentrates stress at its end points with increasing slip. If accommodating structures free up the ends, then the fault responds by rotating, lengthening, and unclamping. This concept is illustrated by a comparison between simple strike-slip faulting and a mid-ocean-ridge model with the same total transform length; calculations show that the more complex system unclapms the transforms and operates at lower energy. In another example, the overlapping San Andreas fault system in the San Francisco Bay region is modeled; this system is complicated by junctions and stepovers. A finite element model indicates that the normal stress along parts of the faults could be reduced to hydrostatic levels after ???60-100 k.y. of system-wide slip. If this process occurs in the earth, then parts of major transform fault zones could appear nearly frictionless.

Geology↗

Dislocation models of interseismic deformation in the western United States

The GPS-derived crustal velocity field of the western United States is used to construct dislocation models in a viscoelastic medium of interseismic crustal deformation. The interseismic velocity field is constrained by 1052 GPS velocity vectors spanning the ???2500-km-long plate boundary zone adjacent to the San Andreas fault and Cascadia subduction zone and extending ???1000 km into the plate interior. The GPS data set is compiled from U.S. Geological Survey campaign data, Plate Boundary Observatory data, and the Western U.S. Cordillera velocity field of Bennett et al. (1999). In the context of viscoelastic cycle models of postearthquake deformation, the interseismic velocity field is modeled with a combination of earthquake sources on ???100 known faults plus broadly distributed sources. Models that best explain the observed interseismic velocity field include the contributions of viscoelastic relaxation from faulting near the major plate margins, viscoelastic relaxation from distributed faulting in the plate interior, as well as lateral variations in depth-averaged rigidity in the elastic lithosphere. Resulting rigidity variations are consistent with reduced effective elastic plate thickness in a zone a few tens of kilometers wide surrounding the San Andreas fault (SAF) system. Primary deformation characteristics are captured along the entire SAF system, Eastern California Shear Zone, Walker Lane, the Mendocino triple junction, the Cascadia margin, and the plate interior up to ???1000 km from the major plate boundaries.

Journal of Geophysical Research B: Solid Earth↗

Analysis of two-color geodimeter measurements of deformation within the Long Valley caldera: June 1983 to October 1985

Line length changes from several baselines in a trilateration network within the Long Valley caldera clearly define a decrease in strain rate from June 1983 through October 1985. The data consist of more than 1600 length measurements on 23 baselines using a two-color geodimeter, which has a precision of 0.2 ppm of the line length. Initial measurements made during the summer of 1983 show extension rates were as high as 6 ppm/yr on some baselines. The rates subsequently decreased steadily, to less than 1 ppm/yr during the summer of 1985. Ih contrast, the strain rates on a few other baselines were less than 0.2 ppm/yr averaged over the 2.3-year interval. A model is constructed using these observations as well as yearly observations of a trilateration network within and near the caldera. The model contains two points of inflation located at 5 and 10 km beneath the resurgent dome of the caldera plus dextral slip on a fault plane within the south moat within an elastic half-space. To mimic the observed time dependence in the two-color geodimeter data, we conclude that the rates of both slip and inflation have steadily decreased over the survey interval. Because the two-color geodimeter network is closer to all three sources than the single-color geodimeter network, these data are important for defining the factor of 2–3 rate decrease of slip and inflation. Comparing our estimates of slip and inflation with similar estimates for the year prior to July 1983, we conclude that the rate decrease is probably a result of postseismic and postinflation relaxation following the earthquake swarm in the south moat in January 1983.

Journal of Geophysical Research Solid Earth↗

Antecedent flow conditions and nitrate concentrations in the Mississippi River basin

The relationship between antecedent flow conditions and nitrate concentrations was explored at eight sites in the 2.9 million square kilometers (km 2 ) Mississippi River basin, USA. Antecedent flow conditions were quantified as the ratio between the mean daily flow of the previous year and the mean daily flow from the period of record (Qratio), and the Qratio was statistically related to nitrate anomalies (the unexplained variability in nitrate concentration after filtering out season, long-term trend, and contemporaneous flow effects) at each site. Nitrate anomaly and Qratio were negatively related at three of the four major tributary sites and upstream in the Mississippi River, indicating that when mean daily streamflow during the previous year was lower than average, nitrate concentrations were higher than expected. The strength of these relationships increased when data were subdivided by contemporaneous flow conditions. Five of the eight sites had significant negative relationships ( p ≤ 0.05) at high or moderately high contemporaneous flows, suggesting nitrate that accumulates in these basins during a drought is flushed during subsequent high flows. At half of the sites, when mean daily flow during the previous year was 50 percent lower than average, nitrate concentration can be from 9 to 27 percent higher than nitrate concentrations that follow a year with average mean daily flow. Conversely, nitrate concentration can be from 8 to 21 percent lower than expected when flow during the previous year was 50 percent higher than average. Previously documented for small, relatively homogenous basins, our results suggest that relationships between antecedent flows and nitrate concentrations are also observable at a regional scale. Relationships were not observed (using all contemporaneous flow data together) for basins larger than 1 million km 2 , suggesting that above this limit the overall size and diversity within these basins may necessitate the use of more complicated statistical approaches or that there may be no discernible basin-wide relationship with antecedent flow. The relationships between nitrate concentration and Qratio identified in this study serve as the basis for future studies that can better define specific hydrologic processes occurring during and after a drought (or high flow period) which influence nitrate concentration, such as the duration or magnitude of low flows, and the timing of low and high flows.

Mississippi River Basin↗

Chronology and provenance of last-glacial (Peoria) loess in western Iowa and paleoclimatic implications

Geologic archives show that the Earth was dustier during the last glacial period. One model suggests that increased gustiness (stronger, more frequent winds) enhanced dustiness. We tested this at Loveland, Iowa, one of the thickest deposits of last-glacial-age (Peoria) loess in the world. Based on K/Rb and Ba/Rb, loess was derived not only from glaciogenic sources of the Missouri River, but also distal loess from non-glacial sources in Nebraska. Optically stimulated luminescence (OSL) ages provide the first detailed chronology of Peoria Loess at Loveland. Deposition began after ~ 27 ka and continued until ~ 17 ka. OSL ages also indicate that mass accumulation rates (MARs) of loess were not constant. MARs were highest and grain size was coarsest during the time of middle Peoria Loess accretion, ~ 23 ka, when ~ 10 m of loess accumulated in no more than ~ 2000 yr and possibly much less. The timing of coarsest grain size and highest MAR, indicating strongest winds, coincides with a summer-insolation minimum at high latitudes in North America and the maximum southward extent of the Laurentide ice sheet. These observations suggest that increased dustiness during the last glacial period was driven largely by enhanced gustiness, forced by a steepened meridional temperature gradient.

Quaternary Research↗

Social-ecological predictors of global invasions and extinctions

Most assessments of resilience have been focused on local conditions. Studies focused on the relationship between humanity and environmental degradation are rare, and are rarely comprehensive. We investigated multiple social-ecological factors for 100 countries around the globe in relation to the percentage of invasions and extinctions within each country. These 100 countries contain approximately 87% of the world’s population, produce 43% of the world’s per capita gross domestic product (GDP), and take up 74% of the earth’s total land area. We used an information theoretic approach to determine which models were most supported by our data, utilizing an a priori set of plausible models that included a combination of 15 social-ecological variables, each social-ecological factor by itself, and selected social-ecological factors grouped into three broad classes. These variables were per capita GDP, export-import ratio, tourism, undernourishment, energy efficiency, agricultural intensity, rainfall, water stress, wilderness protection, total biodiversity, life expectancy, adult literacy, pesticide regulation, political stability, and female participation in government. Our results indicate that as total biodiversity and total land area increase, the percentage of endangered birds also increases. As the independent variables (agricultural intensity, rainfall, water stress, and total biodiversity) in the ecological class model increase, the percentage of endangered mammals in a country increases. The percentage of invasive birds and mammals in a country increases as per capita GDP increases. As life expectancy increases, the percentage of invasive and endangered birds and mammals increases. Although our analysis does not determine mechanisms, the patterns observed in this study provide insight into the dynamics of a complex, global, social-ecological system.

Ecology and Society↗

Modeling of barrier breaching during Hurricanes Sandy and Matthew

Physical processes driving barrier island change during storms are important to understand to mitigate coastal hazards and to evaluate conceptual models for barrier evolution. Spatial variations in barrier island topography, landcover characteristics, and nearshore and back-barrier hydrodynamics can yield complex morphological change that requires models of increasing resolution and physical complexity to predict. Using the Coupled Ocean-Atmosphere-Wave-Sediment Transport (COAWST) modeling system, we investigated two barrier island breaches that occurred on Fire Island, NY during Hurricane Sandy (2012) and at Matanzas, FL during Hurricane Matthew (2016). The model employed a recently implemented infragravity (IG) wave driver to represent the important effects of IG waves on nearshore water levels and sediment transport. The model simulated breaching and other changes with good skill at both locations, resolving differences in the processes and evolution. The breach simulated at Fire Island was 250 m west of the observed breach, whereas the breach simulated at Matanzas was within 100 m of the observed breach. Implementation of the vegetation module of COAWST to allow three-dimensional drag over dune vegetation at Fire Island improved model skill by decreasing flows across the back-barrier, as opposed to varying bottom roughness that did not positively alter model response. Analysis of breach processes at Matanzas indicated that both far-field and local hydrodynamics influenced breach creation and evolution, including remotely generated waves and surge, but also surge propagation through back-barrier waterways. This work underscores the importance of resolving the complexity of nearshore and back-barrier systems when predicting barrier island change during extreme events.

Florida, New York↗

A model integrating satellite-derived shoreline observations for predicting fine-scale shoreline response to waves and sea-level rise across large coastal regions

Satellite-derived shoreline observations combined with dynamic shoreline models enable fine-scale predictions of coastal change across large spatiotemporal scales. Here, we present a satellite-data-assimilated, “littoral-cell”-based, ensemble Kalman-filter shoreline model to predict coastal change and uncertainty due to waves, sea-level rise (SLR), and other natural and anthropogenic processes. We apply the developed ensemble model to the entire California coastline (approximately 1,760 km), much of which is sparsely monitored with traditional survey methods (e.g., Lidar/GPS). Water-level-corrected, satellite-derived shoreline observations (obtained from the CoastSat toolbox) offer a nearly unbiased representation of in situ surveyed shorelines (e.g., mean sea-level elevation contours) at Ocean Beach, San Francisco. We demonstrate that model calibration with satellite observations during a 20-year hindcast period (1995–2015) provides nearly equivalent model forecast accuracy during a validation period (2015–2020) compared to model calibration with monthly in situ observations at Ocean Beach. When comparing model-predicted shoreline positions to satellite-derived observations, the model achieves an accuracy of <10 m RMSE for nearly half of the entire California coastline for the validation period. The calibrated/validated model is then applied for multi-decadal simulations of shoreline change due to projected wave and sea-level conditions, while holding the model parameters fixed. By 2100, the model estimates that 24%–75% of California's beaches may become completely eroded due to SLR scenarios of 1.0–3.0 m, respectively. The satellite-data-assimilated modeling system presented here is generally applicable to a variety of coastal settings around the world owing to the global coverage of satellite imagery.

JGR Earth Surface↗

Petrology, composition, and age of intrusive rocks associated with the Quartz Hill molybdenite deposit, southeastern Alaska

A large porphyry molybdenum deposit (Quartz Hill deposit) was recently discovered in the heart of the Coast Range batholithic complex about 70 km east of Ketchikan, southeastern Alaska. Intrusive rocks associated with the mineral deposit form two composite epizonal to hypabyssal stocks and many dikes in country rocks. The stocks are characterized by a variety of textural rock types varying from equigranular or weakly seriate biotite granite to porphyries with aphanitic or very fine grained and aplitic groundmasses. These rocks contain about equal amounts of quartz, albitic plagioclase, and microperthitic microcline and less than 2.5% biotite. Unaltered rocks contain between 0.2 and 1% CaO, less than 1.7% combined Fe 2 O 3 , FeO, and MgO, and 74.4 to 77.7% SiO 2 . Total alkalis are between 8 and 9%, and K 2 O/Na 2 O is about 1.1. The range of major-oxide variation is small, but it is systematically related to lithology. Many trace-elements, including B, Pb, Sn, and Li have low concentrations. Intrusive rocks associated with the Quartz Hill deposit are more albitic and possibly trace-element depleted compared to some other rocks associated with porphyry molybdenum deposits.All observed metallization and alteration is within the Quartz Hill stock. Molybdenite forms fracture coatings and occurs in veins with quartz. Alteration is widespread and includes development of secondary quartz, pyrite, K-feldspar, biotite, white mica, chlorite, and zeolite. Field relations indicate that the stocks were emplaced after regional uplift and erosion of the Coast Range batholithic complex, and K–Ar data show that intrusion and alteration took place in late Oligocene time, about 27 to 30 Ma ago. Data from the Ketchikan quadrangle indicate that porphyry molybdenum metallization in the Coast Range batholithic complex is associated with regionally extensive but spotty, middle Tertiary or younger, felsic magmatism.

Alaska↗

Probabilistic assessment of postfire debris-flow inundation in response to forecast rainfall

Communities downstream of burned steep lands face increases in debris-flow hazards due to fire effects on soil and vegetation. Rapid postfire hazard assessments have traditionally focused on quantifying spatial variations in debris-flow likelihood and volume in response to design rainstorms. However, a methodology that provides estimates of debris-flow inundation downstream of burned areas based on forecast rainfall would provide decision-makers with information that directly addresses the potential for downstream impacts. We introduce a framework that integrates a 24 h lead-time ensemble precipitation forecast with debris-flow likelihood, volume, and runout models to produce probabilistic maps of debris-flow inundation. We applied this framework to simulate debris-flow inundation associated with the 9 January 2018 debris-flow event in Montecito, California, USA. When the observed debris-flow volumes were used to drive the probabilistic forecast model, analysis of the simulated inundation probabilities demonstrates that the model is both reliable and sharp. In the fully predictive model, however, in which debris-flow likelihood and volume were computed from the atmospheric model ensemble's predictions of peak 15 min rainfall intensity, I 15 , the model generally under-forecasted the inundation area. The observed peak I 15 lies in the upper tail of the atmospheric model ensemble spread; thus a large fraction of ensemble members forecast lower I 15 than observed. Using these I 15 values as input to the inundation model resulted in lower-than-observed flow volumes which translated into under-forecasting of the inundation area. Even so, approximately 94 % of the observed inundated area was forecast to have an inundation probability greater than 1 %, demonstrating that the observed extent of inundation was generally captured within the range of outcomes predicted by the model. Sensitivity analyses indicate that debris-flow volume and two parameters associated with debris-flow mobility exert significant influence on inundation predictions, but reducing uncertainty in postfire debris-flow volume predictions will have the largest impact on reducing inundation outcome uncertainty. This study represents a first step toward a near-real-time hazard assessment product that includes probabilistic estimates of debris-flow inundation and provides guidance for future improvements to this and similar model frameworks by identifying key sources of uncertainty.

California↗

USGS Digital Spectral Library splib06a

Introduction We have assembled a digital reflectance spectral library that covers the wavelength range from the ultraviolet to far infrared along with sample documentation. The library includes samples of minerals, rocks, soils, physically constructed as well as mathematically computed mixtures, plants, vegetation communities, microorganisms, and man-made materials. The samples and spectra collected were assembled for the purpose of using spectral features for the remote detection of these and similar materials. Analysis of spectroscopic data from laboratory, aircraft, and spacecraft instrumentation requires a knowledge base. The spectral library discussed here forms a knowledge base for the spectroscopy of minerals and related materials of importance to a variety of research programs being conducted at the U.S. Geological Survey. Much of this library grew out of the need for spectra to support imaging spectroscopy studies of the Earth and planets. Imaging spectrometers, such as the National Aeronautics and Space Administration (NASA) Airborne Visible/Infra Red Imaging Spectrometer (AVIRIS) or the NASA Cassini Visual and Infrared Mapping Spectrometer (VIMS) which is currently orbiting Saturn, have narrow bandwidths in many contiguous spectral channels that permit accurate definition of absorption features in spectra from a variety of materials. Identification of materials from such data requires a comprehensive spectral library of minerals, vegetation, man-made materials, and other subjects in the scene. Our research involves the use of the spectral library to identify the components in a spectrum of an unknown. Therefore, the quality of the library must be very good. However, the quality required in a spectral library to successfully perform an investigation depends on the scientific questions to be answered and the type of algorithms to be used. For example, to map a mineral using imaging spectroscopy and the mapping algorithm of Clark and others (1990a, 2003b), one simply needs a diagnostic absorption band. The mapping system uses continuum-removed reference spectral features fitted to features in observed spectra. Spectral features for such algorithms can be obtained from a spectrum of a sample containing large amounts of contaminants, including those that add other spectral features, as long as the shape of the diagnostic feature of interest is not modified. If, however, the data are needed for radiative transfer models to derive mineral abundances from reflectance spectra, then completely uncontaminated spectra are required. This library contains spectra that span a range of quality, with purity indicators to flag spectra for (or against) particular uses. Acquiring spectral measurements and performing sample characterizations for this library has taken about 15 person-years of effort. Software to manage the library and provide scientific analysis capability is provided (Clark, 1980, 1993). A personal computer (PC) reader for the library is also available (Livo and others, 1993). The program reads specpr binary files (Clark, 1980, 1993) and plots spectra. Another program that reads the specpr format is written in IDL (Kokaly, 2005). In our view, an ideal spectral library consists of samples covering a very wide range of materials, has large wavelength range with very high precision, and has enough sample analyses and documentation to establish the quality of the spectra. Time and available resources limit what can be achieved. Ideally, for each mineral, the sample analysis would include X-ray diffraction (XRD), electron microprobe (EM) or X-ray fluorescence (XRF), and petrographic microscopic analyses. For some minerals, such as iron oxides, additional analyses such as Mossbauer would be helpful. We have found that to make the basic spectral measurements, provide XRD, EM or XRF analyses, and microscopic analyses, document the results, and complete an entry of one spectral library sample, all takes about

Data Series↗

Formation of natural gas hydrates in marine sediments 1. Conceptual model of gas hydrate growth conditioned by host sediment properties

The stability of submarine gas hydrates is largely dictated by pressure and temperature, gas composition, and pore water salinity. However, the physical properties and surface chemistry of deep marine sediments may also affect the thermodynamic state, growth kinetics, spatial distributions, and growth forms of clathrates. Our conceptual model presumes that gas hydrate behaves in a way analogous to ice in a freezing soil. Hydrate growth is inhibited within fine-grained sediments by a combination of reduced pore water activity in the vicinity of hydrophilic mineral surfaces, and the excess internal energy of small crystals confined in pores. The excess energy can be thought of as a "capillary pressure" in the hydrate crystal, related to the pore size distribution and the state of stress in the sediment framework. The base of gas hydrate stability in a sequence of fine sediments is predicted by our model to occur at a lower temperature (nearer to the seabed) than would be calculated from bulk thermodynamic equilibrium. Capillary effects or a build up of salt in the system can expand the phase boundary between hydrate and free gas into a divariant field extending over a finite depth range dictated by total methane content and pore-size distribution. Hysteresis between the temperatures of crystallization and dissociation of the clathrate is also predicted. Growth forms commonly observed in hydrate samples recovered from marine sediments (nodules, and lenses in muds; cements in sands) can largely be explained by capillary effects, but kinetics of nucleation and growth are also important. The formation of concentrated gas hydrates in a partially closed system with respect to material transport, or where gas can flush through the system, may lead to water depletion in the host sediment. This "freeze-drying" may be detectable through physical changes to the sediment (low water content and overconsolidation) and/or chemical anomalies in the pore waters and metastable presence of free gas within the normal zone of hydrate stability.

Journal of Geophysical Research B: Solid Earth↗

Oligocene–Miocene northward growth of the Tibetan Plateau: Insights from intermontane basins in the West Qinling Belt, NW China

Growth of the Tibetan Plateau, Earth’s broadest and highest elevation collisional system, shapes orographic barriers, reorganizes drainage networks, and influences surface erosion and sediment delivery, whose changes in space and provenance feed back to intracontinental tectonic processes. Studies of interior basins within the northern Tibetan Plateau provide new sediment accumulation, provenance, paleodrainage, and deformation timing data that enable a reconstruction of the far-field tectono-geomorphic evolution of the rising Tibetan Plateau. Along the northern plateau margin, topographic growth in the West Qinling Belt is inferred to have initiated in the Eocene, nearly coeval with the India-Asia collision, as well as in the late Miocene. However, geological knowledge about the intervening period remains at present enigmatic, and the kinematics and dynamics are uncertain. This study presents a multidisciplinary data set from the intermontane Anhua-Huicheng Basin (AHB; Gansu Province, China) to fill this gap. Magnetostratigraphic dating, regional mapping, and sedimentological analysis imply that contractional deformation and thrust-top basin systems formed within the West Qinling Belt in the Oligocene (not later than ca. 24 Ma). A combination of observations including paleocurrent changes, detrital zircon U-Pb age variations, and appearance of growth strata along the Anhua-Huicheng Basin reveal the rapid uplift of the West Qinling Belt at ca. 15 Ma. Sedimentation in the intermontane basins ended after the late Miocene (ca. 8 Ma), when the region experienced intrabasinal deformation, uplift, and erosion with the establishment of an external drainage system. Since the late Miocene, the growth of the West Qinling Belt reached a climax with the lack of substantial contractional deformation in Cenozoic sequences heralding the onset of the modern kinematic regime and attainment of high elevation. Observed transitions in the tectonostratigraphy and paleodrainage define different phases of deformation and plateau-wide shifts in stress reorganization, which led to the northward growth and later lateral expansion of the Tibetan Plateau.

GSA Bulletin↗

The digital archivist: Automating legacy macroseismic data processing using large language models

Macroseismic data are a key resource to investigate shaking and damage from preinstrumental and early instrumental eras. However, data are often stored as inconsistently formatted reports describing observed shaking and damage, making manually parsing and interpreting accounts labor‐intensive. We introduce a novel workflow using Google’s Gemini 2.5 Pro large language model (LLM) to automate the extraction and structuring of macroseismic observations from summary reports. We apply this workflow to the 22 March 1957 M 5.3 Daly City, California, earthquake as a case study. We used Gemini to extract addresses, originally assigned modified Mercalli intensity values, and descriptions from each report. To address coordinate precision limits, addresses were geocoded via Google’s Geocoding application programming interface. This workflow yielded over 2300 geocoded intensity reports for the Daly City earthquake. We use the geocoded accounts, with the original report intensity assignments, to develop a shaking intensity map that in some respects rivals modern Did You Feel It? Maps. We also extract and present data for the 9 February 1971 M L 6.7 Sylmar, California, earthquake. Our results demonstrate the potential of LLMs for reliably extracting and analyzing large, unstructured macroseismic datasets. LLMs offer a scalable solution for rapidly digitizing macroseismic archives, enabling their broader use to constrain ground‐motion models in modern seismic hazard analysis and to improve our understanding of site effects in urban areas. The concepts explored here may also be applied to the handling of other legacy seismological and earth science data.

Seismological Research Letters↗