USGS ScienceSearch

SEARCH · USGS Science

Results for “Journal of Geophysical Research - Oceans”

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 253 records · Page 14Linked to original sources

Source spectral properties of small-to-moderate earthquakes in southern Kansas

The source spectral properties of injection-induced earthquakes give insight into their nucleation, rupture processes, and influence on ground motion. Here we apply a spectral decomposition approach to analyze P -wave spectra and estimate Brune-type stress drop for more than 2000 M L 1.5–5.2 earthquakes occurring in southern Kansas from 2014 to 2016. We find that these earthquakes are characterized by low stress drop values (median ∼0.4MPa) compared to natural seismicity in California. We observe a significant increase in stress drop as a function of depth, but the shallow depth distribution of these events is not by itself sufficient to explain their lower stress drop. Stress drop increases with magnitude from M1.5–M3.5, but this scaling trend may weaken above M4 and also depends on the assumed source model. Although we observe a nonstationary, sequence-specific temporal evolution in stress drop, we find no clear systematic relation with the activity of nearby injection wells.

Kansas

Deformation of the Pacific/North America plate boundary at Queen Charlotte Fault: The possible role of rheology

The Pacific/North America (PA/NA) plate boundary between Vancouver Island and Alaska is similar to the PA/NA boundary in California in its kinematic history and the rate and azimuth of current relative motion, yet their deformation styles are distinct. The California plate boundary shows a broad zone of parallel strike slip and thrust faults and folds, whereas the 49‐mm/yr PA/NA relative plate motion in Canada and Alaska is centered on a single, narrow, continuous ~900‐km‐long fault, the Queen Charlotte Fault (QCF). Using gravity analysis, we propose that this plate boundary is centered on the continent/ocean boundary (COB), an unusual location for continental transform faults because plate boundaries typically localize within the continental lithosphere, which is weaker. Because the COB is a boundary between materials of contrasting elastic properties, once a fault is established there, it will probably remain stable. We propose that deformation progressively shifted to the COB in the wake of Yakutat terrane's northward motion along the margin. Minor convergence across the plate boundary is probably accommodated by fault reactivation on Pacific crust and by an eastward dipping QCF. Underthrusting of Pacific slab under Haida Gwaii occurs at convergence angles >14°–15° and may have been responsible for the emergence of the archipelago. The calculated slab entry dip (5°–8°) suggests that the slab probably does not extend into the asthenosphere. The PA/NA plate boundary at the QCF can serve as a structurally simple site to investigate the impact of rheology and composition on crustal deformation and the initiation of slab underthrusting.

Journal of Geophysical Research B: Solid Earth

Mass addition at Mount St. Helens, Washington, inferred from repeated gravity surveys

Measurements of subtle changes in the Earth’s gravityfield can provideinformation on the addition/loss of mass (e.g., magma or aqueousfluids) beneath a volcano. In this study, wemeasured gravity at Mount St. Helens from 2010 to 2016 to investigate possible mass changes followingthe 2004–2008 dome-forming eruption. The raw gravity measurements were corrected for changes in themass and shape of Crater Glacier and for deformation of the volcanic edifice to obtain residual gravity values.We found positive residual gravity changes that we interpret as evidence for partial recharge of the magmareservoir that fed the 2004–2008 eruption and/or for accumulation of groundwater in one or more shallowaquifers. Most of the gravity signal can be explained by groundwater accumulation; magma recharge is notrequired. However, there is additional evidence for at least partial recharge from seismic, deformation, andgeochemical observations

Washington

Constraints on friction, dilatancy, diffusivity, and effective stress from low-frequency earthquake rates on the deep San Andreas Fault

Families of recurring low-frequency earthquakes (LFEs) within nonvolcanic tremor on the San Andreas Fault in central California are sensitive to tidal stresses. LFEs occur at all levels of the tides, are strongly correlated and in phase with the ~200 Pa shear stresses, and weakly and not systematically correlated with the ~2 kPa tidal normal stresses. We assume that LFEs are small sources that repeatedly fail during shear within a much larger scale, aseismically slipping fault zone and consider two different models of the fault slip: (1) modulation of the fault slip rate by the tidal stresses or (2) episodic slip, triggered by the tides. LFEs are strongly clustered with duration much shorter than the semidiurnal tide; they cannot be significantly modulated on that time scale. The recurrence times of clusters, however, are many times longer than the semidiurnal, leading to an appearance of tidal triggering. In this context we examine the predictions of laboratory-observed triggered frictional (dilatant) fault slip. The undrained end-member model produces no sensitivity to the tidal normal stress, and slip onsets are in phase with the tidal shear stress. The tidal correlation constrains the diffusivity to be less than ~1 × 10 −6 /s and the product of the friction and dilatancy coefficients to be at most 5 × 10 −7 , orders of magnitude smaller than observed at room temperature. In the absence of dilatancy the effective normal stress at failure would be about ~55 kPa. For this model the observations require intrinsic weakness, low dilatancy, and lithostatic pore fluid.

Journal of Geophysical Research

Combining multiphase groundwater flow and slope stability models to assess stratovolcano flank collapse in the Cascade Range

Hydrothermal alteration can create low‐permeability zones, potentially resulting in elevated pore‐fluid pressures, within a volcanic edifice. Strength reduction by rock alteration and high pore‐fluid pressures have been suggested as a mechanism for edifice flank instability. Here we combine numerical models of multiphase heat transport and groundwater flow with a slope‐stability code that incorporates three‐dimensional distributions of strength and pore‐water pressure to address the following questions: (1) What permeability distributions and contrasts produce elevated pore‐fluid pressures in a stratovolcano? (2) What are the effects of these elevated pressures on flank stability? (3) Finally, what are the effects of magma intrusion on potential flank failure in an edifice? Simulation results show that under a range of plausible parameters, water tables in a stratovolcano can be elevated or perched. These elevated water tables result in universally lower stability (lower factor of safety) compared with equivalent dry edifices, indicating a higher likelihood of flank collapse. Low‐permeability (<1 × 10 −17 m 2 ) layers such as altered pyroclastic deposits or breccias can result in locally saturated regions (perched water) and lower factors of safety near the ground surface but may actually reduce liquid water saturation and pore pressures in the core of the edifice and thus may favor small, shallow collapses over larger, deeper collapses. Magma intrusion into the base of the edifice increases pore‐fluid pressures and decreases the factor of safety. However, the shear strength of edifice rocks also exerts a significant control on stability, so both mechanical properties and pore‐fluid pressures are important for stability assessments.

Washington

Seismicity induced by wastewater injection in Washington County, Ohio: Influence of preexisting structure, regional stress regime, and well operations

Recent seismicity in Washington County, Ohio, has been suggested to be induced by wastewater disposal operations despite injection ~2 km above the Precambrian basement. We investigated the relationships between disposal well locations and operational histories, spatiotemporal patterns of seismicity enhanced by waveform correlation, and mapped subsurface structures. We also analyzed proxies for in situ stress conditions to evaluate existing principal stress orientations and magnitudes. Double‐difference relocations of small ( M L ≤ 2.1) seismicity revealed linear clusters defining a NE‐SW orientation that parallels the trend of small‐amplitude folds in Upper Paleozoic rocks within 1–2 km of the injection well. Hypocenter depths of ~3.1 to 4.3 km indicate seismicity occurred along faults below the Silurian injection interval in lower Paleozoic strata and Precambrian crystalline basement. Similar folds/faults imaged in nearby seismic reflection lines were associated with basement‐involved fault systems that cut the injection interval targeted by the injection well, providing permeability pathways for fluid pressure increases that initiated slip. Interpreted faults were also optimally oriented in the regional stress field for reactivation due to injection‐related decreases in effective stress. Well completion records suggest that during the recorded seismicity, the Long Run‐1 well was injecting disposal fluids close to or above the fracture pressure of the targeted Silurian injection interval. A comparison of observed seismicity rates and monthly injection records from the disposal well also indicate a close correspondence between earthquake frequency and injection rate. Taken together, these results provide a physical framework for understanding why wastewater disposal operations induced the Washington County seismic sequence.

Ohio

Investigating the origin of continual radio frequency impulses during explosive volcanic eruptions

Volcanic lightning studies have revealed that there is a relatively long‐lasting source of very high frequency radiation associated with the onset of explosive volcanic eruptions that is distinct from radiation produced by lightning. This very high frequency signal is referred to as “continual radio frequency (CRF)” due to its long‐lasting nature. The discharge mechanism producing this signal was previously hypothesized to be caused by numerous, small (10–100 m) leader‐forming discharges near the vent of the volcano. To test this hypothesis, a multiparametric data set of electrical and volcanic activity occurring during explosive eruptions of Sakurajima Volcano in Japan was collected from May to June 2015. Our observations show that a single CRF impulse has a duration on the order of 160 ns (giving an upper limit on discharge length of 10 m) and is distinct from near‐vent lightning discharges that were on the order of 30 m in length. CRF impulses did not produce discernible electric field changes and occurred in the absence of a net static electric field. Lightning mapping data and infrared video observations of the eruption column showed that CRF impulses originated from the gas thrust region of the column. These observations indicate that CRF impulses are not produced by small, leader‐forming discharges but rather are more similar to a streamer discharge, likely on the order of a few meters in length.

Sakurajima Volcano

Presentation and analysis of a worldwide database of earthquake-induced landslide inventories

Earthquake-induced landslide (EQIL) inventories are essential tools to extend our knowledge of the relationship between earthquakes and the landslides they can trigger. Regrettably, such inventories are difficult to generate and therefore scarce, and the available ones differ in terms of their quality and level of completeness. Moreover, access to existing EQIL inventories is currently difficult because there is no centralized database. To address these issues, we compiled EQIL inventories from around the globe based on an extensive literature study. The database contains information on 363 landslide-triggering earthquakes and includes 66 digital landslide inventories. To make these data openly available, we created a repository to host the digital inventories that we have permission to redistribute through the U.S. Geological Survey ScienceBase platform. It can grow over time as more authors contribute their inventories. We analyze the distribution of EQIL events by time period and location, more specifically breaking down the distribution by continent, country, and mountain region. Additionally, we analyze frequency distributions of EQIL characteristics, such as the approximate area affected by landslides, total number of landslides, maximum distance from fault rupture zone, and distance from epicenter when the fault plane location is unknown. For the available digital EQIL inventories, we examine the underlying characteristics of landslide size, topographic slope, roughness, local relief, distance to streams, peak ground acceleration, peak ground velocity, and Modified Mercalli Intensity. Also, we present an evaluation system to help users assess the suitability of the available inventories for different types of EQIL studies and model development.

Journal of Geophysical Research F: Earth Surface

Compositional signatures in acoustic backscatter over vegetated and unvegetated mixed sand-gravel riverbeds

Multibeam acoustic backscatter has considerable utility for remote characterization of spatially heterogeneous bed sediment composition over vegetated and unvegetated riverbeds of mixed sand and gravel. However, the use of high-frequency, decimeter-resolution acoustic backscatter for sediment classification in shallow water is hampered by significant topographic contamination of the signal. In mixed sand-gravel riverbeds, changes in the abiotic composition of sediment (such as homogeneous sand to homogeneous gravel) tend to occur over larger spatial scales than is characteristic of small-scale bedform topography (ripples, dunes, and bars) or biota (such as vascular plants and periphyton). A two-stage method is proposed to filter out the morphological contributions to acoustic backscatter. First, the residual supragrain-scale topographic effects in acoustic backscatter with small instantaneous insonified areas, caused by ambiguity in the local (beam-to-beam) bed-sonar geometry, are removed. Then, coherent scales between high-resolution topography and backscatter are identified using cospectra, which are used to design a frequency domain filter that decomposes backscatter into the (unwanted) high-pass component associated with bedform topography (ripples, dunes, and sand waves) and vegetation, and the (desired) low-frequency component associated with the composition of sediment patches superimposed on the topography. This process strengthens relationships between backscatter and sediment composition. A probabilistic framework is presented for classifying vegetated and unvegetated substrates based on acoustic backscatter at decimeter resolution. This capability is demonstrated using data collected from diverse settings within a 386 km reach of a canyon river whose bed varies among sand, gravel, cobbles, boulders, and submerged vegetation.

Colorado River, Grand Canyon, Glen Canyon, Marble

Can beaches survive climate change?

Anthropogenic climate change is driving sea level rise, leading to numerous impacts on the coastal zone, such as increased coastal flooding, beach erosion, cliff failure, saltwater intrusion in aquifers, and groundwater inundation. Many beaches around the world are currently experiencing chronic erosion as a result of gradual, present-day rates of sea level rise (about 3 mm/year) and human-driven restrictions in sand supply (e.g., harbor dredging and river damming). Accelerated sea level rise threatens to worsen coastal erosion and challenge the very existence of natural beaches throughout the world. Understanding and predicting the rates of sea level rise and coastal erosion depends on integrating data on natural systems with computer simulations. Although many computer modeling approaches are available to simulate shoreline change, few are capable of making reliable long-term predictions needed for full adaption or to enhance resilience. Recent advancements have allowed convincing decadal to centennial-scale predictions of shoreline evolution. For example, along 500 km of the Southern California coast, a new model featuring data assimilation predicts that up to 67% of beaches may completely erode by 2100 without large-scale human interventions. In spite of recent advancements, coastal evolution models must continue to improve in their theoretical framework, quantification of accuracy and uncertainty, computational efficiency, predictive capability, and integration with observed data, in order to meet the scientific and engineering challenges produced by a changing climate.

Journal of Geophysical Research F: Earth Surface

A framework for modeling connections between hydraulics, water surface roughness, and surface reflectance in open channel flows

This paper introduces a framework for examining connections between the flow field, the texture of the air-water interface, and the reflectance of the water surface and thus evaluating the potential to infer hydraulic information from remotely sensed observations of surface reflectance. We used a spatial correlation model describing water surface topography to illustrate the application of our framework. Nondimensional relations between model parameters and flow intensity were established based on a prior flume study. Expressing the model in the spatial frequency domain allowed us to use an efficient Fourier transform-based algorithm for simulating water surfaces. Realizations for both flume and field settings had water surface slope distributions positively correlated with velocity and water surface roughness. However, most surface facets were gently sloped and thus unlikely to yield strong specular reflections; the model exaggerated the extent of water surface features, leading to underestimation of facet slopes. A ray tracing algorithm indicated that reflectance was greatest when solar and view zenith angles were equal and the sensor scanned toward the Sun to capture specular reflections of the solar beam. Reflected energy was concentrated in a small portion of the sky, but rougher water surfaces reflected rays into a broader range of directions. Our framework facilitates flight planning to avoid surface-reflected radiance while mapping other river attributes, or to maximize this component to exploit relationships between hydraulics and surface reflectance. This initial analysis also highlighted the need for improved models of water surface topography in natural rivers.

Journal of Geophysical Research F: Earth Surface

Using stereo satellite imagery to account for ablation, entrainment, and compaction in volume calculations for rock avalanches on Glaciers: Application to the 2016 Lamplugh Rock Avalanche in Glacier Bay National Park, Alaska

The use of preevent and postevent digital elevation models (DEMs) to estimate the volume of rock avalanches on glaciers is complicated by ablation of ice before and after the rock avalanche, scour of material during rock avalanche emplacement, and postevent ablation and compaction of the rock avalanche deposit. We present a model to account for these processes in volume estimates of rock avalanches on glaciers. We applied our model by calculating the volume of the 28 June 2016 Lamplugh rock avalanche in Glacier Bay National Park, Alaska. We derived preevent and postevent 2‐m resolution DEMs from WorldView satellite stereo imagery. Using data from DEM differencing, we reconstructed the rock avalanche and adjacent surfaces at the time of occurrence by accounting for elevation changes due to ablation and scour of the ice surface, and postevent deposit changes. We accounted for uncertainties in our DEMs through precise coregistration and an assessment of relative elevation accuracy in bedrock control areas. The rock avalanche initially displaced 51.7 ± 1.5 Mm3 of intact rock and then scoured and entrained 13.2 ± 2.2 Mm3 of snow and ice during emplacement. We calculated the total deposit volume to be 69.9 ± 7.9 Mm3. Volume estimates that did not account for topographic changes due to ablation, scour, and compaction underestimated the deposit volume by 31.0–46.8 Mm3. Our model provides an improved framework for estimating uncertainties affecting rock avalanche volume measurements in glacial environments. These improvements can contribute to advances in the understanding of rock avalanche hazards and dynamics.

Alaska

Methane fluxes from tropical coastal lagoons surrounded bymangroves, Yucatán, Mexico

Methane concentrations in the water column and emissions to the atmosphere were determined for three tropical coastal lagoons surrounded by mangrove forests on the Yucatán Peninsula, Mexico. Surface water dissolved methane was sampled at different seasons over a period of 2 years in areas representing a wide range of salinities and anthropogenic impacts. The highest surface water methane concentrations (up to 8378 nM) were measured in a polluted canal associated with Terminos Lagoon. In Chelem Lagoon, methane concentrations were typically lower, except in the polluted harbor area (1796 nM). In the relatively pristine Celestún Lagoon, surface water methane concentrations ranged from 41 to 2551 nM. Methane concentrations were negatively correlated with salinity in Celestún, while in Chelem and Terminos high methane concentrations were associated with areas of known pollution inputs, irrespective of salinity. The diffusive methane flux from surface lagoon water to the atmosphere ranged from 0.0023 to 15 mmol CH 4 m −2 d −1 . Flux chamber measurements revealed that direct methane release as ebullition was up to 3 orders of magnitude greater than measured diffusive flux. Coastal mangrove lagoons may therefore be an important natural source of methane to the atmosphere despite their relatively high salinity. Pollution inputs are likely to substantially enhance this flux. Additional statistically rigorous data collected globally are needed to better consider methane fluxes from mangrove-surrounded coastal areas in response to sea level changes and anthropogenic pollution in order to refine projections of future atmospheric methane budgets.

Yucatam

Dissolved organic matter compositional change and biolability during two storm runoff events in a small sgricultural watershed

Agricultural watersheds are globally pervasive, supporting fundamentally different organic matter source, composition, and concentration profiles in comparison to natural systems. Similar to natural systems, agricultural storm runoff exports large amounts of organic carbon from agricultural land into waterways. But intense management of upper soil layers, waterway channelization, wetland and riparian habitat removal, and postharvest vegetation removal promise to uniquely drive organic matter release to waterways. During a winter first flush and a subsequent storm event, this study investigated the influence of a small agricultural watershed on dissolved organic matter (DOM) source, composition, and biolability. Storm water discharge released strongly terrestrial yet biolabile (23 to 32%) dissolved organic carbon (DOC). Following a 21 day bioassay, a parallel factor analysis identified an 80% reduction in a protein-like (phenylpropyl) component (C2) that was previously correlated to lignin phenol concentration, and a 10% reduction in a humic-like, terrestrially sourced component (C4). Storm-driven releases tripled DOC concentration (from 2.8 to 8.7 mg L −1 ) during the first flush event in comparison to base flow and were terrestrially sourced, with an eightfold increase in vascular plant derived lignin phenols (23.0 to 185 μg L −1 ). As inferred from system hydrology, lignin composition, and nitrate as a groundwater tracer, an initial pulse of dilute water from the upstream watershed caused a counterclockwise DOC hysteresis loop. DOC concentrations peaked after 3.5 days, with the delay between peak discharge and peak DOC attributed to storm water hydrology and a period of initial water repellency of agricultural soils, which delayed DOM leaching.

California

Estimation and uncertainty of recent carbon accumulation and vertical accretion in drained and undrained forested peatlands of the southeastern USA

The purpose of this study was to determine how drainage impacts carbon densities and recent rates (past 50 years) of vertical accretion and carbon accumulation in southeastern forested peatlands. We compared these parameters in drained maple-gum (MAPL), Atlantic white cedar (CDR), and pocosin (POC) communities in the Great Dismal Swamp National Wildlife Refuge (GDS) of Virginia/North Carolina and in an intact (undrained) CDR swamp in the Alligator River National Wildlife Refuge (AR) of North Carolina. Peat cores were analyzed for bulk density, percent organic carbon, and 137 Cs and 210 Pb. An uncertainty analysis of both 137 Cs and 210 Pb approaches was used to constrain error at least partially related to mobility of both radioisotopes. GDS peats had lower porosities (89.6% (SD = 1.71) versus 95.3% (0.18)) and higher carbon densities (0.082 (0.021) versus 0.037 (0.009) g C cm −3 ) than AR. Vertical accretion rates (0.10–0.56 cm yr −1 ) were used to estimate a time period of ~84–362 years for reestablishment of peat lost during the 2011 Lateral West fire at the GDS. Carbon accumulation rates ranged from 51 to 389 g C m −2 yr −1 for all sites. In the drained (GDS) versus intact (AR) CDR sites, carbon accumulation rates were similar with 137 Cs (87 GDS versus 92 AR g C m −2 yr −1 ) and somewhat less at the GDS than AR as determined with 210 Pb (111 GDS versus 159 AR g C m −2 yr −1 ). Heightened productivity and high polyphenol content of peat may be responsible for similar rates of carbon accumulation in both drained and intact CDR peatlands.

Alligator Rivers National Wildlife Refuge, Great D

The impact of late Holocene land-use change, climate variability, and sea-level rise on carbon storage in tidal freshwater wetlands on the southeastern United States Coastal Plain

This study examines Holocene impacts of changes in climate, land use, and sea-level rise (SLR) on sediment accretion, carbon accumulation rates (CAR), and vegetation along a transect of tidal freshwater forested wetlands (TFFW) to oligohaline marsh along the Waccamaw River, South Carolina (4 sites) and along the Savannah River, Georgia (4 sites). We use pollen, plant macrofossils, accretion, and CAR from cores, spanning the last 1500-6000 years to test the hypothesis that TFFW have remained stable throughout the late Holocene and that marshes transitioned from TFFW during elevated SLR during the Medieval Climate Anomaly, with further transformation resulting from colonial land-use change. Results show low and stable accretion and CAR through much of the Holocene, despite moderate changes associated with Holocene paleoclimate. In all records, the largest observed change occurred within the last ~400 years, driven by colonial land clearance, shifting terrigenous sediment into riparian wetlands, resulting in order-of-magnitude increases in accretion and C accumulation. The oligohaline marshes transitioned from TFFW ~300-500 years ago, coincident with colonial land clearance. Post-colonial decreases in CAR and accretion occur because of watershed reforestation over the last century. All sites show evidence of recent (decades to century) swamp forest decline due to increasing salinity and tidal inundation from SLR. This study suggests that allochthonous sediment input during colonialization helped maintain TFFW, but that current SLR rates are too high for TFFW to persist, although higher accretion rates in oligohaline marshes increases the resilience of tidal wetlands as they transition from TFFW to marsh.

Georgia, South Carolina

Origin of methane and sources of high concentrations in Los Angeles groundwater

In 2014, samples from 37 monitoring wells at 17 locations, within or near oil fields, and one site >5 km from oil fields, in the Los Angeles Basin, California, were analyzed for dissolved hydrocarbon gas isotopes and abundances. The wells sample a variety of depths of an aquifer system composed of unconsolidated and semiconsolidated sediments under various conditions of confinement. Concentrations of methane in groundwater samples ranged from 0.002 to 150 mg/L—some of the highest concentrations reported in a densely populated urban area. The δ13C and δ2H of the methane ranged from −80.8 to −45.5 per mil (‰) and −249.8 to −134.9‰, respectively, and, along with oxidation‐reduction processes, helped to identify the origin of methane as microbial methanogenesis and CO2 reduction as its main formation pathway. The distribution of methane concentrations and isotopes is consistent with the high concentrations of methane in Los Angeles Basin groundwater originating from relatively shallow microbial production in anoxic or suboxic conditions. Source of the methane is the aquifer sediments rather than the upward migration or leakage of thermogenic methane associated with oil fields in the basin.

California

Direct and indirect effects of tides on ecosystem-scale CO2 exchange in a brackish tidal marsh in Northern California

We investigated the direct and indirect influence of tides on net ecosystem exchange (NEE) of carbon dioxide (CO 2 ) in a temperate brackish tidal marsh. NEE displayed a tidally driven pattern with obvious characteristics at the multiday scale, with greater net CO 2 uptake during spring tides than neap tides. Based on the relative mutual information between NEE and biophysical variables, this was driven by a combination of higher water table depth (WTD), cooler air temperature, and lower vapor pressure deficit (VPD) during spring tides relative to neap tides, as the fortnightly tidal cycle not only influenced water levels but also strongly modulated water and air temperature and VPD. Tides also influenced NEE at shorter timescales, with a reduction in nighttime fluxes during growing season spring tides when the higher of the two semidiurnal tides caused inundation at the site. WTD significantly influenced ecosystem respiration ( R eco ), with lower R eco during spring tides than neap tides. While WTD did not appear to affect ecosystem photosynthesis (gross ecosystem production, GPP) directly, the impact of tides on temperature and VPD influenced GPP, with higher daily light‐use efficiency and photosynthetic activity during spring tides than neap tides when temperature and VPD were lower. The strong direct and indirect influence of tides on NEE across the diel and multiday timescales has important implications for modeling NEE in tidal wetlands and can help inform the timing and frequency of chamber measurements as annual or seasonal net CO 2 uptake may be underestimated if measurements are only taken during nonflooded periods.

California