Geology Reports⌕ Search

SEARCH · Geology Reports

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

Accuracy of EAARL lidar ground elevations using a bare-earth algorithm in marsh and beach grasses on the Chandeleur Islands, Louisiana

The NASA Experimental Advanced Airborne Lidar (EAARL) is an airborne lidar (light detection and ranging) instrument designed to map coastal topography and bathymetry. The EAARL system has the capability to capture each laser-pulse return over a large signal range and can digitize the full waveform of the backscattered energy. Because of this ability to capture the full waveform, the EAARL system can map features such as coral reefs, beaches, coastal vegetation, and trees, where extreme variations in the laser backscatter are caused by different physical and optical characteristics. Post-processing of the EAARL data is accomplished using the Airborne Lidar Processing System (ALPS) (Nayegandhi and others, 2009). In ALPS, the waveform of the lidar is analyzed and split into first and last returns. The 'first returns' are indicative of vegetation-canopy height, or bare ground in the absence of vegetation, whereas 'last returns' typically represent 'bare-earth' elevations under vegetation. To test the accuracy of the first-return and bare-earth EAARL data, topographic and vegetation height surveys were conducted in the Chandeleur Islands, concurrent with an EAARL lidar survey and an aerial oblique-photographic survey from September 20 to 27, 2006. The Chandeleur Islands are a north-south-oriented chain of low-lying islands located approximately 100 kilometers east of the city of New Orleans, Louisiana. The islands are narrow north-south strips of land with marsh on the landward (west sides) and sandy beaches on their gulfward (east sides). Prior to Hurricane Katrina, which made landfall at Buras, Louisiana, as a Category 3 storm on August 29, 2005, prominent, 3- to 4-meter-high sand dunes were present in the northern Chandeleurs. The storm removed them, leaving post-storm island elevations of generally less than 2 meters above 0.0 NAVD88. This report is part of a study of the impact of Hurricane Katrina on the Chandeleur Islands using pre-storm and post-storm lidar surveys to detect morphological changes. The islands lost over 80 percent of their land area during Hurricane Katrina, and in the first 2 years following Katrina, many of the islands experienced continued shoreline retreat (Sallenger and others, 2007). In addition to land-area losses, the loss of dunes made the islands increasingly vulnerable to future storm impacts. The U.S. Geological Survey, along with partners in the Louisiana Department of Natural Resources and the U.S. Army Corps of Engineers, continues to monitor changes in shoreline position, land area, and elevation in the Chandeleur Islands.

Open-File Report↗

Simulated effects of alternative withdrawal strategies on groundwater flow in the unconfined Kirkwood-Cohansey aquifer system, the Rio Grande water-bearing zone, and the Atlantic City 800-foot sand in the Great Egg Harbor and Mullica River Basins, New Jersey

Groundwater is essential for water supply and plays a critical role in maintaining the environmental health of freshwater and estuarine ecosystems in the Atlantic Coastal basins of New Jersey. The unconfined Kirkwood-Cohansey aquifer system and the confined Atlantic City 800-foot sand are major sources of groundwater in the area, and each faces different water-supply concerns. The U.S. Geological Survey (USGS), in cooperation with the New Jersey Department of Environmental Protection (NJDEP), conducted a study to simulate the effects of withdrawals in the Kirkwood-Cohansey aquifer system, the Atlantic City 800-foot sand, and the Rio Grande water-bearing zone and to evaluate potential scenarios. The study area encompasses Atlantic County and parts of Burlington, Camden, Gloucester, Ocean, Cape May, and Cumberland Counties. The major hydrogeologic units affecting water supply in the study area are the surficial Kirkwood-Cohansey aquifer system, a thick diatomaceous clay confining unit in the upper part of Kirkwood Formation; the Rio Grande water-bearing zone; and the Atlantic City 800-foot sand of the Kirkwood Formation. Hydrogeologic data from 18 aquifer tests and specific capacity data from 230 wells were analyzed to provide horizontal hydraulic conductivity of the aquifers. Groundwater withdrawals are greatest from the Kirkwood-Cohansey aquifer system, and 65 percent of the water is used for public supply. Groundwater withdrawals from the Atlantic City 800-foot sand are about half those from the Kirkwood-Cohansey aquifer system. Ninety-five percent of the withdrawals from the Atlantic City 800-foot sand is used for public supply. Data from six streamgaging stations and 51 low-flow partial record sites were used to estimate base flow in the area. Base flow ranges from 60 to 92 percent of streamflow. A groundwater flow model of the Kirkwood-Cohansey aquifer system, the Rio Grande water-bearing zone, and the Atlantic City 800-foot sand was developed and calibrated using water-level data from 148 wells and base-flow data from 22 gaging or low-flow partial record stations. The Kirkwood-Cohansey aquifer system within the Great Egg Harbor River and the Mullica River Basins was simulated on a monthly basis from 1998 through 2006. An existing regional model of the New Jersey Coastal Plain was revised to provide boundary conditions for the Great Egg Harbor and Mullica River Basin model (referred to as the Great Egg-Mullica model). In the Great Egg-Mullica model, monthly groundwater recharge rates used in the model ranged from 10-15 inches per year in 2001 to 20-25 inches per year in 2005. The mean-absolute error for 10 of the 14 long-term hydrographs used in model calibration was less than 5 ft. Groundwater flow budgets for the Great Egg-Mullica model calibration periods, May 2005 and September 2006, and for the entire model calibration period 1998 to 2006, showed that nearly 70 percent of the water entering the Atlantic City 800-foot sand came from the horizontal connection with the Kirkwood-Cohansey aquifer system in updip areas. The groundwater flow model was used to simulate scenarios under three possible conditions: average 1998 to 2006 withdrawals (Average scenario), full-allocation withdrawals (Full Allocation scenario), and projected 2050-demand withdrawals (2050 Demand scenario). Withdrawals in the Full Allocation scenario are nearly twice the withdrawals from the Average scenario, primarily because of the potential for large agricultural withdrawals if all allocations are used. Withdrawals for the 2050 Demand scenario are about 50 percent greater than those for the Average scenario, primarily due to expected increases in withdrawals for public supply. Monthly base-flow depletion criteria were determined using the Low-Flow Margin method, currently under consideration by NJDEP, to estimate available water on an annual basis at the Hydrologic Unit Code 11 (HUC11) level and to determine whether a water-supply deficit exists. Simulations of various groundwater-withdrawal scenarios were made using the calibrated model, and results were compared with baseline conditions (no withdrawals) to determine where and when base-flow deficits may be occurring and may be expected to occur in the future. Scenarios were simulated to assess base-flow depletion that could occur from different groundwater-withdrawal situations. In the Average scenario, deficits occurred in 7 of the 14 subbasins. In the Full Allocation scenario, deficits occurred in 11 of the subbasins. In the 2050 Demand scenario, deficits occurred in 9 of the 14 subbasins. The largest deficits occurred in the Absecon Creek subbasin because the base-flow depletion criteria for this subbasin is small due to the surface-water diversions that are already occurring there and because existing groundwater withdrawals in the subbasin have resulted in base-flow depletion under current (1998-2006) conditions. Three adjusted scenarios, variations of the Average, Full Allocation, and 2050 Demand scenarios, were simulated; for the adjusted scenarios, the withdrawals were modified in stages with the intent to successively eliminate or minimize the base-flow deficits. Modifications included shifting withdrawals to a deeper part of the Kirkwood-Cohansey aquifer system, implementing seasonal conjunctive use of shallow and deep aquifers, and specifying reductions in withdrawals within a HUC11 subbasin in deficit. The adjusted scenarios are intended to show the relative effectiveness of each of the three approaches in reducing the deficits. Most of the deficits under the Average, Full Allocation, and 2050 Demand scenarios were eliminated by reductions in withdrawals or allocations. Shifting withdrawals to a deeper part of the Kirkwood-Cohansey aquifer system or seasonal conjunctive use did not eliminate deficits for any subbasin. Reductions in withdrawals accounted for more than 95 percent of the total reduction of deficits in all but one subbasin.

New Jersey↗

Comparative analyses of hydrological responses of two adjacent watersheds to climate variability and change using the SWAT model

Water quality problems in the Chesapeake Bay Watershed (CBW) are expected to be exacerbated by climate variability and change. However, climate impacts on agricultural lands and resultant nutrient loads into surface water resources are largely unknown. This study evaluated the impacts of climate variability and change on two adjacent watersheds in the Coastal Plain of the CBW, using the Soil and Water Assessment Tool (SWAT) model. We prepared six climate sensitivity scenarios to assess the individual impacts of variations in CO 2 concentration (590 and 850 ppm), precipitation increase (11 and 21 %), and temperature increase (2.9 and 5.0 °C), based on regional general circulation model (GCM) projections. Further, we considered the ensemble of five GCM projections (2085–2098) under the Representative Concentration Pathway (RCP) 8.5 scenario to evaluate simultaneous changes in CO 2 , precipitation, and temperature. Using SWAT model simulations from 2001 to 2014 as a baseline scenario, predicted hydrologic outputs (water and nitrate budgets) and crop growth were analyzed. Compared to the baseline scenario, a precipitation increase of 21 % and elevated CO 2 concentration of 850 ppm significantly increased streamflow and nitrate loads by 50 and 52 %, respectively, while a temperature increase of 5.0 °C reduced streamflow and nitrate loads by 12 and 13 %, respectively. Crop biomass increased with elevated CO 2 concentrations due to enhanced radiation- and water-use efficiency, while it decreased with precipitation and temperature increases. Over the GCM ensemble mean, annual streamflow and nitrate loads showed an increase of ∼ 70 % relative to the baseline scenario, due to elevated CO 2 concentrations and precipitation increase. Different hydrological responses to climate change were observed from the two watersheds, due to contrasting land use and soil characteristics. The watershed with a larger percent of croplands demonstrated a greater increased rate of 5.2 kg N ha −1 in nitrate yield relative to the watershed with a lower percent of croplands as a result of increased export of nitrate derived from fertilizer. The watershed dominated by poorly drained soils showed increased nitrate removal due do enhanced denitrification compared to the watershed dominated by well-drained soils. Our findings suggest that increased implementation of conservation practices would be necessary for this region to mitigate increased nitrate loads associated with predicted changes in future climate.

Greensboro Watershed, Tuckahoe Creek Watershed↗

Crustal structure of the coastal and marine San Francisco Bay region, California

As of the time of this writing, the San Francisco Bay region is home to about 6.8 million people, ranking fifth among population centers in the United States. Most of these people live on the coastal lands along San Francisco Bay, the Sacramento River delta, and the Pacific coast. The region straddles the tectonic boundary between the Pacific and North American Plates and is crossed by several strands of the San Andreas Fault system. These faults, which are stressed by about 4 cm of relative plate motion each year, pose an obvious seismic hazard. We have many ways to study earthquake faults. Where faults break the land surface, we may learn valuable information needed for hazard assessment, such as cumulative offset, slip rate, and earthquake history. However, many of the major faults in the region are partly submerged beneath San Francisco and Monterey Bays. Although this situation poses problems in gathering observational data for hazard assessment, bay-region waterways provide an opportunity to study faultzone structure by using marine subsurface-imaging techniques, which are easier and cheaper than equivalent studies on land. In 1993, the U.S. Geological Survey (USGS) launched a 5-year project aimed at unearthing the basic science of the submerged San Andreas strike-slip fault system in the San Francisco Bay region with its many interacting strands. Primary project goals were structural, such as to discover how the San Andreas and Hayward Faults are connected or related at depth, to learn how the complex of faults in the San Andreas stepover zone on the Golden Gate platform functions, and to locate previously unknown faults. This volume thus contains mostly structural information about the San Francisco Bay region, much of it gathered through exploratory geophysical experiments. The volume is organized “top down,” from studies in the shallowest crust to the base of the crust. The first three chapters are linked through their use of novel geophysical techniques to study earthquake effects, coseismic slip, and shallow stratigraphy. Kayen and others examine crustal structure at very high resolution and demonstrate the use of ground-penetrating-radar tomography to measure the liquefaction potential of coastal sedimentary deposits. McGann and others use microfossils from drill cores along the San Francisco-Oakland Bay Bridge to determine a more detailed late Pleistocene stratigraphy of San Francisco Bay than was previously available. Geist and Zoback use the historical record of a small local tsunami generated by the great 1906 San Francisco earthquake to model the rupture process of that earthquake. The last four chapters are dedicated to studies of fault related structure of the seismogenic crust in and around the San Andreas Fault system in the San Francisco Bay region. Jachens and others compile an aeromagnetic anomaly map from new high-resolution flights across the bay region. Some of these anomalies mark the positions of offshore faults, and others are offset by faults, providing constraints on cumulative slip. Hart and others concisely summarize the marine seismic data recorded in and around San Francisco Bay, map the coverage, and provide archival information for those interested in acquiring data. The last two chapters present the results of the seismic data that have been analyzed. Bruns and others present their analysis of high-quality intermediate-resolution (~5-km penetration) seismic-reflection data gathered over the complex San Andreas-San Gregorio Fault junction. This junction, which is thought to be where the 1906 San Francisco earthquake originated (see Geist and Zoback, this volume), contains an apparent extensional right stepover in the San Andreas Fault. Finally, Parsons and others review and summarize the results of deep-crustal seismic-reflection experiments and local-earthquake tomographic studies, including previously unpublished data, and provide additional support and discussion for already-published studies. In summary, these studies were carried out in an environment where background information on faults in the San Francisco Bay region was sought. Much of the structural information presented here comes from experiments of a style unlikely to be conducted by the USGS in the near future. Together, the chapters in this volume provide a structural framework for a major part of a complex strike-slip fault system.

California↗

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↗