JOURNAL OF GEOPHYSICAL RESEARCH, VOL. ???, XXXX, DOI:10.1029/, 1 2 Surfzone to inner-shelf exchange estimated from dye tracer balances K. Hally-Rosendahl, Integrative Oceanography Division, Scripps Institution of Oceanography, UCSD, 9500 Gilman Dr., La Jolla CA 92093-0209 (kai@coast.ucsd.edu) F. Feddersen, Integrative Oceanography Division, Scripps Institution of Oceanography, UCSD, 9500 Gilman Dr., La Jolla CA 92093-0209 (falk@coast.ucsd.edu) D. B. Clark, Woods Hole Oceanographic Institution, MS#12, 266 Woods Hole Rd., Woods Hole, MA 02543 (dclark@whoi.edu) R. T. Guza, Integrative Oceanography Division, Scripps Institution of Oceanography, UCSD, 9500 Gilman Dr., La Jolla CA 92093-0209 (rguza@ucsd.edu) D R A F T March 17, 2015, 12:07pm D R A F T X-2 3 4 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE Abstract. Surfzone and inner-shelf tracer dispersion are observed at an approximately 5 alongshore-uniform beach. Fluorescent Rhodamine WT dye tracer, released 6 near the shoreline continuously for 6.5 h, is advected alongshore by break- 7 ing wave- and wind-driven currents, and ejected offshore from the surfzone 8 to the inner-shelf by transient rip currents. Aerial multispectral imaging of 9 dye concentration and in situ measurements of dye, waves, and currents pro- 10 vide tracer transport and dilution observations spanning approximately 400 11 m cross-shore and 3 km alongshore. The alongshore dilution of near-shoreline 12 dye follows weak power-law decay; a 10-fold increase in downstream distance 13 from the source only reduces the concentration by 50%. Surfzone and inner- 14 shelf dye mass balances close, and in 5 h roughly 1/2 of the surfzone-released 15 dye is transported offshore to the inner-shelf. Observed cross-shore transports 16 from the surfzone to the inner-shelf are parameterized well using a bulk ex- 17 change velocity (correlation r2 = 0.85 and best fit slope 0.7). The best fit 18 cross-shore exchange velocity u∗ ≈ 1 × 10−2 m s−1 is similar to u∗ es- 19 timated using temperature observations on another day with similar wave 20 conditions, confirming the dominance of transient rip currents in surfzone 21 to inner-shelf cross-shore exchange during moderate waves at an alongshore- 22 uniform beach. D R A F T March 17, 2015, 12:07pm D R A F T HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE X-3 1. Introduction 23 The nearshore region, consisting of the surfzone (shoreline to xb , the seaward bound- 24 ary of depth-limited wave breaking) and the inner-shelf (xb to approximately 20 m water 25 depth), is vitally important to coastal economies, recreation, and human and ecosystem 26 health. However, nearshore water quality is often compromised by terrestrial runoff and 27 offshore waste disposal [e.g., Koh and Brooks, 1975; Schiff et al., 2000; Halpern et al., 2008]. 28 Globally, microbial pathogen exposure from polluted nearshore water causes an estimated 29 120 million gastrointestinal illnesses and 50 million severe respiratory illnesses annually 30 [Dorfman and Stoner , 2012] with significant economic impacts. Furthermore, excess nutri- 31 ents in polluted runoff can spur rapid growth of harmful algal blooms (HABs), damaging 32 ecosystems and causing serious and even life-threatening human illnesses through direct 33 ocean exposure or consumption of algal-contaminated seafood [Dorfman and Haren, 2013]. 34 Pathogens, nutrients, HABs, and other contaminants all act as nearshore tracers, the 35 transport and dilution of which are governed by surfzone and inner-shelf physical pro- 36 cesses. Yet, despite the detriment of contaminated coastal water to our health and econ- 37 omy, understanding of nearshore tracer transport and mixing remains relatively poor. 38 Several field experiments [e.g., Harris et al., 1963; Inman et al., 1971; Grant et al., 2005] 39 have tracked near-shoreline-released dye, simulating point source tracers of terrestrial or 40 surfzone origin. However, these studies were limited by sparse sampling and/or small spa- 41 tiotemporal extents of quantitative observations. A rapid-sampling, jetski-based dye mea- 42 surement platform [Clark et al., 2009] provided greatly improved observations. Analyses of 43 dye plume evolution at Huntington Beach, California (HB06) showed that surfzone cross- D R A F T March 17, 2015, 12:07pm D R A F T X-4 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE 44 shore tracer dispersion is dominated by horizontal eddies [Clark et al., 2010; Feddersen 45 et al., 2011; Clark et al., 2011] forced by finite crest length wave breaking [Peregrine, 1998; 46 Spydell and Feddersen, 2009; Clark et al., 2012; Feddersen, 2014]. However, HB06 obser- 47 vations were limited to ≤ 2 hours and usually < 400 m downstream of the dye source with 48 analyses specifically restricted to surfzone-contained portions of the dye plumes. While 49 shoreline-source tracers are first transported and mixed within the surfzone, their fate 50 is ultimately determined by exchange with the inner-shelf [Hally-Rosendahl et al., 2014]. 51 An improved understanding of long distance surfzone tracer dilution requires quantitative 52 estimates of net cross-shore exchange between the the surfzone and inner-shelf. 53 The surfzone and inner-shelf are governed by drastically different dynamics. The sur- 54 fzone is dominated by breaking-wave-driven currents [e.g., Thornton and Guza, 1986] 55 and horizontal eddies [e.g., Peregrine, 1998; Clark et al., 2012], whereas the inner-shelf is 56 forced by a combination of wind, tides, buoyancy, and both surface and internal waves 57 [e.g., Lucas et al., 2011; Lentz and Fewings, 2012; Kumar et al., 2014; Sinnett and Fed- 58 dersen, 2014]. Furthermore, the surfzone is vertically well mixed while the inner-shelf 59 can be strongly stratified immediately offshore of the surfzone [Hally-Rosendahl et al., 60 2014]. The intersection of, and exchange between, these dynamically different regions is 61 particularly complex. The fall 2009 IB09 experiment at Imperial Beach, California was 62 designed to observe shoreline-released dye with better resolution and over longer times 63 and distances than preceding studies, and particularly to quantify surfzone and inner-shelf 64 exchange. Hally-Rosendahl et al. [2014] analyzed 29 Sept dye and temperature observa- 65 tions across the surfzone and inner-shelf, spanning 7-12 hours and 700 m alongshore. 66 Downstream dye dilution followed a power-law decay and was weaker than that expected D R A F T March 17, 2015, 12:07pm D R A F T HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE X-5 67 for a Fickian diffusive process as observed and modeled for surfzone-contained portions 68 of dye plumes [Clark et al., 2010, 2011]. Computed using temperature observations, a 69 bulk cross-shore exchange velocity ≈ 1 × 10−2 m s−1 was 10-15× larger than the estimated 70 Stokes drift driven exchange velocity, indicating that the observed transient rip currents 71 dominated surfzone/inner-shelf cross-shore exchange at the alongshore-uniform Imperial 72 Beach [Hally-Rosendahl et al., 2014]. However, the cross-shore dye transport was not mea- 73 sured, as coupled dye and velocity observations at high temporal and spatial resolutions 74 over large alongshore distances are required to resolve the transient and small scale rip 75 current ejection structure. 76 Here, aerial and in situ dye observations on 13 Oct, spanning the surfzone and inner-shelf 77 for even greater alongshore distances (up to 3 km) than 29 Sept, are used to investigate 78 far-downstream dye dilution and cross-shore transport estimated from surfzone and inner- 79 shelf dye mass balances. The experiment site, dye release, instrument platforms, and 80 sampling schemes are described in section 2. Wave and alongshore current conditions 81 are presented in section 3.1. In section 3.2, aerial dye observations are described, and 82 time periods and spatial regions for subsequent analyses are established. Surfzone cross- 83 shore dye structure is described in section 3.3. Alongshore dye dilution and transport are 84 examined in sections 3.4 and 3.5, respectively. Sections 4.1-4.3 present dye mass balances 85 (estimation methods are described in Appendix A). The closure of these balances allows 86 for observational estimates of cross-shore dye transport (section 4.4) which are compared 87 to parameterized estimates in section 5.1. Section 6 is a summary. 2. IB09 Experiment D R A F T March 17, 2015, 12:07pm D R A F T X-6 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE 2.1. Field Site and Coordinate System 88 IB09 field observations were acquired during fall 2009 at Imperial Beach, California 89 (32.6◦ N, 117.1◦ W), a west (269.6◦ ) facing beach with an approximately straight shoreline 90 (Figure 1). In the right-handed coordinate system, cross-shore coordinate x increases 91 negatively seaward (x = 0 m at the mean shoreline), alongshore coordinate y increases 92 positively toward the north (y = 0 m at the dye release location), and vertical coordinate 93 z increases positively upward (z = 0 at mean sea level). The dye release discussed here 94 took place on 13 October, 2009. Bathymetry surveys from 9 October and 19 October 95 were similar, and are averaged to give a representative bathymetry for 13 October that is 96 approximately alongshore-uniform (Figure 1). All times are in PDT. 2.2. Dye Release 97 Fluorescent Rhodamine WT dye (2.1 × 108 parts per billion (ppb)) was released con- 98 tinuously at 2.4 mL s−1 near the shoreline at (x, y) = (−10, 0) m for approximately 6.5 h 99 (10:39-17:07 h). Visual observations suggested rapid vertical mixing, and measured dye 100 concentrations were reduced from O(108 ) ppb to O(102 ) ppb within 10 m of the release. 101 Therefore, the dye specific gravity was quickly reduced from 1.2 to ≈ 1. Rhodamine WT 102 has a photochemical decay e-folding time of approximately 667 h of sunlight [e.g., Smart 103 and Laidlaw , 1977]; decay over the ≈ 9 h of sunlight during this study is negligible. 2.3. In Situ Instrumentation: Surfzone and Inner-Shelf 104 2.3.1. Cross-Shore Array 105 A 125 m-long cross-shore array of six fixed, near-bed instrument frames (denoted f1-f6, 106 onshore to offshore) was deployed from near the shoreline to approximately 4 m water 107 depth (diamonds, Figure 1). The frames held Paros pressure sensors, SonTek acous- D R A F T March 17, 2015, 12:07pm D R A F T HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE X-7 108 tic Doppler velocimeters (ADVs), Yellow Springs Instrument Company thermistors, and 109 WET Labs ECO Triplet fluorometers (hereafter ET) to measure dye concentration D. 110 One frame (f4), located near the seaward edge of the surfzone, held instruments at three 111 different vertical locations (0.2, 0.7, and 1.3 m above the bed). Cross-shore array instru- 112 ments sampled for 51 min each hour, with the remaining 9 min used by the ADVs to 113 estimate bed location [Feddersen, 2012; Spydell et al., 2014]. On 13 Oct, the dye release 114 (y = 0 m) was 248 m south of this cross-shore array (Figure 1). 115 116 2.3.2. Surfzone Near-Shoreline Alongshore Array Four thermistor-equipped ETs were deployed near the shoreline at y = 82, 546, 1069, 1662 m 117 (circles, Figure 1), referred to as SA1-SA4, respectively. For some analyses, ET data from 118 f2 (at yf = 248 m) are used in conjunction with data from SA1-SA4. The ET on f2 119 sampled throughout (and after) the dye release, while SA1-SA4 were deployed after the 120 dye release started. 121 2.3.3. Cross-Shore Jetski Transects 122 Surface dye concentration and temperature were measured with fluorometers and ther- 123 mistors mounted on two GPS-tracked jetskis [Clark et al., 2009] that drove repeated cross- 124 shore transects from x ≈ −300 m to the shoreline (e.g., Figure 1) at various designated 125 alongshore locations between y = 5 m and y ≈ 2 km. The alongshore spacing between 126 transects varied from approximately 20 m (near the release) to 300 m (far downstream of 127 the release). Analyses only include shoreward transects, when jetskis were driven immedi- 128 ately in front of bores to minimize turbidity from bubbles and suspended sand. Seaward 129 transects, sometimes corrupted when jetskis swerved or became airborne jumping over 130 waves, are discarded. D R A F T March 17, 2015, 12:07pm D R A F T X-8 131 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE 2.3.4. Inner-Shelf Alongshore Boat Transects 132 Offshore of the surfzone, the vertical and alongshore structure of dye concentration 133 and temperature were measured with a vertical array of five thermistor-equipped ETs 134 towed alongshore behind a small boat. The vertical array sampled from z = −1 to 135 −3 m at 0.5 m spacing. During 14:06-17:43 h, repeated ≈ 2 km-long alongshore transects 136 (e.g., Figure 1) were driven at roughly 1 m s−1 at a mean cross-shore location nominally 137 twice the surfzone width. The transects were approximately shore-parallel with deviations 138 to avoid large waves. 2.4. Remote Sensing Instrumentation: Inner-Shelf 139 Aerial observations of near-surface dye concentration were obtained from a small plane 140 with a multispectral camera system and coupled global positioning and inertial navigation 141 systems [Clark et al., 2014]. Two cameras captured images near the peak excitation 142 and emission wavelengths of the fluorescent Rhodamine WT. Dye concentrations were 143 determined by calibrating the ratio of emission to excitation radiances with coincident 144 in situ data. Aerial dye concentration errors range from ±1.5 ppb near D = 0 ppb to 145 ±4.5 ppb near D = 20 ppb. The georeferenced aerial images are combined into mosaics 146 and regridded onto a rectangular grid with 2 m × 2 m lateral resolution. See Clark et al. 147 [2014] for details. 148 Between 11:21 and 15:32 h, 23 mosaic images were obtained, each separated by roughly 149 6 min (with a longer gap from 13:08 to 14:56 h). The dye field was imaged from the 150 shoreline to roughly 400 m offshore and from the release to roughly 3 km downstream. 151 Pixels with excitation image brightness above an empirical threshold [Clark et al., 2014] 152 owing to sun glitter or white foam from breaking waves are discarded. The surfzone is D R A F T March 17, 2015, 12:07pm D R A F T HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE X-9 153 therefore often poorly resolved, and quantitative analyses of aerial images are confined to 154 the inner-shelf. 155 2.4.1. Corrections to Measured Dye Fluorescence 156 All aerial and in situ dye observations are corrected for temperature [Smart and Laidlaw , 157 1977], and all in situ dye observations are corrected for turbidity [Clark et al., 2009]. 158 Corrected D typically differ from measured D by less than 5%. 3. Observations 3.1. Wave, Wind, and Alongshore Current Conditions 159 During the dye release, the incident wave field (with peak period Tp = 13 s) is relatively 160 constant, and the tide varies less than 0.7 m (low tide at 12:33 h). The release-averaged 161 significant wave height Hs (x) shoals to a maximum of 0.87 m at f4 (break point xb = 162 −81 m, mean breaking depth hb = 2.1 m, Figures 2a and 2c). The mean alongshore 163 current V (x) is northward (positive) at all f1-f6 locations, with a near-shoreline maximum 164 of 0.40 m s−1 (Figure 2b). Offshore, V decreases to 0.12 m s−1 at the seaward surfzone 165 boundary xb and then increases slightly to 0.17 m s−1 at inner-shelf f6 (x = −135 m, 166 Figure 2b). The mean surfzone (f1-f4) alongshore current is VSZ = 0.22 m s−1 , and the 167 mean inner-shelf (f4-f6) alongshore current is VIS = 0.14 m s−1 . Wind is from the south 168 at 4-7 m s−1 . 3.2. Inner-shelf Surface Dye Evolution 169 Aerial images (e.g., Figures 3a-3f) spanning 0:42-4:53 h after the t0 = 10:39 h start of the 170 dye release are partitioned into three time periods based on temporal gaps in images and 171 the dye plume evolution: period I (early in release, 11:21-11:53 h), period II (mid-release, D R A F T March 17, 2015, 12:07pm D R A F T X - 10 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE 172 12:08-13:01 h), and period III (late in release, 14:56-15:32 h). Approximately 40 min after 173 the release begins (Figure 3a, period I), surfzone dye has advected about 600 m alongshore 174 at ≈ 0.25 m s−1 , consistent with in situ VSZ = 0.22 m s−1 (Figure 2b). Surfzone dye is 175 ejected onto the inner-shelf in narrow (≈ 50 m) alongshore bands (Figure 3a), presumably 176 due to transient rip currents [e.g., Hally-Rosendahl et al., 2014]. As the dye release contin- 177 ues (Figures 3b-3d, period II), the leading portion of inner-shelf dye is alongshore-patchy 178 with length scales ≈ 50 m (as in Figure 3a). Behind the leading edge, slower alongshore 179 advection of inner-shelf dye (e.g., the feature at y ≈ 1250 m and 1500 m in Figures 3e 180 and 3f, respectively, period III) is apparent at a speed of ≈ 0.15 m s−1 , consistent with 181 in situ VIS = 0.14 m s−1 (Figure 2b). At these longer times and downstream distances 182 (Figures 3e and 3f, period III), inner-shelf dye advects alongshore, disperses cross-shore, 183 and moves to larger alongshore length scales. In particular, Figures 3e and 3f reveal a 184 coherent nearshore eddy feature (at y ≈ 1250 m and 1500 m, respectively) with an along- 185 shore length scale ≈ 300 m, roughly six times larger than the length scales of inner-shelf 186 dye patches when recently ejected from the surfzone. 187 In addition to the temporal partitioning into periods I, II, and III (defined above), 188 dye is cross-shore partitioned into the surfzone and inner-shelf (separated by xb = −81 m, 189 section 3.1) and alongshore partitioned into near- and far-field regions A and B, separated 190 by the cross-shore frame array at yf = 248 m (Figure 3a). 191 The leading alongshore edge of the dye plume yp (t) is defined as the northernmost 192 location where aerial-imaged inner-shelf D exceeds 3 ppb within 40 m of xb (green triangles 193 in Figures 3a-3e). The plume leading edge yp (t) increases roughly linearly during each 194 time period, with the fastest advance during period III (Figure 4). The associated mean D R A F T March 17, 2015, 12:07pm D R A F T HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE X - 11 195 alongshore velocity over periods I, II, and III is 0.17 m s−1 (Figure 4), between the surfzone 196 and inner-shelf means VSZ = 0.22 m s−1 and VIS = 0.14 m s−1 . 3.3. Cross-surfzone Mean Dye Profiles 197 Time-averaged surface dye profiles D(x, yj ) from repeated jetski cross-shore transects 198 at designated alongshore locations yj are cross- and alongshore-binned corresponding to 199 where near-shoreline dye is released (y = 0 m) and measured (SA1, f2, SA2, SA3, and 200 SA4, Figure 1). Near the release (y = 14 m), dye concentration is high (≈ 80 ppb) 201 near the shoreline and decays to ≈ 10 ppb near xb (Figure 5). Immediately downstream 202 (y = 87 m), as dye is dispersed offshore, the mean dye cross-surfzone profile begins to 203 flatten. At y ≥ 207 m, dye is well mixed across the surfzone (Figure 5), indicating that 204 surfzone-representative D can be estimated with near-shoreline measurements. 3.4. Near-Shoreline Alongshore Dye Dilution 205 Here, alongshore dye dilution is examined with near-shoreline in situ data from SA1- 206 SA4 and f2 (yellow circles and diamond, respectively, Figure 3). The SA1-SA4 data start 207 times (beginning progressively later with alongshore distance from the release) indicate 208 SA1-SA4 deployment times instead of dye arrival times at those locations (Figure 6). 209 With the exception of SA3 (y = 1069 m), SA instruments were recovered shortly after 210 the dye release ended. 211 When the dye field is roughly stationary, mean near-shoreline dye concentration decays 212 with downstream distance from the release (Figures 6 and 7) because dye is dispersed 213 cross-shore onto the inner-shelf as it is advected alongshore (Figure 3). Near-shoreline 214 dye variability also decreases significantly with y (Figure 6 and vertical bars in Figure 7). D R A F T March 17, 2015, 12:07pm D R A F T X - 12 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE 215 For example, near the release at y = 248 m, D varies between 0-80 ppb with a mean of 216 15 ppb, while at y = 1662 m, D varies between 4-12 ppb with a mean of 8 ppb (Figures 6b, 217 218 6e, and 7). Furthermore, the time scale of dye variability increases with downstream −1/2 2 2 distance from the release (Figures 6a-6e); the characteristic time scale dD /D dt 219 increases monotonically from 71 s at y = 82 m to 584 s at y = 1662 m. Together with 220 the increasing time scale, the reduction of D variability relative to the mean suggests 221 that surfzone alongshore dye length scales also increase with y, consistent with inner-shelf 222 observations (Figures 3a-3f). The near-shoreline mean dye dilutes following a power-law, D = D0 (y/y0 )α , (1) 223 where y0 = 1 m is chosen for simplicity. The best fit constants are D0 = 98 ppb and 224 α = −0.3 (dashed line in Figure 7). Similar downstream dilution structure was ob- 225 served over shorter distances (700 m) during a 29 Sept dye release with α = −0.2 [Hally- 226 Rosendahl et al., 2014]. The downstream dilution power-law exponents α = (−0.3, −0.2) 227 are smaller than that expected for a shoreline Fickian diffusive process (α = −0.5, dye flux 228 proportional to mean concentration gradient). Deviation from α ≈ −0.5 observed and 229 modeled for surfzone-contained dye plumes [Clark et al., 2010, 2011] indicates that the 230 dispersion is not Fickian after dye saturates the surfzone and spreads onto the inner-shelf 231 [Hally-Rosendahl et al., 2014]. Note that α = −0.3 corresponds to relatively weak decay; 232 a 10-fold increase in y only reduces D by 50%. 3.5. Alongshore Dye Transport D R A F T March 17, 2015, 12:07pm D R A F T HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE X - 13 Dye is advected by the northward alongshore current V (Figures 2b and 3) from the release location (x, y) = (−10, 0) m past the cross-shore array at yf = 248 m. The alongshore dye transport T y,A/B from region A to region B through yf (Figure 3, diamonds) is estimated for the surfzone, y,A/B TSZ (t) = Z 0 d(x, t)V (x, t)D(x, t) dx, (2) d(x, t)V (x, t)D(x, t) dx, (3) xb and the inner-shelf, y,A/B TIS (t) = Z xb xf6 233 using in situ, 30 s-averaged total water depth d = h + η, alongshore current V , and D at 234 f1-f6, assuming vertically uniform V (x, t) and D(x, t). For the surfzone, this assumption 235 is a good approximation. The vertical structure of surfzone dye concentration is measured 236 at the seaward surfzone boundary xb in mean water depth hb = 2.1 m, where three ETs on 237 f4 are located 0.2, 0.7, and 1.3 m above the bed (Figure 2c). As for a 29 Sept dye release 238 [Hally-Rosendahl et al., 2014], surfzone D is vertically uniform on 13 Oct (Figure 8). 239 However, inner-shelf D is not necessarily vertically uniform, as thermal stratification can 240 inhibit inner-shelf vertical dye mixing, even immediately offshore of the vertically mixed 241 surfzone (not shown here, similar to 29 Sept dye release in Hally-Rosendahl et al. [2014], 242 Figure 15). Inner-shelf V may also be vertically sheared, likely larger in the upper water 243 column, driven by southerly wind. Lastly, dye at yf sometimes extends offshore of f6 244 (e.g., Figure 3), but V measurements, and therefore the extent of cross-shore integration 245 for (3), are limited to xf6 = −135 m. For these reasons, TIS 246 247 y,A/B is biased low. y,A/B and TIS The surfzone and inner-shelf alongshore dye transports TSZ y,A/B generally vary between approximately 0-1000 ppb m3 s−1 and 0-400 ppb m3 s−1 , respectively (Figures 9a D R A F T March 17, 2015, 12:07pm D R A F T X - 14 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE y,A/B 248 and 9b). Averaged during the release period, T SZ y,A/B = 320 ppb m3 s−1 and T IS = 249 76 ppb m3 s−1 , about 62% and 15% of the dye release rate Q, respectively. Therefore, 77% 250 of the released dye is alongshore-transported past yf = 248 m between the shoreline and 251 xf6 = −135 m (recall TIS y,A/B 252 is biased low). The cumulative (time-integrated) alongshore dye transports at yf for the surfzone and Rt y,A/B T t0 SZ (τ ) dτ and Rt y,A/B TIS (τ ) dτ , respectively, where t0 = 10:39 h is the 253 inner-shelf are 254 dye release start time and TSZ 255 surfzone alongshore transport 256 ure 9c) with small steps corresponding to pulses of TSZ 257 after the last of the dye is advected past the cross-shore array, approximately four times 258 as much dye has been alongshore-transported through yf between the shoreline and f4 as 259 between f4 and f6 (Figure 9c). y,A/B Rt t0 t0 y,A/B and TIS y,A/B TSZ are defined in (2) and (3). The cumulative dτ is roughly linear during the dye release (Figy,A/B (Figure 9a). At t ≈ 17:30 h, 4. Dye Mass Balances and Cross-shore Exchange 260 In sections 4.1-4.3, dye mass balances are shown to close in total, for near-field region 261 A and far-field region B, and for the surfzone and inner-shelf. These results are used in 262 section 4.4 to infer the surfzone to inner-shelf dye tracer exchange. 4.1. Mass Balance: Total and Regions Surfzone and inner-shelf dye masses integrated over the entire alongshore domain (regions A (0 < y = 0 ≥ 248 m) and B (y > 248 m) combined) are A+B MSZ (t) = Z yp (t) 0m D R A F T Z 0m xb Z 0m D(x, y, z, t) dz dx dy, (4) −h March 17, 2015, 12:07pm D R A F T HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE X - 15 and A+B MIS (t) = Z ∞ −∞ Z xb −∞ Z 0m D(x, y, z, t) dz dx dy, (5) −h where yp (t) is the alongshore location of the leading edge of the northward-advecting dye plume (green triangles in Figure 3). Estimation methods for (4) and (5) are described in Appendix A. The total dye mass balance for the surfzone and inner-shelf is A+B MSZ (t) + A+B MIS (t) = Z t Q dτ , (6) t0 263 where Q is the steady dye release rate, t0 = 10:39 h is the dye release start time, and 264 A+B A+B MSZ (t0 ) ≡ MIS (t0 ) ≡ 0 ppb m3 . 265 The total dye mass balance (6) closes well over the 11:21-15:32 h time span of aerial 266 images (Figure 10a, compare red asterisks with line). On average, 88% of the released dye 267 tracer is accounted for using aerial and (relatively sparse) in situ measurements. Note, 268 A+B B MSZ (and therefore MSZ ) may be biased low because yp (t) may be underestimated 269 (Appendix A1). Analyses are broken into time periods I, II, and III based on temporal 270 gaps in aerial data (e.g., Figure 10a) and dye plume evolution (recall section 3.2 and 271 A+B A+B Figure 3). Early in the release during period I, MSZ ≈ MIS (Figure 10a). Starting in 272 A+B period II, as more dye spreads from the surfzone to the inner-shelf, MIS becomes larger 273 A+B A+B A+B . In period III, MIS ≈ 2MSZ (Figure 10a), indicating significant transport than MSZ 274 of the surfzone-released dye to the inner-shelf. 275 The surfzone and inner-shelf are also decomposed into near-field region A and far- 276 A B field region B (Figure 3a). For the surfzone, MSZ ≈ MSZ during period I (Figure 10b), 277 when dye has not advected very far downstream (e.g., Figure 3a). As the dye plume 278 B A advects farther alongshore during period II, MSZ becomes larger than MSZ . Though D R A F T March 17, 2015, 12:07pm D R A F T X - 16 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE 279 dye concentrations are highest near the release (region A) and decrease downstream, the 280 power-law decay is weak (equation (1) and Figure 7), and the larger alongshore extent of 281 B A the dye plume in region B than region A results in period II MSZ ≈ 2MSZ (Figure 10b). 282 B During period III, dye has advected far downstream (e.g., Figures 3e and 3f), and MSZ ≈ 283 A 4MSZ (Figure 10b). Similar trends are observed for the inner-shelf. During period I, 284 A B B A MIS ≈ MIS (Figure 10c). In period II, MIS begins to dominate MIS , and during period 285 B A III, MIS MIS (Figure 10c). 4.2. Mass Balance: Near-field Region A In near-field region A, the total released dye mass ( Rt t0 Q dτ ) must balance the surfzone A A and inner-shelf accumulated dye mass (MSZ + MIS ) and the time-integrated alongshore R t y,A/B y,A/B transport from A to B ( t0 (TSZ + TIS ) dτ ), i.e., A MSZ (t) + A MIS (t) + Z t t0 y,A/B TSZ + y,A/B TIS dτ = Z t Q dτ . (7) t0 286 On average, the sum of the observed mass and cumulative transports account for 76% of 287 the released dye (Figure 11, compare red asterisks with red line). The largest terms of (7) 288 are the cumulative surfzone and inner-shelf alongshore transports, with 289 0.5 290 A A The region A dye masses MSZ and MIS are relatively small, especially during period III 291 (Figure 11, triangles). The 24% of dye unaccounted for in region A is consistent with the 292 low-bias of TIS Rt t0 Q dτ and Rt t0 y,A/B y,A/B TIS dτ ≈ 0.2 Rt t0 Rt t0 y,A/B TSZ dτ ≈ Q dτ (Figure 11, gray and blue curves, respectively). (described in section 3.5). 4.3. Mass Balance: Far-field Region B D R A F T March 17, 2015, 12:07pm D R A F T HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE X - 17 Similar to the near-field region A, the far-field region B accumulated dye mass must balance the A to B alongshore transport, i.e., B MSZ (t) + B MIS (t) = Z t t0 y,A/B y,A/B TSZ + TIS dτ. (8) 293 The observed mass and transport estimates agree well with 18% relative rms error (Fig- 294 ure 12, compare red circles with squares), confirming the consistency among aerial and 295 in situ data and the validation of mass and transport estimation methods. On average, 296 the region B accumulated dye mass (red squares, Figure 12) is slightly larger than the 297 time-integrated alongshore transport (red circles, Figure 12), again consistent with the 298 low bias of TIS y,A/B . 4.4. Exchange Between the Surfzone and Inner-shelf 299 Because the section 4.1-4.3 dye mass balances close, cross-shore surfzone to inner-shelf 300 x,A x,B transport estimates for regions A and B (TSZ/IS and TSZ/IS , respectively) can be inferred 301 from observations. The region A inner-shelf dye mass must balance cross-shore transport 302 input from the region A surfzone and alongshore transport loss to the region B inner-shelf 303 (Figure 13). The region B inner-shelf dye mass must balance cross-shore transport input 304 from the region B surfzone and alongshore transport input from the region A inner-shelf 305 (Figure 13). The corresponding equations are Z t Zt0t t0 D R A F T x,A TSZ/IS dτ = A MIS (t) + x,B B (t) − TSZ/IS dτ = MIS Z t Zt0t t0 y,A/B dτ , (9a) y,A/B dτ . (9b) TIS TIS March 17, 2015, 12:07pm D R A F T X - 18 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE Adding (9a) and (9b), one recovers the expected inner-shelf balance for regions A and B combined: Z t t0 x,A x,B TSZ/IS + TSZ/IS A B dτ = MIS (t) + MIS (t) , (10) 306 stating that the total (A+B) inner-shelf-accumulated dye mass is the time integral of 307 the total cross-shore transport of surfzone-released dye. The time-integrated cross-shore 308 transports (9a), (9b), and (10) are inferred from the observed inner-shelf dye mass and 309 alongshore transport. 310 The inferred time-integrated cross-shore transports Rt t0 x dτ are approximately linTSZ/IS 311 x are estimated ear in each time period, and the associated cross-shore transports TSZ/IS 312 x,A from the slope of each best fit line (Figure 14). The region A cross-shore transport TSZ/IS 313 is similar for periods I and II (137 and 115 ppb m3 s−1 , respectively), consistent with the 314 fixed 248 m alongshore extent of region A, independent of the dye plume advecting far- 315 x,B ther northward with time. In contrast, the region B cross-shore transport TSZ/IS increases 316 significantly among periods I, II, and III (25, 263, and 495 ppb m3 s−1 , respectively) as 317 yp (t) moves northward (e.g., Figures 3 and 4). As a result, the region A and B combined 318 x,A+B cross-shore transport TSZ/IS also increases with time. During period I, when dye has not 319 x,A+B advected far downstream (e.g., Figure 3a), TSZ/IS = 162 ppb m3 s−1 = 0.32Q. During 320 x,A+B period II, when dye has advected farther downstream (e.g., Figures 3b-3d), TSZ/IS = 321 378 ppb m3 s−1 = 0.74Q. During period III, when dye has advected approximately 3 km 322 x,A+B downstream (e.g., Figures 3e and 3f), TSZ/IS = 498 ppb m3 s−1 = 0.97Q, and most of 323 the cross-shore transport occurs in region B (Figure 14, compare red with green best fit 324 slopes). Over the approximate 5 h of aerial observations, roughly 1/2 of the shoreline- D R A F T March 17, 2015, 12:07pm D R A F T X - 19 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE 325 released dye is cross-shore transported to the inner-shelf, i.e., 326 (Figure 14, compare blue symbols with thin red line). Rt t0 x,A+B TSZ/IS dτ ≈ 1 2 Rt t0 Q dτ 5. Discussion 5.1. Parameterizing Cross-Shore Transport Estimates The box-model-based cross-shore tracer flux parameterization used for temperature by Hally-Rosendahl et al. [2014] is tested here for dye with the inferred estimates of surfzone x to inner-shelf cross-shore dye transport (section 4.4). The cross-shore dye flux FˆSZ/IS (units ppb m2 s−1 ) at the surfzone/inner-shelf boundary xb is parameterized as x FˆSZ/IS = hb u∗ ∆D, (11) where hb is the water depth at xb , u∗ is a bulk cross-shore exchange velocity, and ∆D = DSZ − DIS is the surfzone to inner-shelf mean dye concentration difference. Here ∆D is computed separately for region A and region B and for each time period using A A B B period-averaged dye mass estimates M SZ , M IS , M SZ , and M IS (section 4.1) and apA A A ). The surfzone volumes are proximate volumes V of each region (e.g., DSZ = M SZ /VSZ A,B defined by the integration regions for M SZ (see Appendix A1). The inner-shelf volumes are estimated using hdye (Appendix A2), cross-shore width | − 250 m − xb | = 169 m (e.g., Figure 3), and alongshore extents yf = 248 m and y p − yf for regions A and B, respectively, where y p is the mean of yp for each time period . The parameterized surfzone x to inner-shelf cross-shore dye transports TˆSZ/IS (units ppb m3 s−1 ) are then x,A TˆSZ/IS x,B TˆSZ/IS = D R A F T Z yp yf = Z yf 0 x,A FˆSZ/IS x,B FˆSZ/IS dy = dy = Z Z yf A A hb u∗ ∆D dy = hb u∗ ∆D yf , (12a) 0 yp B hb u∗ ∆D dy = hb u∗ ∆D yf March 17, 2015, 12:07pm B y p − yf (12b) D R A F T X - 20 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE for region A and region B, respectively, and x,A+B x,A x,B TˆSZ/IS = TˆSZ/IS + TˆSZ/IS (13) 327 x,A x,B x,A+B for regions A and B combined. The parameterized TˆSZ/IS , TˆSZ/IS , and TˆSZ/IS are each 328 computed for periods I, II, and III. 329 330 x x Parameterized TˆSZ/IS and inferred TSZ/IS are generally similar (Figure 15) with squared correlation r2 = 0.85 and best fit slope 0.7. Minimizing the rms error among parameterized 331 x x TˆSZ/IS and inferred TSZ/IS transports yields the best fit cross-shore exchange velocity 332 u∗ = 1.2 × 10−2 m s−1 . This is consistent with the u∗ = 0.9 × 10−2 m s−1 found using 333 temperature observations on another day with similar wave conditions [Hally-Rosendahl 334 et al., 2014]. 335 x,A The parameterized TˆSZ/IS are comparable between periods I and II and are similar to 336 x,A the inferred TSZ/IS (Figure 15, green circle and square), while the period III parameterized 337 x,A x,A TˆSZ/IS over-estimates the small inferred TSZ/IS (Figure 15, green triangle). The param- 338 x,B eterized TˆSZ/IS increases among periods I, II, and III (Figure 15, vertical coordinates of 339 red symbols) as yp (t) moves farther northward (Figure 4), consistent with the increase 340 x,B of inferred TSZ/IS (Figure 15, horizontal coordinates of red symbols). Similarly, parame- 341 x,A+B x,A+B terized TˆSZ/IS and inferred TSZ/IS are comparable and both increase with time among 342 periods I, II, and III (Figure 15, blue symbols). 343 x,A During period III, the inferred TSZ/IS is small (horizontal coordinate of green triangle 344 in Figure 15), suggesting that most of the region A surfzone dye is transported alongshore 345 to the region B surfzone rather than offshore to the inner-shelf. Consistent with this 346 inference, the mean period III surfzone alongshore transport T SZ y,A/B D R A F T March 17, 2015, 12:07pm = 507 ppb m3 s−1 is D R A F T HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE 347 X - 21 within 1% of the dye release rate Q = 512 ppb m3 s−1 . Combined with small period III d dt A MIS (Figure 10c, triangles), this confirms that the period III cross-shore 348 A MIS and 349 transport within region A is small. 350 The inferred cross-shore exchange velocity u∗ = 1.2 × 10−2 m s−1 incorporates all po- 351 tential surfzone/inner-shelf exchange mechanisms, including both rip currents and Stokes 352 drift driven exchange (onshore near-surface mass flux balanced by undertow at depth). 353 Clark et al. [2010] found the Stokes drift driven exchange mechanism to be negligible for 354 cross-shore dye dispersion within the surfzone due to the surfzone being vertically well- 355 mixed (assumed in Clark et al. [2010], demonstrated here in Figure 8). However, seaward 356 of the surfzone, where dye concentration is not necessarily vertically uniform, the Stokes 357 drift driven exchange mechanism could potentially be important. Here, the u∗ magnitude 358 is compared to the estimated Stokes drift driven exchange velocity offshore of the surfzone. 359 Assuming normally incident, narrow-banded waves with the observed Hs (Figure 2a) and 360 peak period Tp = 13 s, the Lagrangian (Eulerian plus Stokes) velocity profile [e.g., Hally- 361 Rosendahl et al., 2014] at f6 is shoreward in the upper water column and seaward only 362 below z = −1.8 m. Averaging vertically over the seaward portion of the velocity profile 363 yields uLsea = 5.9 × 10−4 m s−1 , which is 20× smaller than the inferred exchange velocity 364 u∗ = 1.2 × 10−2 m s−1 . Observations at x ≈ 2xb (roughly 25 m offshore of f6) show that 365 inner-shelf dye (which has a vertically-mixed surfzone source) is surface-intensified and 366 decreases with depth (Figure 16). This requires seaward dye advection in the upper water 367 column, inconsistent with the Stokes drift driven exchange mechanism (shoreward in the 368 upper water column and seaward at depth). Furthermore, inner-shelf D is distinctively 369 alongshore-patchy (Figures 3 and 13), but Stokes drift driven exchange is expected to be D R A F T March 17, 2015, 12:07pm D R A F T X - 22 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE 370 approximately alongshore-uniform. These discrepancies in magnitude, vertical structure, 371 and alongshore-patchiness suggest that the observed surfzone to inner-shelf cross-shore 372 dye transport is dominated by transient rip ejections on this day with moderate waves, as 373 also concluded using temperature observations from another day with similar wave condi- 374 tions [Hally-Rosendahl et al., 2014]. Farther seaward of 2xb , Stokes drift driven exchange 375 may be important [e.g., Lentz et al., 2008; Suanda and Feddersen, 2015]. 376 x,A The small TSZ/IS inferred for region A during period III and the associated lack of 377 significant transient rip ejections is not surprising, as transient rip currents are sporadic 378 in space and time, and region A is small (< 250 m alongshore) and period III short 379 (< 40 min). The cross-shore transport parameterization (12a) is a bulk quantity that 380 does not resolve the spatial or temporal variability of individual transient rip ejections. 381 Therefore, as ∆DA is comparable among periods I, II, and III, the parameterization (12a) 382 x,A is unable to reproduce the smaller period III inferred TSZ/IS . However, the period III 383 x,A+B parameterized TˆSZ/IS for the entire alongshore domain (regions A and B combined) still 384 x,A+B (Figure 15, blue triangle). agrees well with the inferred TSZ/IS 6. Summary 385 A continuous 6.5 h, near-shoreline release of fluorescent Rhodamine WT dye tracer 386 was observed on 13 October, 2009, at the alongshore-uniform Imperial Beach, California 387 (IB09 experiment). Surfzone and inner-shelf dye concentration was measured in situ with 388 fixed and mobile (jetski- and boat-mounted) fluorometers, and remotely with an aerial 389 multispectral camera system. Waves and currents were measured between the shoreline 390 and roughly 4 m water depth. Dye was advected alongshore by breaking wave- and wind- 391 driven currents, forming a several km-long plume. D R A F T March 17, 2015, 12:07pm D R A F T HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE X - 23 392 Aerial images showed the plume advecting alongshore while ≈ 50 m-wide dye bands 393 were ejected from the surfzone to the inner-shelf by transient rip currents. Near the 394 shoreline, alongshore dye dilution over ≈ 2 km was relatively weak and followed a power- 395 law with exponent −0.3, similar to 29 September dilution over 700 m with exponent −0.2 396 [Hally-Rosendahl et al., 2014]. Both of these dilution rates were smaller than previously 397 observed and modeled for Fickian dilution of surfzone-contained dye plumes (exponent 398 ≈ −0.5, Clark et al. [2010, 2011]). At a cross-shore instrument array 248 m downstream 399 of the release location, ≈ 75% of the released dye was alongshore-advected through the 400 array within 135 m of the shoreline. 401 Using aerial and in situ measurements, 88% of the released dye mass was accounted for 402 over the entire domain (≈ 400 m cross-shore and 3 km alongshore) during a 5 h period after 403 the start of the release. Mass and alongshore transport observations for separate near- 404 and far-field regions were also in agreement (and the small discrepancies were consistent 405 with low-biased inner-shelf alongshore dye transport estimates). The closure of these dye 406 mass balances allowed for inferred estimates of surfzone to inner-shelf cross-shore dye 407 transport, which amounted to roughly 1/2 of the shoreline-released dye during this 5 h 408 period. 409 The cross-shore dye transport was parameterized well using a bulk exchange velocity and 410 surfzone/inner-shelf mean dye concentration difference (inferred and parameterized cross- 411 shore dye transports were correlated (r2 = 0.85) with best fit slope 0.7). The resulting 412 best fit bulk exchange velocity was u∗ = 1.2 × 10−2 m s−1 , consistent with a temperature- 413 derived u∗ = 0.9×10−2 m s−1 from another day of IB09 with similar waves. Together with 414 the observed inner-shelf dye alongshore-patchiness and surface intensification, a 20-fold D R A F T March 17, 2015, 12:07pm D R A F T X - 24 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE 415 difference between u∗ and Stokes drift driven velocity magnitudes indicates that transient 416 rip currents were the dominant surfzone to inner-shelf tracer exchange mechanism at this 417 alongshore-uniform beach during moderate wave conditions. Appendix A: Dye Mass Estimates A1. Surfzone The total surfzone dye mass in regions A and B is defined as A+B MSZ (t) Z = yp (t) 0m Z 0m xb Z 0m D(x, y, z, t) dz dx dy, (A1) −h where h is the water depth, xb is the seaward surfzone boundary, and yp (t) is the leading alongshore edge of the northward-advecting dye plume (Figures 3 and 4). The surfzone dye mass is estimated at times corresponding to the aerial images. The three MSZ integrals (dz, dx, dy) are estimated as follows. Dye is vertically well mixed in the surfzone (Figure 8 and Hally-Rosendahl et al. [2014]), and therefore the vertical integral becomes Z 0m D(x, y, z, t) dz = hD(x, y, t). (A2) −h Cross-shore D profiles do not exist for all y and t. However, cross-shore jetski transects were repeated at various y and are used to compute time-averaged cross-shore dye profiles (see section 3.3 and Figure 5) at alongshore locations near f2 and SA1-SA4, where nearshoreline dye was measured continuously. These mean profiles are used to compute a surfzone dye cross-shore uniformity parameter ξ(y) defined as ξ(y) = R0m xb h(x)D(x, y)dx hSZ LSZ D(x ≈ −10 m, y) , (A3) where hSZ is the mean surfzone water depth, LSZ = |xb | is the surfzone width, and x ≈ −10 m is the location of the shoreward-most observations. By definition 0 < ξ(y) ≤ 1, D R A F T March 17, 2015, 12:07pm D R A F T HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE X - 25 with ξ ≈ 0 corresponding to shoreline-released dye being highly shoreline-concentrated and ξ = 1 corresponding to dye being perfectly cross-shore uniform. Close to the dye release (y = 14 m), the cross-shore uniformity parameter ξ ≈ 0.5. Downstream, as dye mixes across the surfzone, ξ increases to > 0.9 by y = 810 m (Figure 17). The cross-shore integral is estimated as Z 0m h(x) D(x, y, t) dx = hSZ LSZ ξ(y)Dsl (y, t), (A4) xb where Dsl (y, t) is dye measured near the shoreline. Combining (A1)-(A4) yields A+B MSZ (t) = Z yp (t) hSZ LSZ ξ(y)Dsl (y, t) dy. (A5) 0 418 The integral (A5) is alongshore-integrated numerically using the trapezoid rule with 419 Dsl (y, t) at y = 1 m, ySA1 , yf , ySA2 , ySA3 , ySA4 , and yp (t) (Figure 3, yellow symbols 420 and green triangle). Nearest the release, Dsl (y = 1 m, t) = 98 ppb is estimated via the 421 best-fit (1) (see Figure 7) during times that dye is being released (note that all M esti- 422 mates are during this time period). Downstream, Dsl (yp (t), t) is also estimated using the 423 best-fit (1) (see Figures 4 and 7). 424 A+B A B MSZ (t) are decomposed into MSZ (t) (0 < y ≤ yf ) and MSZ (t) (yf < y ≤ yp (t)) 425 using the alongshore boundary yf = 248 m. MSZ estimates are computed for all aerial 426 B image times when in situ near-shoreline dye data are available. Note that MSZ (t) and 427 A+B (t) may be biased low because the inner-shelf yp (t) may be smaller than the actual MSZ 428 alongshore extent of the dye plume within the surfzone (where the alongshore current is 429 fastest (Figure 2b)) and the transient rip ejections from the surfzone to the inner-shelf 430 are episodic in space and time. D R A F T March 17, 2015, 12:07pm D R A F T X - 26 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE A2. Inner-Shelf The total inner-shelf dye mass in regions A and B is defined as A+B MIS (t) = Z Z ∞ Z xb −∞ −∞ 0m D(x, y, z, t) dz dx dy, (A6) −h where h is the water depth, and xb is the surfzone/inner-shelf boundary. Inner-shelf dye mass estimates are calculated using surface dye concentration maps Ds (x, y, t) from the aerial images (e.g., Figure 3) and in situ observations of inner-shelf vertical dye structure (Figure 16) from boat-towed vertical array data (section 2.3.4). These vertical array data resolve inner-shelf D for z = −1 to −3 m (section 2.3.4), thus requiring assumptions for the vertical structure outside this range. As inner-shelf dye comes from the vertically mixed surfzone (Figure 8 and Hally-Rosendahl et al. [2014]), inner-shelf D(x, y, z, t) is assumed vertically uniform in the upper 1 m. For z < −3 m, the best fit of mean D(z) (Figure 16) is extrapolated to the depth where it would vanish. This structure is then vertically integrated, and hdye is computed as the depth that yields an equivalent vertical integral hdye Ds (x, y, t). The inner-shelf vertical dye integral is thus estimated as Z 0 D(x, y, z, t) dz = hdye Ds (x, y, t), (A7) −h where h is the water depth, Ds (x, y, t) is the aerial-measured surface dye concentration, and hdye = min (2.67 m, h). The inner-shelf dye mass estimates are then A+B MIS (t) = Z ∞ −∞ Z xb hdye Ds (x, y, t) dx dy, (A8) −∞ 431 A+B integrated using the trapezoid rule in each lateral direction. MIS (t) are decomposed 432 A B into MIS (t) (−∞ < y ≤ yf ) and MIS (t) (yf < y < ∞) using the alongshore boundary 433 yf = 248 m. D R A F T March 17, 2015, 12:07pm D R A F T HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE X - 27 434 Acknowledgments. 435 IB09 field work and analysis was funded by ONR, NSF, and California Sea Grant. 436 K. Hally-Rosendahl was supported by the National Science Foundation Graduate Re- 437 search Fellowship under Grant No. DGE1144086 and California Sea Grant under Project 438 No. R/CONT-207TR. Staff and students from the Integrative Oceanography Division 439 (B. Woodward, B. Boyd, K. Smith, D. Darnell, R. Grenzeback, A. Gale, M. Spydell, 440 M. Omand, M. Yates, M. Rippy, A. Doria) were instrumental in acquiring the field obser- 441 vations. K. Millikan, D. Ortiz-Suslow, M. Fehlberg, and E. Drury provided field assistance. 442 Imperial Beach lifeguards, supervised by Captain R. Stabenow, helped maintain public 443 safety. The YMCA Surf Camp management generously allowed extensive use of their 444 facility for staging and recuperation. The U. S. Navy provided access to Naval property 445 for data collection. Dr. M. Okihiro coordinated permits and logistics. We thank these 446 people and organizations. This work comprises a portion of K. Hally-Rosendahl’s PhD 447 thesis. As subsequent thesis chapters will use these same IB09 data, they are not shared 448 with this work. However, upon completion of K. Hally-Rosendahl’s thesis, IB09 data will 449 be made available at http://falk.ucsd.edu in accordance with the AGU data policy. References 450 Clark, D. B., F. Feddersen, M. M. Omand, and R. T. Guza (2009), Measuring fluo- 451 rescent dye in the bubbly and sediment laden surfzone, Water Air Soil Poll., doi: 452 10.1007/s11270-009-0030-z. 453 454 Clark, D. B., F. Feddersen, and R. T. Guza (2010), Cross-shore surfzone tracer dispersion in an alongshore current, J. Geophys. Res., 115 (C10035), doi:10.1029/2009JC005683. D R A F T March 17, 2015, 12:07pm D R A F T X - 28 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE 455 Clark, D. B., F. Feddersen, and R. T. Guza (2011), Modeling surfzone tracer plumes, 456 part 2: Transport and dispersion, J. Geophys. Res., Oceans, 116 (C11028), doi: 457 10.1029/2011JC007211. 458 459 Clark, D. B., S. Elgar, and B. Raubenheimer (2012), Vorticity generation by short-crested wave breaking, Geophys. Res. Lett., 39, doi:10.1028/2012GL054034. 460 Clark, D. B., L. Lenain, F. Feddersen, E. Boss, and R. T. Guza (2014), Aerial imaging 461 of fluorescent dye in the nearshore, J. Atmos. and Ocean. Tech., 31, 1410–1421, doi: 462 10.1175/JTECH-D-13-00230.1. 463 464 465 466 467 468 469 470 Dorfman, M., and A. Haren (2013), Testing the Waters: A guide to water quality at vacation beaches, Tech. rep., National Resrouces Defense Council. Dorfman, M., and N. Stoner (2012), Testing the Waters: A guide to water quality at vacation beaches, Tech. rep., National Resrouces Defense Council. Feddersen, F. (2012), Scaling surfzone dissipation, Geophys. Res. Lett., 39 (L18613), doi: 10.1029/2012GL052970. Feddersen, F. (2014), The generation of surfzone eddies in a strong alongshore current, J. Phys. Ocean., 44, 600–617, doi:10.1175/JPO-D-13-051.1. 471 Feddersen, F., D. B. Clark, and R. T. Guza (2011), Boussinesq modeling of surfzone tracer 472 plumes, Part 1: Eulerian wave and current comparisons, J. Geophys. Res., 116 (C11027), 473 doi:10.1029/2011JC007210. 474 Grant, S. B., J. H. Kim, B. H. Jones, S. A. Jenkins, J. Wasyl, and C. Cudaback (2005), Surf 475 zone entrainment, along-shore transport, and human health implications of pollution 476 from tidal outlets, J. Geophys. Res., 110 (C10025), doi:10.1029/2004JC002,401. D R A F T March 17, 2015, 12:07pm D R A F T HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE X - 29 477 Hally-Rosendahl, K., F. Feddersen, and R. T. Guza (2014), Cross-shore tracer exchange 478 between the surfzone and inner-shelf, Journal of Geophysical Research: Oceans, 119 (7), 479 4367–4388, doi:10.1002/2013JC009722. 480 Halpern, B. S., S. Walbridge, K. A. Selkoe, C. V. Kappel, F. Micheli, C. D’Agrosa, J. F. 481 Bruno, K. S. Casey, C. Ebert, H. E. Fox, R. Fujita, D. Heinemann, H. S. Lenihan, 482 E. M. P. Madin, M. T. Perry, E. R. Selig, M. Spalding, R. Steneck, and R. Watson 483 (2008), A global map of human impact on marine ecosystems, SCIENCE, 319 (5865), 484 948–952, doi:10.1126/science.1149345. 485 486 487 488 Harris, T. F. W., J. M. Jordaan, W. R. McMurray, C. J. Verwey, and F. P. Anderson (1963), Mixing in the surf zone, Int. J. Air Wat. Polut., 7, 649–667. Inman, D. L., R. J. Tait, and C. E. Nordstrom (1971), Mixing in the surf zone, J. Geophys. Res., 26, 3493–3514. 489 Koh, R., and N. Brooks (1975), Fluid mechanics of waste-water disposal in the ocean, in 490 Annual review of fluid mechanics. Vol.7, edited by M. Van Dyke and G. Vincenti, pp. 491 187–211. 492 Kumar, N., F. Feddersen, J. C. McWilliams, Y. Uchiyama, and W. O’Reilly (2014), Mid- 493 shelf to surf zone coupled ROMS-SWAN model-data comparison of waves, currents, and 494 temperature: Diagnosis of subtidal forcings and response, J. Phys. Ocean., revised. 495 496 Lentz, S. J., and M. R. Fewings (2012), The Wind- and Wave-Driven Inner-Shelf Circulation, 4, 317–343, doi:10.1146/annurev-marine-120709-142745. 497 Lentz, S. J., M. Fewings, P. Howd, J. Fredericks, and K. Hathaway (2008), Observations 498 and a Model of Undertow over the Inner Continental Shelf, J. Phys. Ocean., 38 (11), 499 2341–2357, doi:10.1175/2008JPO3986.1. D R A F T March 17, 2015, 12:07pm D R A F T X - 30 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE 500 Lucas, A., P. Franks, and C. Dupont (2011), Horizontal internal-tide fluxes support 501 elevated phytoplankton productivity over the inner continental shelf, Limnology and 502 Oceangraphy: Fluids and Environment, 1, 56–74, doi:10.1215/21573698-1258185. 503 Peregrine, D. H. (1998), Surf zone currents, Theor. Comput. Fluid Dyn., 10, 295–309. 504 Schiff, K. C., M. J. Allen, E. Y. Zeng, and S. M. Bay (2000), Southern California, Mar. 505 Polut. Bull., 41, 76–93. 506 Sinnett, G., and F. Feddersen (2014), The surf zone heat budget: The effect of wave 507 heating, Geophysical Research Letters, 41 (20), 7217–7226, doi:10.1002/2014GL061398. 508 Smart, P. L., and I. M. S. Laidlaw (1977), An evaluation of some fluorescent dyes for 509 510 511 water tracing, Water Resources Research, 13 (1), 15–33. Spydell, M. S., and F. Feddersen (2009), Lagrangian drifter dispersion in the surf zone: Directionally spread, normally incident waves, J. Phys. Ocean., 39, 809–830. 512 Spydell, M. S., F. Feddersen, J. H. MacMahan, and R. T. Guza (2014), Relating La- 513 grangian and Eulerian horizontal eddy statistics in the surf zone, J. Geophys. Res., 119, 514 doi:10.1002/2013JC009415, in press. 515 516 517 518 Suanda, S. H., and F. Feddersen (2015), Self-similar scaling for cross-shelf exchange by transient rip currents, submitted. Thornton, E. B., and R. T. Guza (1986), Surf zone longshore currents and randome waves: Field data and models, J. Phys. Ocean., 16 (7), 1165–1178. D R A F T March 17, 2015, 12:07pm D R A F T X - 31 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE S A4 boat 1600 1400 1200 S A3 800 b e ac h y ( m) 1000 je t s ki 600 S A2 400 f6 f4 f2 200 S A1 0 −250 Figure 1. −200 −150 −100 x ( m) −50 0 Planview of IB09 bathymetry contours versus cross-shore coordinate x and alongshore coordinate y. Star indicates dye release location. Diamonds denote the crossshore array of bottom-mounted instrument frames f1-f6 (onshore to offshore). Circles indicate SA1-SA4 fluorometer locations. Vertical dashed line represents an idealized boat alongshore transect driven repeatedly near this cross-shore location. Horizontal dashed line represents an idealized jetski cross-shore surface transect driven repeatedly at various alongshore locations. D R A F T March 17, 2015, 12:07pm D R A F T X - 32 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE H s ( m) 1 ( a) 0.5 V ( m s − 1) 0 0.4 xb (b) 0.2 0 z ( m) 0 xb (c) f3 −2 −4 −150 f2 f1 f4 f5 f6 −100 xb −50 x ( m) 0 Figure 2. Time-averaged (11:00-16:00 h) (a) significant wave height Hs , (b) alongshore current V , and (c) vertical locations of f1-f6 versus cross-shore coordinate x. In (c), the black curve gives the bathymetry h(x). The mean seaward surfzone boundary xb = −81 m is defined as the location of maximum Hs . D R A F T March 17, 2015, 12:07pm D R A F T X - 33 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE Figure 3. Aerial multispectral images of surface dye concentration (see colorbar) versus cross-shore coordinate x and alongshore coordinate y for six times (indicated in each panel). The mean shoreline is at x = 0 m. Green star indicates location of continuous dye release (starting at t0 = 10:39 h). Yellow diamonds indicate cross-shore array f1-f6 locations, and yellow circles indicate SA1-SA4 locations. Light gray indicates regions outside the imaged area, and black indicates unresolved regions due to foam from wave breaking. Vertical dashed cyan line at xb divides the surfzone (SZ) and inner-shelf (IS), and horizontal cyan line divides the near- and far-field regions A and B (see panel (a)). Plume leading edge yp (t) is shown with green triangles at x ≈ −100 m (for panel (f), yp ≈ 3250 m). Panels (a), (b)-(d), and (e)-(f) are in time periods I, II, and III, respectively. D R A F T March 17, 2015, 12:07pm D R A F T X - 34 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE 3500 3000 y p ( m) 2500 2000 1500 1000 500 I 0 11 Figure 4. II 12 III 13 14 t (h) 15 16 17 Alongshore coordinate of dye plume leading edge yp versus time. The determination of yp (t) is described in section 3.2. Magenta bar indicates duration of nearshoreline, continuous dye release at y = 0 m. Blacks bars denote time periods I, II, and III. 100 D (ppb) 80 y y y y y = = = = = 14 m 87 m 207 m 619 m 810 m 60 40 20 0 −80 Figure 5. −60 −40 x ( m) −20 0 Time-averaged, cross- and alongshore-binned surfzone D from jetski surface transects versus cross-shore coordinate x (see legend for alongshore locations y). D R A F T March 17, 2015, 12:07pm D R A F T X - 35 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE 200 ( a) S A1, y = 82 m D (ppb) 150 100 50 0 80 ( b ) f 2, y = 248 m D (ppb) 60 40 20 0 50 ( c ) S A2, y = 546 m D (ppb) 40 30 20 10 0 20 ( d ) S A3, y = 1069 m D (ppb) 15 10 5 0 20 ( e ) S A4, y = 1662 m D (ppb) 15 10 5 0 10 Figure 6. 11 12 13 14 15 t (h) 16 17 18 19 20 Dye concentration D versus time at the near-shoreline f2 and SA1-SA4 (diamond and circles, respectively, in Figures 1 and 3). Alongshore location is indicated in each panel. Magenta bars indicate duration of near-shoreline, continuous dye release at y = 0 m. Vertical axes differ, and each time series begins at the instrument deployment time (not dye arrival time). D R A F T March 17, 2015, 12:07pm D R A F T X - 36 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE 2 D (ppb) 10 1 10 0 10 0 10 Figure 7. 2 10 y ( m) 4 10 Mean (time-averaged) dye concentration versus alongshore coordinate y at the near-shoreline f2 and SA1-SA4 (yellow diamond and circles in Figure 3, respectively). Vertical bars are standard deviations about the means. Best fit line (dashed) is D = D0 (y/y0 )α , where y0 = 1 m is chosen for simplicity, and best fit constants are D0 = 98 ppb and α = −0.3. D R A F T March 17, 2015, 12:07pm D R A F T X - 37 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE 50 D (ppb) f 4, 1. 3 mab 40 f 4, 0. 7 mab 30 f 4, 0. 2 mab 20 10 0 10 Figure 8. 11 12 13 14 t (h) 15 16 17 18 Dye concentration D versus time at three ETs with different vertical ele- vations (mab is meters above bottom) on f4 at the seaward surfzone boundary xb (see legend and Figure 2c). Magenta bar indicates duration of near-shoreline, continuous dye release (Figure 1, star) 248 m south of the cross-shore array (Figure 1, diamonds). Gaps in the time series result from sampling for 51 minutes of each hour. D R A F T March 17, 2015, 12:07pm D R A F T X - 38 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE TSZ A/B , y ( p p b m 3 s − 1) 4000 SZ ( a) IS (b) SZ IS (c) 3000 2000 1000 ( p p b m 3 s − 1) 0 800 600 TI S A/B , y 400 200 6 4 2 Rt t0 T A/B , y 3 d τ (ppb m ) 6 0 x 10 8 0 10 Figure 9. 11 12 13 14 t (h) 15 16 17 18 Time series of alongshore dye transport from region A to B in the (a) y,A/B surfzone (TSZ y,A/B defined in (2)) and (b) inner-shelf (TIS defined in (3)). Vertical scales differ. (c) Time series of cumulative (time-integrated) surfzone and inner-shelf alongshore dye transports (see legend). Magenta bars indicate duration of near-shoreline, continuous dye release (Figure 3a, star) 248 m south of the cross-shore array (Figure 3a, diamonds) dividing regions A and B . The dye release rate Q = 512 ppb m3 s−1 , and the total dye released is 1.19 × 107 ppb m3 . D R A F T March 17, 2015, 12:07pm D R A F T X - 39 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE 6 x 10 10 Rt II III t0 Q d τ MSZ + MI S MSZ MI S 8 M ( p p b m 3) I 6 ( a) 4 2 6 0 x 10 10 (b) Rt (c) t0 Q d τ M SAZ + M SBZ M SAZ M SBZ 8 M ( p p b m 3) Rt 6 4 2 6 0 x 10 10 M ( p p b m 3) 8 6 t0 Q d τ M IAS + M IBS M IAS M IBS 4 2 0 11 Figure 10. 11.5 12 12.5 13 13.5 t (h) 14 14.5 15 15.5 16 Dye mass M versus time. Black bars denote time periods I, II, and III. (a) Surfzone estimates MSZ (gray) are from in situ observations, and inner-shelf estimates MIS (blue) are from aerial images. Red asterisks are MSZ +MIS . Red line shows the timeRt integrated dye mass released since t0 =10:39 h ( t0 Q dτ , where Q is the dye release rate). Estimation methods for MSZ and MIS are described in Appendix A1 and A2, respectively. (b) Surfzone dye mass MSZ versus time for the near-field region A (y ≤ 248 m, triangles) A+B and the far-field region B (y > 248 m, squares). Solid gray diamonds are MSZ . (c) Inner-shelf dye mass MIS versus time for region A (triangles) and region B (squares). A+B Solid blue circles are MIS . D R A F T March 17, 2015, 12:07pm D R A F T X - 40 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE 6 x 10 10 Rt II III t0 Q d τ 2 R t 1 y , A/B y , A/B d τ + M SAZ ( t ) + M IAS ( t ) TSZ + TI S t0 Rt y , A/B TSZ dτ Rtt0 y , A/B dτ t0 T I S A MSZ (t) M IAS ( t ) 8 M ( p p b m 3) I 6 4 2 0 11 11.5 12 12.5 13 13.5 t (h) 14 14.5 15 15.5 16 Dye mass balance terms versus time for near-field region A (0 < y ≤ Figure 11. yf = 248 m). See legend and equation (7). Red line shows the time-integrated dye mass released since t0 =10:39 h. 6 x 10 10 I Rt 1 t0 M ( p p b m 3) 8 6 4 II y , A/B y , A/B + TI S TSZ M SBZ ( t ) + M IBS Rt y , A/B dτ t0 T S Z Rt y , A/B dτ t0 T I S B MSZ (t) M IBS ( t ) III 2 dτ (t) 2 0 11 11.5 12 12.5 13 13.5 t (h) 14 14.5 15 15.5 16 Figure 12. Dye mass balance terms versus time for far-field region B (y > yf = 248 m). See legend and equation (8). D R A F T March 17, 2015, 12:07pm D R A F T X - 41 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE SZ IS B MIS R B R A x,B TSZ/IS d⌧ y,A/B TIS R A MIS d⌧ x,A TSZ/IS d⌧ R Q d⌧ Figure 13. Planview photograph and superposed schematic of dye mass balances (9a), (9b), and (10). Star denotes location of dye released at steady rate Q. D R A F T March 17, 2015, 12:07pm D R A F T X - 42 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE 6 x 10 Rt x t 0 T S Z /I S d τ ( p p b m3 ) 10 8 6 Rt I II III Q dτ 2 Rtt0 1 x , A x, B T + T S Z /I S S Z /I S d τ t R t0 x , A T dτ Rtt0 SxZ, B/I S t 0 T S Z /I S d τ 4 2 0 11 11.5 12 12.5 13 13.5 t (h) 14 14.5 15 15.5 16 Figure 14. Time series of cumulative (time-integrated) cross-shore dye transports from the surfzone to inner-shelf (circles) inferred from inner-shelf dye mass observations MIS y,A/B and alongshore transport measurements TIS . See (9a), (9b), and (10). Line segments are least squares fits for each time period, and line segment slopes yield inferred crossx . Thin red line shows the time-integrated dye mass released shore dye transports TSZ/IS since t0 = 10:39 h. D R A F T March 17, 2015, 12:07pm D R A F T X - 43 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE 500 A+B A B TˆSxZ /I S ( p p b m 3 s − 1 ) 400 I II III 300 200 100 0 0 Figure 15. 100 200 300 400 T SxZ /I S ( p p b m 3 s − 1 ) 500 x versus inferred cross-shore Parameterized cross-shore dye transport TˆSZ/IS x dye transport TSZ/IS for regions A, B, and A+B during periods I, II, and III (see legend). x x are inferred from follow (12a), (12b), and (13). The TSZ/IS The parameterized TˆSZ/IS aerial and in situ observations (see (9a), (9b), and (10)). D R A F T March 17, 2015, 12:07pm D R A F T X - 44 HALLY-ROSENDAHL ET AL.: SURFZONE TO INNER-SHELF TRACER EXCHANGE −1 z ( m) −1.5 −2 −2.5 −3 0 5 D (ppb) 10 Figure 16. Mean inner-shelf dye concentration D versus vertical coordinate z from the boat-towed vertical array for data within inner-shelf dye patches (D(x, y, z = −1 m, t) ≥ 2 ppb). Dashed curves indicate standard deviations about the mean. 1 0.8 ξ ( ) 0.6 0.4 0.2 0 0 Figure 17. 200 400 600 y (m ) 800 1000 1200 Surfzone dye cross-shore uniformity parameter ξ versus alongshore coordi- nate y (ξ is defined in (A3) and described in Appendix A1). D R A F T March 17, 2015, 12:07pm D R A F T
© Copyright 2024