A single MATLAB binary file: CSIDE_SZplume.mat, contains data referenced in manuscript Grimes et al. 2021, Long-distance/time surf-zone tracer evolution affected by inner-shelf tracer retention and recirculation: The MATLAB datenum() time of the dye release spans, 7.362452213425926e5 <= tr <= 7.362453811805556e5 08-Oct-2015 05:18:44 <= tr <= 08-Oct-2015 09:08:54 herafter, time coordinates are hours relative 05:18:44 on 10/08/2015 ------------ NDWI, Unit-Normalized Total Radiance, DYE (Figure 2) -------------- Example Dye, N, R from ~6.8h after dye release start, A: [501x1201 double], alongshore coordinates [m] C: [501x1201 double], cross-shore coordinates [m] DYE: [501x1201 double], dye concentration [ppb] (estimate_DYE.m; eq.1-2) NDWI: [501x1201 double], normalized difference water index (estimate_NDWI.m; eq.3) RAD: [501x1201 double], unit-normalized total radiance (eq.4) ------------ NDWI, Unit-Normalized Total Radiance, DYE -------------- ------------ Experiment Conditions (Figure 3) ----------------------- Tijuana River National Estuary Research Reserve Meteorologial Station (MET) Winds: t_wind: [60x1 double], time [hours since release] v_wind: [60x1 double], alongshore speed u_wind: [60x1 double], cross-shore speed Imperial Beach-West CDIP wave buoy #218: t_wave: [1x667 double], time [hours since release] Hs: [668x1 double], significant wave height [m] Dp: [668x1 double], peak wave direction [deg] Lwave: [668x1 single], deep water wave length SxyONr: [1x668 single], S_xy/rho [m^3/s^2] San Diego Bay NOAA tide station (9410170) t_eta: [43200x1 double] eta: [43200x1 double] Water level estimates based on MASS shoreline evolution, x_sl(t,y): t_MASS: [1x26 double], North-bound MASS time vector [h] eta_MASS: [1x26 double], water-level estimate Surfzone width stats based on normalized difference water index and unit based total radiance derived x_sl(t,y) & x_sz(t,y), respectively: Lsz_MASS_yavg: [1x26 double], y-averaged surf-zone width Lsz_MASS_ystd: [1x26 double], y-standard deviation Depth and 30-minute averaged WireWalker 1 adjacent Aqquadopp currents tWW1: [1x1440 double], time [h] uWW1: [1x1440 double], cross-shore velocity [m/s] vWW1: [1x1440 double], alongshore vel [m/s] Depth and 30-minute averaged RDI WorkHorse currents t12N: [1x16 double], time [h] u12N: [16x1 double], cross-shore vel [m/s] v12N: [16x1 double], alongshore vel [m/s] ------------ Experiment Conditions (Figure 3) ----------------------- ------------ Qualitative Dye Evolution (Figure 4) ----------------------- MASS derived surface dye concentration at 5.1, 6.9, 9.95, & 11.7 h after release start: A: [501x1201 double] C: [501x1201 double] DYEvsTIME: [3-D double], DYE(A,C,t) [ppb] SZmask: [3-D double], mask between x_sl and x_sz ------------ Qualitative Dye Evolution (Figure 4) ----------------------- ------------ Inner-shelf Alongshore Transect (Figure 5) ----------------------- Afternoon inner-shelf alongshore (north/south bound) D/T transect data: Northbound transect: tTnorth: [1x210 double], time [h since release start] aTnorth: [1x210 double], alongshore coord cTnorth: [1x210 double], cross-shore coord Tnorth: [22x210 double], Temperature [C] Dnorth: [22x210 double], Dye [ppb] Northbound transect: tTsouth: [1x195 double], same as above aTsouth: [1x195 double] cTsouth: [1x195 double] Tsouth: [22x195 double] Dsouth: [22x195 double] ------------ Inner-shelf Alongshore Transect (Figure 5) ----------------------- ------------ Surfzone Dye Observations (Figure 6) ----------------------- tsz: [145x1 double], time [h] Asz: [6x1 double], alongshore coord [m] Dsz: [145x6 double], 30 min avg. dye [ppb] Dsz_pm: [3-D double], +/- 30 min std. dev. [ppb] szDmax: [1x6 double], maximum D observed [ppb] szDmax_pm: [6x2 double], szDmax +/- std. dev. [ppb] tSZplume: [6x2 double], arrival time tf and passage tp tSZplume_m5: [6x2 double], ... 5% decrease in threshold tSZplume_p5: [6x2 double], ... 5% increase in threshold ------------ Surfzone Dye Observations (Figure 6) ----------------------- ------------ Surfzone Plume Alongshore Transport, Dilution, Widening (Figure 7) ----------------------- tSxy: [1x24 double], Sxy time vector [h], aSxy: [24x1 single], alongshore front based on Sxy Lsz_mu: [1x1 single], coefficient (Lsz*mu) with Lsz~100m (eq.6) t_MASS: [1x26 double], MASS time vector aISplume: [26x1 double], alongshore front based on MASS inner-shelf dye, aISplume_pm1: [26x2 double], ... +/- 1% change in threshold szDmax: [1x6 double], maximum D observed [ppb] szDmax_pm: [6x2 double], szDmax +/- std. dev. [ppb] tSZplume: [6x2 double], arrival time tf and passage tp tSZplume_m5: [6x2 double], ... 5% decrease in threshold tSZplume_p5: [6x2 double], ... 5% increase in threshold ------------ Surfzone Plume Alongshore Transport, Dilution, Widening (Figure 7) ----------------------- ------------ Coupled SZ/IS Tracer Model Optimization ----------------------- Model constants: Kyy: [1x1 double], surfzone y-dispersion coefficient Lsz: [1x1 double], surfzone width Hsz: [1x1 double], outer surfzone depth, slope 0.002 Vis: [1x1 double], inner-shelf bin advection speed Vsz: [1x1 double], surf-zone bin advection speed Parameter ranges: Ksz0: [1x9 double], 2-way SZ/IS exchange rate Kis0: [1x8 double], 1-way IS exchange rate Lis0: [1x13 double], width of IS bin error: [3-D double], resulting error matrix Optimized: Kis: [1x1 double], minimum error Kis_uncertainty: [1x1 double], 1D curvature of error at minimum Ksz: [1x1 double], ... Ksz_uncertainty: [1x1 double] Kyy: [1x1 double] gamma: [1x1 double], Asz/Ais with min error Lis: [1x1 double], Lis for gamma, assuming well mixed Ais Lis_uncertainty: [1x1 double], gamma and Lis uncertainty see: run_SZ_IS_box_model_v2.m & solve_SZ_IS_dye_adv_diff_v2.m for numerics ------------ Coupled SZ/IS Tracer Model Optimization ----------------------- ------------ Surfzone Model/Obs time series (Figure 9) -------------------- tsz_shifted: [145x6 double], shifted time (t-vsz/y) Dsz: [145x6 double], 30 min avg. dye [ppb] Dsz_pm: [3-D double], +/- 30 min std. dev. [ppb] tmod_shifted: [1860x6 double], model shifted time vector DmodSZ: [1860x6 double], model SZ concentration DmodIS: [1860x6 double], model IS concentration ------------ Surfzone Model/Obs time series (Figure 9) -------------------- ------------ Surfzone Model/Obs Dmax and Width (Figure 10) -------------------- Asz: [6x1 double], alongshore coord [m] szDmax: [1x6 double], maximum D observed [ppb] szDmax_pm: [6x2 double], szDmax +/- std. dev. [ppb] tSZplume: [6x2 double], arrival time tf and passage tp tSZplume_m5: [6x2 double], ... 5% decrease in threshold tSZplume_p5: [6x2 double], ... 5% increase in threshold modDmax: [1x6 double], modeled maximum D [ppb] tMODplume: [6x2 double], tf and tp for model tMODplume_m5: [6x2 double], ... with 5% decrease in threshold tMODplume_p5: [6x2 double], ... with 5% increase in threshold ------------ Surfzone Model/Obs Dmax and Width (Figure 10) -------------------- Ancillary Data: ----------- NOAA DEM Bathymetry --------------- MATLAB Binary File Name: CSIDEbathy_smooth.mat 100mx100m smoothed 2012 NOAA Tsunami DEM: https://www.ncei.noaa.gov/metadata/geoportal/rest/metadata/item/gov.noaa.ngdc.mgg.dem: 3542/html# variables= xi: [1501x1501 double], Imperial Beach x-coordinate; see lltoxy_imperialbeach.m yi: [1501x1501 double], y-coord Zi: [1501x1501 double], vertical coordinate (elevation>0) Zi_smooth: [1501x1501 double], 100m x 100m smoothed Zi_noJ_smooth: [1501x1501 double], same but smoothed after removing sediment retention groins ----------- NOAA DEM Bathymetry --------------- ----------- Cross-shore Dye Correction --------------- To determine and account for the foam/turbidity induced remotely sensed dye errors cross-shore transects of remotely sensed D were extracted from regions known to be dye free. The mean or background cross-shore profile was removed before D analysis: data was exracted from the following array of pass/indices outside dye plume [pass, min(a), max(a)]; input = [1 , -1e3 , -500;... 2 , -700 , -650;... 3 , 4200 , 4300;... 4 , -1e3 , -500;... artifact in mapped D 5 , 4200 , 4400;... 6 , 5e3 , 5500;... 7 , -1e3 , -500;... 8 , 4800 , 5200;... artifact in mapped D 9 , -1e3 , -500;... 12 , 4750 , 5250;... 13 , -1750 ,-1500;... 14 , 5100 , 5300;... 18 , 5200 , 5400;... 21 , 5800 , 6100;... 22 , 5500 , 5600;... 28 , 5400 , 5600;... 32 ,-2500 ,-2250;... ]; MATLAB Binary file name: background_xshore_DYE.mat Variables= Dbar: [501x1 double], mean cross-shore D profile (from "input" above) Nbar: [501x1 double], mean Normalized difference water index profile Rbar: [501x1 double], mean Unit-normalized total radiance profile Rmax: [1x1 double], maximum R Rmin: [1x1 double], minimum R c: [501x1 double], cross-shore coordinate csbar: [1x1 double], mean shoreline location (based on N) cszbar: [1x1 double], mean outer surf-zone location (based on R) ----------- Cross-shore Dye Correction --------------- ----------- Scripts --------------- Also included in the data repository are a few scripts used in manuscript analyses/figures. File Name: CSIDE_SZplume_scripts map_xy_to_ac_CSIDE_v2.m: USAGE: [A,C] = map_xy_to_ac_CSIDE_v2(X,Y); Calls xytoac_v2.m which converts MASS spectral data from Imperial Beach (X,Y) coordinates to MSL quasi-shorenormal coordinates (\tilde{x},\tilde{y}). A = \tilde{y} is the alongshore coordinate and C = \tilde{x} is the local shorenormal coordinate. See also: lltoxy_imperialbeach for lat/lon to (x,y) transformation, or xytoll_imperialbeach.m for (x,y) to lat/lon transformation grid_RAD_CSIDE.m: USAGE: [rad,a,c] = grid_RAD_CSIDE(dat,PASS,A,C,alims,da,dc); Interpolates push-broom scattered pixel data to regular grid within alims = [Amin, Amax, Cmin, Cmax]; default alongshore spacing da=10m, and cross-shore spacing dc=2m. estimate_DYE.m USAGE: [DYE] = estimate_DYE(dat,PAR,xcorrect,name,rblock,cblock); Estimates surface dye concentration from spectral radiance measurements; PAR has MASS header info and dat has radiance data in the form of a matlab.io file or in structures w/ fields: PAR.Wavelength = center of spectral bins; PAR.Wavelencth_width= bin width of band dat.rad = spectral radiance xcorrect = logical; 1 == apply cross-shore correction, *need to specify path to file* name = string if 'rad' is not default variable rblock = row-vector indices on which to estimate D cblock = column-vector indices ... list of scripts to include: xytoll_imperialbeach.m lltoxy_imperialbeach.m xytoac_v2.m map_xy_to_ac_CSIDE_v2.m grid_RAD_CSIDE.m estimate_DYE.m estimate_NDWI.m run_SZ_IS_box_model_v2.m solve_SZ_IS_dye_adv_diff_v2.m