Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Geophysical Journal”

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 559 records · Page 31Linked to original sources

Continental rupture and the creation of new crust in the Salton Trough rift, southern California and northern Mexico: Results from the Salton Seismic Imaging Project

A refraction and wide-angle reflection seismic profile along the axis of the Salton Trough, California and Mexico, was analyzed to constrain crustal and upper mantle seismic velocity structure during active continental rifting. From the northern Salton Sea to the southern Imperial Valley, the crust is 17-18 km thick and approximately one-dimensional. The transition at depth from Colorado River sediment to underlying crystalline rock is gradual and is not a depositional surface. The crystalline rock from ~3 to ~8 km depth is interpreted as sediment metamorphosed by high heat flow. Deeper felsic crystalline rock could be stretched pre-existing crust or higher grade metamorphosed sediment. The lower crust below ~12 km depth is interpreted to be gabbro emplaced by rift-related magmatic intrusion by underplating. Low upper-mantle velocity indicates high temperature and partial melting. Under the Coachella Valley, sediment thins to the north and the underlying crystalline rock is interpreted as granitic basement. Mafic rock does not exist at 12-18 depth as it does to the south, and a weak reflection suggests Moho at ~28 km depth. Structure in adjacent Mexico has slower mid-crustal velocity and rocks with mantle velocity must be much deeper than in the Imperial Valley. Slower velocity and thicker crust in the Coachella and Mexicali valleys define the rift zone between them to be >100 km wide in the direction of plate motion. North American lithosphere in the central Salton Trough has been rifted apart and is being replaced by new crust created by magmatism, sedimentation, and metamorphism.

Journal of Geophysical Research↗

Origin of the pulse-like signature of shallow long-period volcano seismicity

Short-duration, pulse-like long-period (LP) events are a characteristic type of seismicity accompanying eruptive activity at Mount Etna in Italy in 2004 and 2008 and at Turrialba Volcano in Costa Rica and Ubinas Volcano in Peru in 2009. We use the discrete wave number method to compute the free surface response in the near field of a rectangular tensile crack embedded in a homogeneous elastic half space and to gain insights into the origin of the LP pulses. Two source models are considered, including (1) a vertical fluid-driven crack and (2) a unilateral tensile rupture growing at a fixed sub-Rayleigh velocity with constant opening on a vertical crack. We apply cross correlation to the synthetics and data to demonstrate that a fluid-driven crack provides a natural explanation for these data with realistic source sizes and fluid properties. Our modeling points to shallow sources (<1 km depth), whose signatures are representative of the Rayleigh pulse sampled at epicentral distances >∼1 km. While a slow-rupture failure provides another potential model for these events, the synthetics and resulting fits to the data are not optimal in this model compared to a fluid-driven source. We infer that pulse-like LP signatures are parts of the continuum of responses produced by shallow fluid-driven sources in volcanoes.

Journal of Geophysical Research B: Solid Earth↗

Low stress drops observed for aftershocks of the 2011 Mw 5.7 Prague, Oklahoma, earthquake

In November 2011, three M w ≥ 4.8 earthquakes and thousands of aftershocks occurred along the structurally complex Wilzetta fault system near Prague, Oklahoma. Previous studies suggest that wastewater injection induced a M w 4.8 foreshock, which subsequently triggered a M w 5.7 mainshock. We examine source properties of aftershocks with a standard Brune-type spectral model and jointly solve for seismic moment ( M 0 ), corner frequency ( f 0 ), and kappa ( κ ) with an iterative Gauss-Newton global downhill optimization method. We examine 934 earthquakes with initial moment magnitudes ( M w ) between 0.33 and 4.99 based on the pseudospectral acceleration and recover reasonable M 0 , f 0 , and κ for 87 earthquakes with M w 1.83–3.51 determined by spectral fit. We use M 0 and f 0 to estimate the Brune-type stress drop, assuming a circular fault and shear-wave velocity at the hypocentral depth of the event. Our observations suggest that stress drops range between 0.005 and 4.8 MPa with a median of 0.2 MPa (0.03–26.4 MPa with a median of 1.1 MPa for Madariaga-type), which is significantly lower than typical eastern United States intraplate events (>10 MPa). We find that stress drops correlate weakly with hypocentral depth and magnitude. Additionally, we find the stress drops increase with time after the mainshock, although temporal variation in stress drop is difficult to separate from spatial heterogeneity and changing event locations. The overall low median stress drop suggests that the fault segments may have been primed to fail as a result of high pore fluid pressures, likely related to nearby wastewater injection.

Oklahoma↗

The role of crystallization-driven exsolution on the sulfur mass balance in volcanic arc magmas

The release of large amounts of sulfur to the stratosphere during explosive eruptions affects the radiative balance in the atmosphere and consequentially impacts climate for up to several years after the event. Quantitative estimations of the processes that control the mass balance of sulfur between melt, crystals, and vapor bubbles is needed to better understand the potential sulfur yield of individual eruption events and the conditions that favor large sulfur outputs to the atmosphere. The processes that control sulfur partitioning in magmas are (1) exsolution of volatiles (dominantly H 2 O) during decompression (first boiling) and during isobaric crystallization (second boiling), (2) the crystallization and breakdown of sulfide or sulfate phases in the magma, and (3) the transport of sulfur-rich vapor (gas influx) from deeper unerupted regions of the magma reservoir. Vapor exsolution and the formation/breakdown of sulfur-rich phases can all be considered as closed-system processes where mass balance arguments are generally easier to constrain, whereas the contribution of sulfur by vapor transport (open system process) is more difficult to quantify. The ubiquitous “excess sulfur” problem, which refers to the much higher sulfur mass released during eruptions than what can be accounted for by amount of sulfur originally dissolved in erupted melt, as estimated from melt inclusion sulfur concentrations (the “petrologic estimate”), reflects the challenges in closing the sulfur mass balance between crystals, melt, and vapor before and during a volcanic eruption. In this work, we try to quantify the relative importance of closed- and open-system processes for silicic arc volcanoes using kinetic models of sulfur partitioning during exsolution. Our calculations show that crystallization-induced exsolution (second boiling) can generate a significant fraction of the excess sulfur observed in crystal-rich arc magmas. This result does not negate the important role of vapor migration in sulfur mass balance but rather points out that second boiling (in situ exsolution) can provide the necessary yield to drive the excess sulfur to the levels observed for crystal-rich systems. In contrast, in crystal-poor systems, magma recharge that releases sulfur-rich bubbles is necessary and most likely the primary contributor to sulfur mass balance. Finally, we apply our model to account for the effect of sulfur partitioning during second boiling and its impact on sulfur released during the Cerro Galan supereruption in Argentina (2.08 Ma) and show the potential importance of second boiling in releasing a large amount of sulfur to the atmosphere during the eruption of large crystal-rich ignimbrites.

Journal of Geophysical Research B: Solid Earth↗

Reexamination of the magnitudes for the 1906 and 1922 Chilean earthquakes using Japanese tsunami amplitudes: Implications for source depth constraints

Far-field tsunami records from the Japanese tide gauge network allow the reexamination of the moment magnitudes ( M w ) for the 1906 and 1922 Chilean earthquakes, which to date rely on limited information mainly from seismological observations alone. Tide gauges along the Japanese coast provide extensive records of tsunamis triggered by six great ( M w >8) Chilean earthquakes with instrumentally determined moment magnitudes. These tsunami records are used to explore the dependence of tsunami amplitudes in Japan on the parent earthquake magnitude of Chilean origin. Using the resulting regression parameters together with tide gauge amplitudes measured in Japan we estimate apparent moment magnitudes of M w 8.0–8.2 and M w 8.5–8.6 for the 1906 central and 1922 north-central Chile earthquakes. The large discrepancy of the 1906 magnitude estimated from the tsunami observed in Japan as compared with those previously determined from seismic waves ( M s 8.4) suggests a deeper than average source with reduced tsunami excitation. A deep dislocation along the Chilean megathrust would favor uplift of the coast rather than beneath the sea, giving rise to a smaller tsunami and producing effects consistent with those observed in 1906. The 1922 magnitude inferred from far-field tsunami amplitudes appear to better explain the large extent of damage and the destructive tsunami that were locally observed following the earthquake than the lower seismic magnitudes ( M s 8.3) that were likely affected by the well-known saturation effects. Thus, a repeat of the large 1922 earthquake poses seismic and tsunami hazards in a region identified as a mature seismic gap.

Journal of Geophysical Research B: Solid Earth↗

A new strategy for earthquake focal mechanisms using waveform-correlation-derived relative polarities and cluster analysis: Application to the 2014 Long Valley Caldera earthquake swarm

In microseismicity analyses, reliable focal mechanisms can typically be obtained for only a small subset of located events. We address this limitation here, presenting a framework for determining robust focal mechanisms for entire populations of very small events. To achieve this, we resolve relative P and S wave polarities between pairs of waveforms by using their signed correlation coefficients—a by-product of previously performed precise earthquake relocation. We then use cluster analysis to group events with similar patterns of polarities across the network. Finally, we apply a standard mechanism inversion to the grouped data, using either catalog or correlation-derived P wave polarity data sets. This approach has great potential for enhancing analyses of spatially concentrated microseismicity such as earthquake swarms, mainshock-aftershock sequences, and industrial reservoir stimulation or injection-induced seismic sequences. To demonstrate its utility, we apply this technique to the 2014 Long Valley Caldera earthquake swarm. In our analysis, 85% of the events (7212 out of 8494 located by Shelly et al . [2016]) fall within five well-constrained mechanism clusters, more than 12 times the number with network-determined mechanisms. Of the earthquakes we characterize, 3023 (42%) have magnitudes smaller than 0.0. We find that mechanism variations are strongly associated with corresponding hypocentral structure, yet mechanism heterogeneity also occurs where it cannot be resolved by hypocentral patterns, often confined to small-magnitude events. Small (5–20°) rotations between mechanism orientations and earthquake location trends persist when we apply 3-D velocity models and might reflect a geometry of en echelon, interlinked shear, and dilational faulting.

California↗

Ground-rupturing earthquakes on the northern Big Bend of the San Andreas Fault, California, 800 A.D. to Present

Paleoseismic data on the timing of ground-rupturing earthquakes constrain the recurrence behavior of active faults and can provide insight on the rupture history of a fault if earthquakes dated at neighboring sites overlap in age and are considered correlative. This study presents the evidence and ages for 11 earthquakes that occurred along the Big Bend section of the southern San Andreas Fault at the Frazier Mountain paleoseismic site. The most recent earthquake to rupture the site was the M w 7.7–7.9 Fort Tejon earthquake of 1857. We use over 30 trench excavations to document the structural and sedimentological evolution of a small pull-apart basin that has been repeatedly faulted and folded by ground-rupturing earthquakes. A sedimentation rate of 0.4 cm/yr and abundant organic material for radiocarbon dating contribute to a record that is considered complete since 800 A.D. and includes 10 paleoearthquakes. Earthquakes have ruptured this location on average every ~100 years over the last 1200 years, but individual intervals range from ~22 to 186 years. The coefficient of variation of the length of time between earthquakes (0.7) indicates quasiperiodic behavior, similar to other sites along the southern San Andreas Fault. Comparison with the earthquake chronology at neighboring sites along the fault indicates that only one other 1857-size earthquake could have occurred since 1350 A.D., and since 800 A.D., the Big Bend and Mojave sections have ruptured together at most 50% of the time in M w ≥ 7.3 earthquakes.

California↗

Frictional strength of wet and dry montmorillonite

Montmorillonite is a common mineral in fault zones, and its low strength relative to other common gouge minerals is important in many models of fault rheology. However, the coefficient of friction, μ , varies with degree of saturation and is not well constrained in the literature due to the difficulty of establishing fully drained or fully dried states in the laboratory. We measured μ of both saturated and oven-dried montmorillonite at normal stresses up to 700 MPa. Care was taken to shear saturated samples slowly enough to avoid pore fluid overpressure. For saturated samples, μ increased from 0.10 to 0.28 with applied effective normal stress, while for dry samples μ decreased from 0.78 to 0.45. The steady state rate dependence of friction, ( a − b ), was positive, promoting stable sliding. The wide disparity in reported frictional strengths can be attributed to experimental procedures that promote differing degrees of partial saturation or overpressured pore fluid conditions.

Journal of Geophysical Research↗

Sediment gravity flows triggered by remotely generated earthquake waves

Recent great earthquakes and tsunamis around the world have heightened awareness of the inevitability of similar events occurring within the Cascadia Subduction Zone of the Pacific Northwest. We analyzed seafloor temperature, pressure, and seismic signals, and video stills of sediment-enveloped instruments recorded during the 2011–2015 Cascadia Initiative experiment, and seafloor morphology. Our results led us to suggest that thick accretionary prism sediments amplified and extended seismic wave durations from the 11 April 2012 M w 8.6 Indian Ocean earthquake, located more than 13,500 km away. These waves triggered a sequence of small slope failures on the Cascadia margin that led to sediment gravity flows culminating in turbidity currents. Previous studies have related the triggering of sediment-laden gravity flows and turbidite deposition to local earthquakes, but this is the first study in which the originating seismic event is extremely distant (> 10,000 km). The possibility of remotely triggered slope failures that generate sediment-laden gravity flows should be considered in inferences of recurrence intervals of past great Cascadia earthquakes from turbidite sequences. Future similar studies may provide new understanding of submarine slope failures and turbidity currents and the hazards they pose to seafloor infrastructure and tsunami generation in regions both with and without local earthquakes.

Cascadia subduction zone↗

Integrating puffing and explosions in a general scheme for Strombolian-style activity

Strombolian eruptions are among the most common subaerial styles of explosive volcanism worldwide. Distinctive features of each volcano lead to a correspondingly wide range of variations of magnitude and erupted products, but most papers focus on a single type of event at a single volcano. Here, in order to emphasize the common features underlying this diversity of styles, we scrutinize a database from 35 different erupting vents, including 21 thermal infrared videos from Stromboli (Italy), Etna (Italy), Yasur (Vanuatu), and Batu Tara (Indonesia), from puffing, through rapid explosions to normal explosions, with variable ejection parameters and relative abundance of gas, ash, and bombs. Using field observations and high-speed thermal infrared videos processed by a new algorithm, we identify the distinguishing characteristics of each type of activity and how they may relate and interact. In particular, we record that ash-poor normal explosions may be preceded and followed by the onset or the increase of the puffing activity, while ash-rich explosions are emergent, i.e., with inflation of the free surface followed directly by emission of increasingly large gas pockets. Overall, we see that all Strombolian activities form a continuum arising from a common mechanism and are modulated by the combination of two well-established controls: (1) the length of the bursting gas pocket with respect to the vent diameter and (2) the presence and thickness of a high-viscosity layer in the uppermost part of the volcanic conduit.

Journal of Geophysical Research B: Solid Earth↗

Temporal variation of tectonic tremor activity in southern Taiwan around the 2010 ML6.4 Jiashian Earthquake

Deep tectonic tremor, which is extremely sensitive to small stress variations, could be used to monitor fault-zone processes during large earthquake cycles and aseismic processes before large earthquakes. In this study, we develop an algorithm for the automatic detection and location of tectonic tremor beneath the southern Central Range of Taiwan and examine the spatio-temporal relationship between tremor and the 4 March 2010 ML6.4 Jiashian earthquake, located about 20 km from active tremor sources. We find that tremor in this region has a relatively short duration, short recurrence time, and no consistent correlation with surface GPS data. We find a short-term increase in the tremor rate 19 days before the Jiashian mainshock, and around the time when the tremor rate began to rise, one GPS station recorded a flip in its direction of motion. We hypothesize that tremor is driven by a slow-slip event that preceded the occurrence of the shallower Jiashian mainshock, even though the inferred slip is too small to be observed by all GPS stations. Our study shows that tectonic tremor may reflect stress variation during the pre-nucleation process of a nearby earthquake.

Journal of Geophysical Research B: Solid Earth↗

Processes of multibathyal aragonite undersaturation in the Arctic Ocean

During 3 years of study (2010–2012), the western Arctic Ocean was found to have unique aragonite saturation profiles with up to three distinct aragonite undersaturation zones. This complexity is produced as inflow of Atlantic-derived and Pacific-derived water masses mix with Arctic-derived waters, which are further modified by physiochemical and biological processes. The shallowest aragonite undersaturation zone, from the surface to ∼30 m depth is characterized by relatively low alkalinity and other dissolved ions. Besides local influence of biological processes on aragonite undersaturation of shallow coastal waters, the nature of this zone is consistent with dilution by sea-ice melt and invasion of anthropogenic CO 2 from the atmosphere. A second undersaturated zone at ∼90–220 m depth (salinity ∼31.8–35.4) occurs within the Arctic Halocline and is characterized by elevated p CO 2 and nutrients. The nature of this horizon is consistent with remineralization of organic matter on shallow continental shelves bordering the Canada Basin and the input of the nutrients and CO 2 entrained by currents from the Pacific Inlet. Finally, the deepest aragonite undersaturation zone is at greater than 2000 m depth and is controlled by similar processes as deep aragonite saturation horizons in the Atlantic and Pacific Oceans. The comparatively shallow depth of this deepest aragonite saturation horizon in the Arctic is maintained by relatively low temperatures, and stable chemical composition. Understanding the mechanisms controlling the distribution of these aragonite undersaturation zones, and the time scales over which they operate will be crucial to refine predictive models.

Journal of Geophysical Research: Oceans↗

Behavior of a wave-driven buoyant surface jet on a coral reef

A wave-driven surface buoyant jet exiting a coral reef was studied in order to quantify the amount of water re-entrained over the reef crest. Both moored observations and Lagrangian drifters were used to study the fate of the buoyant jet. To investigate in detail the effects of buoyancy and along-shore flow variations, we developed an idealized numerical model of the system. Consistent with previous work, the ratio of along-shore velocity to jet-velocity and the jet internal Froude number were found to be important determinants of the fate of the jet. In the absence of buoyancy, the entrainment of fluid at the reef crest, creates a significant amount of retention, keeping 60% of water in the reef system. However, when the jet is lighter than the ambient ocean-water, the net effect of buoyancy is to enhance the separation of the jet from shore, leading to a greater export of reef water. Matching observations, our modeling predicts that buoyancy limits retention to 30% of the jet flow for conditions existing on the Moorea reef. Overall, the combination of observations and modeling we present here shows that reef-ocean temperature gradients can play an important role in reef-ocean exchanges.

Journal of Geophysical Research C: Oceans↗

Sediment transport in the presence of large reef bottom roughness

The presence of large bottom roughness, such as that formed by benthic organisms on coral reef flats, has important implications for the size, concentration, and transport of suspended sediment in coastal environments. A 3 week field study was conducted in approximately 1.5 m water depth on the reef flat at Ningaloo Reef, Western Australia, to quantify the cross-reef hydrodynamics and suspended sediment dynamics over the large bottom roughness (∼20–40 cm) at the site. A logarithmic mean current profile consistently developed above the height of the roughness; however, the flow was substantially reduced below the height of the roughness (canopy region). Shear velocities inferred from the logarithmic profile and Reynolds stresses measured at the top of the roughness, which are traditionally used in predictive sediment transport formulations, were similar but much larger than that required to suspend the relatively coarse sediment present at the bed. Importantly, these stresses did not represent the stresses imparted on the sediment measured in suspension and are therefore not relevant to the description of suspended sediment transport in systems with large bottom roughness. Estimates of the bed shear stresses that accounted for the reduced near-bed flow in the presence of large roughness vastly improved the relationship between the predicted and observed grain sizes that were in suspension. Thus, the impact of roughness, not only on the overlying flow but also on bed stresses, must be accounted for to accurately estimate suspended sediment transport in regions with large bottom roughness, a common feature of many shallow coastal ecosystems.

Journal of Geophysical Research C: Oceans↗

Identification and classification of very low frequency waves on a coral reef flat

Very low frequency (VLF, 0.001&ndash;0.005 Hz) waves are important drivers of flooding of low-lying coral reef-islands. In particular, VLF wave resonance is known to drive large wave runup and subsequent overwash. Using a 5 month data set of water levels and waves collected along a cross-reef transect on Roi-Namur Island in the Republic of the Marshall Islands, the observed VLF motions were categorized into four different classes: (1) resonant, (2) (nonresonant) standing, (3) progressive-growing, and (4) progressive-dissipative waves. Each VLF class is set by the reef flat water depth and, in the case of resonance, the incident-band offshore wave period. Using an improved method to identify VLF wave resonance, we find that VLF wave resonance caused prolonged (&sim;0.5&ndash;6.0 h), large-amplitude water surface oscillations at the inner reef flat ranging in wave height from 0.14 to 0.83 m. It was induced by relatively long-period, grouped, incident-band waves, and occurred under both storm and nonstorm conditions. Moreover, observed resonant VLF waves had nonlinear, bore-like wave shapes, which likely have a larger impact on the shoreline than regular, sinusoidal waveforms. As an alternative technique to the commonly used Fast Fourier Transformation, we propose the Hilbert-Huang Transformation that is more computationally expensive but can capture the wave shape more accurately. This research demonstrates that understanding VLF waves on reef flats is important for evaluating coastal flooding hazards.

Roi-Namur Island↗

Barrier island breach evolution: Alongshore transport and bay-ocean pressure gradient interactions

Physical processes controlling repeated openings and closures of a barrier island breach between a bay and the open ocean are studied using aerial photographs and atmospheric and hydrodynamic observations. The breach site is located on Pea Island along the Outer Banks, separating Pamlico Sound from the Atlantic Ocean. Wind direction was a major control on the pressure gradients between the bay and the ocean to drive flows that initiate or maintain the breach opening. Alongshore sediment flux was found to be a major contributor to breach closure. During the analysis period from 2011 to 2016, three hurricanes had major impacts on the breach. First, Hurricane Irene opened the breach with wind-driven flow from bay to ocean in August 2011. Hurricane Sandy in October 2012 quadrupled the channel width from pressure gradient flows due to water levels that were first higher on the ocean side and then higher on the bay side. The breach closed sometime in Spring 2013, most likely due to an event associated with strong alongshore sediment flux but minimal ocean-bay pressure gradients. Then, in July 2014, Hurricane Arthur briefly opened the breach again from the bay side, in a similar fashion to Irene. In summary, opening and closure of breaches are shown to follow a dynamic and episodic balance between along-channel pressure gradient driven flows and alongshore sediment fluxes.

Pea Island↗

Rip currents and alongshore flows in single channels dredged in the surf zone

To investigate the dynamics of flows near nonuniform bathymetry, single channels (on average 30 m wide and 1.5 m deep) were dredged across the surf zone at five different times, and the subsequent evolution of currents and morphology was observed for a range of wave and tidal conditions. In addition, circulation was simulated with the numerical modeling system COAWST, initialized with the observed incident waves and channel bathymetry, and with an extended set of wave conditions and channel geometries. The simulated flows are consistent with alongshore flows and rip-current circulation patterns observed in the surf zone. Near the offshore-directed flows that develop in the channel, the dominant terms in modeled momentum balances are wave-breaking accelerations, pressure gradients, advection, and the vortex force. The balances vary spatially, and are sensitive to wave conditions and the channel geometry. The observed and modeled maximum offshore-directed flow speeds are correlated with a parameter based on the alongshore gradient in breaking-wave-driven-setup across the nonuniform bathymetry (a function of wave height and angle, water depths in the channel and on the sandbar, and a breaking threshold) and the breaking-wave-driven alongshore flow speed. The offshore-directed flow speed increases with dissipation on the bar and reaches a maximum (when the surf zone is saturated) set by the vertical scale of the bathymetric variability.

Journal of Geophysical Research: Oceans↗

Observations and a linear model of water level in an interconnected inlet-bay system

A system of barrier islands and back-barrier bays occurs along southern Long Island, New York, and in many coastal areas worldwide. Characterizing the bay physical response to water level fluctuations is needed to understand flooding during extreme events and evaluate their relation to geomorphological changes. Offshore sea level is one of the main drivers of water level fluctuations in semienclosed back-barrier bays. We analyzed observed water levels (October 2007 to November 2015) and developed analytical models to better understand bay water level along southern Long Island. An increase (∼0.02 m change in 0.17 m amplitude) in the dominant M 2 tidal amplitude (containing the largest fraction of the variability) was observed in Great South Bay during mid-2014. The observed changes in both tidal amplitude and bay water level transfer from offshore were related to the dredging of nearby inlets and possibly the changing size of a breach across Fire Island caused by Hurricane Sandy (after December 2012). The bay response was independent of the magnitude of the fluctuations (e.g., storms) at a specific frequency. An analytical model that incorporates bay and inlet dimensions reproduced the observed transfer function in Great South Bay and surrounding areas. The model predicts the transfer function in Moriches and Shinnecock bays where long-term observations were not available. The model is a simplified tool to investigate changes in bay water level and enables the evaluation of future conditions and alternative geomorphological settings.

Journal of Geophysical Research C: Oceans↗