USGS Science⌕ Search

SEARCH · USGS Science

Results for “Coastal Futures”

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

Geological processes and sedimentation rates of wind-tidal flats, Laguna Madre, Texas

Coastal flats worldwide that are periodically exposed to arid climates and periodically flooded by marine waters are unique depositional environments because they receive sediments surficially and interstitially from both land and sea. The wind-tidal flats bordering Laguna Madre, Texas, which fit this unique category, are modified by eolian processes when subaerially exposed, and by wave and current processes when submerged. Floodwater is derived from the lagoon and driven onto the flats by strong and persistent winds during the passage of cold fronts and tropical cyclones. Low surface gradients of the flats prevent rapid drainage and promote seawater evaporation. The depositional products of these processes are interbedded and interlaminated sand, mud, marine shells, algal mats, and evaporites. This assemblage of sediments is geologically diagnostic evidence for intertidal marine deposition and the same assemblage of sediments have been reported for modern marginal-marine flats in the Middle East. The wind-tidal flat surface at Laguna Madre is constantly changing. However, the net effect of natural changes during the past century is either negligible or the changes occur at such a slow rate that they are almost imperceptible. Sediments are repeatedly added to and removed from the surface of the flats in minor increments and in different areas at different times. Preservation potential is enhanced at a particular site by the development of thick mats of blue-green algae. The 14 C ages of buried algal mats yield average long-term (centuries to millennia) sedimentation rates for the wind-tidal flats that range from 0.13 to 0.96 mm/yr and average 0.57 mm/yr. The 210 Pb profiles yield average short-term (150 yr) sedimentation rates that are an order of magnitude higher, ranging from 0.7 to 8.3 mm/yr and averaging 2.9 mm/yr. The minimum present rate of relative sea-level rise in Laguna Madre (3.4 mm/yr) exceeds the historical sedimentation rates for most of the flats. If future sea-level rise is faster than the rates of aggradation, then the wind-tidal flats will progressively become more frequently flooded and will eventually become permanently submerged.

Texas↗

Oyster model inventory: Identifying critical data and modeling approaches to support restoration of oyster reefs in coastal U.S. Gulf of Mexico waters

Executive Summary Along the coast of the U.S. Gulf of Mexico, the eastern oyster ( Crassostrea virginica ) plays important ecological and economic roles. Commercial landings from this region account for more than 50 percent of all U.S. landings; these oyster reefs also provide varied ecosystem services, including nursery habitat for many fish and macroinvertebrate species, shoreline protection, and water-quality maintenance. Declining trends in both total oyster production and functional reef area across this region have spurred investment in restoration of oyster resources, with specific calls for restoration projects to develop a network of reefs and identify broodstock and sanctuary reef restoration sites. Decision making related to restoration and establishment of a network of oyster reefs in the Gulf of Mexico requires information on both the environment and the effects of the environment on the oyster life cycle (including larval movement, survival, oyster recruitment, reproduction, growth, and mortality). Here, we examined the current state of data and model development in this region with the goal of providing an overview of oyster modeling approaches and an inventory of available data and existing oyster models. This report is meant to provide an overview to managers for understanding existing efforts and identify a path forward to most efficiently inform oyster resource management and restoration planning in moving from a single reef management approach to a reef network management approach. Numerous models related to some aspect of the oyster life cycle have been built, calibrated, and validated for various Gulf of Mexico estuaries over the last few decades (over 30 models identified). These models, which could inform site restoration, can be classified into four approaches: (1) oyster Habitat Suitability Index (HSI) models; (2) larval transport models; (3) on-reef oyster models that may include oyster growth, mortality and reproduction, and substrate persistence; and (4) coupled larval transport on-reef metapopulation models that simulate the entire oyster life cycle. The data requirements, model complexity and assumptions, and transferability vary by approach. Specifically, some approaches may offer greater accessibility, flexibility, and transferability spatially or temporally, with minimal data input, but only provide broad information to support site selection. In contrast, other approaches may require significant site-specific data for their construction and validation but may provide more accurate and location-specific data to support site selection for broodstock reefs. Regardless of modeling approach used, data on environmental drivers, such as salinity, water temperature, or water flow impacting oyster metabolism and movement, are required at appropriate spatial and temporal scales. While numerous data collection platforms, environmental models, and research products exist within Gulf of Mexico estuaries to provide important environmental data to use as drivers in the oyster models, significant variability in temporal and spatial coverage of the data, and variation in the availability of future condition models, exists across estuaries. This variation influences the spatial and temporal scales at which oyster models may be developed and impacts the calibration and validation of the oyster models within a given estuary, affecting its potential ability to address specific management or restoration questions. While multiple modeling approaches exist for informing site selection of broodstock or sanctuary oyster reefs, the development, calibration, and validation of a single modeling platform presents the most efficient, transferable, and useful tool for managers across the Gulf of Mexico. The development of a single modeling platform would involve using standardized input variables, governing equations, and assumptions for the modeled oyster processes and outputs, and for standardized calibration and validation procedures that could be applied within each estuary. The differences among estuary applications would require substituting only estuary-specific environmental data, and calibrating and validating the modeling approach with local oyster data. Two modeling approaches likely to be useful include (1) development of a general geospatial HSI modeling framework that could be applied consistently across estuaries and (2) a mechanistic coupled larval transport on-reef metapopulation model requiring only estuarine specific calibration and hydrodynamic models. Both approaches benefit from existing work across multiple Gulf of Mexico estuaries and could provide valuable support for oyster restoration, but may differ in their ability to address specific questions related to oyster restoration. HSI models specifically guide restoration practitioners in determining suitable habitat based on available data. The HSI approach, while currently more widely used and accessible, requires more development of larval suitability and larval input and output components in order to inform reef connectivity. A metapopulation approach considering the full oyster life cycle that simulates both on-reef oyster growth, mortality, reproduction, substrate persistence, and larval transport (ideally with larval growth and mortality) would provide the greatest detail and level of understanding but requires significant up-front investment. The larval oyster model and on-reef oyster model are usually developed independently for systems, although the two approaches can be coupled to represent the entire oyster life cycle in order to characterize and assess a reef metapopulation. This approach may be less accessible and much more data-intensive, however, and it requires some expertise to run and apply to inform oyster resource management. Ultimately, the development of single modeling platforms for each of these approaches would provide flexible tools applicable across all Gulf of Mexico oyster supporting estuaries. By using a single platform for model development, testing, calibrating and validating, and evaluation of modeled future scenarios, oyster restoration scientists and managers would not only be able to examine different scenario outcomes within a single estuary, but could also have comparable modeled results to evaluate potential outcomes, across estuaries and regions, that are not confounded by varying modeled data inputs, governing equations, assumptions, or user judgement.

Alabama, Florida, Louisiana, Mississippi, Texas↗

Hyperspectral (VNIR-SWIR) analysis of roll front uranium host rocks and industrial minerals from Karnes and Live Oak Counties, Texas Coastal Plain

VNIR-SWIR (400–2500 nm) reflectance measurements were made on the surfaces of various cores, cuttings and sample splits of sedimentary rocks from the Tertiary Jackson Group, and Catahoula, Oakville and Goliad Formations. These rocks vary in composition and texture from mudstone and claystone to sandstone and are known host rocks for roll front uranium occurrences in Karnes and Live Oak Counties, Texas. Spectral reflectance profiles, 569 in total, were reduced to 125 representative spectral signatures, which were analyzed using the U.S. Geological Survey's (USGS) Material Identification and Characterization Algorithm (MICA). MICA uses an automated continuum-removal procedure together with a least-squares linear regression to determine the fit of observed sample spectral absorption features to those of reference mineral standards in a spectral library. The reference minerals include various clay, mica, carbonate, ferric and ferrous iron minerals and their mixtures. In addition, absorption feature band-depth analysis was done to identify rock surfaces exhibiting absorption features related to uranium and zeolite minerals, which were not included in the command files used to execute MICA. Rocks from each of the four geologic units produced broadly similar spectral signatures as a result of comparable mineral compositions, but there were some notable differences. For example, Ca- and Na-montmorillonite was matched most frequently to the spectral absorption features in 2-μm (∼2000–2500 nm) wavelengths, while goethite occurred often at 1-μm (∼400–1000 nm) wavelengths. The latter is related to limonitic iron-staining in and around oxidized zones of the uranium roll front as described in previous papers. Rocks of the Jackson Group differed from those of the Catahoula, Oakville and Goliad units in that the former exhibited spectral features we interpret as being due to the presence of lignite-bearing mudstone layers. Goliad rocks exhibit spectral features related to dolomite, gypsum, anhydrite, and an unidentified green clay mineral that is possibly glauconite . Jackson Group rocks also exhibit weak but well-resolved absorption features at 964 and 1157 nm related to either or both zeolite minerals clinoptilolite and heulandite. These zeolite minerals and a few spectra exhibiting hydrous silica absorption features are indicative of alteration of volcanic glass in tuffaceous mudstone and claystone layers. A few sample spectra exhibited strong absorption features at around 1135 nm related to the uranium mineral coffinite. Both the 1135 nm coffinite and 1157 nm zeolite absorption features overlap somewhat, potentially making them difficult to distinguish without additional hyperspectral field, laboratory or remote sensing data. The results of this study were compared to mixtures of minerals described for ore, gangue and alteration minerals in deposit models for sandstone-hosted uranium, sedimentary bentonite and sedimentary zeolite. Use of these spectra can help facilitate mapping of both waste materials from the legacy mining of the above commodities, as well as future exploration and resource assessment activities.

Texas↗

Baseline data for evaluating development trajectory and provision of ecosystem services of created fringing oyster reefs in Vermilion Bay, Louisiana

Understanding the time frame in which ecosystem services (that is, water quality maintenance, shoreline protection, habitat provision) are expected to be provided is important when restoration projects are being designed and implemented. Restoration of three-dimensional shell habitats in coastal Louisiana and elsewhere presents a valuable and potentially self-sustaining approach to providing shoreline protection, enhancing nekton habitat, and providing water quality maintenance. As with most restoration projects, the development of expected different ecosystem services often occurs over varying time frames, with some services provided immediately and others taking longer to develop. This project was designed initially to compare the provision and development of ecosystem services by created fringing shoreline reefs in subtidal and intertidal environments in Vermilion Bay, Louisiana. Specifically, the goal was to test the null hypothesis that over time, the oyster recruitment and development of a sustainable oyster reef community would be similar at both intertidal and subtidal reef bases, and these sustainable reefs would in time provide similar shoreline stabilization, nekton habitat, and water quality services over similar time frames. Because the ecosystem services hypothesized to be provided by oyster reefs reflect long-term processes, fully testing the above-stated null hypothesis requires a longer-time frame than this project allowed. As such, this project was designed to provide the initial data on reef development and provision of ecosystem services, to identify services that may develop immediately, and to provide baseline data to allow for longer-term follow up studies tracking reef development over time. Unfortunately, these initially created reef bases (subtidal, intertidal) were not constructed as planned because of the Deepwater Horizon oil spill in April 2010, which resulted in reef duplicates being created 6 months apart. Further confounding the project were additional construction and restoration projects along the same shorelines which occurred between 2011 and June 2012. Because of constant activity near and around the reefs and continuing construction, development trajectories could not be compared among reef types at this time. This report presents the data collected at the sites over 3 years (2010–2012), describing only conditions and trends. In addition, these data provide an extensive and detailed dataset documenting initial conditions and initial ecosystem changes which will prove valuable in future data collection and analyses of reef development at this site. Data collection characterized the local water quality conditions (salinity, temperature, total suspended sediments, dissolved oxygen, chlorophyll a), adjacent marsh vegetation, soils, and shoreline position along the project shoreline at Vermilion Bay. During the study, marsh vegetation and soil characteristics were similar across the study area and did not change over time. Shoreline movement indicated shoreline loss at all sites, which varied by reefs. Water quality conditions followed expected seasonal patterns for this region, and no significant nonseasonal changes were measured throughout the study period. Despite oyster recruitment in fall 2010 and 2011, few if any oysters survived from the 2010 year class to 2012. At the last sampling of this project, some oysters recruited in fall 2011 survived through 2012, resulting in an on-reef density of 18.3 ± 2.1 individuals per square meter (mean size: 85.6 ± 2.2 millimeters). Because project goals were to compare reef development and provision of ecosystem services over time, as well as many of the processes identified for monitoring reflect long-term processes, results and data are presented only qualitatively, and trends or observations should be interpreted cautiously at this point. Measurable system responses to reef establishment require more time than was available for this study. These data provide a valuable baseline that can be ultimately used to help inform site selections for future restoration projects as well to further investigate the development trajectories of ecosystem provision of created reefs in this region.

Louisiana↗

Simulation of regional groundwater flow and advective transport of per- and polyfluoroalkyl substances, Joint Base McGuire-Dix-Lakehurst and vicinity, New Jersey, 2018

A three-dimensional numerical model of groundwater flow was developed and calibrated for the unconsolidated New Jersey Coastal Plain aquifers underlying Joint Base McGuire-Dix-Lakehurst (JBMDL) and vicinity, New Jersey, to evaluate groundwater flow pathways of per- and polyfluoroalkyl substances (PFAS) contamination associated with use of aqueous film forming foam (AFFF) at the base. The regional subsurface flow model spans an area of approximately 518 square miles around JBMDL and is based on a previously developed hydrogeologic framework of the area. Steady-state flow in the unconsolidated aquifers was simulated using the MODFLOW 6 groundwater flow model, which is able to account for hydrostratigraphic pinchouts and discontinuities in the Coastal Plain aquifers underlying JBMDL. To account for local patterns of fluid flow driving advective subsurface migration of PFAS, the grid was refined using quadtree meshes spanning 21 areas where historical AFFF use was identified, five off-site reconnaissance areas identified by AFCEC as areas in which the occurrence of PFAS is most likely to pose a potential danger to local drinking water supplies, and along streams that behave as drains in the base-flow-dominated Coastal Plain. Following grid refinement, four physical processes known to govern subsurface flow were introduced to the model. These included effective precipitation recharge, discharge to streams and stream-connected wetlands, regional inflows and outflows along the model bottom, and withdrawals from wells, each of which were incorporated into the model as either external or internal boundary conditions. To account for effective precipitation recharge, a specified-flow boundary was assigned along the top of the model. Similarly, regional flows predicted using the modified U.S Geological Survey’s New Jersey Coastal Plain Regional Aquifer System Analysis model were treated as specified-flow boundary conditions along the bottom of the model. Base-flow losses were treated as drains along streams delineated using a 10-foot LiDAR dataset. Drains were also assigned to cells falling within stream-connected National Hydrologic Database wetlands. Finally, well-pumpage data mined from the New Jersey Water Transfer database were added to the model to account for extraction of groundwater through pumping from industrial-supply and drinking-water-supply wells. Along model edges established at groundwater divides, where the net flux of water across the boundary is equal to zero, natural no-flow boundary conditions were imposed. The refined flow model was calibrated using the parameter-estimation (PEST) program, which adjusts model parameters by performing a gradient search over the sum-of-squared-error objective function until the parameter set that produces simulated water levels and base flows most closely matches 544 water levels and 20 estimated base flows and closely adheres to initial parameter estimates. Based on the analysis of calibration residuals, the model did not appear to be affected by significant model structural error. The MODPATH particle-tracking algorithm was used to estimate advective transport paths of PFAS in the vicinity of JBMDL. Forward tracking was used to determine paths of PFAS away from AFFF source areas to streams, wetlands, pumping wells, and geographic areas that PFAS may contaminate. Additionally, reverse tracking was used to determine particle pathlines away from off-site PFAS reconnaissance areas, or areas within which all sources of PFAS might be advectively transported into subsurface drinking-water supplies, to locations at land surface that may indicate a source of PFAS. The coupled and calibrated groundwater flow and particle-tracking transport model provide valuable tools for predicting the relative extent of PFAS contamination from onsite legacy source areas. The calibrated model also provides measures of water-level and base-flow observation influence that can help guide future data-collection efforts related to groundwater and surface water sampling for PFAS.

New Jersey↗

Procedures and rationale for marking sea otters captured and treated during the T/V Exxon Valdez oil spill

Four methods were used for marking sea otters ( Enhydra lutris ) captured and treated during the response to the T/V Exxon Valdez oil spill. Colored and numbered flipper tags were placed on each sea otter that was brought to the otter treatment centers. These tags allowed individual recognition and permitted the tracking of individuals through the treatment and holding processes. Recovery of tagged carcasses may provide a crude measure of the fates of rehabilitated otters. Seven sea otters were instrumented with radio transmitters attached to flipper tags as part of a pilot release program. The results of that study were inconclusive. Forty-five sea otters were implanted with radio transmitters as part of a study to assess the fate of rehabilitated sea otters. Specific objectives of the study include estimating survived rates and monitoring the reproductive success of the sample of rehabilitated sea otters and comparing the results with similar variables in a control population. Transponder chips were injected in the perianal region of all but seven of the instrumented sea otters, and all were tagged with red flipper tags. The advantages and disadvantages of each type of tag are discussed.

Alaska↗

Analysis for agreement of the Northern Gulf of Mexico topobathymetric digital elevation model with 3-Dimensional Elevation Program 1/3 arc-second digital elevation models

Topographical differencing and edge-matching analyses were used to evaluate agreement of the Coastal National Elevation Database Applications Project’s Northern Gulf of Mexico topobathymetric digital elevation model (TBDEM) with The National Map 3-Dimensional Elevation Program (3DEP) 1/3 arc-second digital elevation models (DEMs). In addition to topographic map products provided through the National Geospatial Program, the model integrates bathymetric and topobathymetric datasets for three-dimensional (3D) mapping of rivers, lakes, and bays in the upland and intertidal wetlands to offshore environments in coastal zones from the border between Texas and Louisiana to east of Mobile Bay, Alabama. Contoured elevation differences between the Northern Gulf of Mexico TBDEM and the 3DEP 1/3 arc-second DEMs indicate that 85 percent of elevation data in the Northern Gulf of Mexico TBDEM agree (no difference for contoured elevations) between 95 and 100 percent with 3DEP 1/3 arc-second DEMs. Edge matching differences between adjacent Northern Gulf of Mexico TBDEM source projects or between the TBDEM and 3DEP DEMs indicate most seams between integrated and 3DEP DEMs are smooth. Where seams did not match, most differences were in the range of tenths to hundredths of a meter. Valid differences that are greater than plus or minus 2 meters in areas of bathymetric data are found in the Mississippi River, Atchafalaya River, Lower Atchafalaya River, Wax Lake Pass channel, the Vermilion Bay bathymetric datasets, and where topobathymetric datasets are integrated in the model. Areas with positive or negative outlier difference elevations seem to be a result of site conditions that affect light detection and ranging (lidar) waveform return signals, misclassification of surface features, or possibly because of interpolation required to develop a smooth elevation surface. Results of this analysis provide information to help understand model parameters and agreement of the Northern Gulf of Mexico TBDEM developed using different data types from different sources with The National Map 3DEP DEMs. Inclusion of bathymetric and topobathymetric data types in the 3DEP aligns with the mission to respond to growing needs for a wide range of three-dimensional representations of the Nation and supports the U.S. Geological Survey strategy for developing a National Terrain Model to provide hydrographic and elevation data that extend the elevation surface below water bodies. The 3D Nation Requirements and Benefits Study sponsored by the U.S. Geological Survey and National Oceanic and Atmospheric Administration to assess local to regional Tribal, State, and Federal technical requirements, needs, and benefits for using topographic and bathymetric 3DEP elevation data will be used to help develop and refine future program alternatives for 3D elevation data that include a category for bathymetry and topobathymetry. At the time of this report (2019), 3DEP acquisition is specific to topographic lidar that meets lidar DEM specifications and which requires surface-water feature areas to be hydroflattened. Cataloging bathymetric and topobathymetric DEMs as part of the 3DEP will require new specifications for acoustic, lidar, merged acoustic and lidar, and possibly other bathymetric and topobathymetric survey data types.

Alabama, Florida, Louisiana, Mississippi, Texas↗

Morphodynamic modelling of the wilderness breach, Fire Island, New York. Part I: Model set-up and validation

On October 29, 2012, storm surge and large waves produced by Hurricane 13 Sandy resulted in the formation of a breach in eastern Fire Island, NY. The goals of this study 14 are to gain a better understanding of the physical processes that govern breach behavior and 15 to assess whether process-based models can be used to forecast the evolution of future 16 breaches. The Wilderness Breach grew rapidly in size during the first winter following 17 formation. Growth of the breach was accompanied by the formation of a complex of flood 18 shoals inside Great South Bay, a primary channel that flowed through the eastern part of the 19 flood shoals, and an ebb shoal on the ocean side of the breach. From the summer of 2013 20 through late 2015, the breach continued to change and evolve, albeit at a much slower pace 21 than in the first year after formation. A hybrid combination of Delft3D and XBeach models is 22 used to hindcast the morphodynamic evolution of the Wilderness Breach over the first three 23 years after formation. The formation of the breach during Hurricane Sandy is not part of the 24 simulations. Model simulations are initiated with a post-storm topography in which the 25 breach is already present. The models are capable of hindcasting the main morphodynamic 26 changes of the Wilderness Breach. The spatial patterns, as well as the bulk statistics, such as 27 2 breach geometry and sediment volume changes, are reasonably 28 well reproduced by the model. 29 The model sheds light on previously unknown processes of breach evolution, especially 30 regarding sediment transport and flow regimes within the breach complex.

New York↗

Management strategy evaluation to assess trade-offs associated with invasive Blue Catfish fisheries and predation impacts

Objective Many species are intentionally introduced beyond their native range to provide benefits to humans (e.g., food, recreation, or biocontrol). However, introduced species can become invasive and can harm native species, prompting resource managers to explore options to simultaneously conserve native biota and enhance fishing opportunities. Management of Chesapeake Bay Blue Catfish Ictalurus furcatus is complicated because the species supports a renowned trophy fishery and commercial and recreational harvest but also negatively affects native species. Consequently, there is uncertainty and disagreement on how to balance trade-offs associated with varying stakeholder interests. Methods We used a management strategy evaluation to project Blue Catfish population dynamics into the future under different fishery policies to understand whether fishery yield, trophy fishing opportunities, and predation on the economically and ecologically important blue crab Callinectes sapidus could be optimized in the James River, a Chesapeake Bay subestuary. Results Simulated population trajectories indicated that objectives related to maintaining fisheries and conserving prey populations were in conflict. Policies that increased the yield and abundance of trophy-size Blue Catfish (≥100 cm total length) generally increased predation on the blue crab, which supports a valuable fishery. Intense harvest of smaller length-classes and protection of larger Blue Catfish yielded outcomes in which trophy fish abundance increased and blue crab predation declined compared to baseline conditions. However, these outcomes were generally associated with lower Blue Catfish yields after 25 years. There were zero scenarios in which the fishery yield increased and blue crab predation decreased after 25 years. Policies limiting the harvest of small Blue Catfish resulted in large population abundances, suggesting that reducing the abundance of small fish could be important for reducing impacts on native species. Conclusions This study supports the importance of management planning to develop objectives and performance measures based on an improved understanding of trade-offs associated with harvest management for a nonnative fish with economic value.

Virginia↗

The marine terraces of Santa Cruz Island, California: Implications for glacial isostatic adjustment models of last-interglacial sea-level history

Glacial isostatic adjustment (GIA) models hypothesize that along coastal California, last interglacial (LIG, broadly from ~130 to ~115 ka) sea level could have been as high as +11 m to +13 m, relative to present, substantially higher than the commonly estimated elevation of +6 m. Areas with low uplift rates can test whether such models are valid. Marine terraces on Santa Cruz Island have previously been reported to occur at low (<10 m) elevations, but ages of many such localities are not known. Using lidar imagery as a base, marine terraces on Santa Cruz Island were newly mapped, elevations were measured, fossils were collected for U-series dating (corals), strontium isotope compositions and amino acid geochronology (mollusks), and paleozoogeography (all taxa). Sr isotope compositions of mollusks from the highest of three marine terraces give ages of ~2.5 Ma to 1.9 Ma, along with Pliocene ages, from shells interpreted to be reworked. U-series ages of corals from the western part of the island indicate that low-elevation terraces north of the Santa Cruz Island fault correlate to the LIG. Where corals are lacking, amino acid ratios and faunal aspects support terrace correlation to the LIG high stand of sea. Elevations of most terrace localities north of the east-west trending Santa Cruz Island fault, in both the western and eastern parts of the island, range from 5.75 m to 8 m above sea level, well below the modeled paleo-sea-level range. Subsidence is ruled out as a mechanism for explaining the lower-than-modeled elevations, because higher-elevation terraces are present along much of the Santa Cruz Island coast north of the fault, indicating long-term tectonic uplift. The low elevations of the LIG terrace fragments are, however, consistent with a low rate of uplift derived from the higher, ~2.5–1.9 Ma terrace. A number of other localities on the Pacific Coast, also dated to the LIG, have marine terrace elevations below the modeled level. GIA models may have overestimated last interglacial sea level by a substantial amount and need to be revised if used for forecasts for future sea-level rise.

California↗

Two centuries of southwest Iceland annually-resolved marine temperature reconstructed from Arctica islandica shells

Iceland's exposure to major ocean current pathways of the central North Atlantic makes it a useful location for developing long-term proxy records of past marine climate. Such records provide more detailed understanding of the full range of past variability which is necessary to improve predictions of future changes. We constructed a 225-year (1791–2015 CE) master shell growth chronology from 29 shells of Arctica islandica collected at 100 m water depth in southwest Iceland (Faxaflói). The growth chronology provides a robust age model for shell oxygen isotope (δ 18 O shell ) data produced at annual resolution for 251 years (1765–2015 CE). The temperature reconstruction derived from δ 18 O shell shows coherence with May–October local surface temperature records and sea surface temperatures in the North Atlantic region, suggesting it is a useful proxy indicator of water temperature variability at 100 m depth within Faxaflói. Field correlations between the shell-based records and gridded sea surface temperature data reveal strong positive correlations between the 1-year lagged shell growth and temperatures within the subpolar gyre post-1972, suggesting a delayed influence of subpolar gyre dynamics on ecological indicators in southwest Iceland in recent decades. However, the shell growth chronology and δ 18 O shell record generally show relatively weak and insignificant correlations with larger region climate indices including the Atlantic Multidecadal Variability, North Atlantic Oscillation, and East Atlantic pattern. Therefore the interannual variations in the newly produced shell-based records appear to reflect more local to regional dynamics around southwest Iceland than large-scale modes of climate variability.

Estuarine, Coastal and Shelf Science↗

Use of geoelectrical methods in groundwater pollution surveys in a coastal environment

The pollution of coastal aquifers by old landfills can contaminate valuable and scarce water resources in the freshwater lens utilized seasonably by overcrowded communities. The pollutants will ultimately flow into the sea where they may also cause a coastal water pollution problem. We have detected pollution in the freshwater lens from a sanitary landfill near Provincetown, Cape Cod, using the geoelectrical resistivity method. This survey included Schlumberger geoelectrical depth soundings and a horizontal geoelectrical profile using the Wenner configuration. The geoelectrical survey was conducted at a site along Highway 6 where it passes the coastal town of Provincetown and a sanitary landfill that has been in operation since 1954. The depth soundings suggest the characteristic decrease in resistivity vs. depth from the high resistivity of the unsaturated zone to the low resistivity of the saltwater saturated zone. The freshwater lens is clearly identified by the change in slope of the steeply dipping curve of resistivity versus electrode spacing. Interpretations made using a multilayer program, Geomate, resulted in layer resistivities between 460 and 95 ohm?m for the freshwater lens. A comparison with well water resistivities suggests that a layer resistivity of 230 ohm?m or lower is indicative of pollution in the freshwater lens. The results of the geoelectrical depth soundings were confirmed in the Wenner horizontal profile. Both measurements suggest that the pollutants do not spread evenly as one would expect for a homogeneous and isotropic medium. Instead, a preferred channel for the flow of the pollutants is observed along a path from the landfill toward the shoreline. The depth to the saltwater/freshwater interface or, more specifically, to the low resistivity-high resistivity interface appears to be shallow where the freshwater lens is polluted. This was confirmed by pore water well samples that were highly mineralized. The equilibrium postulated by the Ghyben-Herzberg relation appears to be disturbed in the area of aquifer pollution. This rise in the conductivity boundary is caused by the highly mineralized bottom of the contaminant plume that submerges into the saltwater saturated zone. In the area of high freshwater pollution the groundwater can be subdivided into three layers that show a decrease in resistivity with depth. The formation factor, F, defined as the ratio of bulk aquifer resistivity to pore water resistivity, shows unusually high values between 10 and 12. These high values are unexpected for an unconsolidated sand. Pollution residues are suspected to clog the pores and thus to increase the resistivity. It is possible that iron-oxidizing bacteria and the precipitation of dissolved iron or organic pollutants are the cause of the high values of F. If proven correct, these interesting possibilities could lead to future new applications of the geoelectrical resistivity method in contaminant hydroloy.

Journal of Applied Geophysics↗

Early growth interactions between a mangrove and an herbaceous salt marsh species are not affected by elevated CO2 or drought

Increasing atmospheric carbon dioxide (CO 2 ) concentrations are likely to influence future distributions of plants and plant community structure in many regions of the world through effects on photosynthetic rates. In recent decades the encroachment of woody mangrove species into herbaceous marshes has been documented along the U.S. northern Gulf of Mexico coast. These species shifts have been attributed primarily to rising sea levels and warming winter temperatures, but the role of elevated CO 2 and water availability may become more prominent drivers of species interactions under future climate conditions. Drought has been implicated as a major factor contributing to salt marsh vegetation dieback in this region. In this greenhouse study we examined the effects of CO 2 concentration (∼380 ppm, ∼700 ppm) and water regime (drought, saturated, flooded) on early growth of Avicennia germinans , a C 3 mangrove species, and Spartina alterniflora , a C 4 grass. Plants were grown in monocultures and in a mixed-species assemblage. We found that neither species responded to elevated CO 2 over the 10-month duration of the experiment, and there were few interactions between experimental factors. Two effects of water regime were documented: lower A. germinans pneumatophore biomass under drought conditions, and lower belowground biomass under flooded conditions regardless of planting assemblage. Evidence of interspecific interactions was noted. Competition for aboveground resources (e.g., light) was indicated by lower S. alterniflora stem biomass in mixed-species assemblage compared to biomass in S. alterniflora monocultures. Pneumatophore biomass of A. germinans was reduced when grown in monoculture compared to the mixed-species assemblage, indicating competition for belowground resources. These interactions provide insight into how these species may respond following major disturbance events that lead to vegetation dieback. Site variation in propagule availability and physico-chemical conditions will determine plant community composition and structure following such disturbances when these two species co-occur.

Estuarine, Coastal and Shelf Science↗

Expanding the North American Breeding Bird Survey analysis to include additional species and regions

The North American Breeding Bird Survey (BBS) contains data for >700 bird species, but analyses often focus on a core group of ∼420 species. We analyzed data for 122 species of North American birds for which data exist in the North American Breeding Bird Survey (BBS) database but are not routinely analyzed on the BBS Summary and Analysis Website. Many of these species occur in the northern part of the continent, on routes that fall outside the core survey area presently analyzed in the United States and southern Canada. Other species not historically analyzed occur in the core survey area with very limited data but have large portions of their ranges in Mexico and south. A third group of species not historically analyzed included species thought to be poorly surveyed by the BBS, such as rare, coastal, or nocturnal species. For 56 species found primarily in regions north of the core survey area, we expanded the scope of the analysis, using data from 1993 to 2014 during which ≥3 survey routes had been sampled in 6 northern strata (Bird Conservation regions in Alaska, Yukon, and Newfoundland and Labrador) and fitting log-linear hierarchical models for an augmented BBS survey area that included both the new northern strata and the core survey area. We also applied this model to 168 species historically analyzed in the BBS that had data from these additional northern strata. For both groups of species we calculated survey-wide trends for the both core and augmented survey areas from 1993 to 2014; for species that did not occur in the newly defined strata, we computed trends from 1966 to 2014. We evaluated trend estimates in terms of established credibility criteria for BBS results, screening for imprecise trends, small samples, and low relative abundance. Inclusion of data from the northern strata permitted estimation of trend for 56 species not historically analyzed, but only 4 of these were reasonably monitored and an additional 13 were questionably monitored; 39 of these species were likely poorly monitored because of small numbers of samples or very imprecisely estimated trends. Only 4 of 66 “new” species found in the core survey area were reasonably monitored by the BBS; 20 were questionably monitored; and 42 were likely poorly monitored by the BBS because of inefficiency in precision, abundance, or sample size. The hierarchical analyses we present provide a means for reasonable inclusion of the additional species and strata in a common analysis with data from the core area, a critical step in the evolution of the BBS as a continent-scale survey. We recommend that results be presented both 1) from 1993 to the present using the expanded survey area, and 2) from 1966 to the present for the core survey area. Although most of the “new” species we analyzed were poorly monitored by the BBS during 1993–2014, continued expansion of the BBS will improve the quality of information in future analyses for these species and for the many other species presently monitored by the BBS.

Journal of Fish and Wildlife Management↗

Spatial and temporal assessment of back-barrier erosion on Cumberland Island National Seashore, Georgia, 2011–2013

Much research has been conducted to better understand erosion and accretion processes for the seaward zones of coastal barrier islands; however, at Cumberland Island National Seashore, Georgia, the greater management concern is the effect that erosion is having on the resources of the island’s western shoreline, or the back barrier. Catastrophic slumping and regular rates of erosion greater than 1 meter per year threaten important habitat, historical and pre-historical resources, and modern infrastructure on the island. Prior research has helped National Park Service (NPS) staff identify the most severe and vulnerable areas, but in order to develop effective management actions, information is needed on what forces and conditions cause erosion. To this end, the U.S. Geological Survey, in cooperation with the NPS, conducted two longitudinal surveys, one each at the beginning and end of the approximately year-long monitoring period from late 2011 to early 2013, along five selected segments of the back barrier of the Cumberland Island National Seashore. Monitoring stations were constructed at four of these locations that had previously been identified as erosional hotspots. The magnitude of erosion at each location was quantified to determine the relative influence of causative agents. Results indicate that erosion is, in general, highly variable within and among these segments of the Cumberland Island National Seashore’s back barrier. Observed erosion ranged from a maximum of 2.5 meters of bluff-line retreat to some areas that exhibited no net erosion over the 1-year study period. In terms of timing of erosion, three of the four sites were primarily affected by punctuated erosional events that were coincident with above-average high tides and elevated wind speeds. The fourth site exhibited steady, low-magnitude retreat throughout the study period. While it is difficult to precisely subscribe certain amounts of erosion to specific agents, this study provides insight into the mode of erosion among sites and the interaction among factors that set up conditions that may be leading to punctuated events. Estimates of sea-level rise were incorporated into the results of this study to project conditions that could be in place by the end of the 21st century. When using the erosion rates observed in this study to extrapolate future shoreline position, results indicate an average retreat (across all monitored locations) of 15 meters by 2050 and approximately 37 meters by 2100.

Georgia↗

Micropaleontologic record of late Pliocene and Quaternary paleoenvironments in the northern Albemarle Embayment, North Carolina, U.S.A.

Micropaleontological data provide a strong actualistic basis for detailed interpretations of Quaternary paleoenvironmental change. The 90 m-thick Quaternary record of the Albemarle Embayment in the mid-Atlantic coastal plain of the USA provides an excellent opportunity to use such an approach in a region where the details of Quaternary environmental change are poorly known. The foraminiferal record in nine cores from the northern Outer Banks, east of Albemarle Sound, North Carolina, indicates the deposition of subhorizontal, mostly open-marine early to late Pleistocene units unconformably upon a basement of late Pliocene reduced-oxygen, fine-grained, shelf-basin deposits. Pollen data record several warm–cool fluctuations within the early to mid-Pleistocene deposits. Diatom data indicate that some fresh and brackish-water units occur within the generally open-marine Pleistocene succession. A channel cut by the paleo-Roanoke River during the last glacial sea-level lowstand occurs in the northern part of the study area. Pollen indicates that the basal fluvial valley fill accumulated in cooler than modern climate conditions in the latest Pleistocene. Overlying silts and muds accumulated under cool climatic, estuarine conditions according to diatom and pollen data. Radiocarbon ages from the estuarine deposits indicate that the bulk of these sediments accumulated during the latest Pleistocene. The estuarine channel-fill deposits are overlain by Holocene open-marine sands deposited as the rising sea transgressed into the estuary approximately 8.5 to 9.0 kyr BP. Within the barrier island drill cores of this study, fully marine sedimentation occurred throughout the Holocene. However, immediately west of the present barrier island, mid- to late Holocene estuarine deposits underlie the modern Albemarle Sound. The islands that currently form a continuous barrier across the mouth of Albemarle Sound have a complex history of Holocene construction and destruction and large portions of them may be less than 3 kyr old. The barrier island sands overlie open-marine sands of Colington Shoal in the north and to the south overlie fluvial and marine sand filling paleo-Roanoke tributary valleys. The Pleistocene sediments underlying the northern Outer Banks study area are mainly of open inner to mid-shelf origin. If, as is likely, sea level continues to rise, a return to such environmental conditions is likely in the near future.

North Carolina↗

Preliminary investigation of the distribution and resources of coal in the Kaiparowits Plateau, southern Utah

This report on the coal resources of the Kaiparowits Plateau, Utah is a contribution to the U.S. Geological Survey's (USGS) 'National Coal Resource Assessment' (NCRA), a five year effort to identify and characterize the coal beds and coal zones that could potentially provide the fuel for the Nation's coal-derived energy during the first quarter of the twenty-first century. For purposes of the NCRA study, the Nation is divided into regions. Teams of geoscientists, knowledgeable about each region, are developing the data bases and assessing the coal within each region. The five major coal-producing regions of the United States under investigation are: (1) the Appalachian Basin; (2) the Illinois Basin; (3) the Gulf of Mexico Coastal Plain; (4) the Powder River Basin and the Northern Great Plains; and (5) the Rocky Mountains and the Colorado Plateau. Six areas containing coal deposits in the Rocky Mountain and Colorado Plateau Region have been designated as high priority because of their potential for development. This report on the coal resources of the Kaiparowits Plateau is the first of the six to be completed. The coal quantities reported in this study are entirely 'resources' and represent, as accurately as the data allow, all the coal in the ground in beds greater than one foot thick. These resources are qualified and subdivided by thickness of coal beds, depth to the coal, distance from known data points, and inclination (dip) of the beds. The USGS has not attempted to estimate coal 'reserves' for this region. Reserves are that subset of the resource that could be economically produced at the present time. The coal resources are differentiated into 'identified' and 'hypothetical' following the standard classification system of the USGS (Wood and others, 1983). Identified resources are those within three miles of a measured thickness value, and hypothetical resources are further than three miles from a data point. Coal beds in the Kaiparowits Plateau are laterally discontinuous relative to many other coal bearing regions of the United States. That is, they end more abruptly and are more likely to fragment or split into thinner beds. Because of these characteristics, the data from approximately 160 drill holes and 40 measured sections available for use in this study are not sufficient to determine what proportion of the resources is technologically and economically recoverable. The Kaiparowits Plateau contains an original resource of 62 billion short tons of coal in the ground. Original resource is defined to include all coal beds greater than one foot thick in the area studied. None of the resource is recoverable by surface mining. However, the total resource figure must be regarded with caution because it does not reflect geologic, technological, land-use, and environmental restrictions that may affect the availability and the recoverability of the coal. At least 32 billion tons of coal are unlikely to be mined in the foreseeable future because the coal beds are either too deep, too thin to mine, inclined at more than 12?, or in beds that are too thick to be completely recovered in underground mining. The estimated balance of 30 billion tons of coal resources does not reflect land use or environmental restrictions, does not account for coal that would be bypassed due to mining of adjacent coal beds, does not consider the amount of coal that must remain in the ground for roof support, and does not take into consideration the continuity of beds for mining. Although all of these factors will reduce the amount of coal that could be recovered, there is not sufficient data available to estimate recoverable coal resources. For purposes of comparison, studies of coal resources in the eastern United States have determined that less than 10 percent of the original coal resource, in the areas studied, could be mined economically at today's prices (Rohrbacher and others, 1994).

Open-File Report↗

Detecting sea-level hazards: Simple regression-based methods for calculating the acceleration of sea level

This report documents the development of statistical tools used to quantify the hazard presented by the response of sea-level elevation to natural or anthropogenic changes in climate and ocean circulation. A hazard is a physical process (or processes) that, when combined with vulnerability (or susceptibility to the hazard), results in risk. This study presents the development and comparison of new and existing sea-level analysis methods, exploration of the strengths and weaknesses of the methods using synthetic time series, and when appropriate, synthesis of the application of the method to observed sea-level time series. These reports are intended to enhance material presented in peer-reviewed journal articles where it is not always possible to provide the level of detail that might be necessary to fully support or recreate published results. The purpose of this report is to document and compare three simple methodologies that have previously been used to provide estimates with associated errors of the acceleration of sea-level elevation. These techniques have been used by coastal scientists and planners in assessing coastal risk over a wide range of spatial and temporal scales. Because relative sea-level (SL) elevation time series contain energetic fluctuations at many time scales, extracting what can be relatively small rate and acceleration signals (along with estimates of the error) from much larger &ldquo;noise&rdquo; has proven to be both difficult and controversial. Acceleration is a preferred measure of SL response to recent changes in the Earth&rsquo;s climate because over time scales of 100 years or less slow vertical land motions (such as glacial isostatic adjustment) contribute only to the linear signal and not to acceleration, thus reducing the complexity of the analysis. Hence acceleration is useful if the goal of a study is to characterize and quantify the hazard associated with the changing relative elevation of water with respect to land on decadal time scales. Although in some cases it may be necessary to determine the cause of relative sea level rise, as a first step, it is important to accurately estimate the magnitude of the threat. Most researchers agree that global sea level (GSL) rose persistently through much of the 20th century at about 1.5&ndash;2.0 millimeters per year (mm/yr). There is far less agreement about whether the rate of sea-level rise (SLR) is increasing (that is, an acceleration in SL). Recent studies, and most of their predecessors, use tide gage data to quantify SL acceleration, A SL (t). In the current study, three techniques were used to calculate acceleration from tide gage data, and of those examined, it was determined that the two techniques based on sliding a regression window through the time series are more robust compared to the technique that fits a single quadratic form to the entire time series, particularly if there is temporal variation in the magnitude of the acceleration. The single-fit quadratic regression method has been the most commonly used technique in determining acceleration in tide gage data. The inability of the single-fit method to account for time-varying acceleration may explain some of the inconsistent findings between investigators. Properly quantifying A SL (t) from field measurements is of particular importance in evaluating numerical models of past, present, and future SLR resulting from anticipated climate change.

Open-File Report↗