USGS Science⌕ Search

SEARCH · USGS Science

Results for “Ocean Modelling”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,153 records · Page 64Linked to original sources
USGS2016

Forcing and variability of nonstationary rip currents

Surface wave transformation and the resulting nearshore circulation along a section of coast with strong alongshore bathymetric gradients outside the surf zone are modeled for a consecutive 4 week time period. The modeled hydrodynamics are compared to in situ measurements of waves and currents collected during the Nearshore Canyon Experiment and indicate that for the entire range of observed conditions, the model performance is similar to other studies along this stretch of coast. Strong alongshore wave height gradients generate rip currents that are observed by remote sensing data and predicted qualitatively well by the numerical model. Previous studies at this site have used idealized scenarios to link the rip current locations to undulations in the offshore bathymetry but do not explain the dichotomy between permanent offshore bathymetric features and intermittent rip current development. Model results from the month‐long simulation are used to track the formation and location of rip currents using hourly statistics, and results show that the direction of the incoming wave energy strongly controls whether rip currents form. In particular, most of the offshore wave spectra were bimodal and we find that the ratio of energy contained in each mode dictates rip current development, and the alongshore rip current position is controlled by the incident wave period. Additionally, model simulations performed with and without updating the nearshore morphology yield no significant change in the accuracy of the predicted surf zone hydrodyanmics indicating that the large‐scale offshore features (e.g., submarine canyon) predominately control the nearshore wave‐circulation system.

Journal of Geophysical Research C: Oceans↗
USGS1993

Oxygen buffering of Kilauea volcanic gases and the oxygen fugacity of Kilauea basalt

Volcanic gases collected during episode 1 of the Puu Oo eruption along the east rift zone of Kilauea Volcano, Hawaii, have uniform C-O-H-S-Cl-F compositions that are sharply depleted in CO 2 . The CO 2 -poor gases are typical of Type II volcanic gases ( gerlach and Graeber , 1985) and were emitted from evolved magma stored for a prolonged period of time in the east rift zone after releasing CO 2 -rich gases during an earlier period of temporary residence in the summit magma chamber. The samples are remarkably free of contamination by atmospheric gases and meteoric water. Thermodynamic evaluation of the analytical data shows that the episode 1 gases have equilibrium compositions appropriate for temperatures between 935 and 1032°C. Open- and closed-system equilibrium models of species distributions for the episode 1 gases show unequivocally that coexisting lavas buffered the gas oxygen fugacities during cooling. These models indicate that the f O 2 buffering process occurs by transfer of oxygen from the major species in the gas phase (H 2 O, CO 2 , SO 2 ) to the lava during cooling and that the transfer of oxygen also controls the fugacities of several minor and trace species (H 2 , CO, H 2 S, S 2 , Cl 2 , F 2 ), in addition to O 2 during cooling. Gas/lava exchanges of other components are apparently insignificant and exert little influence, compared to oxygen exchange, during cooling. Oxygen transfer during cooling is variable, presumably reflecting short-term fluctuations in gas flow rates. Higher flow rates restrict the time available for gas/lava oxygen transfer and result in gases with higher equilibrium temperatures. Lower flow rates favor f O 2 -constrained equilibration by oxygen transfer down to lower temperatures. Thus, the chemical equilibrium preserved in these gases is a heterogeneous equilibrium constrained by oxygen fugacity, and the equilibrium temperatures implied by the compositions of the gases reflect the temperatures at which gas/lava oxygen exchange ceased. This conclusion challenges the common assumption that volcanic gases are released from lava in a state of chemical equilibrium and then continue equilibrating homogeneously with falling temperature until reaction rates are unable to keep pace with cooling. No evidence is found, moreover, that certain gas species are kinetically more responsive and able to equilibrate down to lower temperatures than those of the last gas/lava oxygen exchange. Homogeneous reaction rates in the gas phase are apparently slow compared to the time it took for the gases to move from the last site of gas/lava equilibration to the site of collection. An earlier set of data for higher temperature CO 2 -rich Type I volcanic gases, which come from sustained summit lava lake eruptions supplied by magma that experienced substantially shorter periods of crustal storage, shows f O 2 buffering by oxygen transfer up to 1185°C. Oxygen fugacity measurements in drill holes into ponded lava flows suggest that buffering by oxygen transfer may control the f O 2 of residual gases down to several hundred degrees below the solidus in the early stages of cooling. Although the details of the f O 2 buffering mechanisms for oxygen transfer are unknown, the fact that f O 2 buffering is effective from molten to subsolidus conditions suggests that the reaction mechanisms must change with cooling as the reactants change from predominantly melt, to melt plus crystals, to glass plus crystals. Mass balance calculations suggest that redox reactions between the gas and ferrous/ferric iron in the lava are plausible mechanisms for the oxygen transfer and that the f O 2 of the gases is buffered by sliding ferrous/ferric equilibria in the erupting lavas. Contrary to expectations based on models predicting the oxidation of basalt by H 2 and CO escape during crustal storage, CO 2 -rich Type I gases and CO 2 -poor Type II gases have identical oxygen fugacities despite greatly different crustal storage and degassing histories. Volcanic gas data give a tightly constrained log f O 2 of NNO − 0.5 (±0.05) for subaerially erupted Kilauea basalt from liquidus to solidus temperatures, consistent with recent f O 2 determinations for the mantle source regions of ocean island basalts. Because the oxygen fugacities of volcanic gases emitted by subaerial lavas imply that the f O 2 of Kilauea basalt is unchanged during crustal storage, Kilauea basalt either arrives in the crust with an oxygen fugacity between NNO and FMQ, or it develops an oxygen fugacity in this range immediately upon arrival in the summit chamber.

Geochimica et Cosmochimica Acta↗
USGS2007

Multibeam observations of mine burial near Clearwater, FL, including comparisons to predictions of wave-induced burial

A Kongsberg Simrad EM 3000 multibeam sonar (Kongsberg Simrad, Kongsberg, Norway) was used to conduct a set of six repeat high-resolution bathymetric surveys west of Indian Rocks Beach (IRB), just to the south of Clearwater, FL, between January and March 2003, to observe in situ scour and burial of instrumented inert mines and mine-like cylinders. Three closely located study sites were chosen: two fine-sand sites, a shallow one located in ??? 13 m of water depth and a deep site located in ???14 m of water depth; and a coarse-sand site in ???13 m. Results from these surveys indicate that mines deployed in fine sand are nearly buried within two months of deployment (i.e., they sunk 74.5% or more below the ambient seafloor depth). Mines deployed in coarse sand showed a lesser amount of scour, burying until they present roughly the same hydrodynamic roughness as the surrounding rippled bedforms. These data were also used to test the validity of the Virginia Institute of Marine Science (VIMS, Gloucester Point, VA) 2-D burial model. The model worked well in areas of fine sand, sufficiently predicting burial over the course of the experiment. In the area of coarse sand, the model greatly overpredicted the amount of burial. This is believed to be due to the presence of rippled bedforms around the mines, which affect local bottom morphodynamics and are not accounted for in the model, an issue currently being addressed by the modelers. This paper focuses specifically on two instrumented mines: an acoustic mine located in fine sand and an optical instrumented mine located in coarse sand. ?? 2007 IEEE.

IEEE Journal of Oceanic Engineering↗
USGS2024

Estimating migration timing and abundance in partial migratory systems by integrating continuous antenna detections with physical captures

Many populations migrate between two different habitats (e.g. wintering/foraging to breeding area, mainstem–tributary, river–lake, river–ocean, river–side channel) as part of their life history. Detection technologies, such as passive integrated transponder (PIT) antennas or sonic receivers, can be placed at boundaries between habitats (e.g. near the confluence of rivers) to detect migratory movements of marked animals. Often, these detection systems have high detection probabilities and detect many individuals but are limited in their ability to make inferences about abundance because only marked individuals can be detected. Here, we introduce a mark–recapture modelling approach that uses detections from a double-array PIT antenna system to imply movement directionality from arrays and estimate migration timing. Additionally, when combined with physical captures, the model can be used to estimate abundances for both migratory and non-migratory groups and help quantify partial migration. We first test our approach using simulation, and results indicate our approach displayed negligible bias for total abundance (less than ±1%) and slight biases for state-specific abundance estimates (±1%–6%). We fit our model to array detections and physical captures of three native fishes (humpback chub [ Gila cypha ], flannelmouth sucker [ Catostomus latipinnis ] and bluehead sucker [ Catostomus discobolus ]) in the Little Colorado River (LCR) in Grand Canyon, AZ, a system that exhibits partial migration (i.e. includes residents and migrants). Abundance estimates from our model confirm that, for all three species, migratory individuals are much more numerous than residents. There was little difference in movement timing between 2021 (a year without preceding winter/spring floods) and 2022 (a year with a small flood occurring in early April). In both years, flannelmouth sucker arrived in mid-March whereas humpback chub and bluehead sucker arrivals occurred early- to mid-April. With humpback chub and flannelmouth sucker, movement timing was influenced by body size so that large individuals were more likely to arrive early compared to smaller individuals. With more years of data, this model framework could be used to evaluate ecological questions pertaining to flow cues and movement timing or intensity, relative trends in migrants versus residents and ecological drivers of skipped spawning.

Journal of Animal Ecology↗
USGS2025

Alaskan hydrology in transition: Changing precipitation and evapotranspiration patterns are projected to reshape seasonal streamflow and water temperature by midcentury (2035-2064)

High spatial and temporal resolution models are essential for understanding future climate impacts and developing effective climate resilience plans. However, existing regional and global river models often lack the resolution needed to accurately capture local conditions. This study uses a series of high-resolution models, including the Regional Arctic System Model, mizuRoute, and the river basin model, to analyze Arctic and sub-Arctic Alaskan hydrology. We compare a historical baseline (1991–2020) with six midcentury (2035–64) futures: two pseudo–global warming scenarios based on historical meteorology and four direct dynamically downscaled global climate models. The six futures reveal significant uncertainty in future annual discharge and peak flows, although a widespread increase in discharge during April (+63%) and October (+31%) is consistently shown across models. Projected increases in rain and shifting weather patterns lead to a transition from snow to rain in spring and autumn, reducing the fraction of snowmelt contributing to river discharge. Rising evapotranspiration moderates discharge changes, particularly in autumn, by offsetting precipitation increases. Average summer river temperatures are projected to increase by approximately 1.5°C, doubling the number of river segments that experience 18°C days, a critical threshold for salmon survival, and intensifying the heat flux to the ocean adding an average of 3.3 × 10 12 MJ yr −1 . These changes in the hydrologic cycle could profoundly impact riverine and oceanic ecosystems, posing substantial challenges to communities reliant on these environments.

Alaska, Yukon↗
USGS2011

Mountain Glaciers and Ice Caps

In addition to the Greenland Ice Sheet, the Arctic contains a diverse array of smaller glaciers ranging from small cirque glaciers to large ice caps with areas up to 20 000 km 2 . Together, these glaciers cover an area of more than 400 000 km 2 , over half the global area of mountain glaciers and ice caps. Their total volume is sufficient to raise global sea level by an average of about 0.41 m if they were to melt completely. These glaciers exist in a range of different climatic regimes, from the maritime environments of southern Alaska, Iceland, western Scandinavia, and Svalbard, to the polar desert of the Canadian Arctic. Glaciers in all regions of the Arctic have decreased in area and mass as a result of the warming that has occurred since the 1920s (in two pulses – from the 1920s to the 1940s and since the mid-1980s). A new phase of accelerated mass loss began in the mid-1990s, and has been most marked in Alaska, the Canadian Arctic, and probably Greenland. Current rates of mass loss are estimated to be in the range 150 to 300 Gt/y; comparable to current mass loss rates from the Greenland Ice Sheet. This implies that the Arctic is now the largest regional source of glacier contributions to global sea-level rise. Most of the current mass loss is probably attributable to a change in surface mass balance (the balance between annual mass addition, primarily by snowfall, and annual mass loss by surface melting and meltwater runoff). Iceberg calving is also a significant source of mass loss in areas such as coastal Alaska, Arctic Canada, Svalbard, and the Russian Arctic. However, neither the current rate of calving loss nor its temporal variability have been well quantified in many regions, so this is a significant source of uncertainty in estimates of the total rate of mass loss. It is, however, clear that the larger Arctic ice caps have similar variability in ice dynamics to that of the Greenland Ice Sheet. That is to say, areas of relatively slow glacier flow (which terminate mainly on land) are separated by faster-flowing outlet glaciers (which terminate mainly in the ocean). Several of these outlet glaciers exhibit surge-type behavior, while others have exhibited substantial velocity changes on seasonal and longer timescales. It is very likely that these changes in ice dynamics affect the rate of mass loss by calving both from individual glaciers and the total ice cover. Projections of future rates of mass loss from mountain glaciers and ice caps in the Arctic focus primarily on projections of changes in the surface mass balance. Current models are not yet capable of making realistic forecasts of changes in losses by calving. Surface mass balance models are forced with downscaled output from climate models driven by forcing scenarios that make assumptions about the future rate of growth of atmospheric greenhouse gas concentrations. Thus, mass loss projections vary considerably, depending on the forcing scenario used and the climate model from which climate projections are derived. A new study in which a surface mass balance model is driven by output from ten general circulation models (GCMs) forced by the IPCC (Intergovernmental Panel on Climate Change) A1B emissions scenario yields estimates of total mass loss of between 51 and 136 mm sea-level equivalent (SLE) (or 13% to 36% of current glacier volume) by 2100. This implies that there will still be substantial glacier mass in the Arctic in 2100 and that Arctic mountain glaciers and ice caps will continue to influence global sea-level change well into the 22nd century.

Book chapter↗
USGS2017

Groundwater-derived nutrient and trace element transport to a nearshore Kona coral ecosystem: Experimental mixing model results

Study region The groundwater influenced coastal waters along the arid Kona coast of the Big Island, Hawai’i. Study focus A salinity-and phase partitioning-based mixing experiment was constructed using contrasting groundwater endmembers along the arid Konacoast of the Big Island, Hawai’i and local open seawater to better understand biogeochemical and physicochemical processes that influence the fate of submarine groundwater discharge (SGD)-derived nutrients and trace elements. New Hydrological Insights for the Region Treated wastewater effluent was the main source for nutrient enrichment downstream at the Honokōhau Harbor site. Conservative mixing for some constituents, such as nitrate + nitrite, illustrate the effectiveness of physical mixing to maintain oceanic concentrations in the colloid (0.02–0.45 μm) and truly dissolved (

Hawaii↗
USGS1989

Application of the Graphic Correlation method to Pliocene marine sequences

Biostratigraphy — the use of paleontological evidence to establish relative chronologies, forms the cornerstone of many sedimentary geological investigations. Several different approaches to biochronology are available. Traditional interval zones, defined on lowest and/or highest occurrences of selected taxa, are used to place bodies of rock in a relative chronological framework. Fossil datum levels, which are more numerous than zones, are often used as chronohorizons for correlation purposs. The Graphic Correlation method, like interval zonations, synthesizes information from a number of different taxa but does not assume synchrony of any one taxon. A magnetobiostratigraphic model for deep-sea Pliocene sequences has been constructed by graphic correlation of Deep Sea Drilling project cores from the North Atlantic (606), Caribbean Sea (502), South Atlantic (516), Tasman Sea (590), Equatorial Pacific (573) and North Pacific (577). All cores are hydraulic piston cores which contain abundant planktonic foraminifers, calcareous nannofossils and which record many of the magnetic reversals expected in the Pliocene. The model is based on internally consistent paleontologic data gathered by the author. This study demonstrates the advantages of graphic correlation over conventional biostratigraphic procedures. Accurate inter-regional correlations can be made between core sites without resorting to multiple microfossil zonations and without invoking synchrony of fossil events. Important results of this study are: (1) many Pliocene planktonic foraminifer and calcareous nannofossil events are diachronous by more than 0.20 m.y., (2) Globorotalia truncatulinoides first occurs in the Southwest Pacific Ocean, approximately 0.50 m.y. earlier than previously reported, (3) a previously undetected hiatus of short duration (0.38 m.y.) exists just above the Cochiti subchron at DSDP 577A.

Marine Micropaleontology↗
USGS2011

Biological and geochemical controls on diel dissolved inorganic carbon cycling in a low-order agricultural stream: Implications for reach scales and beyond

Movement of dissolved inorganic carbon (DIC) through the hydrologic cycle is an important component of global carbon budgets, but there is considerable uncertainty about the controls of DIC transmission from landscapes to streams, and through river networks to the oceans. In this study, diel measurements of DIC, d13C-DIC, dissolved oxygen (O2), d18O-O2, alkalinity, pH, and other parameters were used to assess the relative magnitudes of biological and geochemical controls on DIC cycling and flux in a nutrient-rich, net autotrophic stream. Rates of photosynthesis (P), respiration (R), groundwater discharge, air–water exchange of CO2, and carbonate precipitation/dissolution were quantified through a time-stepping chemical/isotope (12C and 13C, 16O and 18O) mass balance model. Groundwater was the major source of DIC to the stream. Primary production and carbonate precipitation were equally important sinks for DIC removed from the water column. The stream was always super-saturated with respect to carbonate minerals, but carbonate precipitation occurred mainly during the day when P increased pH. We estimated more than half (possibly 90%) of the carbonate precipitated during the day was retained in the reach under steady baseflow conditions. The amount of DIC removed from the overlying water through carbonate precipitation was similar to the amount of DIC generated from R. Air–water exchange of CO2 was always from the stream to the atmosphere, but was the smallest component of the DIC budget. Overall, the in-stream DIC reactions reduced the amount of CO2 evasion and the downstream flux of groundwater-derived DIC by about half relative to a hypothetical scenario with groundwater discharge only. Other streams with similar characteristics are widely distributed in the major river basins of North America. Data from USGS water quality monitoring networks from the 1960s to the 1990s indicated that 40% of 652 stream monitoring stations in the contiguous USA were at or above the equilibrium saturation state for calcite, and 77% of all stations exhibited apparent increases in saturation state from the 1960/70s to the 1980/90s. Diel processes including partially irreversible carbonate precipitation may affect net carbon fluxes from many such watersheds.

Chemical Geology↗
USGS2000

Estimating formation properties from early-time oscillatory water levels in a pumped well

Hydrologists often attempt to estimate formation properties from aquifer tests for which only the hydraulic responses in a pumped well are available. Borehole storage, turbulent head losses, and borehole skin, however, can mask the hydraulic behavior of the formation inferred from the water level in the pumped well. Also, in highly permeable formations or in formations at significant depth below land surface, where there is a long column of water in the well casing, oscillatory water levels may arise during the onset of pumping to further mask formation responses in the pumped well. Usually borehole phenomena are confined to the early stages of pumping or recovery, and late-time hydraulic data can be used to estimate formation properties. In many instances, however, early-time hydraulic data provide valuable information about the formation, especially if there are interferences in the late-time data. A mathematical model and its Laplace transform solution that account for inertial influences and turbulent head losses during pumping is developed for the coupled response between the pumped borehole and the formation. The formation is assumed to be homogeneous, isotropic, of infinite areal extent, and uniform thickness, with leakage from an overlying aquifer, and the screened or open interval of the pumped well is assumed to fully penetrate the pumped aquifer. Other mathematical models of aquifer flow can also be coupled with the equations describing turbulent head losses and the inertial effects on the water column in the pumped well. The mathematical model developed in this paper is sufficiently general to consider both underdamped conditions for which oscillations arise, and overdamped conditions for which there are no oscillations. Through numerical inversion of the Laplace transform solution, type curves from the mathematical model are developed to estimate formation properties through comparison with the measured hydraulic response in the pumped well. The mathematical model is applied to estimate formation properties from a singlewell test conducted near Waialua, Oahu, Hawaii. At this site, both the drawdown and recovery showed oscillatory water levels in the pumped well, and a step-drawdown test showed that approximately 86% of the drawdown is attributed to turbulent head losses. Analyses at this site using late-time drawdown data were confounded by the noise present in the measured water levels due primarily to nearby irrigation wells and ocean tides. By analyzing the early-time oscillatory recovery data at the Waialua site, upper and lower bounds were placed on the transmissivity, T, storage coefficient, S, and the leakance of the confining unit, K′/B′. The upper and lower bounds on T differ by a factor of 2. Upper and lower bounds on S and K′/B′ are much larger, because drawdown stabilized relatively quickly after the onset of pumping.

Journal of Hydrology↗
USGS2005

Potential for calibration of geostationary meteorological satellite imagers using the Moon

Solar-band imagery from geostationary meteorological satellites has been utilized in a number of important applications in Earth Science that require radiometric calibration. Because these satellite systems typically lack on-board calibrators, various techniques have been employed to establish "ground truth", including observations of stable ground sites and oceans, and cross-calibrating with coincident observations made by instruments with on-board calibration systems. The Moon appears regularly in the margins and corners of full-disk operational images of the Earth acquired by meteorological instruments with a rectangular field of regard, typically several times each month, which provides an excellent opportunity for radiometric calibration. The USGS RObotic Lunar Observatory (ROLO) project has developed the capability for on-orbit calibration using the Moon via a model for lunar spectral irradiance that accommodates the geometries of illumination and viewing by a spacecraft. The ROLO model has been used to determine on-orbit response characteristics for several NASA EOS instruments in low Earth orbit. Relative response trending with precision approaching 0.1% per year has been achieved for SeaWiFS as a result of the long time-series of lunar observations collected by that instrument. The method has a demonstrated capability for cross-calibration of different instruments that have viewed the Moon. The Moon appears skewed in high-resolution meteorological images, primarily due to satellite orbital motion during acquisition; however, the geometric correction for this is straightforward. By integrating the lunar disk image to an equivalent irradiance, and using knowledge of the sensor's spectral response, a calibration can be developed through comparison against the ROLO lunar model. The inherent stability of the lunar surface means that lunar calibration can be applied to observations made at any time, including retroactively. Archived geostationary imager data that contains the Moon can be used to develop response histories for these instruments, regardless of their current operational status.

Conference Paper↗
USGS2009

The performance of nearshore dredge disposal at Ocean Beach, San Francisco, California, 2005-2007

Ocean Beach, California, contains an erosion hot spot in the shadow of the San Francisco ebb tidal delta that threatens valuable public infrastructure as well as the safe recreational use of the beach. In an effort to reduce the erosion at this location a new plan for the management of sediment dredged annually from the main shipping channel at the mouth of San Francisco Bay was implemented in May 2005 by the United States Army Corps of Engineers, San Francisco District (USACE). The USACE designated a temporary nearshore dredge disposal site for the annual disposal of about 230,000 m3 (300,000 yd3) of sand about 750 m offshore and slightly south of the erosion hot spot, in depths between approximately 9 and 14 m. The site has now been used three times for a total sediment disposal of about 690,000 m3 (about 900,000 yds3). The disposal site was chosen because it is in a location where strong tidal currents and open-ocean waves can potentially feed sediment toward the littoral zone in the reach of the beach that is experiencing critical erosion, as well as prevent further scour on an exposed outfall pipe. The onshore migration of sediment from the target disposal location might feed the primary longshore bar or the nearshore zone, and provide a buffer to erosion that peaks during winter months when large waves impact the region. The United States Geological Survey (USGS) has been monitoring and modeling the bathymetric evolution of the test dredge disposal site and the adjacent coastal region since inception in May 2005. This paper reports on the first 2.5 years of this monitoring program effort (May 2005 to December 2007) and assesses the short-term coastal response. Here are the key findings of this report: *Approximately half of the sediment that has been placed in the nearshore dredge-disposal site during the 2.5 years of this study remains within the dredge focus area. *In the winter of 2006-7, large waves transported the dredge-mound material onshore. *High rates of seasonal cross-shore sediment transport mask any potential profile change in the Coastal Profiling System data due to dredge placement. *Pockets of accretion have been recorded by topographic surveying adjacent to the dredge site, but it is unclear if the accretion is linked to the nourishment. *Cross-shore profile modeling suggests that dredge material must be placed in water depths no greater than 5 m to drive a positive shoreline response. *Area modeling demonstrates that the new dredge site increases wave dissipation and modifies local sediment-transport patterns, although the effect on the nearshore morphology is largely negligible. *Any increase in beach width or wave energy-dissipation related to the nourishment is likely to be realized only in the vicinity directly onshore of the nourishment site, which is several hundred meters south of the area of critical erosion. *Larger waves from the northwest and smaller waves from the west or southwest contribute most to the sediment transport from the dredge mound onshore.

California↗
USGS1987

Geohydrology and simulated effects of withdrawals on the Miocene aquifer system in the Mississippi Gulf Coast area

Intense development of the Miocene aquifer system for water supplies along the Mississippi Gulf Coast has resulted in large water level declines that have altered the groundwater flow pattern in the area. Water levels in some Miocene aquifers have declined about 2 ft/year since 1940; declines exceed 100 ft (80 ft sea level) in large areas along the coast. Water levels in the surficial aquifer system, generally less than 20 ft below land surface, have not declined. The Miocene and younger interbedded and lenticular sands and clays crop out in southern Mississippi and dip to the south and southwest. These sediments have large vertical variations in head and locally respond to stresses as separate aquifers. Freshwater recharge to the Miocene aquifer system primarily is from rainfall on the surficial aquifers. The water generally moves to the south and southeast along the bedding planes toward the Mississippi Gulf Coast where the water is either withdrawn by wells, discharges to the ocean, or gradually percolates upward into overlying aquifers. Drawdowns caused by large groundwater withdrawals along the coast probably have resulted in the gradual movement of the saltwater toward the pumping centers. In parts of the Miocene aquifer system commonly used for water supplies, the water generally is a sodium bicarbonate type. Increasing chloride concentrations in a few wells indicate that saline water is migrating into parts of all layers in the Pascagoula area. A quasi three-dimensional numerical model of the groundwater flow system was constructed and calibrated on the basis of the both pre- and post-development conditions. The effects of an expected 1.5% annual increase in groundwater withdrawals during the period 1985-2005 were evaluated by the flow model. Additional water level declines expected by the year 2005 in response to estimated pumpage are as follows: Gulfport, 135 ft in layer 4; Biloxi-Gulfport area, 100 ft in layer 5 and 50 ft in layer 3; Pascagoula area, 40 ft in layer 6 and 30 ft in layer 4. The most serious threats of saltwater encroachment occur in layers 4, 5, and 6 (the 800-, 600- and 400-ft sands) in the Pascagoula area where contamination of the southern edges of the production areas is expected to occur in less than 10 years. (Author 's abstract)

Water-Resources Investigations Report↗
USGS2017

Developing multi-model ensemble projections of ecologically relevant climate variables for Puerto Rico and the US Caribbean

The global increases in surface air temperature are the most widespread and direct consequence of anthropogenic climate change. However, while 21st century temperatures are projected to increase in the Caribbean, the low variability and high average temperatures suggest that impacts on ecosystems and water resources are more likely through changes to the availability, timing, and pattern of moisture. The lack of local-scale climate model information that can resolve the complex topography and small scale climate features hinders the development of robust adaptation strategies. The goal of this project was to develop a suite of local-scale climate projections using dynamic downscaling to aid the development of adaptation strategies in Puerto Rico and the U.S. Virgin Islands (USVI). This project began by engaging the ecologists, hydrologists, and conservation biologists in the region to determine the most valuable types of information to aid research and decision making. The final product provides projections of future climate at a 2km horizontal resolution based on three global climate models and two regional climate models for a scenario with high greenhouse gas emissions. Results from the projections suggest that for Puerto Rico, annual temperature would increase between 1°C and 1.3°C by mid-century with larger temperature increases located in the interior portion of the island. Precipitation totals decrease for much of the island with island average decline between 12% and 19%, with some potentially large localized decreases exceeding 30%. The projected changes for the USVI are dominated by the surrounding ocean environment. The resulting projections will be provided to stakeholders in the region via the USGS and the CLCC.

Puerto Rico↗
USGS2009

Coastal ocean transport patterns in the central Southern California Bight

In the past decade, several large programs that monitor currents and transport patterns for periods from a few months to a few years were conducted by a consortium of university, federal, state, and municipal agencies in the central Southern California Bight, a heavily urbanized section of the coastal ocean off the west coast of the United States encompassing Santa Monica Bay, San Pedro Bay, and the Palos Verdes shelf. These programs were designed in part to determine how alongshelf and cross-shelf currents move sediments, pollutants, and suspended material through the region. Analysis of the data sets showed that the current patterns in this portion of the Bight have distinct changes in frequency and amplitude with location, in part because the topography of the shelf and upper slope varies rapidly over small spatial scales. However, because the mean, subtidal, and tidal-current patterns in any particular location were reasonably stable with time, one could determine a regional pattern for these current fields in the central Southern California Bight even though measurements at the various locations were obtained at different times. In particular, because the mean near-surface flows over the San Pedro and Palos Verdes shelves are divergent, near-surface waters from the upper slope tend to carry suspended material onto the shelf in the northwestern portion of San Pedro Bay. Water and suspended material are also carried off the shelf by the mean and subtidal flow fields in places where the orientation of the shelf break changes abruptly. The barotropic tidal currents in the central Southern California Bight flow primarily alongshore, but they have pronounced amplitude variations over relatively small changes in alongshelf location that are not totally predicted by numerical tidal models. Nonlinear internal tides and internal bores at tidal frequencies are oriented more across the shelf. They do not have a uniform transport direction, since they move fine sediment from the shelf to the slope in Santa Monica Bay, but carry suspended material from the mid-shelf to the beach in San Pedro Bay. It is clear that there are a large variety of processes that transport sediments and contaminants along and across the shelf in the central Southern California Bight. However, because these processes have a variety of frequencies and relatively small spatial scales, the dominant transport processes tend to be localized and have dissimilar characteristics even in adjacent regions of this small part of the coastal ocean. ?? 2009 The Geological Society of America.

Special Paper of the Geological Society of America↗
USGS2000

Simulated influences of Lake Agassiz on the climate of central North America 11,000 years ago

Eleven thousand years ago, large lakes existed in central and eastern North America along the margin of the Laurentide Ice Sheet. The large-scale North American climate at this time has been simulated with atmospheric general circulation models, but these relatively coarse global models do not resolve potentially important features of the mesoscale circulation that arise from interactions among the atmosphere, ice sheet, and proglacial lakes. Here we present simulations of the climate of central and eastern North America 11,000 years ago with a high-resolution, regional climate model nested within a general circulation model. The simulated climate is in general agreement with that inferred from palaeoecological evidence. Our experiments indicate that through mesoscale atmospheric feedbacks, the annual delivery of moisture to the Laurentide Ice Sheet was diminished at times of a large, cold Lake Agassiz relative to periods of lower lake stands. The resulting changes in the mass balance of the ice sheet may have contributed to fluctuations of the ice margin, thus affecting the routing of fresh water to the North Atlantic Ocean. A retreating ice margin during periods of high lake level may have opened an outlet for discharge of Lake Agassiz into the North Atlantic. A subsequent advance of the ice margin due to greater moisture delivery associated with a low lake level could have dammed the outlet, thereby reducing discharge to the North Atlantic. These variations may have been decisive in causing the Younger Dryas cold even.

Nature↗
USGS2022

One shell of a problem: Cumulative threat analysis of male sea turtles indicates high anthropogenic threat for migratory individuals and Gulf of Mexico residents

Human use of oceans has dramatically increased in the 21st century. Sea turtles are vulnerable to anthropogenic stressors in the marine environment because of lengthy migrations between foraging and breeding sites, often along coastal migration corridors. Little is known about how movement and threat interact specifically for male sea turtles. To better understand male sea turtle movement and the threats they encounter, we satellite-tagged 40 adult male sea turtles of four different species. We calculated movement patterns using state-space modeling (SSM), and quantified threats in seven unique categories; shipping, fishing, light pollution, oil rigs, proximity to coast, marine protected area (MPA) status, and location within or outside of the U.S. Exclusive Economic Zone (EEZ). We found significantly higher threat severity in northern and southern latitudes for green turtles ( Chelonia mydas ) and Kemp’s ridleys ( Lepidochelys kempii ) in our study area. Those threats were pervasive, with only 35.9% of SSM points encountering no high threat exposure, of which 47% belong to just two individuals. Kemp’s ridleys were most exposed to high threats among tested species. Lastly, turtles within MPA boundaries face significantly lower threat exposure, indicating MPAs could be a useful conservation tool.

Remote Sensing↗
USGS1989

Geologic framework of the offshore region adjacent to Delaware

Several multichannel, common depth point (CDP) seismic reflection profiles concentrated in the area of the entrance to Delaware Bay provide a tie between the known onshore geology of the Coastal Plain of Delaware and the offshore geology of the Baltimore Canyon Trough. The data provide a basis for understanding the geologic framework and petroleum resource potential of the area immediately offshore Delaware. Our research has focused on buried early Mesozoic rift basins and their geologic history. Assuming that the buried basins are analogous to the exposed Newark Supergroup basins of Late Triassic-Early Jurassic age, the most likely possibility for occurrence of hydrocarbon source beds in the area of the landward margin of the Baltimore Canyon Trough is presumed to be lacustrine, organic-rich shales probably present in the basins. Although buried basins mapped offshore Delaware are within reach of drilling, no holes have been drilled to date; therefore, direct knowledge of source, reservoir, and sealing beds is absent. Buried rift basins offshore Delaware show axial trends ranging from NW-SE to NNE-SSW. Seismic reflection profiles are too widely spaced to delineate basin boundaries accurately. Isopleths of two-way travel time representing basin fill suggest that, structurally, the basins are grabens and half-grabens. As shown on seismic reflection profiles, bounding faults of the basins intersect or merge with low-angle fault surfaces that cut the pre-Mesozoic basement. The rift basins appear to have formed by Mesozoic extension that resulted in reverse motion on reactivated basement thrust faults that originated from compressional tectonics during the Paleozoic. Computer-plotted structure contour maps derived from analysis of seismic reflection profiles provide information on the burial history of the rift basins. The postrift unconformity bevels the rift basins and, in the offshore area mapped, ranges from 2000 to 12,000 m below present sea level. The oldest postrift sediments that cover the more deeply buried rift basins are estimated to be of Middle Jurassic age (Bajocian-Bathonian), the probable time of opening of the Atlantic Ocean basin and onset of continental drift about 175–180 m.y. ago. By late Oxfordian-early Kimmeridgian time, the less deeply buried basins nearshore Delaware had been covered. A time-temperature index of maturity plot of one of the basins indicates that only dry gas would be present in reservoirs in synrift rocks buried by more than 6000 m of postrift sediments and in the oldest (Bathonian?-Callovian?) postrift rocks. Less deeply buried synrift rocks landward of the basin modeled might still be within the oil generation window.

Marine Geology↗