USGS Science⌕ Search

SEARCH · USGS Science

Results for “Geosciences”

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 721 records · Page 40Linked to original sources

Perspectives on premetamorphic stratabound tourmalinites

Stratabound tourmalinites are metallogenically important rocks that locally show a close spatial association with diverse types of mineralization, especially volcanogenic massive sulfides (VMS) and clastic-dominated (CD) Zn-Pb deposits. These tourmalinite occurrences pan the geologic record from Eoarchean to Jurassic. Host lithologies are dominated by clastic metasedimentary rocks but in some areas include metavolcanic rocks, marble, or metaevaporites. Stratabound and stratiform (conformable) tourmalinites commonly display sedimentary structures such as graded beds, cross-beds, and rip-up clasts. In most cases, field and microtextural relationships are consistent with a synsedimentary to the early diagenetic introduction of boron as a precursor to tourmaline formation. Whole-rock geochemical data for major, trace, and rare earth elements (REE) provide valuable insights into tourmalinite origins. Al-normalized values relative to those for least-altered host metasedimentary rocks suggest that tourmalinites in proximal settings at or near hydrothermal vent sites characterized by high fluid/rock regimes (e.g., Sullivan Pb-Zn-Ag deposit, Canada) have very different signatures than those in low fluid/rock, distal settings (e.g., Broken Hill Pb-Zn-Ag deposit, Australia). The high fluid/rock regimes at Sullivan show large mass changes of +60 % for Mg and +180 % for Mn, as well as large variations in abundances of light and middle REE. In contrast, tourmalinite formation in low fluid/rock regimes yields minimal Al-normalized changes in major elements, trace elements, and REE. Boron isotope values of tourmalinite-hosted tourmaline vary widely from -26.1 to +27.5 ‰, and are attributed mainly to boron sources (e.g., sediments, evaporites) with generally minor influence from processes such as formational temperature, fluid/rock ratio, and secular variation in seawater δ 11 B values. Laterally extensive stratiform tourmalinites formed mainly by syngenetic or early diagenetic processes on or beneath the seafloor. The syngenetic process is attributed to the interaction of vented B-rich brines with aluminous minerals in sediments, whereas the diagenetic process involves the selective replacement of aluminous sediments by B-rich fluids. Modern examples of tourmalinites, as yet undiscovered, may exist in metalliferous sediments of the Red Sea and the eastern Pacific Ocean, in altered volcaniclastic sediments within active seafloor-hydrothermal systems of the South Pacific, and in hydrothermal mounds and vents associated with mafic sill complexes in extensional basins as in the North Sea and South China Sea. Stratabound tourmalinites that contain base-metal sulfides, high Mn concentrations (>1 wt. % MnO), or positive Eu anomalies can be valuable exploration guides for base-metal sulfide deposits in sedimentary and volcanic terranes.

Journal of Geosciences↗

Metallogeny, exploitation and environmental impact of the Mt. Amiata mercury ore district (Southern Tuscany, Italy)

The Mt. Amiata mining district (Southern Tuscany, Italy) is a world class Hg district, with a cumulate production of more than 100,000 tonnes of Hg, mostly occurring between 1870 and 1980. The Hg mineralization at Mt. Amiata is younger than 0.3 Ma, and is directly related to shallow hydrothermal systems similar to present-day geothermal fields of the region. There is likely a continuum of Hg deposition to present day, because Hg emission from geothermal power plants is on-going. In this sense, the Mt. Amiata deposits present some analogies with “hot-spring type” deposits of western USA, although an ore deposit model for the district has not been established. Specifically, the source of Hg remains highly speculative. The mineralizing hydrothermal fluids are of low temperature, and of essentially meteoric origin. Recent results by our research group indicate that, 30 years after mine closure, the environmental effects of Hg contamination related to mining are still recorded by the ecosystem, namely on waterways of the Paglia and Tiber River basins. In particular, the close spatial connection between the town of Abbadia San Salvatore, the Hg mine within its immediate neighborhood, and the drainage catchment of the Paglia River has an influence also on Hg speciation, transported mainly in the particulate form by the river system. The extent of Hg contamination has been identified at least 100 km from Abbadia San Salvatore along the Paglia-Tiber River system. Estimated annual Hg mass loads transported by the Paglia River to the Tiber River were about 11 kg yr −1 . However, there is evidence that flood events may enhance Hg mobilization in the Paglia River basin, increasing Hg concentrations in stream sediment. The high methyl-Hg/Hg ratio in water in this area is an additional factor of great concern due to the potential harmful effects on human and wildlife health. Results of our studies indicate that the Mt. Amiata region is at present a source of Hg of remarkable environmental concern at the local, regional (Tiber River), and Mediterranean scales. Ongoing studies are aimed to a more detailed quantification of the Hg mass load input to the Mediterranean Sea, and to unravel the processes concerning Hg transport and fluid dynamics.

Tuscany↗

Colony-forming unit spreadplate assay versus liquid culture enrichment-polymerase chain reaction assay for the detection of Bacillus Endospores in soils

A liquid culture enrichment-polymerase chain reaction (E-PCR) assay was investigated as a potential tool to overcome inhibition by chemical component, debris, and background biological impurities in soil that were affecting detection assay performance for soil samples containing Bacillus atrophaeus subsp. globigii (a surrogate for B. anthracis). To evaluate this assay, 9 g of matched sets of three different soil types (loamy sand [sand], sandy loam [loam] and clay) was spiked with 0, ~4.5, 45, 225, 675 and 1350 endospores. One matched set was evaluated using a previously published endospore concentration and colony-forming unit spreadplate (CFU-S) assay and the other matched set was evaluated using an E-PCR assay to investigate differences in limits of detection between the two assays. Data illustrated that detection using the CFU-S assay at the 45-endospore spike level started to become sporadic whereas the E-PCR assay produced repeatable detection at the ~4.5-endospore spike concentration. The E-PCR produced an ~2-log increase in sensitivity and required slightly less time to complete than the CFU-S assay. This study also investigated differences in recovery among pure and blended sand and clay soils and found potential activation of B. anthracis in predominately clay-based soils.

Geosciences↗

Using fission-track radiography coupled with scanning electron microscopy for efficient identification of solid-phase uranium mineralogy at a former uranium pilot mill (Grand Junction, Colorado)

At a former uranium pilot mill in Grand Junction, Colorado, mine tailings and some subpile sediments were excavated to various depths to meet surface radiological standards, but residual solid-phase uranium below these excavation depths still occurs at concentrations above background. The combination of fission-track radiography and scanning electron microscope energy-dispersive X-ray spectroscopy (SEM-EDS) provides a uniquely efficient and quantitative way of determining mineralogic associations of uranium that can influence uranium mobility. After the creation of sample thin sections, a mica sheet is placed on those thin sections and irradiated in a nuclear research reactor. Decay of the irradiated uranium creates fission tracks that can be viewed with a microscope. The fission-track radiography images indicate thin section sample areas with elevated uranium that are focus areas for SEM-EDS work. EDS spectra provide quantitative elemental data that indicate the mineralogy of individual grains or grain coatings associated with the fission-track identification of elevated uranium. For the site in this study, the results indicated that uranium occurred (1) with coatings of aluminum–silicon (Al/Si) gel and gypsum, (2) dispersed in the unsaturated zone associated with evaporite-type salts, and (3) sorbed onto organic carbon. The Al/Si gel likely formed when low-pH waters were precipitated during calcite buffering, which in turn retained or precipitated trace amounts of Fe, As, U, V, Ca, and S. Understanding these mechanisms can help guide future laboratory and field-scale efforts in determining long-term uranium release rates to groundwater.

Colorado↗

The Yorktown Formation: Improved stratigraphy, chronology and paleoclimate interpretations from the U.S. mid-Atlantic Coastal Plain

The Yorktown Formation records paleoclimate conditions along the mid-Atlantic Coastal Plain during the mid-Piacenzian Warm Period (3.264 to 3.025 Ma), a climate interval of the Pliocene in some ways analogous to near future climate projections. To gain insight into potential near future changes, we investigated Yorktown Formation outcrops and cores in southeastern Virginia, refining the stratigraphic framework. We analyzed 485 samples for alkenone-based sea surface temperature (SST) and productivity estimates from the Holland and Dory cores, an outcrop at Morgarts Beach, Virginia, and the lectostratotype of the Yorktown Formation at Rushmere, Virginia, and analyzed planktonic foraminferal assemblage data from the type section. Using the structure of the SST record, we improved the chronology of the Yorktown Formation by establishing the maximum age ranges of the Rushmere (3.3–3.2 Ma) and Morgarts Beach (3.2–3.15 Ma) Members. SST values for these members average ~26 °C, corroborating existing sclerochronological data. Increasing planktonic foraminifer abundance, productivity, and species diversity parallel increasing SST over the MIS M2/M1 transition. These records constitute the greatest temporal concentration of paleoecological estimates within the Yorktown Formation, aiding our understanding of western North Atlantic temperature patterns, seasonality and ocean circulation during this interval. We provide a chronologic framework for future studies analyzing ecological responses to profound climate change.

Virginia↗

Electrical properties and anisotropy of schists and fault rocks from New Zealand’s Southern Alps under confining pressure

Magnetotelluric models spanning the Pacific–Australian Plate boundary in New Zealand’s South Island indicate a localized zone of low electrical resistivity that is spatially coincident with theductile mid-crustal part of the Alpine Fault Zone (AFZ). We explored the source of this anomaly bymeasuring the electrical properties of samples collected from surface outcrops approaching the AFZthat have accommodated a gradient of systematic strain and deformation conditions. We investigatedthe effects of tectonite fabric, fluid saturated pore/fracture networks and surface conductivity on the bulk electrical response and the anisotropy of resistivity measured under increasing confining pressures up to 200 MPa. We find that porosity and resistivity increase while porosity and the change in anisotropy of resistivity with confining pressure (δ (ρk/ρ⊥)/δ (peff)) decreases approaching theAFZ, indicating the electrical response is controlled by pore fluid conductivity and modified during progressive metamorphism. Conversely, Alpine mylonites exhibit relatively low resistivities at low porosities, and lower δ (ρk/ρ⊥)/δ (peff) than the schists. These findings indicate a transition in both the porosity distribution and electrical charge transport processes in rocks that have experienced progressive grain size reduction and mixing of phases during development of mylonitic fabrics due to creep shear strain within the AFZ.

Southern Alps↗

The over-prediction of seismically induced soil liquefaction during the 2016 Kumamoto, Japan earthquake sequence

Following the M7.0 strike-slip earthquake near Kumamoto, Japan, in April of 2016, most geotechnical engineering experts believed that there would be significant soil liquefaction and liquefaction-induced infrastructure damage observed in the densely populated city of Kumamoto during the post-event engineering reconnaissance. This belief was driven by several factors including the young geologic environment, alluvially deposited soils, a predominance of loose sandy soils documented in publicly available boring logs throughout the region, and the high intensity ground motions observed from the earthquake. To the surprise of many of the researchers, soil liquefaction occurred both less frequently and less severely than expected. This paper summarizes findings from our field, laboratory, and simplified analytical studies common to engineering practice to assess the lower occurrence of liquefaction. Measured in situ SPT and CPT resistance values were evaluated with current liquefaction triggering procedures. Minimally disturbed samples were subjected to cyclic triaxial testing. Furthermore, an extensive literature review on Kumamoto volcanic soils was performed. Our findings suggest that current liquefaction triggering procedures over-predict liquefaction frequency and effects in alluvially deposited volcanic soils. Volcanic soils were found to possess properties of soil crushability, high fines content, moderate plasticity, and unanticipated organic constituents. Cyclic triaxial tests confirm the high liquefaction resistance of these soils. Moving forward, geotechnical engineers should holistically consider the soil’s mineralogy and geology before relying solely on simplified liquefaction triggering procedures when evaluating volcanic soils for liquefaction.

Mount Aso↗

The relative stability of planktic foraminifer thermal preferences over the past 3 million years

Stationarity of species’ ecological tolerances is a first-order assumption of paleoenvironmental reconstruction based upon analog methods. To test this and other assumptions used in quantitative analysis of foraminiferal faunas for paleoceanographic reconstruction, we analyzed paired alkenone unsaturation ratio ( U 37 K ′ )   "> U K ′ 37 ) 37′) sea surface temperature (SST) estimates and relative abundances of planktic foraminifera within Late Pliocene assemblages. We established Pliocene temperature preferences for nine species in the North Atlantic: Dentoglobigerina altispira, Globorotalia menardii, Globoconella puncticulata, Neogloboquadrina atlantica, Neogloboquadrina incompta, Neogloboquadrina pachyderma, Trilobatus sacculifer, Globigerinita glutinata, and Globigerina bulloides. We compared these to the temperature preferences of the same extant species, and in the three cases where the species are now extinct ( Dentoglobigerina altispira, Neogloboquadrina atlantica, and Globoconella puncticulata ), comparisons were made to either the descendant species or other modern species commonly used as analogs. In general, the taxa tested show similar temperature responses in both Late Pliocene and present-day (core-top) distributions. The data from these comparisons are mostly encouraging, supporting past paleoceanographic conclusions, and are otherwise valuable for testing previous taxonomic grouping decisions that are often necessary for interpreting the paleoenvironment based upon Pliocene foraminiferal assemblages.

Geosciences↗

Using an open-source tool to develop a three-dimensional hydrogeologic framework of the Kobo Valley, Ethiopia

Groundwater resource management requires understanding the groundwater basin’s hydrogeology and would be improved with the development of a three-dimensional hydrogeologic framework model (HFM). A wide range of methods and software exist to quantify the extent, structure, and properties of geologic systems. However, most geologic software is proprietary and cost-prohibitive for use in developing countries. GemPy is a Python-based, open-source (no-cost) tool for generating three-dimensional geological models. This study uses available data and GemPy to develop the Kobo Valley Hydrogeologic Framework Model (KV-HFM), a three-dimensional HFM for Kobo Valley in northern Ethiopia, which is part of the East African Rift System. The KV-HFM is a conceptual model that comprises the hydrostratigraphy, structural features, and hydraulic properties of the Kobo Valley groundwater system. The limited data described the extent and altitude of the hydrostratigraphic units using the GemPy implicit potential–field interpolation. The KV-HFM showed the existence of an east-to-west, structural-based groundwater divide composed of volcanic rock and clay. This divide splits the catchment into two groundwater systems with limited interconnected flow. This study illustrates the use of open-source software for developing an HFM using sparse, existing geologic data.

Afar Depression, Kobo Valley catchment↗

The mineral chemistry networks of tin and tungsten reflect metallogenic eras of the Mesozoic

Continental remobilization is a crucial driver of metallogenesis and the formation of ore deposits. Some of the world’s largest mineral deposits of the economically valuable elements tin (Sn), tungsten (W), gold (Au), copper (Cu), lead (Pb), and zinc (Zn) formed during the Mesozoic Era. Additionally, the chemistry and distribution of the elements Sn and W have been investigated in previous studies to understand planetary formation and differentiation processes. These two elements are largely co-located during certain South China Mesozoic metallogenic events but are not co-located during other time periods in the same regions. Here, we investigated the mineral chemistry network similarities and dissimilarities of Sn and W to understand their mineral formation and distribution during the Mesozoic Era and throughout Earth history. Mineral chemistry network community detection analysis and electronegativity associations among mineral constituent elements of Sn minerals and W minerals indicate that the elements have similar chemistry among their oxide minerals. However, Sn forms a much wider range of minerals that also contain S compared to W, which occurs in a limited number of S-containing minerals. The divergent constituent element interactions among S-containing Sn minerals and W minerals reflect the redox sensitivity and importance of oxygen (O) fugacity in Sn mineral formation. Conversely, extensive W mineral deposits are known to form at both high and low O fugacities. The similarities and differences between the mineral chemistry networks of Sn and W reflect the mineral distribution of the two elements in the Sn-W mineralization event from 160 to 139 Ma vs. the Sn–uranium (U) mineralization event from 125 to 98 million years ago (Ma). The mineral chemistry and distribution of Mesozoic Sn and W deposits illustrate the contrasting importance of redox and O fugacity on the mineral formation of different elements, and the dynamic crustal evolution that took place during this period of Earth history.

Geosciences↗

Structural evolution of the east Sierra Valley system (Owens Valley and vicinity), California: a geologic and geophysical synthesis

The tectonically active East Sierra Valley System (ESVS), which comprises the westernmost part of the Walker Lane-Eastern California Shear Zone, marks the boundary between the highly extended Basin and Range Province and the largely coherent Sierra Nevada-Great Valley microplate (SN-GVm), which is moving relatively NW. The recent history of the ESVS is characterized by oblique extension partitioned between NNW-striking normal and strike-slip faults oriented at an angle to the more northwesterly relative motion of the SN-GVm. Spatially variable extension and right-lateral shear have resulted in a longitudinally segmented valley system composed of diverse geomorphic and structural elements, including a discontinuous series of deep basins detected through analysis of isostatic gravity anomalies. Extension in the ESVS probably began in the middle Miocene in response to initial westward movement of the SN-GVm relative to the Colorado Plateau. At ca. 3-3.5 Ma, the SN-GVm became structurally separated from blocks directly to the east, resulting in significant basin-forming deformation in the ESVS. We propose a structural model that links high-angle normal faulting in the ESVS with coeval low-angle detachment faulting in adjacent areas to the east.

California↗

Anthrax and the geochemistry of soils in the contiguous United States

Soil geochemical data from sample sites in counties that reported occurrences of anthrax in wildlife and livestock since 2000 were evaluated against counties within the same states (MN, MT, ND, NV, OR, SD and TX) that did not report occurrences. These data identified the elements, calcium (Ca), manganese (Mn), phosphorus (P) and strontium (Sr), as having statistically significant differences in concentrations between county type (anthrax occurrence versus no occurrence). Tentative threshold values of the lowest concentrations of each of these elements (Ca = 0.43 wt %, Mn = 142 mg/kg, P = 180 mg/kg and Sr = 51 mg/kg) and average concentrations (Ca = 1.3 wt %, Mn = 463 mg/kg, P = 580 mg/kg and Sr = 170 mg/kg) were identified from anthrax-positive counties as prospective investigative tools in determining whether an outbreak had “potential” or was “likely” at any given geographic location in the contiguous United States.

Alabama, California, Colorado, Minnesota, Montana,↗

Climate change and future fire regimes: Examples from California

Climate and weather have long been noted as playing key roles in wildfire activity, and global warming is expected to exacerbate fire impacts on natural and urban ecosystems. Predicting future fire regimes requires an understanding of how temperature and precipitation interact to control fire activity. Inevitably this requires historical analyses that relate annual burning to climate variation. Fuel structure plays a critical role in determining which climatic parameters are most influential on fire activity, and here, by focusing on the diversity of ecosystems in California, we illustrate some principles that need to be recognized in predicting future fire regimes. Spatial scale of analysis is important in that large heterogeneous landscapes may not fully capture accurate relationships between climate and fires. Within climatically homogeneous subregions, montane forested landscapes show strong relationships between annual fluctuations in temperature and precipitation with area burned; however, this is strongly seasonal dependent; e.g., winter temperatures have very little or no effect but spring and summer temperatures are critical. Climate models that predict future seasonal temperature changes are needed to improve fire regime projections. Climate does not appear to be a major determinant of fire activity on all landscapes. Lower elevations and lower latitudes show little or no increase in fire activity with hotter and drier conditions. On these landscapes climate is not usually limiting to fires but these vegetation types are ignition-limited. Moreover, because they are closely juxtaposed with human habitations, fire regimes are more strongly controlled by other direct anthropogenic impacts. Predicting future fire regimes is not rocket science; it is far more complicated than that. Climate change is not relevant to some landscapes, but where climate is relevant, the relationship will change due to direct climate effects on vegetation trajectories, as well as by feedback processes of fire effects on vegetation distribution, plus policy changes in how we manage ecosystems.

Geoscience Canada↗

Spatio-temporal mapping of plate boundary faults in California using geodetic imaging

The Pacific–North American plate boundary in California is composed of a 400-km-wide network of faults and zones of distributed deformation. Earthquakes, even large ones, can occur along individual or combinations of faults within the larger plate boundary system. While research often focuses on the primary and secondary faults, holistic study of the plate boundary is required to answer several fundamental questions. How do plate boundary motions partition across California faults? How do faults within the plate boundary interact during earthquakes? What fraction of strain accumulation is relieved aseismically and does this provide limits on fault rupture propagation? Geodetic imaging, broadly defined as measurement of crustal deformation and topography of the Earth’s surface, enables assessment of topographic characteristics and the spatio-temporal behavior of the Earth’s crust. We focus here on crustal deformation observed with continuous Global Positioning System (GPS) data and Interferometric Synthetic Aperture Radar (InSAR) from NASA’s airborne UAVSAR platform, and on high-resolution topography acquired from lidar and Structure from Motion (SfM) methods. Combined, these measurements are used to identify active structures, past ruptures, transient motions, and distribution of deformation. The observations inform estimates of the mechanical and geometric properties of faults. We discuss five areas in California as examples of different fault behavior, fault maturity and times within the earthquake cycle: the M6.0 2014 South Napa earthquake rupture, the San Jacinto fault, the creeping and locked Carrizo sections of the San Andreas fault, the Landers rupture in the Eastern California Shear Zone, and the convergence of the Eastern California Shear Zone and San Andreas fault in southern California. These examples indicate that distribution of crustal deformation can be measured using interferometric synthetic aperture radar (InSAR), Global Navigation Satellite System (GNSS), and high-resolution topography and can improve our understanding of tectonic deformation and rupture characteristics within the broad plate boundary zone.

California↗

Probabilistic substrate classification with multispectral acoustic backscatter: A comparison of discriminative and generative models

We propose a probabilistic graphical model for discriminative substrate characterization, to support geological and biological habitat mapping in aquatic environments. The model, called a fully-connected conditional random field (CRF), is demonstrated using multispectral and monospectral acoustic backscatter from heterogeneous seafloors in Patricia Bay, British Columbia, and Bedford Basin, Nova Scotia. Unlike previously proposed discriminative algorithms, the CRF model considers both the relative backscatter magnitudes of different substrates and their relative proximities. The model therefore combines the statistical flexibility of a machine learning algorithm with an inherently spatial treatment of the substrate. The CRF model predicts substrates such that nearby locations with similar backscattering characteristics are likely to be in the same substrate class. The degree of allowable proximity and backscatter similarity are controlled by parameters that are learned from the data. CRF model results were evaluated against a popular generative model known as a Gaussian Mixture model (GMM) that doesn’t include spatial dependencies, only covariance between substrate backscattering response over different frequencies. Both models are used in conjunction with sparse bed observations/samples in a supervised classification. A detailed accuracy assessment, including a leave-one-out cross-validation analysis, was performed using both models. Using multispectral backscatter, the GMM model trained on 50% of the bed observations resulted in a 75% and 89% average accuracies in Patricia Bay and Bedford Basin, respectively. The same metrics for the CRF model were 78% and 95%. Further, the CRF model resulted in a 91% mean cross-validation accuracy across four substrate classes at Patricia Bay, and a 99.5% mean accuracy across three substrate classes at Bedford Basin, which suggest that the CRF model generalizes extremely well to new data. This analysis also showed that the CRF model was much less sensitive to the specific number and locations of bed observations than the generative model, owing to its ability to incorporate spatial autocorrelation in substrates. The CRF therefore may prove to be a powerful ‘spatially aware’ alternative to other discriminative classifiers.

Geosciences↗

Geometric versus anemometric surface roughness for a shallow accumulating snowpack

When applied to a snow-covered surface, aerodynamic roughness length, z 0 , is typically considered as a static parameter within energy balance equations. However, field observations show that z 0 changes spatially and temporally, and thus z 0 incorporated as a dynamic parameter may greatly improve models. To evaluate methods for characterizing snow surface roughness, we compared concurrent estimates of z 0 based on (1) terrestrial light detection and ranging derived surface geometry of the snowpack surface (geometric, z 0 G ) and (2) vertical wind profile measurements (anemometric, z 0 A ). The value of z 0 G was computed from Lettau’s equation and underestimated z 0 A but compared well when scaled by a factor of 2.34. The Counihan method for computing z 0 G was found to be unsuitable for estimating z 0 on a snow surface. During snowpack accumulation in early winter, z 0 varied as a function of the snow-covered area (SCA). Our results show that as the SCA increases, z 0 decreases, indicating there is a topographic influence on this relation.

Geosciences↗

Optimizing an inner-continental shelf geologic framework investigation through data repurposing and machine learning

The U.S. Geological Survey (USGS) and the National Oceanic Atmospheric Administration (NOAA) have collected approximately 5,400 km2 of geophysical and hydrographic data on the Atlantic continental shelf between Delaware and Virginia over the past decade and a half. Although originally acquired for different objectives, the comprehensive coverage and variety of data (bathymetry, backscatter, imagery and physical samples) presents an opportunity to merge collections and create high-resolution, broad-scale geologic maps of the seafloor. This compilation of data repurposes hydrographic data, expands the area of geologic investigation, highlights the versatility of mapping data, and creates new geologic products that would not have been independently possible. The data are classified using a variety of machine learning algorithms, including unsupervised and supervised methods. Four unique classes were targeted for classification, and source data include bathymetry, backscatter, slope, curvature, and shaded-relief. A random forest classifier used on all five source data layers was found to be the most accurate method for these data. Geomorphologic and sediment texture maps are derived from the classified acoustic data using over 200 groundtruth samples. The geologic data products can be used to identify sediment sources, inform resource management, link seafloor environments to sediment texture, improve our understanding of seafloor structure and sediment pathways, and demonstrate how ocean mapping resources can be useful beyond their original intent to maximize the footprint and scientific impact of a study.

Delaware, Maryland, Virginia↗

Bathymetry and geomorphology of Shelikof Strait and the western Gulf of Alaska

We defined the bathymetry of Shelikof Strait and the western Gulf of Alaska (WGOA) from the edges of the land masses down to about 7000 m deep in the Aleutian Trench. This map was produced by combining soundings from historical National Ocean Service (NOS) smooth sheets (2.7 million soundings); shallow multibeam and LIDAR (light detection and ranging) data sets from the NOS and others (subsampled to 2.6 million soundings); and deep multibeam (subsampled to 3.3 million soundings), single-beam, and underway files from fisheries research cruises (9.1 million soundings). These legacy smooth sheet data, some over a century old, were the best descriptor of much of the shallower and inshore areas, but they are superseded by the newer multibeam and LIDAR, where available. Much of the offshore area is only mapped by non-hydrographic single-beam and underway files. We combined these disparate data sets by proofing them against their source files, where possible, in an attempt to preserve seafloor features for research purposes. We also attempted to minimize bathymetric data errors so that they would not create artificial seafloor features that might impact such analyses. The main result of the bathymetry compilation is that we observe abundant features related to glaciation of the shelf of Alaska during the Last Glacial Maximum including abundant end moraines, some medial moraines, glacial lineations, eskers, iceberg ploughmarks, and two types of pockmarks. We developed an integrated onshore–offshore geomorphic map of the region that includes glacial flow directions, moraines, and iceberg ploughmarks to better define the form and flow of former ice masses.

Alaska↗