USGS Science⌕ Search

SEARCH · USGS Science

Results for “Earth”

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,279 records · Page 71Linked to original sources

A novel automatic phenology learning (APL) method of training sample selection using multiple datasets for time-series land cover mapping

The long record of Landsat imagery, which is the cornerstone of Earth observation, provides an opportunity to monitor land use and land cover (LULC) change and understand the interactions between the climate and earth system through time. A few change detection algorithms such as Continuous Change Detection and Classification (CCDC) have been developed to utilize all available Landsat images for change detection and characterization at local or global scales. However, the reliable, rapid, and reproducible collection of training samples have become a challenge for time series land cover classification at a large scale. To meet the challenge, we proposed an automatic phenology learning (APL) method with the assumption that the temporal profiles of samples within the same land cover type are the same or similar at a local scale to generate evenly distributed training samples automatically. We designed the method to build land cover patterns for each category based on consensus samples derived from multiple existing scientific datasets including LANDFIRE's (LF) Existing Vegetation Type (EVT), USGS National Land Cover Database (NLCD), National Agricultural Statistics Service (NASS) Cropland Data Layer (CDL), and National Wetlands Inventory (NWI). Then we calculated the Time-Weighted Dynamic Time Warping (twDTW) distance between any undefined samples and land cover patterns in the same geographical region as prior knowledge. Finally, we selected the optimal land cover category for each undefined sample from the land cover products based on the designed criteria iteratively using the twDTW distance as an indicator. The method was applied in the footprint of 10 selected Landsat Analysis Ready Data (ARD) tiles in the eastern and western conterminous United States (CONUS) to produce annual land cover maps from 1985 to 2017. The accuracy assessment and visual comparison revealed that the APL method can generate reliable training samples without any manual interpretation, producing better land cover results especially for the grass/shrub and wetland land cover classes. Applying the APL method, the overall accuracy of the annual land cover maps was improved by 2% over the accuracy of Land Change Monitoring, Assessment, and Projection (LCMAP) Collection 1.0 Science Products in the research regions. Our results also indicate that the APL method provides an approach for best use of different land cover products and meets the requirement of intensive sampling for training data collection.

Remote Sensing of Environment↗

An accuracy assessment of the surface reflectance product from the EMIT imaging spectrometer

The Earth surface Mineral dust source InvesTigation (EMIT) is an imaging spectrometer launched to the International Space Station in July 2022 to measure the mineral composition of Earth’s dust-producing regions. We present a systematic accuracy assessment of the EMIT surface reflectance product in two parts. First, we characterize the surface reflectance product’s overall performance using multiple independent vicarious calibration field experiments with hand-held and automated field spectrometers. We find that the EMIT surface reflectance product has a standard error of ±1.0 % in absolute reflectance units for temporally coincident observations. Discrepancies rise to ±2.7 % for spectra acquired at different dates and times of day, which we attribute mainly to changes in solar geometry. Second, we develop an error budget that explains the differences between EMIT and in-situ field spectrometer data. We find that uncertainties in spatial footprints, field spectroscopy, and the EMIT-reported measurement were sufficient to explain discrepancies in most cases. Our approach did not detect any systematic calibration or reflectance errors in the timespan considered. Together, these findings demonstrate that a space-based imaging spectrometer can acquire high-quality spectra across a wide range of observational and atmospheric conditions.

Nevada↗

Enhancements to the USGS Landsat Level 2 surface temperature and emissivity product for Collection 3 reprocessing

The Landsat program provides the longest continuous global record of thermal infrared observations of the Earth's surface, underpinning critical applications in climate monitoring, water resources, ecosystem dynamics, urban heat analysis, and natural hazard assessment. The release of a global inventory of Landsat Collection 2 Level 2 surface temperature products by the U.S. Geological Survey (USGS) marked a major milestone in operational provision of Landsat thermal infrared analysis-ready data. Ongoing validations and community uses of Collection 2 have identified opportunities to further improve accuracy, uncertainty characterization, and emissivity correction across diverse atmospheric and surface conditions. In preparation for the planned Landsat Collection 3 reprocessing of the Landsat data record in the late 2020s, the USGS is implementing a coordinated set of enhancements to the Level 2 surface temperature products. These include revised emissivity estimation that leverages external datasets, improved atmospheric characterization and uncertainty propagation, expanded dynamic range for high temperature targets, consideration of split window atmospheric correction algorithm for Landsat 8 and 9, and decoupling of thermal infrared processing from visible to shortwave infrared constraints to enable surface temperature retrievals under low or no solar illumination conditions. These changes are designed to improve product quality and consistency across the Landsat record. Beyond near-term performance gains, the Collection 3 design establishes a scalable processing architecture to accommodate the expanded spectral and radiometric measurement capabilities of the forthcoming Landsat 10 mission. By preserving continuity across the Landsat 4–9 record while enabling future algorithm evolution, Landsat Collection 3 will provide a foundation for long-term, multi-decadal Earth system thermal infrared observations.

Remote Sensing of Environment↗

Ammonia and aquatic ecosystems – A review of global sources, biogeochemical cycling, and effects on fish

The purpose of this review is to better understand the full life cycle and influence of ammonia from an aquatic biology perspective. While ammonia has toxic properties in water and air, it also plays a central role in the biogeochemical nitrogen (N) cycle and regulates mechanisms of normal and abnormal fish physiology. Additionally, as the second most synthesized chemical on Earth, ammonia contributes economic value to many sectors, particularly fertilizers, energy storage, explosives, refrigerants, and plastics. But, with so many uses, industrial N 2 -fixation effectively doubles natural reactive N concentrations in the environment. The consequence is global, with excess fixed nitrogen driving degradation of soils, water, and air; intensifying eutrophication, biodiversity loss, and climate change; and creating health risks for humans, wildlife, and fisheries. Thus, the need for ammonia research in aquatic systems is growing. In response, we prepared this review to better understand the complexities and connectedness of environmental ammonia. Even the term “ammonia” has multiple meanings. So, we have clarified the nomenclature, identified units of measurement, and summarized methods to measure ammonia in water. We then discuss ammonia in the context of the N-cycle, review its role in fish physiology and mechanisms of toxicity, and integrate the effects of human N-fixation, which continuously expands ammonia's sources and uses. Ammonia is being developed as a carbon-free energy carrier with potential to increase reactive nitrogen in the environment. With this in mind, we review the global impacts of excess reactive nitrogen and consider the current monitoring and regulatory frameworks for ammonia. The presented synthesis illustrates the complex and interactive dynamics of ammonia as a plant nutrient, energy molecule, feedstock, waste product, contaminant, N-cycle participant, regulator of animal physiology, toxicant, and agent of environmental change. Few molecules are as influential as ammonia in the management and resilience of Earth's resources.

Science of the Total Environment↗

iBluff: An open-source R package for geomorphic analysis of coastal bluffs/cliffs

The R package iBluff is designed for coastal bluffs/bluffs morphological analysis and offers an automatic and reproducible alternative to identify bluff edges using a bare earth digital elevation model (DEM) instead of hand digitizing. This package extracts elevation profiles along automatically identified transects on the bluff-face, bluff top, toe, secondary inflections, relative concavity/convexity of bluff-face, and beach dunes (crests and troughs). The package requires at a minimum a bare earth DEM as a raster and a generalized line shapefile (shoreline) approximately parallel with the bluff-face. Both files should be in the same projected coordinate system. The iBluff package was developed to expand and generalize studies of high-relief coastal areas, investigate erosion and seasonality, and could be extended to use three-dimensional (3D) point-cloud data instead of a DEM.

SoftwareX↗

Simple equations guide high-frequency surface-wave investigation techniques

We discuss five useful equations related to high-frequency surface-wave techniques and their implications in practice. These equations are theoretical results from published literature regarding source selection, data-acquisition parameters, resolution of a dispersion curve image in the frequency-velocity domain, and the cut-off frequency of high modes. The first equation suggests Rayleigh waves appear in the shortest offset when a source is located on the ground surface, which supports our observations that surface impact sources are the best source for surface-wave techniques. The second and third equations, based on the layered earth model, reveal a relationship between the optimal nearest offset in Rayleigh-wave data acquisition and seismic setting - the observed maximum and minimum phase velocities, and the maximum wavelength. Comparison among data acquired with different offsets at one test site confirms the better data were acquired with the suggested optimal nearest offset. The fourth equation illustrates that resolution of a dispersion curve image at a given frequency is directly proportional to the product of a length of a geophone array and the frequency. We used real-world data to verify the fourth equation. The last equation shows that the cut-off frequency of high modes of Love waves for a two-layer model is determined by shear-wave velocities and the thickness of the top layer. We applied this equation to Rayleigh waves and multi-layer models with the average velocity and obtained encouraging results. This equation not only endows with a criterion to distinguish high modes from numerical artifacts but also provides a straightforward means to resolve the depth to the half space of a layered earth model. ?? 2005 Elsevier Ltd. All rights reserved.

Soil Dynamics and Earthquake Engineering↗

Numerical simulation of the paleohydrology of glacial Lake Oshkosh, eastern Wisconsin, USA

Proglacial lakes, formed during retreat of the Laurentide ice sheet, evolved quickly as outlets became ice-free and the earth deformed through glacial isostatic adjustment. With high-resolution digital elevation models (DEMs) and GIS methods, it is possible to reconstruct the evolution of surface hydrology. When a DEM deforms through time as predicted by our model of viscoelastic earth relaxation, the entire surface hydrologic system with its lakes, outlets, shorelines and rivers also evolves without requiring assumptions of outlet position. The method is applied to proglacial Lake Oshkosh in Wisconsin (13,600 to 12,900??cal yr BP). Comparison of predicted to observed shoreline tilt indicates the ice sheet was about 400??m thick over the Great Lakes region. During ice sheet recession, each of the five outlets are predicted to uplift more than 100??m and then subside approximately 30??m. At its maximum extent, Lake Oshkosh covered 6600??km2 with a volume of 111??km3. Using the Hydrologic Engineering Center-River Analysis System model, flow velocities during glacial outburst floods up to 9??m/s and peak discharge of 140,000??m3/s are predicted, which could drain 33.5??km3 of lake water in 10??days and transport boulders up to 3??m in diameter. ?? 2007 University of Washington.

Quaternary Research↗

Loess studies in central United States: Evolution of concepts

Few words in the realm of earth science have caused more debate than "loess". It is a common term that was first used as a name of a silt deposit before it was defined in a scientific sense. Because this "loose" deposit is easily distinguished from other more coherent deposits, it was recognized as a matter of practical concern and later became the object of much scientific scrutiny. Loess was first recognized along the Rhine Valley in Germany in the 1830s and was first noted in the United States in 1846 along the lower Mississippi River where it later became the center of attention. The use of the name eventually spread around the world, but its use has not been consistently applied. Over the years some interpretations and stratigraphic correlations have been validated, but others have been hotly contested on conceptual grounds and semantic issues. The concept of loess evolved into a complex issue as loess and loess-like deposits were discovered in different parts of the US. The evolution of concepts in the central US developed in four indefinite stages: the eras of (1) discovery and development of hypotheses, (2) conditional acceptance of the eolian origin of loess, (3) "bandwagon" popularity of loess research, and (4) analytical inquiry on the nature of loess. Toward the end of the first era around 1900, the popular opinion on the meaning of the term loess shifted from a lithological sense of loose silt to a lithogenetic sense of eolian silt. However, the dual use of the term fostered a lingering skepticism during the second era that ended in 1944 with an explosion of interest that lasted for more than a decade. In 1944, R.J. Russell proposed and H.N. Fisk defended a new non-eolian, property-based, concept of loess. The eolian advocates reacted with surprise and enthusiasm. Each side used constrained arguments to show their view of the problem, but did not examine the fundamental problem, which was not in the proofs of their hypothesis, but in the definition of the term. Between 1944 and about 1950, the debates about loess reached a maximum level of complexity. The main semantic problem was submersed in peripheral arguments about physical properties and genetic interpretations. The scholarly treatment of the subject by Fisk and Russell stimulated quality responses from a diversity of earth scientists interested in academic and applied studies, particularly geo-history, pedology, soil mechanics and stratigraphy. The long-lasting popularity of loess studies during the bandwagon era lasted to about 1970. By that time, the analytical and technical interests had attracted the mainstream into the fourth era with a focus beyond the old arguments. Although Fisk and Russell found themselves defending an unpopular theory, they stimulated a scientific interest in the late Quaternary history of the Mississippi Valley that may never be exceeded.

Engineering Geology↗

Global occurrence of tellurium-rich ferromanganese crusts and a model for the enrichment of tellurium

Hydrogenetic ferromanganese oxyhydroxide crusts (Fe-Mn crusts) precipitate out of cold ambient ocean water onto hard-rock surfaces (seamounts, plateaus, ridges) at water depths of about 400 to 4000 m throughout the ocean basins. The slow-growing (mm/Ma) Fe-Mn crusts concentrate most elements above their mean concentration in the Earth’s crust. Tellurium is enriched more than any other element (up to about 50,000 times) relative to its Earth’s crustal mean of about 1 ppb, compared with 250 times for the next most enriched element. We analyzed the Te contents for a suite of 105 bulk hydrogenetic crusts and 140 individual crust layers from the global ocean. For comparison, we analyzed 10 hydrothermal stratabound Mn-oxide samples collected from a variety of tectonic environments in the Pacific. In the Fe-Mn crust samples, Te varies from 3 to 205 ppm, with mean contents for Pacific and Atlantic samples of about 50 ppm and a mean of 39 ppm for Indian crust samples. Hydrothermal Mn samples have Te contents that range from 0.06 to 1 ppm. Continental margin Fe-Mn crusts have lower Te contents than open-ocean crusts, which is the result of dilution by detrital phases and differences in growth rates of the hydrogenetic phases. Correlation coefficient matrices show that for hydrothermal deposits, Te has positive correlations with elements characteristic of detrital minerals. In contrast, Te in open-ocean Fe-Mn crusts usually correlates with elements characteristic of the MnO 2 , carbonate fluorapatite, and residual biogenic phases. In continental margin crusts, Te also correlates with FeOOH associated elements. In addition, Te is negatively correlated with water depth of occurrence and positively correlated with crust thickness. Q-mode factor analyses support these relationships. However, sequential leaching results show that most of the Te is associated with FeOOH in Fe-Mn crusts and ≤10% is leached with the MnO 2 . Thermodynamic calculations indicate that Te occurs predominantly as H 5 TeO 6 − in ocean water. The speciation of Te in ocean water and charge balance considerations indicate that Te should be scavenged by FeOOH, which is in agreement with our leaching results. The thermodynamically more stable Te(IV) is less abundant by factors of 2 to 3.5 than Te(VI) in ocean water. This can be explained by preferential (not exclusive) scavenging of Te(IV) by FeOOH at the Fe-Mn crust surface and by Fe-Mn colloids in the water column. We propose a model in which the extreme enrichment of Te in Fe-Mn crusts is likely the result of an oxidation reaction on the surface of FeOOH. A similar oxidation process has been confirmed for Co, Ce, and Tl at the surface of MnO 2 in crusts, but has not been suggested previously to occur in association with FeOOH in Fe-Mn crusts. Mass-balance considerations indicate that ocean floor Fe-Mn deposits are the major sink for Te in the oceans. The concentration and redox chemistry of Te in the global ocean are likely controlled by scavenging on Fe-Mn colloids in the water column and Fe-Mn deposits on the ocean floor, as is also the case for Ce.

Geochimica et Cosmochimica Acta↗

The Landsat 7 mission: terrestrial research and applications for the 21st century

The Landsat Earth observation approach introduced in 1972 created a new way of monitoring land cover and land use globally. The Landsat 7 mission, successfully launched on April 15, 1999, continues those observations and demonstrates significant progress in precise numerical radiometry, spectral differentiation, and seasonally repetitive monitoring. Substantial improvements in calibration procedures, both prior to launch and during normal operations, have also been made to ensure long-term stability in the acquired spectral radiometry. Landsat 7 data acquisitions are being driven by a long-term data acquisition plan that was designed to ensure that substantially cloud-free, seasonal coverage would be recorded and archived in the US for all land areas of the globe. NASA competitively selected a Landsat Science Team, consisting of representatives from US universities and government agencies, to exploit the Landsat 7 record for global change research. This team is addressing the technical and analytical means to process and analyze the core of this observation record, and for the first time in the history of the Landsat mission, the technical and operational aspects of the mission are being driven by the goals of the US science community. The expected outcome of these efforts is a rapid improvement in understanding the Earth system, as well as conceptual knowledge that will underpin significant advancements in the application of this technology for commercial, operational, educational, and research purposes. Pathways to achieve effective Landsat continuity in the early decades of the 21st century are also being given careful attention, and there is no question that the lessons learned from the Landsat 7 mission will strongly influence these next-generation sensor systems.

Remote Sensing of Environment↗

35 - Strength and energetics of active fault zones

The strength of active fault zones, i.e., the shear stress level required to cause fault slip, is fundamental to understanding the physics of earthquakes and to assessing earthquake hazard. Although many researchers have concluded that fault zones are weak (shear stresses 10 MPa or less averaged between 0 and ˜20 km depth), others maintain that faults are strong (˜100 MPa average of an approximate linear increase with depth). Thus, despite 30 y of dedicated research, relevant data remain inconclusive and fault strength remains uncertain by an order of magnitude. In part, this is because the main source of energy release in earthquakes is at depths greater than 5 km, inaccessible to direct instrumental observation. Very large earthquakes rupture to the Earth's surface where direct observation of the shallow rupture process is possible. However, the rupture characteristics at shallow depth may differ from those at seismogenic depths. To date no great ( M ≥ 7.75) earthquakes have occurred within a network of modern strong-motion instruments, but the large, well-recorded 1999 earthquakes in Turkey and Taiwan, both M = 7.6, show that this data gap is rapidly being filled. Furthermore, shear stress changes at the earthquake source (and the resulting seismic waves) are nearly linear perturbations of the absolute stress field. Thus, an unknown absolute background stress does not greatly affect the basic characteristics of the observed low frequency seismic waves and the observed geodetic deformation. Other, less-direct data must then be used to infer the physical state and ambient stress levels on active faults. This chapter discusses available evidence and current ideas about fault zone strength and energetics. In our review we begin by outlining the general physical conditions prevailing in the Earth around active faults and summarize the generally agreed upon observational features of faulting and earthquake occurrence ( Section 2 ). In doing so we make as few assumptions as possible, because making uncertain assumptions may lead to logical inconsistencies and apparent paradoxes. For example, the absence of a measurable, frictionally generated heat flow anomaly near active faults, the so-called stress-heat flow paradox, rests upon assumptions that must somehow be incorrect. We wish to avoid such inconsistencies and to begin we introduce the observations and briefly state what they imply about fault zones. Subsequently, in Sections 3–8 , we place these observations in an interpretative context and show how they have been used to infer various measures of fault shear stress. Section 9 summarizes our assessment of the average shear stress state near faults, argues for the general importance of stress heterogeneity in faulting processes, and discusses its implications.

International Geophysics↗

Seismic velocity structure of the continental lithosphere from controlled source data

The purpose of this chapter is to provide a summary of the seismic velocity structure of the continental lithosphere, i.e., the crust and uppermost mantle. We define the crust as the outer layer of the Earth that is separated from the underlying mantle by the Mohorovičić discontinuity (Moho). We adopted the usual convention of defining the seismic Moho as the level in the Earth where the seismic compressional-wave ( P -wave) velocity increases rapidly or gradually to a value greater than or equal to 7.6 km sec −1 (Steinhart, 1967), defined in the data by the so-called “ Pn ” phase ( P-n ormal). Here we use the term uppermost mantle to refer to the 50–200+ km thick lithospheric mantle that forms the root of the continents and that is attached to the crust (i.e., moves with the continental plates). This summary has been preceded by 90 y of intense scientific activity. Mohorovičić (1910) was the first to publish an estimate for crustal thickness (54 km near Zagreb, Croatia), and to describe the seismically defined boundary between the crust and mantle that now bears his name (often shortened to “the Moho”). The fact that the oceanic crust is significantly thinner than continental crust (5 km versus about 40 km) was documented 40 y later (e.g., Hersey et al., 1952). Numerous later studies demonstrated that the continental crust varies in thickness from about 15 km to greater than 70 km beneath the Tibet plateau. Jarchow and Thompson (1989) provide a useful summary of early crustal studies, and Table 1 provides additional references. We emphasize results from active-source seismic refraction profiles that provide detailed P -wave velocity information. The shear-wave ( S -wave) structure of the crust and uppermost mantle, as determined by surface waves, is discussed, for example, by Ekström et al. (1997) and Ritzwoller et al. (1998).

International Geophysics↗

Lanthanide, yttrium, and zirconium anomalies in the Fire Clay coal bed, Eastern Kentucky

The Fire Clay coal bed in the Central Appalachian basin region contains a laterally-persistent tonstein that is found in the coal throughout most of its areal extent. The tonstein contains an array of minerals, including sanidine, ??-quartz, anatase and euhedral zircon, thhat constitutes strong evidence for a volcanic origin of the parting. For this study, five samples of the tonstein and four sets of coal samples underlying the tonstein were collected from five sites in eastern Kentucky. Inductively coupled plasma-mass spectroscopy (ICP-MS) analysis of the tonstein and underlying coal collected from four sites in eastern Kentucky show that although Zr concentrations are high in the tonstein (570-1820 ppm on a coal-ash basis (cab)), they are highest in the coal directly underlying the tonstein (2870-4540 ppm (cab)). A similar enrichment pattern is observed in the concentration of Y plus the sum of the rare earth elements (Y + ??REE): total Y + ??REE concentrations in the five tonstein samples range from 511 to 565 ppm (cab). However, Y + ??REE contents are highest in the coals directly underlying the tonsteins: values range from 1965 to 4198 ppm (cab). Scanning electron microscopy of samples from coal which directly underlies two of the tonstein samples show that REE-rich phosphate, tentatively identified as monazite, commonly infills cracks in clays and cells in clarain and vitrain. Zircon is rare and commonly subhedral. On the basis of coal chemistry and grain morphology, we suggest that volcanic components in the tonstein were leached by ground water. The leachate, rich in Y and REE precipitated as authigenic mineral phases in the underlying coal.The Fire Clay coal bed in the Central Appalachian basin region contains a laterally-persistent tonstein that is found in the coal throughout most of its areal extent. The tonstein contains an array of minerals, including sanidine, ??-quartz, anatase and euhedral zircon, that constitutes strong evidence for a volcanic origin of the parting. For this study, five samples of the tonstein and four sets of coal samples underlying the tonstein were collected from five sites in eastern Kentucky. Inductively coupled plasma-mass spectroscopy (ICP-MS) analysis of the tonstein and underlying coal collected from four sites in eastern Kentucky show that although Zr concentrations are high in the tonstein (570-1820 ppm on a coal-ash basis (cab)), they are highest in the coal directly underlying the tonstein (2870-4540 ppm (cab)). A similar enrichment pattern is observed in the concentration of Y plus the sum of the rare earth elements (Y+???REE): total Y+???REE concentrations in the five tonstein samples range from 511 to 565 ppm (cab). However, Y+???REE contents are highest in the coals directly underlying the tonsteins: values range from 1965 to 4198 ppm (cab). Scanning electron microscopy of samples from coal which directly underlies two of the tonstein samples show that REE-rich phosphate, tentatively identified as monazite, commonly infills cracks in clays and cells in clarain and vitrain. Zircon is rare and commonly subhedral. On the basis of coal chemistry and grain morphology, we suggest that volcanic components in the tonstein were leached by ground water. The leachate, rich in Y and REE precipitated as authigenic mineral phases in the underlying coal.

International Journal of Coal Geology↗

International collaboration: The cornerstone of satellite land remote sensing in the 21st century

Satellite land remotely sensed data are used by scientists and resource managers world-wide to study similar multidisciplinary earth science problems. Most of their information requirements can be met by a small number of satellite sensor types. Moderate-resolution resource satellites and low-resolution environmental satellites are the most prominent of these, and they are the focus of this paper. Building, launching, and operating satellite systems are very expensive endeavors. Consequently, nations should change the current pattern of independently launching and operating similar, largely redundant resource and environmental satellite systems in favor of true and full collaboration in developing, launching, operating, and sharing the data from such systems of the future. The past decade has seen encouraging signs of increasing international collaboration in earth remote sensing, but full collaboration has not yet been attempted. A general strategy to achieve such international collaboration is presented here, including discussion of potential obstacles, ideas for organizing and overseeing the long-term process toward collaboration, and short-term objectives whereby early successes critical to accomplishing long-term goals can be achieved.

Space Policy↗

Nano-metrology and terrain modelling - convergent practice in surface characterisation

The quantification of magnetic-tape and disk topography has a macro-scale counterpart in the Earth sciences - terrain modelling, the numerical representation of relief and pattern of the ground surface. The two practices arose independently and continue to function separately. This methodological paper introduces terrain modelling, discusses its similarities to and differences from industrial surface metrology, and raises the possibility of a unified discipline of quantitative surface characterisation. A brief discussion of an Earth-science problem, subdividing a heterogeneous terrain surface from a set of sample measurements, exemplifies a multivariate statistical procedure that may transfer to tribological applications of 3-D metrological height data.

Tribology International↗

Implications of crustal permeability for fluid movement between terrestrial fluid reservoirs

A classic paper by Rubey [Geol. Soc. Amer. Bull 62 (1951) 1111] examined various hypotheses regarding the origin of sea water and concluded that the most likely hypothesis was volcanic outgassing, a view that was generally accepted by Earth scientists for the next several decades. More recent work suggests that the rate of subduction of water is much larger than the volcanic outgassing rate, lending support to hypotheses that either ocean volume has decreased with time, or that the imbalance is offset by continuous replenishment of water by cometary impacts. These alternatives are required in the absence of additional mechanisms for the return of water from subducting lithosphere to the Earth's surface. Our recent work on crustal permeability suggests a large capacity for water upflow through tectonically active continental crust, resulting in a heretofore-unrecognized degassing pathway that can accommodate the water-subduction rate. Escape of recycled water via delivery from the mantle through zones of active metamorphism eliminates the mass-balance argument for the loss of ocean volume or extraterrestrial sources. ?? 2003 Elsevier Science B.V. All rights reserved.

Conference Paper↗

Geochemical and isotopic (Nd-Pb-Sr-O) variations bearing on the genesis of volcanic rocks from Vesuvius, Italy

Alkaline volcanism produced by Monte Somma-Vesuvius volcano includes explosive plinian and subplinian activity in addition to effusive lava flows. Pumice, scoria, and lava (150 samples) exhibit major- and trace-element gradients as a function of SiO2 (58.9-47.2 wt%) and MgO (0-7.8 wt%); Mg value are ???50. Internally gradational chemical groups or cycles are distinguished by age: (1) 25 000 to 14 000 yr B.P.; (2) 8000 yr B.P. to A.D. 79; and (3) A.D. 79 to 1944. A small number of lavas, dikes and scora were also analysed from the Somma formation (~ 35 000 to 25 000 yr B.P.). Within each group, contents of Na2O + K2O increas with decreasing MgO along distinct rocks. Nb/Y values are variable from 0.66 to 3.14 (at SiO2 ??? 50 wt%) generally in the range of alkaline and ultra-alkaline rocks. Variations in contents of some majro elements (e.g., P and Ti), and trace elements (e.g., Th, Nb, Ta, Zr, Hf, Pb, La, and Sc), as well as contrasting trends in ratios of various elements (e.g., Ta/Yb, Hf/U, Th/Ta, Th/Hf, Th/Yb, etc.) are also generally consistent with the group subdivisions. For example, Th/Hf increases from ??? 5 to ??? 10 with decreasing age for the Vesuvius system as a whole, yielding similar compositions in the least evolved rocks (low-silica, high-MgO, imcompatible element-poor) erupted at the end of each cycle. Internal variations within individual eruptions also systematically changed generally towards a common mafic composition at the end of each cycle, thus reflecting the dominanit volume in the magma chamber. At the start of a new eruptive cycle, the rocks are relatively enriched in incompatible elements; younger groups also contain higher abundances than other groups. N-MORB-normalized multielement diagrams exhibit selective enrichments of Sr, K, Rb, Th, and the light rare-earth elements; deep Nb and Ta negative anomalies commonly seen in rocks generated at orogenic margins are absent in the light rare-earth elements; deep Nb and Ta netgative anomalies commonly seen in rocks generated at orogenic margins are absent in our samples. Sr isotopic compositions are known to be variable within some of the units, in agreement with our data (87Sr/86Sr ~ 0.70699 to 0.70803) and with contributions from several isotopic components. Isotopic compositions for ??18O (7.3 to 10.2%), Pb for mineral separates and whole rocks (206Pb/204Pb ~ 18.947 to 19.178, 207Pb/204/Pb ~ 15.617 to 15.769, 208Pb/204Pb ~38.915 to 39.345), and Nd (143Nd ~ 0.51228 to 0.51251) also show variability. Oxygen isotope data show that pumices have higher ??18O values than cogenetic lavas, and that ??18O values and SiO2 are correlated. Radiogenic and stable isotope data plot within range of isotopic compositions for the Roman comagmatic province. Fractional crystallization cannot account for the radiogenic isotopic compositions of the Vesuvius magmas. We favor instead the combined effects of heterogeneous magma sources, together with isotopic exchange near the roof of the magma chamber. We suggest that metasomatized continental mantle lithosphere is the principal source of the magmas. This kind of enriched mantle was melted and reactivated in an area of continental extension (incipient rift setting) without direct reliance on contemporaneous subduction processes but possibly with input from mantle sources that resemble those that produce ocean island basalts.

Journal of Volcanology and Geothermal Research↗

Chapter 21 Western phosphate field - Depositional and economic deposit models

The Western Phosphate Field (WPF), composed of Permian marine sedimentary strata that cover over 300,000 km 2 in the middle Rocky Mountains of Idaho, Montana, Utah, and Wyoming in the United States, contains vast resources of phosphate mined for fertilizer and a range of other industrial applications. The richest deposits of phosphate in the WPF occur in the Meade Peak Phosphatic Shale Member of the Phosphoria Formation in southeast Idaho. Phosphate is an essential and even limiting nutrient of algal production, which occurs at the bottom of the marine food web in the oceanic photic zone. The high concentrations of phosphate and trace elements in the Phosphoria Formation reflect a low accumulation rate of diluting phases, such as terrigenous siliciclastic debris and carbonate, rather than an unusually high level of primary productivity at the time of deposition. Indeed, the mean rate of accumulation of PO 4 3 required a continuous flux of PO 4 3 into the basin and the photic zone of the water column, but only at a moderate level. This flux was maintained by upwelling of nutrient-rich seawater, imported at depth from the open ocean. Although only a fraction of the organic matter that hosted the PO 4 3 and other nutrients (NO 3 , Cd, Cu, Mo, Ni, and Zn) actually escaped oxidation in the water column, their rate of accu- mulation on the sea floor defined the basin hydrography. Rates of accumulation of Cr, U, V, and rare-earth elements by precipitation and adsorp- tion reactions identify redox conditions in the bottom water as having been denitrifying, maintained by a balance between the rate of oxidation of organic matter settling through the water column and the flux of open-ocean seawater at depth. Atmospheric mixing maintained oxygen respiration in the uppermost several tens of meters of the water column. This hydrography and seawater chemistry are present in several sedimentary envi- ronments in the ocean today. In the WPF, there is an estimated surface mineable reserve base and subeconomic resource of 7.6 billion mt, at an average grade of 24% P 2 O 5 ; a subeconomic underground- mineable resource of 17 billion mt, at a grade of 28%; and 507 billion mt of subresource- grade phosphatic material that underlie the WPF at a depth greater than 305 m. The relationship between phosphate-ore specifications and weathering suggests that significant changes in processing, with associated cost increases, will be required to extend recovery of ore below the relatively strongly weathered zone near the surface. Four open pit mines currently extract phosphate from two moderately to steeply dipping ore zones that typically contain between 20% and 35% P 2 O 5 . Although the shales are enriched in trace elements, especially As, Cd, Cr, Cu, Mo, Se, U, V, Zn, and rare-earth elements, the relative concentration of organic carbon and selected major element oxides determines the suitability of phosphate-rich rock for feed to processing plants and its other applications. Selected specifications from the four operating mines include the following: minimum P 2 O 5 of 18-20% and average of 26-27%; maximum A1 2 O 3 of 1.6-5.0%; maximum MgO of 0.3-0.6%; a CaO/P 2 O 5 ratio of 1.5-1.6; and total carbon content of 4%-5%. Weathering to a depth of as much as 100 m significantly enhances ore quality by decreasing the proportions of calcite, dolomite, and organic matter relative to carbonate fluorapatite, the primary ore mineral

Colorado, Idaho, Montana, Nebraska, North Dakota, ↗