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

Sediment deposition, erosion, and bathymetric change in San Francisco Bay, California, 1971–1990 and 1999–2020

Bathymetric change analyses document historical patterns of sediment deposition and erosion, providing valuable insight into the sediment dynamics of coastal systems, including pathways of sediment and sediment-bound contaminants. In 2014 and 2015, the Office for Coastal Management, in partnership with the National Oceanic and Atmospheric Administration (NOAA) Office of Coastal Management, provided funding for new bathymetric surveys of large portions of San Francisco Bay. A total of 93 bathymetric surveys were conducted during this 2-year period, using a combination of interferometric sidescan and multibeam sonar systems. These data, along with recent NOAA, U.S. Geological Survey (USGS), U.S. Army Corps of Engineers, and private contractor surveys collected from 1999 to 2020 (hereinafter referred to as 2010s), were used to create the most comprehensive bathymetric digital elevation models (DEMs) of San Francisco Bay since the 1980s. Comparing DEMs created from these 2010s surveys with USGS DEMs created from NOAA’s 1971–1990 (hereinafter referred to as 1980s) surveys provides information on the quantities and patterns of erosion and deposition in San Francisco Bay during the 9 to 47 years between surveys. This analysis reveals that in the areas surveyed in both the 1980s and 2010s, the bay floor lost about 34 million cubic meters of sediment since the 1980s. Results from this study can be used to assess how San Francisco Bay has responded to changes in the system, such as sea-level rise and variation in sediment supply from the Sacramento-San Joaquin Delta and local tributaries, and supports the creation of a new, system-wide sediment budget. This report provides data on the quantities and patterns of sediment volume change in San Francisco Bay for ecosystem managers that are pertinent to various sediment-related issues, including restoration of tidal marshes, exposure of legacy contaminated sediment, and strategies for the beneficial use of dredged sediment.

California↗

Rheologic and structural controls on the deformation of Okmok volcano, Alaska: FEMs, InSAR, and ambient noise tomography

Interferometric synthetic aperture radar (InSAR) data indicate that the caldera of Okmok volcano, Alaska, subsided more than a meter during its eruption in 1997. The large deformation suggests a relatively shallow magma reservoir beneath Okmok. Seismic tomography using ambient ocean noise reveals two low‐velocity zones (LVZs). The shallow LVZ corresponds to a region of weak, fluid‐saturated materials within the caldera and extends from the caldera surface to a depth of 2 km. The deep LVZ clearly indicates the presence of the magma reservoir beneath Okmok that is significantly deeper (>4 km depth) compared to previous geodetic‐based estimates (3 km depth). The deep LVZ associated with the magma reservoir suggests magma remains in a molten state between eruptions. We construct finite element models (FEMs) to simulate deformation caused by mass extraction from a magma reservoir that is surrounded by a viscoelastic rind of country rock embedded in an elastic domain that is partitioned to account for the weak caldera materials observed with tomography. This configuration allows us to reduce the estimated magma reservoir depressurization to within lithostatic constraints, while simultaneously maintaining the magnitude of deformation required to predict the InSAR data. More precisely, the InSAR data are best predicted by an FEM simulating a rind viscosity of 7.5 × 10 16 Pa s and a mass flux of −4.2 × 10 9 kg/d from the magma reservoir. The shallow weak layer within the caldera provides a coeruption stress regime and neutral buoyancy horizon that support lateral magma propagation from the central magma reservoir to extrusion near the rim of the caldera.

Alaska↗

Many atolls may be uninhabitable within decades due to climate change

Observations show global sea level is rising due to climate change, with the highest rates in the tropical Pacific Ocean where many of the world’s low-lying atolls are located. Sea-level rise is particularly critical for low-lying carbonate reef-lined atoll islands; these islands have limited land and water available for human habitation, water and food sources, and ecosystems that are vulnerable to inundation from sea-level rise. Here we demonstrate that sea-level rise will result in larger waves and higher wave-driven water levels along atoll islands’ shorelines than at present. Numerical model results reveal waves will synergistically interact with sea-level rise, causing twice as much land forecast to be flooded for a given value of sea-level rise than currently predicted by current models that do not take wave-driven water levels into account. Atolls with islands close to the shallow reef crest are more likely to be subjected to greater wave-induced run-up and flooding due to sea-level rise than those with deeper reef crests farther from the islands’ shorelines. It appears that many atoll islands will be flooded annually, salinizing the limited freshwater resources and thus likely forcing inhabitants to abandon their islands in decades, not centuries, as previously thought.

Hawaii↗

Across borders: External factors and prior behaviour influence North Pacific albatross associations with fishing vessels

Understanding encounters between marine predators and fisheries across national borders and outside national jurisdictions offers new perspectives on unwanted interactions to inform ocean management and predator conservation. Although seabird–fisheries overlap has been documented at many scales, remote identification of vessel encounters has lagged because vessel movement data often are lacking. Here, we reveal albatrosses–fisheries associations throughout the North Pacific Ocean. We identified commercial fishing operations using Global Fishing Watch data and algorithms to detect fishing vessels. We compiled GPS tracks of adult black-footed Phoebastria nigripes and Laysan Phoebastria immutabilis albatrosses, and juvenile short-tailed albatrosses Phoebastria albatrus . We quantified albatrosses-vessel encounters based on the assumed distance that birds perceive a vessel (≤30 km), and associations when birds approached vessels (≤3 km). For each event we quantified bird behaviour, environmental conditions and vessel characteristics and then applied Boosted Regression Tree models to identify drivers and the duration of these associations. In regions of greater fishing effort short-tailed and Laysan albatrosses associated with fishing vessels more frequently. However, fishing method (e.g. longline, trawl) and flag nation did not influence association prevalence nor the duration short-tailed albatrosses attended fishing vessels. Laysan albatrosses were more likely to approach longer vessels. Black-footed albatrosses were the most likely to approach vessels (61.9%), but limited vessel encounters ( n = 21) prevented evaluation of meaningful explanatory models for this species of high bycatch concern. Temporal variables (time of day and month) and bird behavioural state helped explain when short-tailed albatrosses were in close proximity to a vessel, but environmental conditions were more important for explaining interaction duration. Laysan albatrosses were more likely to associate with vessels while searching and during the last 60% (by time) of their trips. Our results provide specific species–fisheries insight regarding contributing factors of high-risk associations that could lead to bycatch of albatrosses within national waters and on the high seas. Policy implications . Given the availability of Global Fishing Watch data, our analysis can be applied to other marine predators—if tracking data are available—to identify spatio-temporal patterns, vessel specific attributes and predator behaviours associated with fishing vessel associations, thus enabling predictive modelling and targeted mitigation measures.

Journal of Applied Ecology↗

Decadal-scale phenology and seasonal climate drivers of migratory baleen whales in a rapidly warming marine ecosystem

Species' response to rapid climate change can be measured through shifts in timing of recurring biological events, known as phenology. The Gulf of Maine is one of the most rapidly warming regions of the ocean, and thus an ideal system to study phenological and biological responses to climate change. A better understanding of climate-induced changes in phenology is needed to effectively and adaptively manage human-wildlife conflicts. Using data from a 20+ year marine mammal observation program, we tested the hypothesis that the phenology of large whale habitat use in Cape Cod Bay has changed and is related to regional-scale shifts in the thermal onset of spring. We used a multi-season occupancy model to measure phenological shifts and evaluate trends in the date of peak habitat use for North Atlantic right ( Eubalaena glacialis ), humpback ( Megaptera novaeangliae ), and fin ( Balaenoptera physalus ) whales. The date of peak habitat use shifted by +18.1 days (0.90 days/year) for right whales and +19.1 days (0.96 days/year) for humpback whales. We then evaluated interannual variability in peak habitat use relative to thermal spring transition dates (STD), and hypothesized that right whales, as planktivorous specialist feeders, would exhibit a stronger response to thermal phenology than fin and humpback whales, which are more generalist piscivorous feeders. There was a significant negative effect of western region STD on right whale habitat use, and a significant positive effect of eastern region STD on fin whale habitat use indicating differential responses to spatial seasonal conditions. Protections for threatened and endangered whales have been designed to align with expected phenology of habitat use. Our results show that whales are becoming mismatched with static seasonal management measures through shifts in their timing of habitat use, and they suggest that effective management strategies may need to alter protections as species adapt to climate change.

Massachusetts↗

Seasonal electrical resistivity surveys of a coastal bluff, Barter Island, North Slope Alaska

Select coastal regions of the North Slope of Alaska are experiencing high erosion rates that can be attributed in part to recent warming trends and associated increased storm intensity and frequency. The upper sediment column of the coastal North Slope of Alaska can be described as continuous permafrost underlying a thin (typically less than 1–2 m) active layer that responds variably to seasonal thaw cycles. Assessing the temporal and spatial variability of the active layer and underlying permafrost is essential to better constrain how heightened erosion may impact material fluxes to the atmosphere and the coastal ocean, and how enhanced thaw cycles may impact the stability of the coastal bluffs. In this study, multi-channel electrical resistivity tomography (ERT) was used to image shallow subsurface features of a coastal bluff west of Kaktovik, on Barter Island, northeast Alaska. A comparison of a suite of paired resistivity surveys conducted in early and late summer 2014 provided detailed information on how the active layer and permafrost are impacted during the short Arctic summer. Such results are useful in the development of coastal resilience models that tie together fluvial, terrestrial, climatic, geologic, and oceanographic forcings on shoreline stability.

Alaska↗

Probabilistic assessment of erosion and flooding risk in the northern Gulf of Mexico

We assess erosion and flooding risk in the northern Gulf of Mexico by identifying interdependencies among oceanographic drivers and probabilistically modeling the resulting potential for coastal change. Wave and water level observations are used to determine relationships between six hydrodynamic parameters that influence total water level and therefore erosion and flooding, through consideration of a wide range of univariate distribution functions and multivariate elliptical copulas. Using these relationships, we explore how different our interpretation of the present-day erosion/flooding risk could be if we had seen more or fewer extreme realizations of individual and combinations of parameters in the past by simulating 10,000 physically and statistically consistent sea-storm time series. We find that seasonal total water levels associated with the 100 year return period could be up to 3 m higher in summer and 0.6 m higher in winter relative to our best estimate based on the observational records. Impact hours of collision and overwash—where total water levels exceed the dune toe or dune crest elevations—could be on average 70% (collision) and 100% (overwash) larger than inferred from the observations. Our model accounts for non-stationarity in a straightforward, non-parametric way that can be applied (with little adjustments) to many other coastlines. The probabilistic model presented here, which accounts for observational uncertainty, can be applied to other coastlines where short record lengths limit the ability to identify the full range of possible wave and water level conditions that coastal mangers and planners must consider to develop sustainable management strategies.

Journal of Geophysical Research C: Oceans↗

Influence of welded boundaries in anelastic media on energy flow, and characteristics of P, S-I, and S-II waves: Observational evidence for inhomogeneous body waves in low-loss solids

A general computer code, developed to calculate anelastic reflection-refraction coefficients, energy flow, and the physical characteristics for general P , S -I, and S -II waves, quantitatively describes physical characteristics for wave fields in anelastic media that do not exist in elastic media. Consideration of wave fields incident on boundaries between anelastic media shows that scattered wave fields experience reductions in phase and energy speeds, increases in maximum attenuation and Q −1 , and directions of maximum energy flow distinct from phase propagation. Each of these changes in physical characteristics are shown to vary with angle of incidence. Finite relaxation times for anelastic media result in energy flow due to interaction of superimposed radiation fields and contribute to energy flow across anelastic boundaries for all angles of incidence. Agreement of theoretical and numerical results with laboratory measurements argues for the validity of the theoretical and numerical formulations incorporating inhomogeneous wave fields. The agreement attests to the applicability of the model and helps confirm the existence of inhomogeneous body waves and their associated set of distinct physical characteristics in the earth. The existence of such body waves in layered, low-loss anelastic solids implies the need to reformulate some seismological models of the earth. The exact anelastic formulation for a liquid-solid interface with no low-loss approximations predicts the existence of a range of angles of incidence or an anelastic Rayleigh window, through which significant amounts of energy are transmitted across the boundary. The window accounts for the discrepancy apparent between measured reflection data presented in early textbooks and predictions based on classical elasticity theory. Characteristics of the anelastic Rayleigh window are expected to be evident in certain sets of wide-angle, ocean-bottom reflection data and to be useful in estimating Q −1 for some ocean bottom reflectors.

Journal of Geophysical Research B: Solid Earth↗

Surface temperatures of the Mid-Pliocene North Atlantic Ocean: Implications for future climate

The Mid-Pliocene is the most recent interval in the Earth's history to have experienced warming of the magnitude predicted for the second half of the twenty-first century and is, therefore, a possible analogue for future climate conditions. With continents basically in their current positions and atmospheric CO 2 similar to early twenty-first century values, the cause of Mid-Pliocene warmth remains elusive. Understanding the behaviour of the North Atlantic Ocean during the Mid-Pliocene is integral to evaluating future climate scenarios owing to its role in deep water formation and its sensitivity to climate change. Under the framework of the Pliocene Research, Interpretation and Synoptic Mapping (PRISM) sea surface reconstruction, we synthesize Mid-Pliocene North Atlantic studies by PRISM members and others, describing each region of the North Atlantic in terms of palaeoceanography. We then relate Mid-Pliocene sea surface conditions to expectations of future warming. The results of the data and climate model comparisons suggest that the North Atlantic is more sensitive to climate change than is suggested by climate model simulations, raising the concern that estimates of future climate change are conservative.

Philosophical Transactions of the Royal Society A:↗

Biomarker chemistry and flux quantification methods for natural petroleum seeps and produced oils, offshore southern California

Sustained, natural oil seepage from the seafloor is common off southern California, and is of great interest to resource managers, who are tasked with distinguishing natural from anthropogenic oil sources. The major purpose of this study was to build upon the work previously funded by the Bureau of Ocean Energy Management (BOEM) and the U.S. Geological Survey (USGS) that has refined the oil-fingerprinting process to enable differentiation of the highly similar Monterey Formation oils from Outer Continental Shelf (OCS) production and adjacent natural seeps. In these initial studies, biomarker and stable carbon isotope ratios were used to infer the age, lithology, organic-matter input, and depositional environment of the source rocks for 388 samples of produced crude oil, seep oil, and tarballs mainly from coastal California. The analysis resulted in a predictive model of oil source families that could be applied to samples of unknown origin.

California↗

Design and performance of a horizontal mooring for upper-ocean research

This paper describes the design and performance of a two-dimensional moored array for sampling horizontal variability in the upper ocean. The mooring was deployed in Massachusetts Bay in a water depth of 84 m for the purpose of measuring the horizontal structure of internal waves. The mooring was instrumented with three acoustic current meters (ACMs) spaced along a 170-m horizontal cable that was stretched between two subsurface buoys 20 m below the sea surface. Five 25-m-long vertical instrument strings were suspended from the horizontal cable. A bottom-mounted acoustic Doppler current profiler (ADCP) was deployed nearby to measure the current velocity throughout the water column. Pressure sensors mounted on the subsurface buoys and the vertical instrument strings were used to measure the vertical displacements of the array in response to the currents. Measurements from the ACMs and the ADCP were used to construct time-dependent, two-dimensional current fields. The current fields were used as input to a numerical model that calculated the deformation of the array with respect to the nominal zero-current configuration. Comparison of the calculated vertical offsets of the downstream subsurface buoy and downstream vertical instrument string with the pressure measurements were used to verify the numerical code. These results were then used to estimate total deformation of the array due to the passage of the internal waves. Based on the analysis of the three internal wave events with the highest measured vertical offsets, it is concluded that the geometry of the main structure (horizontal cable and anchor legs) was kept to within ±2.0 m, and the geometry of the vertical instrument strings was kept to within ±4.0 m except for one instance when the current velocity reached 0.88 m s −1 .

Massachusetts↗

The oligocene Lund Tuff, Great Basin, USA: A very large volume monotonous intermediate

Unusual monotonous intermediate ignimbrites consist of phenocryst-rich dacite that occurs as very large volume (> 1000 km3) deposits that lack systematic compositional zonation, comagmatic rhyolite precursors, and underlying plinian beds. They are distinct from countless, usually smaller volume, zoned rhyolite-dacite-andesite deposits that are conventionally believed to have erupted from magma chambers in which thermal and compositional gradients were established because of sidewall crystallization and associated convective fractionation. Despite their great volume, or because of it, monotonous intermediates have received little attention. Documentation of the stratigraphy, composition, and geologic setting of the Lund Tuff - one of four monotonous intermediate tuffs in the middle-Tertiary Great Basin ignimbrite province - provides insight into its unusual origin and, by implication, the origin of other similar monotonous intermediates. The Lund Tuff is a single cooling unit with normal magnetic polarity whose volume likely exceeded 3000 km3. It was emplaced 29.02 ?? 0.04 Ma in and around the coeval White Rock caldera which has an unextended north-south diameter of about 50 km. The tuff is monotonous in that its phenocryst assemblage is virtually uniform throughout the deposit: plagioclase > quartz ??? hornblende > biotite > Fe-Ti oxides ??? sanidine > titanite, zircon, and apatite. However, ratios of phenocrysts vary by as much as an order of magnitude in a manner consistent with progressive crystallization in the pre-eruption chamber. A significant range in whole-rock chemical composition (e.g., 63-71 wt% SiO2) is poorly correlated with phenocryst abundance. These compositional attributes cannot have been caused wholly by winnowing of glass from phenocrysts during eruption, as has been suggested for the monotonous intermediate Fish Canyon Tuff. Pumice fragments are also crystal-rich, and chemically and mineralogically indistinguishable from bulk tuff. We postulate that convective mixing in a sill-like magma chamber precluded development of a zoned chamber with a rhyolitic top or of a zoned pyroclastic deposit. Chemical variations in the Lund Tuff are consistent with equilibrium crystallization of a parental dacitic magma followed by eruptive mixing of compositionally diverse crystals and high-silica rhyolite vitroclasts during evacuation and emplacement. This model contrasts with the more systematic withdrawal from a bottle-shaped chamber in which sidewall crystallization creates a marked vertical compositional gradient and a substantial volume of capping-evolved rhyolite magma. Eruption at exceptionally high discharge rates precluded development of an underlying plinian deposit. The generation of the monotonous intermediate Lund magma and others like it in the middle Tertiary of the western USA reflects an unusually high flux of mantle-derived mafic magma into unusually thick and warm crust above a subducting slab of oceanic lithosphere. ?? 2002 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

Channelized lava flows at the East Pacific Rise crest 9°-10°N: the importance of off-axis lava transport in developing the architecture of young oceanic crust

Submarine lava flows are the building blocks of young oceanic crust. Lava erupted at the ridge axis is transported across the ridge crest in a manner dictated by the rheology of the lava, the characteristics of the eruption, and the topography it encounters. The resulting lava flows can vary dramatically in form and consequently in their impact on the physical characteristics of the seafloor and the architecture of the upper 50–500 m of the oceanic crust. We have mapped and measured numerous submarine channelized lava flows at the East Pacific Rise (EPR) crest 9°–10°N that reflect the high-effusion-rate and high-flow-velocity end-member of lava eruption and transport at mid-ocean ridges. Channel systems composed of identifiable segments 50–1000 m in length extend up to 3 km from the axial summit trough (AST) and have widths of 10–50 m and depths of 2–3 m. Samples collected within the channels are N-MORB with Mg# indicating eruption from the AST. We produce detailed maps of lava surface morphology across the channel surface from mosaics of digital images that show lineated or flat sheets at the channel center bounded by brecciated lava at the channel margins. Modeled velocity profiles across the channel surface allow us to determine flux through the channels from 0.4 to 4.7 × 10 3 m 3 /s, and modeled shear rates help explain the surface morphology variation. We suggest that channelized lava flows are a primary mechanism by which lava accumulates in the off-axis region (1–3 km) and produces the layer 2A thickening that is observed at fast and superfast spreading ridges. In addition, the rapid, high-volume-flux eruptions necessary to produce channelized flows may act as an indicator of the local magma budget along the EPR. We find that high concentrations of channelized lava flows correlate with local, across-axis ridge morphology indicative of an elevated magma budget. Additionally, in locations where channelized flows are located dominantly to the east or west of the AST, the ridge crest is asymmetric, and layer 2A appears to thicken over a greater distance from the AST toward the side of the ridge crest where the channels are located.

Geochemistry, Geophysics, Geosystems↗

Widespread terrestrial ecosystem disruption at the onset of the Paleocene–Eocene Thermal Maximum

The Paleocene–Eocene Thermal Maximum (PETM, ~56 Mya) interval was marked by massive 13 C-depleted carbon emissions into the ocean/atmosphere system, manifested as a negative carbon isotope excursion (CIE) in sedimentary components, and ~5 °C global average warming. Episodes of hydrological perturbations and soil-erosion have been widely documented for the PETM but their link with vegetation- and carbon cycle changes remain poorly constrained. Here, we present organic microfossil evidence showing a strong increase in fern-dominated pioneer vegetation that replaced coniferous forests on the margin of the Norwegian Sea during the first millennia of the CIE. With the present stratigraphic constraints, the “fern spike” occurred simultaneously in terrestrial settings along the North Sea, Arctic Ocean, the US east coast and in southern Australia, indicating that pioneer vegetation persisted for several millennia following a partial collapse of previously stable terrestrial ecosystems. Both the ferns and influx of microcharcoal imply recurrent physical disturbance, including soil destabilization and erosion, potentially linked to droughts, wildfires, and strong hydrological forcing resulting from extreme climate change. Together with evidence for reworked clay minerals and ancient organic matter (kerogen), these findings show that highly disturbed terrestrial ecosystems were widespread across mid- and high-latitude regions globally. Carbon cycle model simulations suggest that a substantial loss of standing and buried biomass, along with oxidation of soil organic matter, acted as important positive feedbacks during the onset of the CIE. Additionally, enhanced kerogen weathering likely contributed as another major positive feedback throughout both the onset and main phase of the CIE.

Proceedings of the National Academy of Sciences↗

Volcanic aquifers of Hawai‘i—Hydrogeology, water budgets, and conceptual models

Hawai‘i’s aquifers have limited capacity to store fresh groundwater because each island is small and surrounded by saltwater. Saltwater also underlies much of the fresh groundwater. Fresh groundwater resources are, therefore, particularly vulnerable to human activity, short-term climate cycles, and long-term climate change. Availability of fresh groundwater for human use is constrained by the degree to which the impacts of withdrawal—such as lowering of the water table, saltwater intrusion, and reduction in the natural discharge to springs, streams, wetlands, and submarine seeps—are deemed acceptable. This report describes the hydrogeologic framework, groundwater budgets (inflows and outflows), conceptual models of groundwater occurrence and movement, and the factors limiting groundwater availability for the largest and most populated of the Hawaiian Islands—Kaua‘i, O‘ahu, Maui, and Hawai‘i Island. The bulk of each of Hawai‘i’s islands is built of many thin lava flows erupted from shield volcanoes; the great piles of lava flows form highly permeable aquifers. In some areas, low-permeability dikes cutting across the lava flows, or low-permeability ash and soil horizons interlayered with the lava flows, can substantially alter groundwater flow. On some islands, sedimentary rocks form thick semiconfining coastal-plain deposits, locally known as caprock, that impede natural groundwater discharge to the ocean. In some regions, thick lava flows that ponded in preexisting depressions form aquifers that are much less permeable than aquifers formed by thin lava flows. Fresh groundwater inflow to Hawai‘i’s aquifers comes from recharge. For predevelopment conditions (1870), estimates of groundwater recharge from this study are 871, 675, 1,279, and 5,291 million gallons per day (Mgal/d) for Kaua‘i, O‘ahu, Maui, and Hawai‘i Island, respectively. Estimates of recharge for recent conditions (2010 land cover and 1978–2007 rainfall for Kaua‘i, O‘ahu, and Maui; 2008 land cover and 1916–1983 rainfall for Hawai‘i Island) are 875, 660, 1,308, and 6,595 Mgal/d for Kaua‘i, O‘ahu, Maui, and Hawai‘i Island, respectively. Recent recharge values differ from predevelopment recharge values by only a few percent for all islands except Hawai‘i Island, where changes in forest cover affected recharge. Spatial distribution of recharge mimics the orographic rainfall pattern—recharge is high on windward slopes and mountain peaks below the top of the trade-wind inversion. Human activity such as irrigation also contributes to recharge in some areas. Outflows from Hawai‘i’s aquifers include withdrawals from wells and natural groundwater discharge to springs, streams, wetlands, and submarine seeps. Under predevelopment conditions, groundwater withdrawal is assumed to be negligible and natural groundwater discharge probably was equal, or close, to recharge. Under recent conditions (2000–2010), groundwater withdrawal averaged 19, 209, 104, and 103 Mgal/d on Kaua‘i, O‘ahu, Maui, and Hawai‘i Island, respectively. If recent withdrawal and recharge rates are maintained until steady state is achieved, natural groundwater discharge will be reduced by an amount equal to the withdrawal rate. Total recent withdrawal for the four islands is only about 5 percent of total recharge, but about half of the withdrawal comes from O‘ahu, whereas O‘ahu receives only 7 percent of the total recharge. Effects of high withdrawals on O‘ahu cannot be mitigated by the lower withdrawals on other islands because no freshwater flows between islands. Even within an island, high withdrawals from one area cannot be completely mitigated by recharge in another area. Water-level, saltwater/freshwater-transition-zone, spring, and stream base-flow data indicate an overall reduction in storage for most areas where groundwater has been developed. Groundwater occurrence and movement in Hawai‘i’s volcanic aquifers can be described in terms of four conceptual models: (1) fresh groundwater lenses in high-permeability lava-flow aquifers, (2) aquifers with groundwater impounded by dikes, (3) thickly saturated low-permeability aquifers, and (4) perched aquifers. In Hawai‘i, most fresh groundwater withdrawn for human use comes from freshwater lenses in the dike-free high-permeability lava-flow aquifers where the principal limiting factor to groundwater availability is saltwater intrusion, but impacts of reduced natural groundwater discharge may also limit availability. Dike-impounded groundwater is common near the center of Hawaiian shield volcanoes, where water moves and is stored in permeable lava flows between the dikes; groundwater availability in these aquifers is primarily limited by storage depletion and reduction of flow to adjacent aquifers and natural groundwater discharge. Thickly saturated low-permeability aquifers have been identified on Kaua‘i and Maui; groundwater availability is primarily limited by streamflow depletion and water-table decline. Perched groundwater is postulated to exist in some areas of Hawai‘i, but store much less water than other modes ofgroundwater occurrence. Limits on groundwater availability in perched aquifers include the potential of reducing inflow to other groundwater settings and reducing natural discharge and stream seepage. Some groundwater bodies in Hawai‘i are enigmatic; consequences of groundwater development in these bodies and their relation to groundwater availability are not completely understood.

Hawai'i↗

The airborne lava-seawater interaction plume at Kilauea Volcano, Hawai'i

Lava flows into the sea at Kīlauea Volcano, Hawaiʻi, and generates an airborne gas and aerosol plume. Water (H 2 O), hydrogen chloride (HCl), carbon dioxide (CO 2 ), nitrogen dioxide (NO 2 ) and sulphur dioxide (SO 2 ) gases were quantified in the plume in 2004–2005, using Open Path Fourier Transform infra-red Spectroscopy. The molar abundances of these species and thermodynamic modelling are used to discuss their generation. The range in molar HCl / H 2 O confirms that HCl is generated when seawater is boiled dry and magnesium salts are hydrolysed (as proposed by [T.M. Gerlach, J.L. Krumhansl, R.O. Fournier, J. Kjargaard, Acid rain from the heating and evaporation of seawater by molten lava: a new volcanic hazard, EOS (Trans. Am. Geophys. Un.) 70 (1989) 1421–1422]), in contrast to models of Na-metasomatism. Airborne droplets of boiled seawater brine form nucleii for subsequent H 2 O and HCl condensation, which acidifies the droplets and liberates CO 2 gas from bicarbonate and carbonate. NO 2 is derived from the thermal decomposition of nitrates in coastal seawater, which takes place as the lava heats droplets of boiled seawater brine to 350–400 °C. SO 2 is derived from the degassing of subaerial lava flows on the coastal plain. The calculated mass flux of HCl from a moderate-sized ocean entry significantly increases the total HCl emission at Kīlauea (including magmatic sources) and is comparable to industrial HCl emitters in the United States. For larger lava ocean entries, the flux of HCl will cause intense local environmental hazards, such as high localised HCl concentrations and acid rain.

Hawaii↗

Hydroclimatology of the Missouri River basin

Despite the importance of the Missouri River for navigation, recreation, habitat, hydroelectric power, and agriculture, relatively little is known about the basic hydroclimatology of the Missouri River basin (MRB). This is of particular concern given the droughts and floods that have occurred over the past several decades and the potential future exacerbation of these extremes by climate change. Here, observed and modeled hydroclimatic data and estimated natural flow records in the MRB are used to 1) assess the major source regions of MRB flow, 2) describe the climatic controls on streamflow in the upper and lower basins , and 3) investigate trends over the instrumental period. Analyses indicate that 72% of MRB runoff is generated by the headwaters in the upper basin and by the lowest portion of the basin near the mouth. Spring precipitation and temperature and winter precipitation impacted by changes in zonal versus meridional flow from the Pacific Ocean play key roles in surface water supply variability in the upper basin. Lower basin flow is significantly correlated with precipitation in late spring and early summer, indicative of Atlantic-influenced circulation variability affecting the flow of moisture from the Gulf of Mexico. Although increases in precipitation in the lower basin are currently overriding the effects of warming temperatures on total MRB flow, the upper basin’s long-term trend toward decreasing flows, reduction in snow versus rain fraction, and warming spring temperatures suggest that the upper basin may less often provide important flow supplements to the lower basin in the future.

Missouri River basin↗

Widespread kelp-derived carbon in pelagic and benthic nearshore fishes

Kelp forests provide habitat for diverse and abundant fish assemblages, but the extent to which kelp provides a source of energy to fish and other predators is unclear. To examine the use of kelp-derived energy by fishes we estimated the contribution of kelp- and phytoplankton-derived carbon using carbon (δ 13 C) and nitrogen (δ 15 N) isotopes measured in muscle tissue. Benthic-foraging kelp greenling (Hexagrammos decagrammus) and pelagic-foraging black rockfish (Sebastes melanops) were collected at eight sites spanning ∼35 to 60°N from the California Current (upwelling) to Alaska Coastal Current (downwelling) in the northeast Pacific Ocean. Muscle δ 13 C values were expected to be higher for fish tissue primarily derived from kelp, a benthic macroalgae, and lower for tissue primarily derived from phytoplankton, pelagic microalgae. Muscle δ 13 C values were higher in benthic-feeding kelp greenling than in pelagic-feeding black rockfish at seven of eight sites, indicating more kelp-derived carbon in greenling as expected. Estimates of kelp carbon contributions ranged from 36 to 89% in kelp greenling and 32 to 65% in black rockfish using carbon isotope mixing models. Isotopic evidence suggests that these two nearshore fishes routinely derive energy from kelp and phytoplankton, across coastal upwelling and downwelling systems. Thus, the foraging mode of nearshore predators has a small influence on their ultimate energy source as energy produced by benthic macroalgae and pelagic microalgae were incorporated in fish tissue regardless of feeding mode and suggest strong and widespread benthic-pelagic coupling. Widespread kelp contributions to benthic- and pelagic-feeding fishes suggests that kelp energy provides a benefit to nearshore fishes and highlights the potential for kelp and fish production to be linked.

Estuarine, Coastal and Shelf Science↗