SEARCH · Geology Reports

Results for “Estuaries”

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

SourceResearch connectionFrom yearThrough year

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

At least 919 records · Page 51Linked to original sources

USGS1981

Simulation model of Skeletonema costatum population dynamics in northern San Francisco Bay, California

A pseudo-two-dimensional model is developed to simulate population dynamics of one dominant phytoplankton species ( Skeletonema costatum ) in northern San Francisco Bay. The model is formulated around a conceptualization of this estuary as two distinct but coupled subsystems—a deep (10–20 m) central channel and lateral areas with shallow (<2 m) water and slow circulation. Algal growth rates are governed by solar irradiation, temperature and salinity, while population losses are assumed to result from grazing bycalanoid copepods. Consequences of estuarine gravitational circulation are approximated simply by reducing convective-dispersive transport in that section of the channel (null zone) where residual bottom currents are near zero, and lateral mixing is treated as a bulkexchange process between the channel and the shoals. Model output is consistent with the hypothesis that, because planktonic algae are light-limited, shallow areas are the sites of active population growth. Seasonal variation in the location of the null zone (a response to variable river discharge) is responsible for maintaining the spring bloom of neritic diatoms in the seaward reaches of the estuary (San Pablo Bay) and the summer bloom upstream (Suisun Bay). Model output suggests that these spring and summer blooms result from the same general process—establishment of populations over the shoals, where growth rates are rapid, coupled with reduced particulate transport due to estuarine gravitational circulation. It also suggests, however, that the relative importance of physical and biological processes to phytoplankton dynamics is different in San Pablo and Suisun Bays. Finally, the model has helped us determine those processes having sufficient importance to merit further refinement in the next generation of models, and it has given new direction to field studies.

California↗

USGS1999

Sediment chronology in San Francisco Bay, California, defined by 210 Pb, 234 Th, 137 Cs, and 239,340 Pu

Sediment chronologies based on radioisotope depth profiles were developed at two sites in the San Francisco Bay estuary to provide a framework for interpreting historical trends in organic compound and metal contaminant inputs. At Richardson Bay near the estuary mouth, sediments are highly mixed by biological and/or physical processes. Excess penetration ranged from 2 to more than 10 cm at eight coring sites, yielding surface sediment mixing coefficients ranging from 12 to 170 cm 2 /year. At the site chosen for contaminant analyses, excess activity was essentially constant over the upper 25 cm of the core with an exponential decrease below to the supported activity between 70 and 90 cm. Both and penetrated to 57-cm depth and have broad subsurface maxima between 33 and 41 cm. The best fit of the excess profile to a steady state sediment accumulation and mixing model yielded an accumulation rate of 0.825 g/cm 2 /year (0.89 cm/year at sediment surface), surface mixing coefficient of 71 cm 2 /year, and 33-cm mixed zone with a half-Gaussian depth dependence parameter of 9 cm. Simulations of and profiles using these parameters successfully predicted the maximum depth of penetration and the depth of maximum and activity. Profiles of successive 1-year hypothetical contaminant pulses were generated using this parameter set to determine the age distribution of sediments at any depth horizon. Because of mixing, sediment particles with a wide range of deposition dates occur at each depth. A sediment chronology was derived from this age distribution to assign the minimum age of deposition and a date of maximum deposition to a depth horizon. The minimum age of sediments in a given horizon is used to estimate the date of first appearance of a contaminant from its maximum depth of penetration. The date of maximum deposition is used to estimate the peak year of input for a contaminant from the depth interval with the highest concentration of that contaminant. Because of the extensive mixing, sediment-bound constituents are rapidly diluted with older material after deposition. In addition, contaminants persist in the mixed zone for many years after deposition. More than 75 years are required to bury 90% of a deposited contaminant below the mixed zone. Reconstructing contaminant inputs is limited to changes occurring on a 20-year time scale. In contrast, mixing is much lower relative to accumulation at a site in San Pablo Bay. Instead, periods of rapid deposition and/or erosion occurred as indicated by frequent sand-silt laminae in the X-radiograph. , , and excess activity all penetrated to about 120 cm. The distinct maxima in the fallout radionuclides at 105–110 cm yielded overall linear sedimentation rates of 3.9 to 4.1 cm/year, which are comparable to a rate of 4.5±1.5 cm/year derived from the excess profile.

California↗

USGS2007

Sensitivity of estuarine turbidity maximum to settling velocity, tidal mixing, and sediment supply

Estuarine turbidity maximum, numerical modeling, settling velocity, stratification The spatial and temporal distribution of suspended material in an Estuarine Turbidity Maxima (ETM) is primarily controlled by particle settling velocity, tidal mixing, shear-stress thresholds for resuspension, and sediment supply. We vary these parameters in numerical experiments of an idealized two-dimensional (x-z) estuary to demonstrate their affects on the development and retention of particles in an ETM. Parameters varied are the settling velocity (0.01, 0.1, and 0.5 mm/s), tidal amplitude (0.4 m 12 hour tide and 0.3 to 0.6 m 14 day spring neap cycle), and sediment availability (spatial supply limited or unlimited; and temporal supply as a riverine pulse during spring vs. neap tide). Results identify that particles with a low settling velocity are advected out of the estuary and particles with a high settling velocity provide little material transport to an ETM. Particles with an intermediate settling velocity develop an ETM with the greatest amount of material retained. For an unlimited supply of sediment the ETM and limit of salt intrusion co-vary during the spring neap cycle. The ETM migrates landward of the salt intrusion during spring tides and seaward during neap tides. For limited sediment supply the ETM does not respond as an erodible pool of sediment that advects landward and seaward with the salt front. The ETM is maintained seaward of the salt intrusion and controlled by the locus of sediment convergence in the bed. For temporal variability of sediment supplied from a riverine pulse, the ETM traps more sediment if the pulse encounters the salt intrusion at neap tides than during spring tides. ?? 2007 Elsevier B.V. All rights reserved.

Proceedings in Marine Science↗

USGS2021

Nutrient improvements in Chesapeake Bay: Direct effect of load reductions and implications for coastal management

In Chesapeake Bay in the United States, decades of management efforts have resulted in modest reductions of nutrient loads from the watershed, but corresponding improvements in estuarine water quality have not clearly materialized. Generalized additive models were used to directly link river flows and nutrient loads from the watershed to nutrient trends in the estuary on a station-by-station basis, which allowed for identification of exactly when and where responses are happening. Results show that Chesapeake Bay total nitrogen and total phosphorus conditions are mostly improving after accounting for variation in freshwater flow. Almost all of these improving nutrient concentrations in the estuary can be explained by reductions in watershed loads entering through 16 rivers and 145 nearby point sources. These two major types of loads from multiple locations across the watershed are together necessary and responsible for improving estuarine nutrient conditions, a finding that is highly relevant to managing valuable estuarine resources worldwide.

Maryland, Virginia↗

USGS2017

Mercury bioaccumulation in estuarine fishes: Novel insights from sulfur stable isotopes

Estuaries are transitional habitats characterized by complex biogeochemical and ecological gradients that result in substantial variation in fish total mercury concentrations (THg). We leveraged these gradients and used carbon (δ 13 C), nitrogen (δ 15 N), and sulfur (δ 34 S) stable isotopes to examine the ecological and biogeochemical processes underlying THg bioaccumulation in fishes from the San Francisco Bay Estuary. We employed a tiered approach that first examined processes influencing variation in fish THg among wetlands, and subsequently examined the roles of habitat and within-wetland processes in generating larger-scale patterns in fish THg. We found that δ 34 S, an indicator of sulfate reduction and habitat specific-foraging, was correlated with fish THg at all three spatial scales. Over the observed ranges of δ 34 S, THg concentrations in fish increased by up to 860% within wetlands, 560% among wetlands, and 291% within specific impounded wetland habitats. In contrast, δ 13 C and δ 15 N were not correlated with THg among wetlands and were only important in low salinity impounded wetlands, possibly reflecting more diverse food webs in this habitat. Together, our results highlight the key roles of sulfur biogeochemistry and ecology in influencing estuarine fish THg, as well as the importance of fish ecology and habitat in modulating the relationships between biogeochemical processes and Hg bioaccumulation.

California↗

USGS2002

Contamination of estuarine water, biota, and sediment by halogenated organic compounds: A field study

Studies conducted in the vicinity of an industrial outfall in the Calcasieu River estuary, Louisiana, have shown that water, bottom and suspended sediment, and four different species of biota are contaminated with halogenated organic compounds (HOC) including haloarenes. A "salting-out" effect in the estuary moderately enhanced the partitioning tendency of the contaminants into biota and sediments. Contaminant concentrations in water, suspended sediments, and biota were found to be far below the values predicted on the basis of the assumption of phase equilibria with respect to concentrations in bottom sediment. Relative concentration factors of HOC between biota (catfish) and bottom sediment increased with increasing octanol/estuarine water partition coefficients (Kow*), maximizing at log Kow* of about 5, although these ratios were considerably less than equilibrium values. In contrast, contaminant concentrations in water, biota, and suspended sediments were much closer to equilibrium values. Bioconcentration factors of HOC determined on the basis of lipid content for four different biotic species correlated reasonably well with equilibrium triolein/water partition coefficients (Ktw).

Louisiana↗

USGS2001

Contaminant profiles of two species of shorebirds foraging together at two neighboring sites in south San Francisco Bay, California

The San Francisco Bay estuary isused by over one million shorebirds during springmigration and is home to several hundred thousandduring the winter. Most shorebird use occurs in thesouthern reach of the estuary (South Bay). Thereduced water circulation and discharge fromindustrial sources in the South Bay are responsiblefor the highest levels of some trace elements in theestuary. Wintering shorebirds have been found to havestrong site fidelity to areas as small as a fewkilometers in the South Bay, which may increase theirexposure to contaminants near local point sources. Inaddition, different shorebird species foraging at thesame site have been shown to have differentcontaminant burdens. Thus, our objectives were totest whether contaminant burdens differed by species,or whether contaminant burdens differed in shorebirdscollected at adjacent sites. We examined thecontaminant profiles of two species of shorebirds,long-billed dowitchers ( Limnodromus scolopaceus )and western sandpipers ( Calidris mauri ) thatforage together at two sites, Hayward and Newark,separated by 8 km in the South Bay. We usedmultivariate analysis of variance tests to compare thecomposition of 14 elemental analytes in their livertissues and estimated their molar ratios of Hg and Se. Composite samples were used for contaminant analysesbecause of the small body size of the shorebirds. Seven elemental analytes (Ag, Ba, Be, Cr, Ni, Pb, V)were below detection limits in a majority of thesamples so statistical analyses were precluded. Inthe measurable analytes (As, Cd, Cu, Hg, Mn, Se, Zn),we found no significant intra-site differences ofcontaminant profiles for the two species. We pooledthe samples to examine inter-site differences andfound significant differences in contaminant profilesbetween shorebirds at the neighboring sites ( P = 0.03). Shorebirds at Newark had higher ( P < 0.05) concentrations of As, Cd, and Se than those at Hayward. Dowitchers at Newark had concentrations of Hg and Se which were highly correlated ( P < 0.003) in a mean molar ratio of 1:19, similar tothat reported in other birds. In the larger dowitcherspecies, we also examined exposure to 20organochlorine compounds. Organic analyses showedthat the dowitchers had been exposed to DDE, PCBs,dieldrin and trans-nonachlor, but with no significantdifferences in concentrations between Hayward andNewark ( P > 0.05).

Environmental Monitoring and Assessment↗

USGS2010

Biological communities in San Francisco Bay track large‐scale climate forcing over the North Pacific

Long‐term observations show that fish and plankton populations in the ocean fluctuate in synchrony with large‐scale climate patterns, but similar evidence is lacking for estuaries because of shorter observational records. Marine fish and invertebrates have been sampled in San Francisco Bay since 1980 and exhibit large, unexplained population changes including record‐high abundances of common species after 1999. Our analysis shows that populations of demersal fish, crabs and shrimp covary with the Pacific Decadal Oscillation (PDO) and North Pacific Gyre Oscillation (NPGO), both of which reversed signs in 1999. A time series model forced by the atmospheric driver of NPGO accounts for two‐thirds of the variability in the first principal component of species abundances, and generalized linear models forced by PDO and NPGO account for most of the annual variability of individual species. We infer that synchronous shifts in climate patterns and community variability in San Francisco Bay are related to changes in oceanic wind forcing that modify coastal currents, upwelling intensity, surface temperature, and their influence on recruitment of marine species that utilize estuaries as nursery habitat. Ecological forecasts of estuarine responses to climate change must therefore consider how altered patterns of atmospheric forcing across ocean basins influence coastal oceanography as well as watershed hydrology.

California↗

USGS2011

Wave-current interaction in Willapa Bay

This paper describes the importance of wave-current interaction in an inlet-estuary system. The three-dimensional, fully coupled, Coupled Ocean-Atmosphere-Wave-Sediment Transport (COAWST) modeling system was applied in Willapa Bay (Washington State) from 22 to 29 October 1998 that included a large storm event. To represent the interaction between waves and currents, the vortex-force method was used. Model results were compared with water elevations, currents, and wave measurements obtained by the U.S. Army Corp of Engineers. In general, a good agreement between field data and computed results was achieved, although some discrepancies were also observed in regard to wave peak directions in the most upstream station. Several numerical experiments that considered different forcing terms were run in order to identify the effects of each wind, tide, and wave-current interaction process. Comparison of the horizontal momentum balances results identified that wave-breaking-induced acceleration is one of the leading terms in the inlet area. The enhancement of the apparent bed roughness caused by waves also affected the values and distribution of the bottom shear stress. The pressure gradient showed significant changes with respect to the pure tidal case. During storm conditions the momentum balance in the inlet shares the characteristics of tidal-dominated and wave-dominated surf zone environments. The changes in the momentum balance caused by waves were manifested both in water level and current variations. The most relevant effect on hydrodynamics was a wave-induced setup in the inner part of the estuary.

Washington↗

USGS2024

Testing megathrust rupture models using tsunami deposits

The 26 January 1700 CE Cascadia subduction zone earthquake ruptured much of the plate boundary and generated a tsunami that deposited sand in coastal marshes from northern California to Vancouver Island. Although the depositional record of tsunami inundation is extensive in some of these marshes, few sites have been investigated in enough detail to map the inland extent of sand deposition and depict variability in tsunami deposit thickness and grain size. We collected 129 cores in marshes of the Salmon River estuary in Oregon and reanalyzed 114 core logs from a 1987–88 study that mapped the inland extent of circa 1700 CE sandy tsunami deposits. The ca. 1700 CE tsunami deposit in the Salmon River estuary is easily recognized in cores ≤1 m deep in which a buried marsh peat is overlain by a well sorted sand bed with a sharp lower contact that thins and fines inland. We use tsunami deposit data and models of sandy tsunami sediment transport (using Delft3D-FLOW) to test 15 rupture models that could represent a ca. 1700 CE earthquake. At least 12–16 m of slip offshore of the Salmon River, which results in 0.8–1.0 m of coastal coseismic subsidence, is required to match the ca. 1700 CE sand deposit's inland extent, which is consistent with models of heterogeneous megathrust slip in ca. 1700 CE. Our methods of detailed tsunami deposit mapping, combined with sediment transport modeling, can be used to test models of megathrust ruptures and their tsunamis to potentially improve earthquake and tsunami hazard assessments.

JGR Earth Surface↗

USGS2023

Biophysical drivers of coastal treeline elevation

Sea level rise is leading to the rapid migration of marshes into coastal forests and other terrestrial ecosystems. Although complex biophysical interactions likely govern these ecosystem transitions, projections of sea level driven land conversion commonly rely on a simplified “threshold elevation” that represents the elevation of the marsh-upland boundary based on tidal datums alone. To determine the influence of biophysical drivers on threshold elevations, and their implication for land conversion, we examined almost 100,000 high-resolution marsh-forest boundary elevation points, determined independently from tidal datums, alongside hydrologic, ecologic, and geomorphic data in the Chesapeake Bay, the largest estuary in the U.S. located along the mid-Atlantic coast. We find five-fold variations in threshold elevation across the entire estuary, driven not only by tidal range, but also salinity and slope. However, more than half of the variability is unexplained by these variables, which we attribute largely to uncaptured local factors including groundwater discharge, microtopography, and anthropogenic impacts. In the Chesapeake Bay, observed threshold elevations deviate from predicted elevations used to determine sea level driven land conversion by as much as the amount of projected regional sea level rise by 2050. These results suggest that local drivers strongly mediate coastal ecosystem transitions, and that predictions based on elevation and tidal datums alone may misrepresent future land conversion.

Delaware, Maryland, Virginia↗

USGS2026

Extreme Potomac floods at Washington D.C. during the past 500 years

Washington D.C. faces one of the highest 100-year flood risks of any major city along the U.S. East Coast. In addition to storm-surge inundation during hurricanes and nor'easters, water-level observations for Washington are strongly skewed by major floods on the Potomac River. Using geologic and historic records we find new evidence for ice-jam flooding at Georgetown during the Little Ice Age, as recently as 1784, that was up to ∼2x the magnitude of the largest events of the past hundred years (1936, 1942). Over the 19th century (a) human modifications to the Potomac estuary as well as (b) increasingly heavy rainfall and (c) land-clearance in the watershed may have contributed to increasingly frequent large floods at Washington. Early surveys of the U.S. Capitol Building and other local landmarks also suggest sea level on the Potomac estuary at Washington has risen by upwards of 0.7 m (2.2 ft) since the 1790s.

Maryland, Virginia↗

USGS2026

Spatial heterogeneity of salt marsh vulnerability to sea-level rise: Dual controls of hydrological setting and salinity regime

Salt marsh vulnerability to sea-level rise (SLR) is typically assessed using point measurements of vertical accretion, neglecting three-dimensionality of geomorphic evolution and spatial variability. Recent studies suggest links between vertical and horizontal vulnerability, with differences between oligohaline and polyhaline marshes, yet these relationships remain untested in estuary-marsh systems. Here we combine geospatial analysis with hydrodynamic modeling to evaluate how unvegetated/vegetated marsh ratio (UVVR), a metric of marsh degradation, relates to elevation across hydrological regions and salinity regimes in the Albemarle-Pamlico Estuarine System, the largest lagoonal estuary in U.S. We show that at given normalized elevation, UVVR decreases across hydrological regions and salinity regimes from offshore to inland. UVVR-elevation relationship varies systematically with both hydrological setting and salinity regime, with hydrology exerting stronger influence. These findings challenge the assumption of a universal marsh deterioration trajectory and underscore the need to account for spatial heterogeneity when predicting responses to SLR.

North Carolina↗

USGS2008

Nitrogen and phosphorus speciation and flux in a large Florida river wetland System

Hydrologic measurements and analyses of various nitrogen and phosphorus species were made on the Apalachicola River system in northern Florida in 1979 and 1980. Annual outflows of total nitrogen (TN) and total phosphorus (TP) were not substantially different from annual inflows. However, there was significant net import of ammonia and soluble reactive phosphorus and net export of some particulate and organic species. The TN: TP ratio ranged from 12 to 15, but the specific ratio of dissolved inorganic nitrogen: soluble reactive phosphorus was much higher (up to 40) and increased in a downstream direction; this contributed to a phosphorus-limiting situation in Apalachicola estuary. Processes within the flood plain ecosystem accounted for much of the release of organic and particulate species and retention of inorganic species. This flood plain function is probably critical for maintaining a nutrient pool in the estuary which supports secondary productivity and a detrital-based food web.

Florida↗

USGS2010

On tide-induced Lagrangian residual current and residual transport: 1. Lagrangian residual current

Residual currents in tidal estuaries and coastal embayments have been recognized as fundamental factors which affect the long-term transport processes. It has been pointed out by previous studies that it is more relevant to use a Lagrangian mean velocity than an Eulerian mean velocity to determine the movements of water masses. Under weakly nonlinear approximation, the parameter k , which is the ratio of the net displacement of a labeled water mass in one tidal cycle to the tidal excursion, is assumed to be small. Solutions for tides, tidal current, and residual current have been considered for two-dimensional, barotropic estuaries and coastal seas. Particular attention has been paid to the distinction between the Lagrangian and Eulerian residual currents. When k is small, the first-order Lagrangian residual is shown to be the sum of the Eulerian residual current and the Stokes drift. The Lagrangian residual drift velocity or the second-order Lagrangian residual current has been shown to be dependent on the phase of tidal current. The Lagrangian drift velocity is induced by nonlinear interactions between tides, tidal currents, and the first-order residual currents, and it takes the form of an ellipse on a hodograph plane. Several examples are given to further demonstrate the unique properties of the Lagrangian residual current.

Water Resources Research↗

USGS2010

On tide-induced Lagrangian residual current and residual transport: 2. Residual transport with application in south San Francisco Bay, California

The transports of solutes and other tracers are fundamental to estuarine processes. The apparent transport mechanisms are convection by tidal current and current-induced shear effect dispersion for processes which take place in a time period of the order of a tidal cycle. However, as emphasis is shifted toward the effects of intertidal processes, the net transport is mainly determined by tide-induced residual circulation and by residual circulation due to other processes. The commonly used intertidal conservation equation takes the form of a convection-dispersion equation in which the convective velocity is the Eulerian residual current, and the dispersion terms are often referred to as the phase effect dispersion or, sometimes, as the “tidal dispersion.” The presence of these dispersion terms is merely the result of a Fickian type hypothesis. Since the actual processes are not Fickian, thus a Fickian hypothesis obscures the physical significance of this equation. Recent research results on residual circulation have suggested that long-term transport phenomena are closely related to the Lagrangian residual current or the Lagrangian residual transport. In this paper a new formulation of an intertidal conservation equation is presented and examined in detail. In a weakly nonlinear tidal estuary the resultant intertidal transport equation also takes the form of a convection-dispersion equation without the ad hoc introduction of phase effect dispersion in a form of dispersion tensor. The convective velocity in the resultant equation is the first-order Lagrangian residual current (the sum of the Eulerian residual current and the Stokes drift). The remaining dispersion terms are important only in higher-order solutions; they are due to shear effect dispersion and turbulent mixing. There exists a dispersion boundary layer adjacent to shoreline boundaries. An order of magnitude estimate of the properties in the dispersion boundary layer is given. The present treatment of intertidal transport processes is illustrated by an analytical solution for an amphidromic system and by a numerical application in South San Francisco Bay, California. The present formulation reveals that the mechanism for long-term transport of solutes is mainly convection due to the Lagrangian residual current in the interior of a tidal estuary. This result also points out the weakness in the tidal dispersion formulation, and explains the large variability of the observed values for tidal dispersion coefficients. Further research on properties of the dispersion boundary layer is needed.

California↗

USGS2024

Top-predator recovery abates geomorphic decline of a coastal ecosystem

The recovery of top predators is thought to have cascading effects on vegetated ecosystems and their geomorphology 1,2 , but the evidence for this remains correlational and intensely debated 3,4 . Here we combine observational and experimental data to reveal that recolonization of sea otters in a US estuary generates a trophic cascade that facilitates coastal wetland plant biomass and suppresses the erosion of marsh edges—a process that otherwise leads to the severe loss of habitats and ecosystem services 5,6 . Monitoring of the Elkhorn Slough estuary over several decades suggested top-down control in the system, because the erosion of salt marsh edges has generally slowed with increasing sea otter abundance, despite the consistently increasing physical stress in the system (that is, nutrient loading, sea-level rise and tidal scour 7,8,9 ). Predator-exclusion experiments in five marsh creeks revealed that sea otters suppress the abundance of burrowing crabs, a top-down effect that cascades to both increase marsh edge strength and reduce marsh erosion. Multi-creek surveys comparing marsh creeks pre- and post-sea otter colonization confirmed the presence of an interaction between the keystone sea otter, burrowing crabs and marsh creeks, demonstrating the spatial generality of predator control of ecosystem edge processes: densities of burrowing crabs and edge erosion have declined markedly in creeks that have high levels of sea otter recolonization. These results show that trophic downgrading could be a strong but underappreciated contributor to the loss of coastal wetlands, and suggest that restoring top predators can help to re-establish geomorphic stability.

California↗

USGS2017

Water quality measurements in San Francisco Bay by the U.S. Geological Survey, 1969–2015

The U.S. Geological Survey (USGS) maintains a place-based research program in San Francisco Bay (USA) that began in 1969 and continues, providing one of the longest records of water-quality measurements in a North American estuary. Constituents include salinity, temperature, light extinction coefficient, and concentrations of chlorophyll- a , dissolved oxygen, suspended particulate matter, nitrate, nitrite, ammonium, silicate, and phosphate. We describe the sampling program, analytical methods, structure of the data record, and how to access all measurements made from 1969 through 2015. We provide a summary of how these data have been used by USGS and other researchers to deepen understanding of how estuaries are structured and function differently from the river and ocean ecosystems they bridge.

California↗