Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Ocean Modelling”

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 793 records · Page 44Linked to original sources

Stress rotation across the Cascadia megathrust requires a weak subduction plate boundary at seismogenic depths

The Mendocino Triple Junction region is the most seismically active part of the Cascadia Subduction Zone. The northward moving Pacific plate collides with the subducting Gorda plate causing intense internal deformation within it. Here we show that the stress field rotates rapidly with depth across the thrust interface from a strike-slip regime within the subducting plate, reflecting the Pacific plate collision, to a thrust regime in the overriding plate. We utilize a dense focal mechanism dataset, including observations from the Cascadia Initiative ocean bottom seismograph experiment, to constrain the stress orientations. To quantify the implications of this rotation for the strength of the plate boundary, we designed an inversion that solves for the absolute stress tensors in a three-layer model subject to assumptions about the strength of the subducting mantle. Our results indicate that the shear stress on the plate boundary fault is likely no more than about ∼50 MPa at ∼20 km depth. Regardless of the assumed mantle strength, we infer a relatively weak megathrust fault with an effective friction coefficient of ∼0 to 0.2 at seismogenic depths. Such a low value for the effective friction coefficient requires a combination of high fluid pressures and/or fault-zone minerals with low inherent friction in the region where a great earthquake is expected in Cascadia.

California↗

Constraints on the Pleistocene chronology of sediments from the Lomonosov Ridge

Despite its importance in the global climate system, age-calibrated marine geologic records reflecting the evolultion of glacial cycles through the Pleistocene are largely absent from the central Arctic Ocean. This is especially true for sediments older than 200 ka. Three sites cored during the Integrated Ocean Drilling Program's Expedition 302, the Arctic Coring Expedition (ACEX), provide a 27 m continuous sedimentary section from the Lomonosov Ridge in the central Arctic Ocean. Two key biostratigraphic datums and constraints from the magnetic inclination data are used to anchor the chronology of these sediments back to the base of the Cobb Mountain subchron (1215 ka). Beyond 1215 ka, two best fitting geomagnetic models are used to investigate the nature of cyclostratigraphic change. Within this chronology we show that bulk and mineral magnetic properties of the sediments vary on predicted Milankovitch frequencies. These cyclic variations record "glacial" and "interglacial" modes of sediment deposition on the Lomonosov Ridge as evident in studies of ice-rafted debris and stable isotopic and faunal assemblages for the last two glacial cycles and were used to tune the age model. Potential errors, which largely arise from uncertainties in the nature of downhole paleomagnetic variability, and the choice of a tuning target are handled by defining an error envelope that is based on the best fitting cyclostratigraphic and geomagnetic solutions. Copyright 2008 by the American Geophysical Union.

Paleoceanography↗

Advances in morphodynamic modeling of coastal barriers: A review

As scientific understanding of barrier morphodynamics has improved, so has the ability to reproduce observed phenomena and predict future barrier states using mathematical models. To use existing models effectively and improve them, it is important to understand the current state of morphodynamic modeling and the progress that has been made in the field. This manuscript offers a review of the literature regarding advancements in morphodynamic modeling of coastal barrier systems and summarizes current modeling abilities and limitations. Broadly, this review covers both event-scale and long-term morphodynamics. Each of these sections begins with an overview of commonly modeled phenomena and processes, followed by a review of modeling developments. After summarizing the advancements toward the stated modeling goals, we identify research gaps and suggestions for future research under the broad categories of improving our abilities to acquire and access data, furthering our scientific understanding of relevant processes, and advancing our modeling frameworks and approaches.

Virginia↗

Quantifying regional effects of best management practices on nutrient losses from agricultural lands

Nitrogen (N) and phosphorus (P) losses from agricultural areas have degraded the water quality of downstream rivers, lakes, and oceans. As a result, investment in the adoption of agricultural best management practices (BMPs) has grown, but assessments of their effectiveness at large spatial scales have lagged. This study applies regional Spatially Referenced Regression On Watershed-attributes (SPARROW) models developed for the Midwest, Northeast, and Southeast United States to quantify potential regional effects of BMPs on nutrient losses from agricultural lands. These models were used because they account for specific BMPs in the prediction of instream nutrient loads. The BMPs included in the models were cover crops, no-till, and conservation tillage. Sensitivity testing for the BMPs on agricultural nutrient loads was done using simulations that varied the intensity of BMPs specified in each region. When the BMP intensity was increased 50% relative to the 2012 intensity, the predicted agricultural load of total P decreased across all regions (4% to 14%). The predicted reduction in average P yields in the Midwest, Northeast, and Southeast was 706, 544, and 26 kg km–2, respectively. Increasing BMPs by 50% decreased predicted agricultural total N loads by 3.5% in the Southeast but increased predicted N loads in the Midwest and Northeast by 4.7% and 1.8%, respectively. Model-predicted average N yields increased by 402 kg km–2 and 302 kg km–2 in the Midwest and Northeast, respectively, and decreased in the Southeast by 329 kg km–2. In model simulations, cover crops were more effective at reducing N and P loads than the tillage BMPs despite lower intensity of implementation in 2012. However, at the regional scale of this investigation, implementation of BMPs result in only moderate predicted effects on agricultural nutrient loads

Journal of Soil and Water Conservation↗

Bathymetry of the Wilderness breach at Fire Island, New York, June 2013

The U.S. Geological Survey (USGS) St. Petersburg Coastal and Marine Science Center in St. Petersburg, Florida, collaborated with the U.S. Army Corps of Engineers Field Research Facility in Duck, North Carolina, to collect shallow water bathymetric data of the Wilderness breach on Fire Island, New York, in June 2013. The breach formed in October 2012 during Hurricane Sandy, and the USGS is involved in a post-Sandy effort to map, monitor, and model the morphologic evolution of the breach as part of Hurricane Sandy Supplemental Project GS2-2B: Linking Coastal Vulnerability and Process, Fire Island. This publication includes a bathymetric dataset of the breach and the adjacent nearshore on the ocean side of the island. The objective of the data collection and analysis is to map the bathymetry of the primary breach channel, ebb shoal, and nearshore bar system.

New York↗

Fifty-Year Flood-Inundation Maps for Santa Rosa de Aguan, Honduras

After the devastating floods caused by Hurricane Mitch in 1998, maps of the areas and depths of the 50-year-flood inundation at 15 municipalities in Honduras were prepared as a tool for agencies involved in reconstruction and planning. This report, which is one in a series of 15, presents maps of areas in the coastal municipality of Santa Rosa de Aguan that are prone to oceanic storm-surge flooding and wave action. The 50-year flood on the Rio Aguan (4,270 cubic meters per second), would inundate most of the area surveyed for this municipality and beyond. Therefore a detailed numerical hydraulic model was not developed for this municipality as it was for the others. The 50-year storm surge would likely produce higher water levels than the 50-year flood on the river during normal astronomical tides. The elevation of the 50-year storm surge was estimated to be 4.35 meters above normal sea level, based on hurricane probabilities and published storm-surge elevations associated with various hurricane categories. Flood-inundation maps, including areas of wave-action hazard and a color-shaded elevation map, were created from the available data and the estimated 50-year storm tide. Geographic Information System (GIS) coverages of the hazard areas are available on a computer in the municipality of Santa Rosa de Aguan as part of the Municipal GIS project and on the Internet at the Flood Hazard Mapping Data Web page (http://mitchnts1.cr.usgs.gov/projects/floodhazard.html). These coverages allow users to view the maps in much more detail than is possible using the maps in this report.

Santa Rosa de Aguan↗

Concentrations and annual fluxes for selected water-quality constituents from the USGS National Stream Quality Accounting Network (NASQAN) 1996-2000

This report contains concentrations and annual mass fluxes (loadings) for a broad range of water-quality constituents measured during 1996-2000 as part of the U.S. Geological Survey National Stream Quality Accounting Network (NASQAN). During this period, NASQAN operated a network of 40-42 stations in four of the largest river basins of the USA: the Colorado, the Columbia, the Mississippi (including the Missouri and Ohio), and the Rio Grande. The report contains surface-water quality data, streamflow data, field measurements (e.g. water temperature and pH), sediment-chemistry data, and quality-assurance data; interpretive products include annual and average loads, regression parameters for models used to estimate loads, sub-basin yield maps, maps depicting percent detections for censored constituents, and diagrams depicting flow-weighted average concentrations. Where possible, a regression model relating concentration to discharge and season was used for flux estimation. The interpretive context provided by annual loads includes identifying source and sink areas for constituents and estimating the loadings to receiving waters, such as reservoirs or the ocean.

Water-Resources Investigations Report↗

Preliminary-assessment and upgrade of a groundwater flow model of the Seacoast Bedrock Aquifer, New Hampshire

In 2003 and 2004, the U.S. Geological Survey investigated the availability of groundwater resources in a 160-square mile area of coastal New Hampshire (Figure 1) using a regional groundwater flow model (Mack, 2009). At that time, population growth and increasing water demand prompted concern for the sustainability of the region’s groundwater resources in a fractured-crystalline bedrock-aquifer with little storage. The groundwater flow model developed for the previous study incorporated detailed water-use information for 2003-4 and simulated the effects of projected increases in water use. However, poor stream representation may reduce the effectiveness of the original model head simulations. Improvements to the model, made by incorporating the USGS’s MODLFOW-2005 Newton formulation (MODFLOW-NWT, Niswonger and others, 2011) and by more accurately representing stream characteristics, are presented in an example simulating approximate changes in water use. Groundwater heads in an area of relatively larger population change, near the center of the Seacoast’s fractured bedrock aquifer, were simulated with the upgraded model using published 2004, and approximated 2015, water use rates. This area is situated at a local topographic high point and near the junction of three towns, where drainages flow westward, toward Great Bay, and eastward, toward the Atlantic Ocean (Figure 1).

New Hampshire↗

Dissolution rates of pure methane hydrate and carbon-dioxide hydrate in undersaturated seawater at 1000-m depth

To help constrain models involving the chemical stability and lifetime of gas clathrate hydrates exposed at the seafloor, dissolution rates of pure methane and carbon-dioxide hydrates were measured directly on the seafloor within the nominal pressure-temperature (P/T) range of the gas hydrate stability zone. Other natural boundary conditions included variable flow velocity and undersaturation of seawater with respect to the hydrate-forming species. Four cylindrical test specimens of pure, polycrystalline CH4 and CO2 hydrate were grown and fully compacted in the laboratory, then transferred by pressure vessel to the seafloor (1028 m depth), exposed to the deep ocean environment, and monitored for 27 hours using time-lapse and HDTV cameras. Video analysis showed diameter reductions at rates between 0.94 and 1.20 ??m/s and between 9.0 and 10.6 ?? 10-2 ??m/s for the CO2 and CH4 hydrates, respectively, corresponding to dissolution rates of 4.15 ?? 0.5 mmol CO2/m2s and 0.37 ?? 0.03 mmol CH4/m2s. The ratio of the dissolution rates fits a diffusive boundary layer model that incorporates relative gas solubilities appropriate to the field site, which implies that the kinetics of the dissolution of both hydrates is diffusion-controlled. The observed dissolution of several mm (CH4) or tens of mm (CO2) of hydrate from the sample surfaces per day has major implications for estimating the longevity of natural gas hydrate outcrops as well as for the possible roles of CO2 hydrates in marine carbon sequestration strategies. ?? 2003 Elsevier Ltd.

Geochimica et Cosmochimica Acta↗

DEM, tide and velocity over sulzberger ice shelf, West Antarctica

Arctic and Antarctic ice sheets preserve more than 77% of the global fresh water and could raise global sea level by several meters if completely melted. Ocean tides near and under ice shelves shifts the grounding line position significantly and are one of current limitations to study glacier dynamics and mass balance. The Sulzberger ice shelf is an area of ice mass flux change in West Antarctica and has not yet been well studied. In this study, we use repeat-pass synthetic aperture radar (SAR) interferometry data from the ERS-1 and ERS-2 tandem missions for generation of a high-resolution (60-m) Digital Elevation Model (DEM) including tidal deformation detection and ice stream velocity of the Sulzberger Ice Shelf. Other satellite data such as laser altimeter measurements with fine foot-prints (70-m) from NASA's ICESat are used for validation and analyses. The resulting DEM has an accuracy of-0.57??5.88 m and is demonstrated to be useful for grounding line detection and ice mass balance studies. The deformation observed by InSAR is found to be primarily due to ocean tides and atmospheric pressure. The 2-D ice stream velocities computed agree qualitatively with previous methods on part of the Ice Shelf from passive microwave remote-sensing data (i.e., LANDSAT). ?? 2005 IEEE.

Conference Paper↗

Slab melting and magma formation beneath the southern Cascade arc

The processes that drive magma formation beneath the Cascade arc and other warm-slab subduction zones have been debated because young oceanic crust is predicted to largely dehydrate beneath the forearc during subduction. In addition, geochemical variability along strike in the Cascades has led to contrasting interpretations about the role of volatiles in magma generation. Here, we focus on the Lassen segment of the Cascade arc, where previous work has demonstrated across-arc geochemical variations related to subduction enrichment, and H-isotope data suggest that H 2 O in basaltic magmas is derived from the final breakdown of chlorite in the mantle portion of the slab. We use naturally glassy, olivine-hosted melt inclusions (MI) from the tephra deposits of eight primitive ( MgO>7 wt% ) basaltic cinder cones to quantify the pre-eruptive volatile contents of mantle-derived melts in this region. The melt inclusions have B concentrations and isotope ratios that are similar to mid-ocean ridge basalt (MORB), suggesting extensive dehydration of the downgoing plate prior to reaching sub-arc depths and little input of slab-derived B into the mantle wedge. However, correlations of volatile and trace element ratios (H 2 O/Ce, Cl/Nb, Sr/Nd) in the melt inclusions demonstrate that geochemical variability is the result of variable addition of a hydrous subduction component to the mantle wedge. Furthermore, correlations between subduction component tracers and radiogenic isotope ratios show that the subduction component has less radiogenic Sr and Pb than the Lassen sub-arc mantle, which can be explained by melting of subducted Gorda MORB beneath the arc. Agreement between pMELTS melting models and melt inclusion volatile, major, and trace element data suggests that hydrous slab melt addition to the mantle wedge can produce the range in primitive compositions erupted in the Lassen region. Our results provide further evidence that chlorite-derived fluids from the mantle portion of the slab ( ∼7–9 km below the slab top) cause flux melting of the subducted oceanic crust, producing hydrous slab melts that migrate into the overlying mantle, where they react with peridotite to induce further melting.

Cascades↗

The distribution of meteoric 36Cl/Cl in the United States: A comparison of models

The natural distribution of 36Cl/Cl in groundwater across the continental United States has recently been reported by Davis et al. (2003). In this paper, the large-scale processes and atmospheric sources of 36Cl and chloride responsible for controlling the observed 36Cl/Cl distribution are discussed. The dominant process that affects 36Cl/Cl in meteoric groundwater at the continental scale is the fallout of stable chloride from the atmosphere, which is mainly derived from oceanic sources. Atmospheric circulation transports marine chloride to the continental interior, where distance from the coast, topography, and wind patterns define the chloride distribution. The only major deviation from this pattern is observed in northern Utah and southern Idaho where it is inferred that a continental source of chloride exists in the Bonneville Salt Flats, Utah. In contrast to previous studies, the atmospheric flux of 36Cl to the land surface was found to be approximately constant over the United States, without a strong correlation between local 36Cl fallout and annual precipitation. However, the correlation between these variables was significantly improved (R 2=0.15 to R 2=0.55) when data from the southeastern USA, which presumably have lower than average atmospheric 36Cl concentrations, were excluded. The total mean flux of 36Cl over the continental United States and total global mean flux of 36Cl are calculated to be 30.5??7.0 and 19.6??4.5 atoms m-2 s-1, respectively. The 36Cl/Cl distribution calculated by Bentley et al. (1996) underestimates the magnitude and variability observed for the measured 36Cl/Cl distribution across the continental United States. The model proposed by Hainsworth (1994) provides the best overall fit to the observed 36Cl/Cl distribution in this study. A process-oriented model by Phillips (2000) generally overestimates 36Cl/Cl in most parts of the country and has several significant local departures from the empirical data.

Hydrogeology Journal↗

Development of the benethic nepheloid layer on the south Texas continental shelf, western Gulf of Mexico

A monitoring study of suspended sediment on the South Texas Continental Shelf indicates that a turbid benthic nepheloid layer is regionally persistent. A sequence of quasi-synoptic measurements of the water column obtained during six cruises in an 18-month period indicates substantial spatial and temporal variability in nepheloidlayer characteristics. Regionally, the thickness of the shelf nepheloid layer increases both seaward and in a convergent alongshelf direction. Greatest thicknesses occur over a muddy substrate, indicating a causal relationship; maximum observed local thickness is 35 m which occurs along the southern shelf break. Analyses of suspended particulate matter in shelf bottom waters indicate mean concentrations ranging from 49 · 10 4 to 111 · 10 4 particle counts/cc; concentrations persistently increase shoreward throughout the region. Bottom particulate matter is predominantly composed of inorganic detritus. Admixtures of organic skeletal particles, primarily diatoms, are generally present but average less than 10% of the total particulate composition. Texturally, the particulate matter in bottom waters is predominantly poorly sorted sediment composed of very fine silt (3.9–7.8 μm). The variability in nepheloid-layer characteristics indicates a highly dynamic shelf feature. The relationship of nepheloid-layer characteristics to hydrographic and substrate conditions suggests a conceptual model whereby nepheloid-layer development and maintenance are the results of the resuspension of sea-floor sediment. Bottom turbulence is attributed primarily to vertical shear and shoaling progressive internal waves generated by migrating shelf-water masses, especially oceanic frontal systems, and secondarily to shoaling surface gravity waves.

Texas↗

Rapid cooling rates at an active mid-ocean ridge from zircon thermochronology

Oceanic spreading ridges are Earth's most productive crust generating environment, but mechanisms and rates of crustal accretion and heat loss are debated. Existing observations on cooling rates are ambiguous regarding the prevalence of conductive vs. convective cooling of lower oceanic crust. Here, we report the discovery and dating of zircon in mid-ocean ridge dacite lavas that constrain magmatic differentiation and cooling rates at an active spreading center. Dacitic lavas erupted on the southern Cleft segment of the Juan de Fuca ridge, an intermediate-rate spreading center, near the intersection with the Blanco transform fault. Their U–Th zircon crystallization ages (29.3 − 4.6 + 4.8 ka; 1σ standard error s.e.) overlap with the (U–Th)/He zircon eruption age (32.7 ± 1.6 ka) within uncertainty. Based on similar 238 U− 230 Th disequilibria between southern Cleft dacite glass separates and young mid-ocean ridge basalt (MORB) erupted nearby, differentiation must have occurred rapidly, within ~ 10–20 ka at most. Ti-in-zircon thermometry indicates crystallization at 850–900 °C and pressures > 70–150 MPa are calculated from H 2 O solubility models. These time-temperature constraints translate into a magma cooling rate of ~ 2 × 10 − 2 °C/a. This rate is at least one order-of-magnitude faster than those calculated for zircon-bearing plutonic rocks from slow spreading ridges. Such short intervals for differentiation and cooling can only be resolved through uranium-series ( 238 U– 230 Th) decay in young lavas, and are best explained by dissipating heat convectively at high crustal permeability.

Juan de Fuca ridge and the Blanco transform fault ↗

Restoration of Circum-Arctic Upper Jurassic source rock paleolatitude based on crude oil geochemistry

Tectonic geochemical paleolatitude (TGP) models were developed to predict the paleolatitude of petroleum source rock from the geochemical composition of crude oil. The results validate studies designed to reconstruct ancient source rock depositional environments using oil chemistry and tectonic reconstruction of paleogeography from coordinates of the present day collection site. TGP models can also be used to corroborate tectonic paleolatitude in cases where the predicted paleogeography conflicts with the depositional setting predicted by the oil chemistry, or to predict paleolatitude when the present day collection locality is far removed from the source rock, as might occur due to long distance subsurface migration or transport of tarballs by ocean currents. Biomarker and stable carbon isotope ratios were measured for 496 crude oil samples inferred to originate from Upper Jurassic source rock in West Siberia, the North Sea and offshore Labrador. First, a unique, multi-tiered chemometric (multivariate statistics) decision tree was used to classify these samples into seven oil families and infer the type of organic matter, lithology and depositional environment of each organofacies of source rock [Peters, K.E., Ramos, L.S., Zumberge, J.E., Valin, Z.C., Scotese, C.R., Gautier, D.L., 2007. Circum-Arctic petroleum systems identified using decision-tree chemometrics. American Association of Petroleum Geologists Bulletin 91, 877-913]. Second, present day geographic locations for each sample were used to restore the tectonic paleolatitude of the source rock during Late Jurassic time (???150 Ma). Third, partial least squares regression (PLSR) was used to construct linear TGP models that relate tectonic and geochemical paleolatitude, where the latter is based on 19 source-related biomarker and isotope ratios for each oil family. The TGP models were calibrated using 70% of the samples in each family and the remaining 30% of samples were used for model validation. Positive relationships exist between tectonic and geochemical paleolatitude for each family. Standard error of prediction for geochemical paleolatitude ranges from 0.9?? to 2.6?? of tectonic paleolatitude, which translates to a relative standard error of prediction in the range 1.5-4.8%. The results suggest that the observed effect of source rock paleolatitude on crude oil composition is caused by (i) stable carbon isotope fractionation during photosynthetic fixation of carbon and (ii) species diversity at different latitudes during Late Jurassic time. ?? 2008 Elsevier Ltd. All rights reserved.

Organic Geochemistry↗

Warming of the Indian Ocean threatens eastern and southern African food security but could be mitigated by agricultural development

Since 1980, the number of undernourished people in eastern and southern Africa has more than doubled. Rural development stalled and rural poverty expanded during the 1990s. Population growth remains very high, and declining per-capita agricultural capacity retards progress toward Millennium Development goals. Analyses of in situ station data and satellite observations of precipitation have identified another problematic trend: main growing-season rainfall receipts have diminished by ???15% in food-insecure countries clustered along the western rim of the Indian Ocean. Occurring during the main growing seasons in poor countries dependent on rain-fed agriculture, these declines are societally dangerous. Will they persist or intensify? Tracing moisture deficits upstream to an anthropogenically warming Indian Ocean leads us to conclude that further rainfall declines are likely. We present analyses suggesting that warming in the central Indian Ocean disrupts onshore moisture transports, reducing continental rainfall. Thus, late 20th-century anthropogenic Indian Ocean warming has probably already produced societally dangerous climate change by creating drought and social disruption in some of the world's most fragile food economies. We quantify the potential impacts of the observed precipitation and agricultural capacity trends by modeling 'millions of undernourished people' as a function of rainfall, population, cultivated area, seed, and fertilizer use. Persistence of current tendencies may result in a 50% increase in undernourished people by 2030. On the other hand, modest increases in per-capita agricultural productivity could more than offset the observed precipitation declines. Investing in agricultural development can help mitigate climate change while decreasing rural poverty and vulnerability. ?? 2008 by The National Academy of Sciences of the USA.

Proceedings of the National Academy of Sciences of↗

Developing a probabilistic tsunami hazard assessment framework for Pacific sources: USGS Powell Center meeting summary

Multi-organizational principal investigators formed a U.S. Geological Survey (USGS) Powell Center Working Group (WG), Tsunami Source Standardization for Hazards Mitigation in the United States, to develop a comprehensive series of sources capable of generating tsunamis that could impact U.S. state and territory coastal areas using probabilistic tsunami hazard analysis (PTHA). PTHA results are commonly used to provide consistent tsunami hazard information for use in engineering and risk assessment and, to a lesser extent, hazard response planning. Following an initial weeklong planning meeting in April 2018, designed to establish the WG’s scope, a series of weeklong meetings devoted to aspects of tsunami hazards placed emphasis on assessment of various tsunami sources, including subduction zones in Alaska, the Atlantic and Caribbean, Cascadia, and the Pacific Basin, as well as landslides in Alaska, the Atlantic, and the Caribbean. The final meeting in the series discussed tsunami sources from crustal faults. These meetings, each with a regional geographic focus, were designed to incorporate reviews and feedback from subject matter experts (SMEs) in academia, private industry, and federal, state, and local governmental organizations. Incorporating consensus from SMEs is important because the results derived from the tsunami source models will be used to inform the public about potential hazards from tsunamis related to safety concerns. This paper describes the USGS Powell Center meeting in March 2023, devoted specifically to developing a PTHA for tsunami sources in the Pacific Ocean Basin other than the Alaska–Aleutian and Cascadia subduction zones that were addressed during previous WG meetings.

Environmental & Engineering Geoscience↗

Marine benthic habitat mapping of Muir Inlet, Glacier Bay National Park and Preserve, Alaska, with an evaluation of the Coastal and Marine Ecological Classification Standard III

Seafloor geology and potential benthic habitats were mapped in Muir Inlet, Glacier Bay National Park and Preserve, Alaska, using multibeam sonar, ground-truth information, and geological interpretations. Muir Inlet is a recently deglaciated fjord that is under the influence of glacial and paraglacial marine processes. High glacially derived sediment and meltwater fluxes, slope instabilities, and variable bathymetry result in a highly dynamic estuarine environment and benthic ecosystem. We characterize the fjord seafloor and potential benthic habitats using the Coastal and Marine Ecological Classification Standard (CMECS) recently developed by the National Oceanic and Atmospheric Administration (NOAA) and NatureServe. Substrates within Muir Inlet are dominated by mud, derived from the high glacial debris flux. Water-column characteristics are derived from a combination of conductivity temperature depth (CTD) measurements and circulation-model results. We also present modern glaciomarine sediment accumulation data from quantitative differential bathymetry. These data show Muir Inlet is divided into two contrasting environments: a dynamic upper fjord and a relatively static lower fjord. The accompanying maps represent the first publicly available high-resolution bathymetric surveys of Muir Inlet. The results of these analyses serve as a test of the CMECS and as a baseline for continued mapping and correlations among seafloor substrate, benthic habitats, and glaciomarine processes.

Scientific Investigations Map↗