Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Ocean Modeling”

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 883 records · Page 49Linked to original sources

Scanning electron microscopy investigations of laboratory-grown gas clathrate hydrates formed from melting ice, and comparison to natural hydrates

Scanning electron microscopy (SEM) was used to investigate grain texture and pore structure development within various compositions of pure sI and sII gas hydrates synthesized in the laboratory, as well as in natural samples retrieved from marine (Gulf of Mexico) and permafrost (NW Canada) settings. Several samples of methane hydrate were also quenched after various extents of partial reaction for assessment of mid-synthesis textural progression. All laboratory-synthesized hydrates were grown under relatively high-temperature and high-pressure conditions from rounded ice grains with geometrically simple pore shapes, yet all resulting samples displayed extensive recrystallization with complex pore geometry. Growth fronts of mesoporous methane hydrate advancing into dense ice reactant were prevalent in those samples quenched after limited reaction below and at the ice point. As temperatures transgress the ice point, grain surfaces continue to develop a discrete "rind" of hydrate, typically 5 to 30 ??m thick. The cores then commonly melt, with rind microfracturing allowing migration of the melt to adjacent grain boundaries where it also forms hydrate. As the reaction continues under progressively warmer conditions, the hydrate product anneals to form dense and relatively pore-free regions of hydrate grains, in which grain size is typically several tens of micrometers. The prevalence of hollow, spheroidal shells of hydrate, coupled with extensive redistribution of reactant and product phases throughout reaction, implies that a diffusion-controlled shrinking-core model is an inappropriate description of sustained hydrate growth from melting ice. Completion of reaction at peak synthesis conditions then produces exceptional faceting and euhedral crystal growth along exposed pore walls. Further recrystallization or regrowth can then accompany even short-term exposure of synthetic hydrates to natural ocean-floor conditions, such that the final textures may closely mimic those observed in natural samples of marine origin. Of particular note, both the mesoporous and highly faceted textures seen at different stages during synthetic hydrate growth were notably absent from all examined hydrates recovered from a natural marine-environment setting.

American Mineralogist↗

Simulation of the effects of rainfall and groundwater use on historical lake water levels, groundwater levels, and spring flows in central Florida

The urbanization of central Florida has progressed substantially in recent decades, and the total population in Lake, Orange, Osceola, Polk, and Seminole Counties more than quadrupled from 1960 to 2010. The Floridan aquifer system is the primary source of water for potable, industrial, and agricultural purposes in central Florida. Despite increases in groundwater withdrawals to meet the demand of population growth, recharge derived by infiltration of rainfall in the well-drained karst terrain of central Florida is the largest component of the long-term water balance of the Floridan aquifer system. To complement existing physics-based groundwater flow models, artificial neural networks and other data-mining techniques were used to simulate historical lake water level, groundwater level, and spring flow at sites throughout the area. Historical data were examined using descriptive statistics, cluster analysis, and other exploratory analysis techniques to assess their suitability for more intensive data-mining analysis. Linear trend analyses of meteorological data collected by the National Oceanic and Atmospheric Administration at 21 sites indicate 67 percent of sites exhibited upward trends in air temperature over at least a 45-year period of record, whereas 76 percent exhibited downward trends in rainfall over at least a 95-year period of record. Likewise, linear trend analyses of hydrologic response data, which have varied periods of record ranging in length from 10 to 79 years, indicate that water levels in lakes (307 sites) were about evenly split between upward and downward trends, whereas water levels in 69 percent of wells (out of 455 sites) and flows in 68 percent of springs (out of 19 sites) exhibited downward trends. Total groundwater use in the study area increased from about 250 million gallons per day (Mgal/d) in 1958 to about 590 Mgal/d in 1980 and remained relatively stable from 1981 to 2008, with a minimum of 559 Mgal/d in 1994 and a maximum of 773 Mgal/d in 2000. The change in groundwater-use trend in the early 1980s and the following period of relatively slight trend is attributable to the concomitant effects of increasing public-supply withdrawals and decreasing use of water by the phosphate industry and agriculture. On the basis of available historical data and exploratory analyses, empirical lake water-level, groundwater-level, and spring-flow models were developed for 22 lakes, 23 wells, and 6 springs. Input time series consisting of various frequencies and frequency-band components of daily rainfall (1942 to 2008) and monthly total groundwater use (1957 to 2008) resulted in hybrid signal-decomposition artificial neural network models. The final models explained much of the variability in observed hydrologic data, with 43 of the 51 sites having coefficients of determination exceeding 0.6, and the models matched the magnitude of the observed data reasonably well, such that models for 32 of the 51 sites had root-mean-square errors less than 10 percent of the measured range of the data. The Central Florida Artificial Neural Network Decision Support System was developed to integrate historical databases and the 102 site-specific artificial neural network models, model controls, and model output into a spreadsheet application with a graphical user interface that allows the user to simulate scenarios of interest. Overall, the data-mining analyses indicate that the Floridan aquifer system in central Florida is a highly conductive, dynamic, open system that is strongly influenced by external forcing. The most important external forcing appears to be rainfall, which explains much of the multiyear cyclic variability and long-term downward trends observed in lake water levels, groundwater levels, and spring flows. For most sites, groundwater use explains less of the observed variability in water levels and flows than rainfall. Relative groundwater-use impacts are greater during droughts, however, and long-term trends in water levels and flows were identified that are consistent with historical groundwater-use patterns. The sensitivity of the hydrologic system to rainfall is expected, owing to the well-drained karst terrain and relatively thin confinement of the Floridan aquifer system in much of central Florida. These characteristics facilitate the relatively rapid transmission of infiltrating water from rainfall to the water table and contribute to downward leakage of water to the Floridan aquifer system. The areally distributed nature of rainfall, as opposed to the site-specific nature of groundwater use, and the generally high transmissivity and low storativity properties of the semiconfined Floridan aquifer system contribute to the prevalence of water-level and flow patterns that mimic rainfall patterns. In general, the data-mining analyses demonstrate that the hydrologic system in central Florida is affected by groundwater use differently during wet periods, when little or no system storage is available (high water levels), compared to dry periods, when there is excess system storage (low water levels). Thus, by driving the overall behavior of the system, rainfall indirectly influences the degree to which groundwater use will effect persistent trends in water levels and flows, with groundwater-use impacts more prevalent during periods of low water levels and spring flows caused by low rainfall and less prevalent during periods of high water levels and spring flows caused by high rainfall. Differences in the magnitudes of rainfall and groundwater use during wet and dry periods also are important determinants of hydrologic response. An important implication of the data-mining analyses is that rainfall variability at subannual to multidecadal timescales must be considered in combination with groundwater use to provide robust system-response predictions that enhance sustainable resource management in an open karst aquifer system. The data-driven approach was limited, however, by the confounding effects of correlation between rainfall and groundwater use, the quality and completeness of the historical databases, and the spatial variations in groundwater use. The data-mining analyses indicate that available historical data when used alone do not contain sufficient information to definitively quantify the related individual effects of rainfall and groundwater use on hydrologic response. The knowledge gained from data-driven modeling and the results from physics-based modeling, when compared and used in combination, can yield a more comprehensive assessment and a more robust understanding of the hydrologic system than either of the approaches used separately.

Florida↗

Estimating hydrodynamic roughness in a wave-dominated environment with a high-resolution acoustic Doppler profiler

Hydrodynamic roughness is a critical parameter for characterizing bottom drag in boundary layers, and it varies both spatially and temporally due to variation in grain size, bedforms, and saltating sediment. In this paper we investigate temporal variability in hydrodynamic roughness using velocity profiles in the bottom boundary layer measured with a high-resolution acoustic Doppler profiler (PCADP). The data were collected on the ebb-tidal delta off Grays Harbor, Washington, in a mean water depth of 9 m. Significant wave height ranged from 0.5 to 3 m. Bottom roughness has rarely been determined from hydrodynamic measurements under conditions such as these, where energetic waves and medium-to-fine sand produce small bedforms. Friction velocity due to current u*c and apparent bottom roughness z0a were determined from the PCADP burst mean velocity profiles using the law of the wall. Bottom roughness kB was estimated by applying the Grant-Madsen model for wave-current interaction iteratively until the model u*c converged with values determined from the data. The resulting kB values ranged over 3 orders of magnitude (10-1 to 10-4 m) and varied inversely with wave orbital diameter. This range of kB influences predicted bottom shear stress considerably, suggesting that the use of time-varying bottom roughness could significantly improve the accuracy of sediment transport models. Bedform height was estimated from kB and is consistent with both ripple heights predicted by empirical models and bedforms in sonar images collected during the experiment. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗

Relationship between water and aragonite barium concentrations in aquaria reared juvenile corals

Coral barium to calcium (Ba/Ca) ratios have been used to reconstruct records of upwelling, river and groundwater discharge, and sediment and dust input to the coastal ocean. However, this proxy has not yet been explicitly tested to determine if Ba inclusion in the coral skeleton is directly proportional to seawater Ba concentration and to further determine how additional factors such as temperature and calcification rate control coral Ba/Ca ratios. We measured the inclusion of Ba within aquaria reared juvenile corals ( Favia fragum ) at three temperatures (∼27.7, 24.6 and 22.5 °C) and three seawater Ba concentrations (73, 230 and 450 nmol kg −1 ). Coral polyps were settled on tiles conditioned with encrusting coralline algae, which complicated chemical analysis of the coral skeletal material grown during the aquaria experiments. We utilized Sr/Ca ratios of encrusting coralline algae (as low as 3.4 mmol mol −1 ) to correct coral Ba/Ca for this contamination, which was determined to be 26 ± 11% using a two end member mixing model. Notably, there was a large range in Ba/Ca across all treatments, however, we found that Ba inclusion was linear across the full concentration range. The temperature sensitivity of the distribution coefficient is within the range of previously reported values. Finally, calcification rate, which displayed large variability, was not correlated to the distribution coefficient. The observed temperature dependence predicts a change in coral Ba/Ca ratios of 1.1 μmol mol −1 from 20 to 28 °C for typical coastal ocean Ba concentrations of 50 nmol kg −1 . Given the linear uptake of Ba by corals observed in this study, coral proxy records that demonstrate peaks of 10–25 μmol mol −1 would require coastal seawater Ba of between 60 and 145 nmol kg −1 . Further validation of the coral Ba/Ca proxy requires evaluation of changes in seawater chemistry associated with the environmental perturbation recorded by the coral as well as verification of these results for Porites species, which are widely used in paleo reconstructions.

Geochimica et Cosmochimica Acta↗

Juvenile continental crust evolution in a modern oceanic arc setting: Petrogenesis of Cenozoic felsic plutons in Fiji, SW Pacific

Viti Levu, Fiji, provides one of the best exposed Phanerozoic analogues for the formation of juvenile continental crust in an intra-oceanic setting. Tonalites and trondhjemites are present in several large (75–150 km 2 ) adjacent, mid-Cenozoic plutons. We report major and trace element data including rare earth element (REE) and high-precision high field strength element (HFSE) compositions, new Hf-Nd-Sr-Pb isotope data, and zircon U/Pb-ages, O-Hf isotopes, and trace elements, from five different plutons. The Eocene Yavuna pluton and the Miocene Colo plutons are mainly composed of tonalites and trondhjemites and represent the exposed middle crust of the former Vitiaz island arc. The plutons can be divided into three suites. One suite is light REE (LREE) depleted with some trace element ratios lower than average normal mid-ocean ridge basalts (N-MORB). A second suite has flat REE patterns similar to local island arc basalts. Both suites occur near the coast of Viti Levu, include a wide compositional spectrum from gabbro to tonalite, and can be produced mostly by fractional crystallization of mafic precursor melts. The third suite is characterized by LREE enrichments with higher La N /Yb N (2.3–4.9), higher Zr/Y (4.3–7.1), and lower Nb/Ta (9.6–12.4). They occur closer to the center of the island and are bimodal trondhjemite-gabbro intrusions. These characteristics are consistent with formation mostly by partial melting of mafic crust. Trace element modeling shows that the trace element ratios of the third suite can be produced by 10–20 % melting of the mafic crust in the presence of residual amphibole, resulting in the retention of the medium REE (MREE) and diagnostic trace element ratios including low Nb/Ta and high Zr/Y. Geochemical similarities of the LREE enriched suite to typical “low”-pressure Archean tonalites-trondhjemites-granodiorites (TTGs) imply a common petrogenetic origin and similar mechanisms for the generation of juvenile Archean and modern differentiated crust by partial melting of mafic crust with residual amphibole. In modern oceanic arcs, genetically unrelated felsic plutonic as well as volcanic rocks co-exist, and in this regard, the Fijian plutons accompany major tectonic disruptions to arc processes.

Geochimica et Cosmochimica Acta↗

Passive acoustic monitoring and convolutional neural networks facilitate high-resolution and broadscale monitoring of a threatened species

Population monitoring is an essential component of biodiversity conservation and management, but low detection probabilities for rare and/or cryptic species makes estimating abundance and occupancy challenging. Passive acoustic monitoring combined with machine learning algorithms represents a potential path forward to effectively and efficiently monitor the occurrence of rare vocalizing species across entire forest landscapes. Our objectives were to develop and implement a convolutional neural network (PNW-Cnet) to identify vocalizations of a rare and threatened forest nesting bird species – the marbled murrelet ( Brachyramphus marmoratus ) – in the Pacific Northwest, U.S.A., 2018–2021. We used PNW-Cnet predictions from broadscale passive acoustic monitoring data to examine spatiotemporal patterns in the distribution of murrelets. PNW-Cnet showed sufficiently high prediction accuracy (overall precision > 0.9) to enable broadscale population monitoring. Spatiotemporal analysis showed that annual peak murrelet call abundance occurs in ordinal weeks 28–32 (late July–Mid August) but this varied by study area. The greatest number of detections typically occurred in the Olympic Peninsula and Oregon Coast Range where late-successional forest dominates and nearer to ocean habitats. We demonstrate that passive acoustic monitoring can be used to understand intensity of use across broad scales for a rare and cryptic species in addition to the typical detection/non-detection data that are often collected. Passive acoustic monitoring combined with PNW-Cnet offers considerable promise for species distribution modeling and long-term population monitoring for rare species.

Oregon, Washington↗

Stress magnitudes in the crust: constraints from stress orientation and relative magnitude data

The World Stress Map Project is a global cooperative effort to compile and interpret data on the orientation and relative magnitudes of the contemporary in situ tectonic stress field in the Earth's lithosphere. The intraplate stress field in both the oceans and continents is largely compressional with one or both of the horizontal stresses greater than the vertical stress. The regionally uniform horizontal intraplate stress orientations are generally consistent with either relative or absolute plate motions indicating that plate-boundary forces dominate the stress distribution within the plates. Current models of stresses due to whole mantle flow inferred from seismic topography models predict a general compressional stress state within continents but do not match the broad-scale horizontal stress orientations. The broad regionally uniform intraplate stress orientations are best correlated with compressional plate-boundary forces and the geometry of the plate boundaries. -from Authors

Philosophical Transactions of the Royal Society A:↗

Modeling the morphological response of a barrier island to Hurricane Matthew

Surge and wave forcing from Hurricane Matthew caused a breach south of Matanzas Inlet (FL, USA) on a complex barrier island, including sandy dunes, hard structures (residential buildings and a highway), wetlands, and the US Intracoastal Waterway. In this paper, the skill of the XBeach model to predict hurricane-induced barrier island overwash, dune erosion, and breaching is demonstrated. The location of the breach is predicted correctly if bottom roughness based on land cover is used to calculate bed shear stresses. While the dunes are initially lowered by wave attack and surge from the ocean side, the main driver for breach formation is the water level difference between the back-barrier and nearshore, causing an ocean-directed outflow of water after the peak of the storm.

Florida↗

The geobiology of sediment-hosted mineral deposits

The role of biological processes in the formation of sediment-hosted ore deposits has long been recognized. In this review, we focus on the biogeochemical cycling of C, Mn, Fe, and S as they relate to the formation of sediment-hosted Mn and Fe deposits, metalliferous black shales, clastic-dominated (CD) Pb-Zn deposits, and phosphorites. Biological mediation of ore-forming processes occurs over large spans of space and time. The most important step is oxygenic photosynthesis, a biological innovation dating from the Archean Eon that releases free oxygen into the surface oceans and atmosphere and delivers chemical potential, in the form of reduced carbon, to the seafloor. Photosynthetic oxygen is available to precipitate dissolved Fe 2 + and Mn 2 +, and therefore it augments the formation of sedimentary Mn and Fe deposits, and drives oxidative weathering of exposed crust, thereby delivering sulfate and transition metals to the ocean. Where reduced carbon accumulates in the deep oceans and on the seafloor, bacterial sulfate reduction produces hydrogen sulfide thereby facilitating the formation of metalliferous black shales, sediment-hosted Pb and Zn sulfide deposits, and phosphorites. Thus, an understanding of major biogeochemical processes and how they have evolved over time is required in order to refine genetic models for sediment-hosted ore deposits and to guide future mineral exploration. A close secular relationship between deposit formation and trends in major biogeochemical cycles provides a potentially powerful tool for mineral resource assessment. Sedimentary basins that formed during a time that is known to lack deposits of a particular metal can be eliminated during exploration programs, whereas others of permissive ages should be considered priorities. For example, sedimentary basins older than ca. 1.8 Ga are unlikely to contain large CD Pb-Zn deposits, and basins that formed between 1.6 and 0.6 Ga are not prospective for phosphorites. Recent technological advances in the application of nanometer-, micron-, and bulk-scale analytical techniques allow for imaging of complex biological structures and have provided new insights into the role of bacteria, not only in direct formation of mineral deposits, but also in leaching of metals from ore and mineralized rocks. Future exploration for, and exploitation of, mineral deposits may include offshore or land-based, low-grade, high-tonnage targets; understanding the role of bacteria in mineral growth, mineral dissolution, and redox transformations will aid in predicting where such deposits exist, and how metal extraction from ores can be enhanced.

Book chapter↗

Isotopic composition of low-latitude paleoprecipitation during the Early Cretaceous

The response of the hydrologic cycle in global greenhouse conditions is important to our understanding of future climate change and to the calibration of global climate models. Past greenhouse conditions, such as those of the Cretaceous, can be used to provide empirical data with which to evaluate climate models. Recent empirical studies have utilized pedogenic carbonates to estimate the isotopic composition of meteoric waters and calculate precipitation rates for the AptianAlbian. These studies were limited to data from mid(35??N) to high (75??N) paleolatitudes, and thus future improvements in accuracy will require more estimates of meteoric water compositions from numerous localities around the globe. This study provides data for tropical latitudes (18.5??N paleolatitude) from the Tlayua Formation, Puebla, Mexico. In addition, the study confirms a shallow nearshore depositional environment for the Tlayua Formation. Petrographic observations of fenestral fabrics, gypsum crystal molds, stromatolitic structures, and pedogenic matrix birefringence fabric support the interpretation that the strata represent deposition in a tidal flat environment. Carbonate isotopic data from limestones of the Tlayua Formation provide evidence of early meteoric diagenesis in the form of meteoric calcite lines. These trends in ?? 18 O versus ?? 13 C were used to calculate the mean ?? 18 O value of meteoric water, which is estimated at -5.46 ?? 0.56??? (Vienna Standard Mean Ocean Water [VSMOW]). Positive linear covariant trends in oxygen and carbon isotopic values from some horizons were used to estimate evaporative losses of vadose groundwater from tropical exposure surfaces during the Albian, and the resulting values range from 8% to 12%. However, the presence of evaporative mineral molds indicates more extensive evaporation. The added tropical data improve latitudinal coverage of paleoprecipitation ??18O estimates. The data presented here imply that earlier isotope mass balance models most likely underestimated tropical to subtropical precipitation and evaporation fluxes. The limited latitudinal constraints for earlier isotope mass balance modeling of the Albian hydrologic cycle of the Northern Hemisphere Americas resulted in extrapolated low-latitude precipitation ?? 18 O values that were much heavier (up to 3???) than the values observed in this study. The lighter values identified in this study indicate a more pronounced rainout effect for tropical regions and quite possibly a more vigorous evaporation effect. These and additional low-latitude data are required to better constrain changes in the hydrologic cycle during the Cretaceous greenhouse period, and to reduce the uncertainties resulting from limited geographic coverage of proxy data. ?? 2009 Geological Society of America.

Geological Society of America Bulletin↗

Enhancing assessments of coastal wetland migration potential with sea-level rise: Accounting for uncertainty in elevation data, tidal data, and future water levels

Sea-level rise rates are predicted to surpass rates of wetland vertical adjustment in the coming decades in many areas, increasing the potential for wetland submergence. Information on where wetland migration is possible can help natural resource managers for planning land acquisition or enhancing habitat connectivity to bolster adaptation of coastal wetlands to rising seas. Elevation-based models of wetland migration are often hampered with uncertainty associated with ground surface elevation, current water levels (i.e., tides and extreme water levels), and future water levels from sea-level rise. Here, we developed an approach that involved digital elevation model error reduction and the use of Monte Carlo simulations that utilize uncertainty assumptions regarding elevation error, contemporary water levels, and future sea levels to identify potential wetland migration areas. Our analyses were developed for Duvall and Nassau Counties in northeastern Florida (USA). We focus on the migration of regularly oceanic-flooded wetlands (i.e., flooded by oceanic water daily) and irregularly oceanic-flooded wetlands (i.e., flooded by oceanic water less frequently than daily). For two relative sea-level rise scenarios based on the 0.5 m and the 1.5 m global mean sea-level rise scenarios, we quantified migration by wetland flooding frequency class and identified land cover and land use types that are vulnerable to future exposure to oceanic waters. The variability in total coverage and relative coverage of wetland migration from our results highlights how topography and accelerated sea-level rise interact. Our wetland migration results communicate uncertainty by showing flooding frequency class as probabilistic outputs.

Florida↗

Geological evidence for solid-state convection in Europa's ice shell

The ice-rich surface of the jovian satellite Europa is sparsely cratered, suggesting that this moon might be geologically active today. Moreover, models of the satellite's interior indicate that tidal interactions with Jupiter might produce enough heat to maintain a subsurface liquid water layer. But the mechanisms of interior heat loss and resurfacing are currently unclear, as is the question of whether Europa has (or had at one time) a liquid water ocean. Here we report on the morphology and geological interpretation of distinct surface features-pits, domes and spots-discovered in high-resolution images of Europa obtained by the Galileo spacecraft. The features are interpreted as the surface manifestation of diapirs, relatively warm localized ice masses that have risen buoyantly through the subsurface. We find that the formation of the features can be explained by thermally induced solid-state convection within an ice shell, possibly overlying a liquid water layer. Our results are consistent with the possibility that Europa has a liquid water ocean beneath a surface layer of ice, but further tests and observations are needed to demonstrate this conclusively.

Nature↗

U-Pb zircon geochronology and geological evolution of the Halaban- Al Amar region of the eastern Arabian Shield, Kingdom of Saudi Arabia.

U-Pb zircon model ages for eleven major units from the Halaban-Al Amar region of the eastern Arabian Shield indicate three stages of evolution: (1) plate convergence, (2) plate collision, and (3) post-orogenic intracratonic activity. Convergence occurred between the western Afif and eastern Ar Rayn plates that were separated by oceanic crust. Remnants of oceanic crust now comprise the ophiolitic complexes of the Urd group. The oldest plutonic unit in the study is from one of these complexes and gave an age of 694 ± 8 Ma. Detrital zircons from the sedimentary Abt formation of the Urd group, which is intercalated with the ophiolitic rocks, were derived from source rocks with a mean age of 710 Ma. The Abt formation may be an accretionary wedge on the western margin of the Ar Rayn plate. Plate convergence was terminated by collision of the Afif and Ar Rayn plates during the Al Amar orogeny which began about 670 Ma. During collision, the Urd group rocks were deformed and in part obducted on to one or both plates. Synorogenic leucogranitoid rocks were intruded from 670 to 640 Ma. From about 640 to 630 Ma, widespread unfoliated dioritic plutons were emplaced in the Ar Rayn block, and represent the end of orogenesis related to collision. There is no definitive evidence for a significantly older basement beneath the study region.

Journal of the Geological Society↗

Coupled effects of vertical mixing and benthic grazing on phytoplankton populations in shallow, turbid estuaries

Coastal ocean waters tend to have very different patterns of phytoplankton biomass variability from the open ocean, and the connections between physical variability and phytoplankton bloom dynamics are less well established for these shallow systems. Predictions of biological responses to physical variability in these environments is inherently difficult because the recurrent seasonal patterns of mixing are complicated by aperiodic fluctuations in river discharge and the high-frequency components of tidal variability. We might expect, then, less predictable and more complex bloom dynamics in these shallow coastal systems compared with the open ocean. Given this complex and dynamic physical environment, can we develop a quantitative framework to define the physical regimes necessary for bloom inception, and can we identify the important mechanisms of physical-biological coupling that lead to the initiation and termination of blooms in estuaries and shallow coastal waters? Numerical modeling provides one approach to address these questions. Here we present results of simulation experiments with a refined version of Cloern's (1991) model in which mixing processes are treated more realistically to reflect the dynamic nature of turbulence generation in estuaries. We investigated several simple models for the turbulent mixing coefficient. We found that the addition of diurnal tidal variation to Cloern's model greatly reduces biomass growth indicating that variations of mixing on the time scale of hours are crucial. Furthermore, we found that for conditions representative of South San Francisco Bay, numerical simulations only allowed for bloom development when the water column was stratified and when minimal mixing was prescribed in the upper layer. Stratification, however, itself is not sufficient to ensure that a bloom will develop: minimal wind stirring is a further prerequisite to bloom development in shallow turbid estuaries with abundant populations of benthic suspension feeders.

Journal of Marine Research↗

Assessment of water quality and discharge in the Herring River, Wellfleet, Massachusetts, November 2015 to September 2017

The U.S. Geological Survey, Cape Cod National Seashore of the National Park Service, and Friends of Herring River cooperated from 2015 to 2017 to assess nutrient concentrations and fluxes across the ocean-estuary boundary at a dike on the Herring River in Wellfleet, Massachusetts. The purpose of this assessment was to characterize environmental conditions prior to a future removal of the dike, which has restricted saltwater inputs into the Herring River watershed for more than 100 years. Water temperature, dissolved oxygen, pH, and specific conductance were monitored continuously, and flow-weighted composite samples were collected approximately twice per month at the ocean-estuary boundary. Bidirectional discharge was computed for the U.S. Geological Survey Herring River at Chequessett Neck Road at Wellfleet, Massachusetts, streamgage (011058798) by using a stage-area rating and index-velocity ratings developed with acoustic Doppler current profile measurements made upstream and downstream from the dike. LOADEST regression modeling software was used to estimate nutrient fluxes (loads) from composite, paired nutrient concentration and discharge data in conjunction with continuous discharge data. Temperature, dissolved oxygen, pH, and specific conductance were also monitored continuously on two tributaries to the Herring River, Pole Dike Creek and Bound Brook, from late-May 2016 to mid-June 2017. Composite or discrete water samples were collected from the tributaries approximately twice per month in most months from late-May 2016 to mid-June 2017 and analyzed for total nitrogen, total phosphorus, and dissolved organic carbon. Flow-weighted concentrations of ammonium, nitrate, and total nitrogen on the Herring River at the dike on the ebb tide generally varied between 0.01 and 0.1, 0.003 and 0.03, and 0.3 and 0.7 milligram per liter as nitrogen, respectively. Flow-weighted concentrations of orthophosphate, total dissolved phosphorus, and total phosphorus generally varied between 0.002 and 0.02, 0.003 and 0.06, and 0.03 and 0.1 milligram per liter as phosphorus, respectively, on the ebb tide. Flow-weighted concentrations of silicate and dissolved organic carbon on the ebb tide generally varied between 0.08 and 3.0 milligrams per liter of silica (silicon dioxide), and 1.7 and 5.6 milligrams per liter of carbon, respectively. Ebb tide concentrations of nitrate were highest in winter and lowest in summer. By contrast, ebb tide concentrations of phosphorus species were highest in late summer and early fall and lowest in winter. Silica and dissolved organic carbon did not exhibit systematic variation in seasonal concentrations. There was uncertainty in estimates of nutrient fluxes, but the LOADEST-estimated fluxes indicated that annual (and in almost all cases seasonal) exports (ebb tides) exceeded inputs (flood tides). Ebb tide concentrations of ammonium, nitrate, total nitrogen, and silica were positively correlated with antecedent cumulative 7-day precipitation.

Massachusetts↗

Topobathymetric model of Mobile Bay, Alabama

Topobathymetric Digital Elevation Models (DEMs) are a merged rendering of both topography (land elevation) and bathymetry (water depth) that provides a seamless elevation product useful for inundation mapping, as well as for other earth science applications, such as the development of sediment-transport, sea-level rise, and storm-surge models. This 1/9-arc-second (approximately 3 meters) resolution model of Mobile Bay, Alabama was developed using multiple topographic and bathymetric datasets, collected on different dates. The topographic data were obtained primarily from the U.S. Geological Survey (USGS) National Elevation Dataset (NED) (http://ned.usgs.gov/) at 1/9-arc-second resolution; USGS Experimental Advanced Airborne Research Lidar (EAARL) data (2 meters) (http://pubs.usgs.gov/ds/400/); and topographic lidar data (2 meters) and Compact Hydrographic Airborne Rapid Total Survey (CHARTS) lidar data (2 meters) from the U.S. Army Corps of Engineers (USACE) (http://www.csc.noaa.gov/digitalcoast/data/coastallidar/). Bathymetry was derived from digital soundings obtained from the National Oceanic and Atmospheric Administration’s (NOAA) National Geophysical Data Center (NGDC) (http://www.ngdc.noaa.gov/mgg/geodas/geodas.html) and from water-penetrating lidar sources, such as EAARL and CHARTS. Mobile Bay is ecologically important as it is the fourth largest estuary in the United States. The Mobile and Tensaw Rivers drain into the bay at the northern end with the bay emptying into the Gulf of Mexico at the southern end. Dauphin Island (a barrier island) and the Fort Morgan Peninsula form the mouth of Mobile Bay. Mobile Bay is 31 miles (50 kilometers) long by a maximum width of 24 miles (39 kilometers) with a total area of 413 square miles (1,070 square kilometers). The vertical datum of the Mobile Bay topobathymetric model is the North American Vertical Datum of 1988 (NAVD 88). All the topographic datasets were originally referenced to NAVD 88 and no transformations were made to these input data. The NGDC hydrographic, multibeam, and trackline surveys were transformed from mean low water (MLW) or mean lower low water (MLLW) to NAVD 88 using VDatum (http://vdatum.noaa.gov). VDatum is a tool developed by the National Geodetic Survey (NGS) that performs transformations among tidal, ellipsoid-based, geoid-based, and orthometric datums using calibrated hydrodynamic models. The vertical accuracy of the input topographic data varied depending on the input source. Because the input elevation data were derived primarily from lidar, the vertical accuracy ranges from 6 to 20 centimeters in root mean square error (RMSE). he horizontal datum of the Mobile Bay topobathymetric model is the North American Datum of 1983 (NAD 83), geographic coordinates. All the topographic and bathymetric datasets were originally referenced to NAD 83, and no transformations were made to the input data. The bathymetric surveys were downloaded referenced to NAD 83 geographic, and therefore no horizontal transformations were required. The topbathymetric model of Mobile Bay and detailed metadata can be obtained from the USGS Web sites: http://nationalmap.gov/.

Alabama↗

Geophysical modeling of the northern Appalachian Brompton-Cameron, Central Maine, and Avalon terranes under the New Jersey Coastal Plain

A regional terrane map of the New Jersey Coastal Plain basement was constructed using seismic, drilling, gravity and magnetic data. The Brompton-Cameron and Central Maine terranes were coalesced as one volcanic island arc terrane before obducting onto Laurentian, Grenville age, continental crust in the Taconian orogeny [Rankin, D.W., 1994. Continental margin of the eastern United States: past and present. In: Speed, R.C., (Ed.), Phanerozoic Evolution of North American Continent-Ocean Transitions. DNAG Continent-Ocean Transect Volume. Geological Society of America, Boulder, Colorado, pp. 129-218]. Volcanic island-arc rocks of the Avalon terrane are in contact with Central Maine terrane rocks in southern Connecticut where the latter are overthrust onto the Brompton-Cameron terrane, which is thrust over Laurentian basement. Similarities of these allochthonous island arc terranes (Brompton-Cameron, Central Maine, Avalon) in lithology, fauna and age suggest that they are faulted segments of the margin of one major late Precambrian to early Paleozoic, high latitude peri-Gondwana island arc designated as "Avalonia", which collided with Laurentia in the early to middle Paleozoic. The Brompton Cameron, Central Maine, and Avalon terranes are projected as the basement under the eastern New Jersey Coastal Plain based on drill core samples of metamorphic rocks of active margin/magmatic arc origin. A seismic reflection profile across the New York Bight traces the gentle dipping (approximately 20 degrees) Cameron's Line Taconian suture southeast beneath allochthonous Avalon and other terranes to a 4 sec TWTT depth (approximately 9 km) where the Avalonian rocks are over Laurentian crust. Gentle up-plunge (approximately 5 degrees) projections to the southwest bring the Laurentian Grenville age basement and the drift-stage early Paleozoic cover rocks to windows in Burlington Co. at approximately 1 km depth and Cape May Co. at approximately 2 km depths. The antiformal Shellburne Falls and Chester domes and Chain Lakes-Pelham dome-Bronson Hill structural trends, and the synformal Connecticut Valley-Gaspe structural trend can be traced southwest into the New Jersey Coastal Plain basement. A Mesozoic rift basin, the "Sandy Hook basin", and associated eastern boundary fault is identified, based upon gravity modeling, in the vicinity of Sandy Hook, New Jersey. The thickness of the rift-basin sedimentary rocks contained within the "Sandy Hook basin" is approximately 4.7 km, with the basin extending offshore to the east of the New Jersey coast. Gravity modeling indicates a deep rift basin and the magnetic data indicates a shallow magnetic basement caused by magnetic diabase sills and/or basalt flows contained within the rift-basin sedimentary rocks. The igneous sills and/or flows may be the eastward continuation of the Watchung and Palisades bodies. ?? 2004 Elsevier Ltd. All rights reserved.

Journal of Geodynamics↗

Structural control on megathrust rupture and slip behavior: Insights from the 2016 Mw 7.8 Pedernales Ecuador earthquake

The heterogeneous seafloor topography of the Nazca Plate as it enters the Ecuador subduction zone provides an opportunity to document the influence of seafloor roughness on slip behavior and megathrust rupture. The 2016 Mw 7.8 Pedernales Ecuador earthquake was followed by a rich and active postseismic sequence. An internationally coordinated rapid response effort installed a temporary seismic network to densify coastal stations of the permanent Ecuadorian national seismic network. A combination of 82 onshore short and intermediate period and broadband seismic stations and six ocean bottom seismometers recorded the postseismic Pedernales sequence for over a year after the mainshock. A robust earthquake catalog combined with calibrated relocations for a subset of magnitude ≥4 earthquakes shows pronounced spatial and temporal clustering. A range of slip behavior accommodates postseismic deformation including earthquakes, slow slip events, and earthquake swarms. Models of plate coupling and the consistency of earthquake clustering and slip behavior through multiple seismic cycles reveal a segmented subduction zone primarily controlled by subducted seafloor topography, accreted terranes, and inherited structure. The 2016 Pedernales mainshock triggered moderate to strong earthquakes (5 ≤ M ≤ 7) and earthquake swarms north of the mainshock rupture close to the epicenter of the 1906 Mw 8.8 earthquake and in the segment of the subduction zone that ruptured in 1958 in a Mw 7.7 earthquake.

Journal of Geophysical Research: Earth Surface↗