Search USGS⌕ Search

SEARCH · Search USGS

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 73 records · Page 4Linked to original sources
USGS2013

Hydrodynamics of spur and groove formations on a coral reef

Spur and groove (SAG) formations are found on the fore reefs of many coral reefs worldwide. Although these formations are primarily present in wave‐dominated environments, their effect on wave‐driven hydrodynamics is not well understood. A two‐dimensional, depth‐averaged, phase‐resolving nonlinear Boussinesq model ( funwaveC ) was used to model hydrodynamics on a simplified SAG system. The modeling results show that the SAG formations together with shoaling waves induce a nearshore Lagrangian circulation pattern of counter‐rotating circulation cells. The mechanism driving the modeled flow is an alongshore imbalance between the pressure gradient (PG) and nonlinear wave (NLW) terms in the momentum balance. Variations in model parameters suggest the strongest factors affecting circulation include spur‐normal waves, increased wave height, weak alongshore currents, increased spur height, and decreased bottom drag. The modeled circulation is consistent with a simple scaling analysis based on the dynamical balance of NLW, PG, and bottom stress terms. Model results indicate that the SAG formations efficiently drive circulation cells when the alongshore SAG wavelength allows for the effects of diffraction to create alongshore differences in wave height without changing the mean wave angle.

Journal of Geophysical Research C: Oceans↗
USGS2013

Seasonal circulation over the Catalan inner-shelf (northwest Mediterranean Sea)

This study characterizes the seasonal cycle of the Catalan inner-shelf circulation using observations and complementary numerical results. The relation between seasonal circulation and forcing mechanisms is explored through the depth-averaged momentum balance, for the period between May 2010 and April 2011, when velocity observations were partially available. The monthly-mean along-shelf flow is mainly controlled by the along-shelf pressure gradient and by surface and bottom stresses. During summer, fall, and winter, the along-shelf momentum balance is dominated by the barotropic pressure gradient and local winds. During spring, both wind stress and pressure gradient act in the same direction and are compensated by bottom stress. In the cross-shelf direction the dominant forces are in geostrophic balance, consistent with dynamic altimetry data.

Journal of Geophysical Research C: Oceans↗
USGS1999

Measurements of Reynolds stress profiles in unstratified tidal flow

In this paper we present a method for measuring profiles of turbulence quantities using a broadband acoustic doppler current profiler (ADCP). The method follows previous work on the continental shelf and extends the analysis to develop estimates of the errors associated with the estimation methods. ADCP data was collected in an unstratified channel and the results of the analysis are compared to theory. This comparison shows that the method provides an estimate of the Reynolds stresses, which is unbiased by Doppler noise, and an estimate of the turbulent kinetic energy (TKE) which is biased by an amount proportional to the Doppler noise. The noise in each of these quantities as well as the bias in the TKE match well with the theoretical values produced by the error analysis. The quantification of profiles of Reynolds stresses simultaneous with the measurement of mean velocity profiles allows for extensive analysis of the turbulence of the flow. In this paper, we examine the relation between the turbulence and the mean flow through the calculation of u * , the friction velocity, and C d the coefficient of drag. Finally, we calculate quantities of particular interest in turbulence modeling and analysis, the characteristic lengthscales, including a lengthscale which represents the stream-wise scale of the eddies which dominate the Reynolds stresses.

Journal of Geophysical Research C: Oceans↗
USGS1999

Estimates of bottom roughness length and bottom shear stress in South San Francisco Bay, California

A field investigation of the hydrodynamics and the resuspension and transport of particulate matter in a bottom boundary layer was carried out in South San Francisco Bay (South Bay), California, during March-April 1995. Using broadband acoustic Doppler current profilers, detailed measurements of turbulent mean velocity distribution within 1.5 m above bed have been obtained. A global method of data analysis was used for estimating bottom roughness length z o and bottom shear stress (or friction velocities u * ). Field data have been examined by dividing the time series of velocity profiles into 24-hour periods and independently analyzing the velocity profile time series by flooding and ebbing periods. The global method of solution gives consistent properties of bottom roughness length z o and bottom shear stress values (or friction velocities u * ) in South Bay. Estimated mean values of z o and u * for flooding and ebbing cycles are different. The differences in mean z o and u * are shown to be caused by tidal current flood-ebb inequality, rather than the flooding or ebbing of tidal currents. The bed shear stress correlates well with a reference velocity; the slope of the correlation defines a drag coefficient. Forty-three days of field data in South Bay show two regimes of z o (and drag coefficient) as a function of a reference velocity. When the mean velocity is >25–30 cm s −1 , the ln z o (and thus the drag coefficient) is inversely proportional to the reference velocity. The cause for the reduction of roughness length is hypothesized as sediment erosion due to intensifying tidal currents thereby reducing bed roughness. When the mean velocity is <25–30 cm s −1 , the correlation between z o and the reference velocity is less clear. A plausible explanation of scattered values of z o under this condition may be sediment deposition. Measured sediment data were inadequate to support this hypothesis, but the proposed hypothesis warrants further field investigation.

Journal of Geophysical Research C: Oceans↗
USGS2000

Estimation of wave phase speed and nearshore bathymetry from video imagery

A new remote sensing technique based on video image processing has been developed for the estimation of nearshore bathymetry. The shoreward propagation of waves is measured using pixel intensity time series collected at a cross-shore array of locations using remotely operated video cameras. The incident band is identified, and the cross-spectral matrix is calculated for this band. The cross-shore component of wavenumber is found as the gradient in phase of the first complex empirical orthogonal function of this matrix. Water depth is then inferred from linear wave theory's dispersion relationship. Full bathymetry maps may be measured by collecting data in a large array composed of both cross-shore and longshore lines. Data are collected hourly throughout the day, and a stable, daily estimate of bathymetry is calculated from the median of the hourly estimates. The technique was tested using 30 days of hourly data collected at the SandyDuck experiment in Duck, North Carolina, in October 1997. Errors calculated as the difference between estimated depth and ground truth data show a mean bias of −35 cm (rms error = 91 cm). Expressed as a fraction of the true water depth, the mean percent error was 13% (rms error = 34%). Excluding the region of known wave nonlinearities over the bar crest, the accuracy of the technique improved, and the mean (rms) error was −20 cm (75 cm). Additionally, under low-amplitude swells (wave height H ≤ 1 m), the performance of the technique across the entire profile improved to 6% (29%) of the true water depth with a mean (rms) error of −12 cm (71 cm).

North Carolina↗
USGS2001

Creation of residual flows in a partially stratified estuary

The creation of residual flows in estuaries is examined using acoustic Doppler current profiler data sets from northern San Francisco Bay. The data sets are analyzed using principal component analysis to examine the temporal variability of the flows which create the residual circulation. It is seen that in this periodically and partially stratified estuary the residual flows are created through a series of pulses with strong variability at the 24‐hour timescale, through the interaction of shear, stratification and mixing. This interaction is captured through the use of a dimensionless number, the horizontal Richardson number ( Ri x ), which is developed to examine the local balance between the stratifying and destratifying forces at the tidal timescale. It is seen that Ri x is a valuable parameter in predicting the onset of the residual‐creating events, with a threshold value of ≈3 on ebb tides. This critical value is argued to be a threshold, above which the stratification and shear flow create a feedback effect, each further intensifying the other. This feedback results in a highly variable exchange flow which creates the estuarine residual in intermittent pulses rather than as a steady flow. Although typically attributed to baroclinic forcing, an argument is made that these pulses of residual‐creating exchange flow could be created by barotropic forcing in the presence of variable stratification which is asymmetric between flood and ebb tides. This result poses a great challenge for turbulence modeling, as the timing and magnitude of stratification and shear must be correctly simulated on the tidal timescale in order to reproduce the effects seen in the data sets presented.

Journal of Geophysical Research C: Oceans↗
USGS2002

Across-shelf sediment transport: Interactions between suspended and bed sediment

We use a two-dimensional, time-dependent sediment-transport model to quantify across-shelf transport, deposition, and sorting during wave-driven resuspension events characteristic of those that dominate sediment transport on many continental shelves. Decreases in wave-orbital velocities as water depth increases, and the resulting cross-shelf gradient in bed shear stress favor a net offshore transport of sediment. On wide, flat shelves (slopes ∼0.1%percnt;), these gradients are low, and the depth to which the seabed is reworked depends mainly on bottom shear stress and local sediment availability. On narrow, steep shelves (slopes ∼0.5%percnt;), however, the gradient in bottom stress generates significant cross-shelf suspended sediment flux gradients that create regions of net erosion and deposition. While the magnitude of waves generally determines the water depth to which sediment can be resuspended, erosional and depositional patterns on narrow shelves are sensitive to cross-shelf gradients in wave energy, nonlocal sediment availability, and the direction and magnitude of the cross-shelf current. During energetic waves, cross-shelf divergence of suspended sediment flux can create a coarsened, erosional area on the inner shelf that abuts a region of fine-grained sediment deposition on the mid-to-outer shelf. If currents are strongly shoreward, however, flux divergence leads to erosion over the entire shelf.

Journal of Geophysical Research C: Oceans↗
USGS2000

Numerical model of frazil ice and suspended sediment concentrations and formation of sediment laden ice in the Kara Sea

A one-dimensional (vertical) numerical model of currents, mixing, frazil ice concentration, and suspended sediment concentration has been developed and applied in the shallow southeastern Kara Sea. The objective of the calculations is to determine whether conditions suitable for turbid ice formation can occur during times of rapid cooling and wind- and wave-induced sediment resuspension. Although the model uses a simplistic approach to ice particles and neglects ice-sediment interactions, the results for low-stratification, shallow (∼20-m) freeze-up conditions indicate that the coconcentrations of frazil ice and suspended sediment in the water column are similar to observed concentrations of sediment in turbid ice. This suggests that wave-induced sediment resuspension is a viable mechanism for turbid ice formation, and enrichment mechanisms proposed to explain the high concentrations of sediment in turbid ice relative to sediment concentrations in underlying water may not be necessary in energetic conditions. However, salinity stratification found near the Ob' and Yenisey Rivers damps mixing between ice-laden surface water and sediment-laden bottom water and probably limits incorporation of resuspended sediment into turbid ice until prolonged or repeated wind events mix away the stratification. Sensitivity analyses indicate that shallow (≤20 m), unstratified waters with fine bottom sediment (settling speeds of ∼1 mm s −1 or less) and long open water fetches (>25 km) are ideal conditions for resuspension.

Kara Sea↗
USGS2003

Interaction of lateral baroclinic forcing and turbulence in an estuary

Observations of density and velocity in a channel in northern San Francisco Bay show that the onset of vertical density stratification during flood tides is controlled by the balance between the cross-channel baroclinic pressure gradient and vertical mixing due to turbulence. Profiles of velocity, salinity, temperature, and suspended sediment concentration were measured in transects across Suisun Cutoff, in northern San Francisco Bay, on two days over the 12.5-hour tidal cycle. During flood tides an axial density front developed between fresher water flowing from the shallows of Grizzly Bay into the northern side of Suisun Cutoff and saltier water flowing up the channel. North of the front, transverse currents were driven by the lateral salinity gradient, with a top-to-bottom velocity difference greater than 30 cm/s. South of the front, the secondary circulation was weak, and along-channel velocities were greater than to the north. The gradient Richardson number shows that stratification was stable north of the front, while the water column was turbulently mixed south of the front. Time-series measurements of velocity and salinity demonstrate that the front develops during each tidal cycle. In estuaries, longitudinal dynamics predict less stratification during flood than ebb tides. These data show that stratification can develop during flood tides due to a lateral baroclinic pressure gradient in estuaries with complex bathymetry.

Journal of Geophysical Research C: Oceans↗
USGS2005

Estimating hydrodynamic roughness in a wave-dominated environment with a high-resolution acoustic Doppler profiler

Hydrodynamic roughness is a critical parameter for characterizing bottom drag in boundary layers, and it varies both spatially and temporally due to variation in grain size, bedforms, and saltating sediment. In this paper we investigate temporal variability in hydrodynamic roughness using velocity profiles in the bottom boundary layer measured with a high-resolution acoustic Doppler profiler (PCADP). The data were collected on the ebb-tidal delta off Grays Harbor, Washington, in a mean water depth of 9 m. Significant wave height ranged from 0.5 to 3 m. Bottom roughness has rarely been determined from hydrodynamic measurements under conditions such as these, where energetic waves and medium-to-fine sand produce small bedforms. Friction velocity due to current u*c and apparent bottom roughness z0a were determined from the PCADP burst mean velocity profiles using the law of the wall. Bottom roughness kB was estimated by applying the Grant-Madsen model for wave-current interaction iteratively until the model u*c converged with values determined from the data. The resulting kB values ranged over 3 orders of magnitude (10-1 to 10-4 m) and varied inversely with wave orbital diameter. This range of kB influences predicted bottom shear stress considerably, suggesting that the use of time-varying bottom roughness could significantly improve the accuracy of sediment transport models. Bedform height was estimated from kB and is consistent with both ripple heights predicted by empirical models and bedforms in sonar images collected during the experiment. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗
USGS2004

Suspended sediment and hydrodynamics above mildly sloped long wave ripples

We investigate the spatial and temporal distribution of suspended sediment and the associated hydrodynamics over mildly sloped long wave ripples on the inner shelf. These bedforms had wavelengths of approximately 1 m and heights of approximately 5 cm, in a mean water depth of 4 m. The vertical and temporal structures of the suspended sediment concentration (SSC) are consistent with the entrainment of sediment on the offshore flank of the ripple, and rapid vertical mixing at the time of flow reversal, followed by advection onshore by the onshore fluid motion. This work confirms that the mechanism for sediment suspension above low-amplitude, long wave ripples is similar to the vortex formation process expected over steeper vortex ripples. Numerical simulations of the flow using the Dune2d model indicate that a separated rotational flow structure is generated at the flank of the ripple on the seaward side of the ripple crest, near the time of flow reversal. The simulations indicate that only one vortex is formed during each wave period, in agreement with the field observations. This asymmetry is due mainly to the presence of an offshore mean near-bed current of approximately 6 to 8 cm/s. The SSC is calculated by the model and compared to the field observations. A hydraulic bed roughness of 10 to 15 median grain diameters ( d 50 ) was used in order to match the model prediction to the observed SSC approximately 1 cm above the seabed (cab). However, the modeled SSC and turbulent kinetic energy were significantly lower than the field observations at elevations exceeding approximately 2 cab.

North Carolina↗
USGS2005

Wave spectral energy variability in the northeast Pacific

The dominant characteristics of wave energy variability in the eastern North Pacific are described from NOAA National Data Buoy Center (NDBC) buoy data collected from 1981 to 2003. Ten buoys at distributed locations were selected for comparison based on record duration and data continuity. Long‐period (LP) [ T > 12] s, intermediate‐period [6 ≤ T ≤ 12] s, and short‐period [ T < 6] s wave spectral energy components are considered separately. Empirical orthogonal function (EOF) analyses of monthly wave energy anomalies reveal that all three wave energy components exhibit similar patterns of spatial variability. The dominant mode represents coherent heightened (or diminished) wave energy along the West Coast from Alaska to southern California, as indicated by composites of the 700 hPa height field. The second EOF mode reveals a distinct El Niño‐Southern Oscillation (ENSO)‐associated spatial distribution of wave energy, which occurs when the North Pacific storm track is extended unusually far south or has receded to the north. Monthly means and principal components (PCs) of wave energy levels indicate that the 1997–1998 El Niño winter had the highest basin‐wide wave energy within this record, substantially higher than the 1982–1983 El Niño. An increasing trend in the dominant PC of LP wave energy suggests that storminess has increased in the northeast Pacific since 1980. This trend is emphasized at central eastern North Pacific locations. Patterns of storminess variability are consistent with increasing activity in the central North Pacific as well as the tendency for more extreme waves in the south during El Niño episodes and in the north during La Niña.

Journal of Geophysical Research C: Oceans↗
USGS2006

February 2003 marine atmospheric conditions and the bora over the northern Adriatic

A winter oceanographic field experiment provided an opportunity to examine the atmospheric marine conditions over the northern Adriatic. Mean February winds are from a northeasterly direction over most of the Adriatic and a more northerly direction along the western coast. Wind speeds are fastest in jets over the NE coast during bora events and weakest in the mid-northwestern Adriatic. Diurnal air temperature cycles are smallest on the NE coast and largest in the midwestern Adriatic. The maximum sea-air difference is +10??C on the eastern coast and near zero on the midwestern Adriatic. Boras are northeasterly (from) wind events that sweep off Croatia and Slovenia, bringing slightly colder and drier air over the northern Adriatic. The main bora season is December to March. Winter 2002-2003 was normal for bora events. Synoptic-scale temporal variations are correlated over the northern Adriatic. Fastest Bora winds and highest wind stress over the northern Adriatic is concentrated in four topographically controlled jets. The strongest is the Senj Jet, while the Trieste Jet extends across the entire northern Adriatic. Between each two jets is a weak wind zone. The greatest mean net heat loss is in bora jets in the NE Adriatic, where it was -438 W m-2 and is weakest in the midwestern northern Adriatic, where it was near zero. Wind stress is concentrated over the NE half of Adriatic in four bora jets, while wind stress is weak in the NW Adriatic. There is significant variation in wind stress mean and standard deviation structure over the northern Adriatic with each bora event. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗
USGS2006

Surface drifter derived circulation in the northern and middle Adriatic Sea: Response to wind regime and season

More than 120 satellite-tracked drifters were deployed in the northern and middle Adriatic (NMA) Sea between September 2002 and November 2003, with the purpose of studying the surface circulation at mesoscale to seasonal scale in relation to wind forcing, river runoff, and bottom topography. Pseudo-Eulerian and Lagrangian statistics were calculated from the low-pass-filtered drifter velocity data between September 2002 and December 2003. The structure of the mean circulation is determined with unprecedented high horizontal resolution by the new data. In particular, mean currents, velocity variance, and kinetic energy levels are shown to be maximal in the Western Adriatic Current (WAC). Separating data into seasons, we found that the mean kinetic energy is maximal in fall, with high values also in winter, while it is significantly weaker in summer. High-resolution Local Area Model Italy winds were used to relate the drifter velocities to the wind fields. The surface currents appear to be significantly influenced by the winds. The mean flow during the northeasterly bora regime shows an intensification of the across-basin recirculating currents. In addition, the WAC is strongly intensified both in intensity and in its offshore lateral extension. In the southeasterly sirocco regime, northward flow without recirculation dominates in the eastern half of the basin, while during northwesterly maestro the WAC is enhanced. Separating the data into low and high Po River discharge rates for low-wind conditions shows that the WAC and the velocity fluctuations in front of the Po delta are stronger for high Po River runoff. Lagrangian covariance, diffusivity, and integral time and space scales are larger in the along-basin direction and are maximal in the southern portion of the WAC. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗
USGS2007

Measurements of storm and nonstorm circulation in the northern Adriatic: October 2002 Through April 2003

Fifteen bottom-mounted Acoustic Doppler Current Profilers were deployed from October 2002 through April 2003 in the northern Adriatic Sea. Average transport from the portion of the Western Adriatic Current (WAC) along the Italian slope was 0.1470 ?? 0.0043 Sv, punctuated by bursts of more than twice that amount during storm events. Monthly means were calculated with times of strong wind-driven circulation excluded. These suggest a 2002/2003 seasonal separation consisting of October, December through February, and March through April. An extreme Po River flood influenced November conditions making seasonal categorization difficult. October generally had more kinetic energy and more vertical structure than other months, and near-inertial waves were more frequent in April and October. The Eastern Adriatic Current (EAC)/WAC (i.e. inflow/outflow) system was clearly present in the means for all months. The cyclonic gyre north of the Po River was present October through February. Generally, in the WAC, over 50% of kinetic energy came from vertically uniform monthly mean flows. Elsewhere, eddy kinetic energy was stronger than mean kinetic energy with 10-40% contributions for vertically uniform monthly mean flows, 40-60% for vertically uniform monthly varying flows, and 10-30% for vertically varying monthly varying flows. Mean currents for bora storms indicate enhancement of the EAC/WAC and the cyclonic northern gyre, a shift toward Kvarner Bay in EAC direction, a circulation null point south of the Po, and double-gyre bifurcation of flow at Istria. Strengthening of both the EAC and WAC also occurs during sirocco storms. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗
USGS2007

Bed forms created by simulated waves and currents in a large flume

The morphology and evolution of bed forms created by combinations of waves and currents were investigated using an oscillating plate in a 4-m-wide flume. Current speed ranged from 0 to 30 cm/s, maximum oscillatory velocity ranged from 20 to 48 cm/s, oscillation period was 8 s (except for one run with 12 s period), and the median grain size was 0.27 mm. The angle between oscillations and current was 90°, 60°, or 45°. At the end of each run the sand bed was photographed and ripple dimensions were measured. Ripple wavelength was also determined from sonar images collected throughout the runs. Increasing the ratio of current to wave (i.e., oscillatory) velocity decreased ripple height and wavelength, in part because of the increased fluid excursion during the wave period. Increasing the ratio of current to waves, or decreasing the angle between current and waves, increased the three-dimensionality of bed forms. During the runs, ripple wavelength increased by a factor of about 2. The average number of wave periods for evolution of ripple wavelength to 90% of its final value was 184 for two-dimensional ripples starting from a flat bed. Bed form orientations at the end of each run were compared to four potential controlling factors: the directions of waves, current, maximum instantaneous bed shear stress, and maximum gross bed form normal transport (MGBNT). The directions of waves and of MGBNT were equally good predictors of bed form orientations, and were significantly better than the other two factors.

Journal of Geophysical Research C: Oceans↗
USGS2007

Tidal asymmetry and residual circulation over linear sandbanks and their implication on sediment transport: a process-oriented numerical study

A series of process-oriented numerical simulations is carried out in order to evaluate the relative role of locally generated residual flow and overtides on net sediment transport over linear sandbanks. The idealized bathymetry and forcing are similar to those present in the Norfolk Sandbanks, North Sea. The importance of bottom drag parameterization and bank orientation with respect to the ambient flow is examined in terms of residual flow and overtide generation, and subsequent sediment transport implications are discussed. The results show that although the magnitudes of residual flow and overtides are sensitive to bottom roughness parameterization and bank orientation, the magnitude of the generated residual flow is always larger than that of the locally generated overtides. Also, net sediment transport is always dominated by the nonlinear interaction of the residual flow and the semidiurnal tidal currents, although cross-bank sediment transport can occur even in the absence of a cross-shore residual flow. On the other hand, net sediment divergence/convergence increases as the bottom drag decreases and as bank orientation increases. The sediment erosion/deposition is not symmetric about the crest of the bank, suggesting that originally symmetric banks would have the tendency to become asymmetric.

Journal of Geophysical Research C: Oceans↗
USGS2007

Role of aquifer heterogeneity in fresh groundwater discharge and seawater recycling: An example from the Carmel coast, Israel

A case study is shown in which the pattern of submarine groundwater discharge and of seawater recycling is controlled by local hydrogeological variability. The coastal aquifer in Dor Bay is composed of two units: a partly confined calcaranitic sandstone (Kurkar) and an overlying loose sand. Groundwater in the Kurkar has elevated activities of 222 Rn (∼390 dpm/L) and relatively low 224 Ra/ 223 Ra activity ratios (3–4), while the sand groundwater is significantly less radiogenic (6–90 dpm/L) and shows higher 224 Ra/ 223 Ra ratios. Groundwater discharging from sand-covered areas of the bay has salinities of 16–31 and an average 222 Rn activity of 168 dpm/L, which lies on a mixing line between Rn-rich Kurkar fresh water and Rn-poor seawater. Another key observation is that seawater infiltrates to some extent into onshore sand groundwater, while the fresh water within the submarine Kurkar can be traced up to 40 m offshore. This implies that while fresh water mainly discharges from the Kurkar unit, seawater recycling is limited to the loose sand, and that the discharge from sand-covered areas is a mixture of Kurkar water with recycled seawater. Advection rates from the bay floor were calculated from Rn time series and found to vary between 0 and 36 cm/d, correlating negatively with bay water depth. The average flux was 8.1 cm/d, and it did not seem to change much during March, May, and July 2006. The average amount of fresh water discharging to the bay was 5.0 m 3 /d per meter of shoreline. Radon activity in the sand groundwater also fluctuates due to influx of Kurkar-type groundwater.

Carmel Coast, Dor Bay↗