1.1Introduction The morphology of sandy coasts, including the shoreline position and the bathymetry of the surf and shoaling zones, quite often displays complex and intriguing patterns. These patterns are sometimes nearly periodic alongshoreoratleastshowingsomesortofregularitywithanalongshorerecurrencelength,L,andarethenknownas rhythmiccoastalmorphologies .Sometypeshavebeendefinedintheliterature(Short,1999;Ribasetal.,2015)but theextremecomplexityofbeachdynamicsandtheincreasingcapacityandfrequencyofbeachmonitoringandfieldobservationsoftenchallengetheirtraditionalclassification(Guillénetal.,2017). Transversesandbarsystemsisonetypeofthosepatterns.Theterm transversebar (TB)isgenericallyappliedtosandbarsextendingperpendicularlytothecoastorwithanobliqueorientation(Shepard,1952).Theyusually occurinpatchesofafewofthemuptotens,theyareseparatedbytroughsandtheyaretypicallyattachedtotheshore.Thealongshorespacing,L,isdefinedasthedistancebetweensuccessivebarcrests.Withadominantlongshore currenttheytendtomigratedowndriftwithmigrationratesupto40 m/d(RibasandKroon,2007;Pellónetal.,2014).Theysometimesshowanasymmetryoftheircross-section(thedown-currentflankbeingsteeperthantheupcurrentflank,Pellónetal.,2014). Severaltypesoftransversebarshavebeenreportedintheliterature(see,e.g.,Pellónetal.,2014;Ribasetal.,2015).Thefirstonesare TBRbars ,whichareassociatedtotheTransverseBarandRip(TBR)stateinthestandard beachstateclassification(WrightandShort,1984).Theyaretypicallywideandshort-crestedandtheiroriginisthemergingofthehornsofacrescenticbarintothebeach.Thesecondtypeare mediumenergyfingerbars ,MEFB, whicharesometimesobservedinopenmicrotidalbeachesundermedium-energyconditions(KonickiandHolman,2000;Khabidov,2001;RibasandKroon,2007)andtheyalwayscoexistwithshore-parallel(orcrescentic)bars.They arethinandelongatedincontrastwiththewiderandshorterTBRbars.MEFBareephemeral(residencetimefrom1 dayto1 month),attachedtothelow-tideshorelineor,occasionally,totheshore-parallelbar.Theyarelinkedtothe Rhythmicmorphologyinamicrotidallow-energybeach AnnaMujal-Colilles1a,⁎
[email protected] ManelGrifoll1a AlbertFalqués2b 1aLaboratorid’'EnginyeriaMarítima,DepartmentofCivilandEnvironmentalEngineering,UPC-BarcelonaTech,C/JordiGirona,1-3,08034Barcelona,Spain 2bDepartmentofPhysics,UPC-BarcelonaTech,C/JordiGirona,1-3,08034Barcelona ⁎Correspondingauthor. Abstract ObservationsofrythmicfeaturesalongtheinnersideoftheTrabucadorbarrierbeacharecoupledtotwonumericalmodelstounravelthemechanismsofitsformation.TheTrabucadorisalong(6Km)narrow(125 m) barrierandmicrotidalbeachattheSWsideoftheEbrodelta(Catalonia).Itsinnersideisalowenergybeachwithasandyshallowterracefeaturinganintricatealongshorerhythmicmorphology.Sixteenaerialorthophotos from1946to2014havebeenanalyzedandcomplementedwithfieldobservationsfrom1986topresent.Thismorphologyisdynamicbutitisusuallycharacterizedby:a)longfingertransversebars(LFTB)andb)largescale shorelineundulations(LSSU).TheLFTBarethinandelongatedwithalengthoftheorderoftheirspacing.Theyareintertidalandtypicallyattachtotheshorelinebyamegacusp,commonlyopeningananti-clockwiseangleof 10°–40°withtheshorenormal.Therecanbemany,upto90,withboththemeanandthemostfrequentalongshorespacingintherange15‐–25 m.Spectralanalysisalwaysshowspeaksinthisrangeandsometimesadditional peaksintherange30‐–65 mthatcorrespondtothespacingbetweenthelargestbarswithsmallerbarsinbetween.TheLSSUtypicallyhavewavelengthsintherange150‐–250 m.Theirapexessometimescoincidewiththe shoreattachmentofthelargestbarsbutnotalways.NumericalmodellingshowsthatbothfeaturescouldemergeoutoffeedbacksbetweenhydrodynamicsandmorphologyduringtheSWwindeventsinvolvinga)deflectionof thelongshorecurrentbythebarscombinedwiththerefractivewavefocusingandb)gradientsintotalalongshoresedimenttransportratetriggeringthehigh-anglewaveinstability. Keywords:Longfingertransversebars;sShorelineundulations;sSpectralanalysis;AlfacsBay;lLowenergybeach;hHigh-anglewaveinstability;sSurfzonemorphodynamicinstabilities
presence of alongshore wavedriven current and theyare up-current oriented, i.e., their distal tip is shifted up-current withrespect tothe shore-attachment. Theirspacing is in therange L ≈ 15 − 200 m. Finally, the long finger transversebars type,LFTB,whichgroupsthe'‘large-scalefingerbars'’andthe'‘low-energyfingerbars'’oftheclassificationofPellónetal.(2014).Theyarecharacterizedbylongcrests,whicharetypicallylargerthanthealongshore spacingwhichmayvaryintherange L ≈ 10 − 500m.Theyaregenerallyobservedtobepersistentfeaturesinlowtomediumwaveenergysystemwhoseforeshoreisaveryflatterrace(Evans,1938;NiederodaandTanner,1970; BrunerandSmosna,1989;Falqués,1989;GelfenbaumandBrooks,2003;).Thewavefocusingcausedbytopographicrefractionbythebarsseemstobeanessentialprocesstothem.Althoughtheyaremostoftenobservedonmicrotidal beaches,theymayalsoexiston meso andmacrotidalcoasts(Levoyetal.,2013;Pellónetal.,2014).Theytendtobeorientedalmostperpendiculartotheshoreorwithasmallangle. Inadditiontosurf/shoalingzonebathymetry,theshorelineitselfcanalsodisplayanalongshorerhythmicityintheformofundulationsorcuspateshapes.Someofthemarelinkedtotransversebarsystems,sothattheirapexes developattheshoreattachmentsofthebarsandtheembaymentsinbetweencorrespondtothetroughsinbetweenbars.Theyareknownas megacusps .Shorelinesmayalsodisplayundulationsatascalewhichislargerthansurfzone rhythmicbars,i.e., L ≫ Xb ,where Xb isthewidthofthesurfzone.Theselargescaleundulationshavetypicalalongshorewavelengths>1Km(onopenoceanbeaches)andarelinkedtosimilarundulationsinthedepthcontourswell offshorethesurfzone.Theyhavebeencalled Km-scaleshorelinesandwaves (IdierandFalqués,2014). Theoriginanddrivingmechanismsofnearshorealongshorerhythmicpatternshavebeenpuzzlingscientistsforlong.Theirstrikingandrelativelyregularmorphologyalongwithitsclearcharacteristicalongshorespacingin spite of the multi-scale high complexity of coastal dynamics suggests that something important occurs in the physics at that particular length scale. A common assumption in the past was that a pre-existing template in the hydrodynamics(currentsand/orwaves)imprintsitsparticularpatternonthemorphology(HolmanandBowen,1982).Sincethefeedbackfromthedevelopingmorphologyintothehydrodynamicsisignoredinthisapproachthesetype ofexplanationshavebeenincreasinglydiscarded(CocoandMurray,2007).Alternatively,ithasbeenfoundthatrhythmicmorphologiescanemergefrominternalinstabilitiesofthecouplingbetweenmorphologyandhydrodynamics throughthesedimenttransport(CocoandMurray,2007;Ribasetal.,2015).Thishasbeenshownforcrescenticbars(Deigaardetal.,1999;Falquésetal.,2000;Calveteetal.,2005;Garnieretal.,2008),transversebars(Garnieretal., 2006;Ribasetal.,2003,2012),beachcusps(Cocoetal.,2000;Doddetal.,2008)andKm-scaleshorelineundulations(Ashtonetal.,2001;FalquésandCalvete,2005;vandenBergetal.,2012;KaergaardandFredsoe,2013),andis knownasthe self-organizationtheory . Checkingtheself-organizationexplanationforalongshorerhythmicfeaturesneedsfrequentanddetailedbathymetricsurveysalongwithmeasurementsofcurrentsandwavesduringthetimewherethefeaturesareincipient. Sincethistypeofinformationisveryscarce,providingdetailedandhighqualitydataonrhythmicmorphologydevelopmentanddynamicsisstronglyvaluable.ThesouthwesternflankoftheEbrodelta(Catalonia)featuresalongspitor barrierbeachcalledTrabucadorbeach,whichseparatestheAlfacsBayfromtheopenMediterraneanSea.TheinnersideoftheTrabucadorbeachcommonlydisplaystransversebarsandmegacuspswithatypicalalongshorespacingof tensofm.Atthesametime,aerialphotos(availablesince1946)alsoshowthepresenceoflargerscaleshorelineundulations,upto~250 m.TheTrabucadorbeachcanbeconsideredas"“low-energy"”sinceitisshelteredfromthe openseawaves.ApreliminarydescriptionofthissystemwasdonebyFalqués(1989). There are in the literature manytheoretical studies onrhythmic surfzone bars andonlarge scaleshorelineundulations. However,theapplication of such studiestospecific features observed in natureismuch scarcer. Moreover,therearenostudiesofrhythmicsurfzonebarscoexistingwithlargescaleundulations.Therefore,theTrabucadorbeachcasehasprovidedauniqueopportunityofstudying,forthefirsttime,bothtypesoffeaturesin coexistenceandfrombothviews:observationsandmodelling. Theprimaryobjectiveofthepresentpaperisamorethoroughdescriptionandcharacterizationoftherhythmicmorphologyinthisbeach.Wemainlyusespectralanalysisoftheshorelinesfromaerialorthophotosduringthe period1946‐–2014 alongwith occasional field observationsbyone ofthe authorsduring the period1986‐–2014. Wehypothesize that their morphology and thelength scales are largely self-organized and thesecond objective is investigatingthephysicalmechanismsdrivingitsformationanddynamicsbymeansofnumericalmodelling. 2.2StudyAreaandmorphology 2.1.2.1Geographicalsettingandmeteo-oceanographicconditions TheEbrodeltaislocatedinthewesternMediterraneansea,NEoftheIberianpeninsula,inCatalonia(seeFigure.1).Asaresultoftheinteractionoftheriverandmarinedynamics(Jiménezetal.,1997),thepresentdaydelta extendsaround25 kmoffshoreandhasanareaofabout325Km2.Onthelateralmarginsitfeaturestwospits,onetrendingtotheNW(ElFangar),anothertrendingtotheSWandconsistingoftheTrabucadorbarrierbeachandLa Banyapeninsula.Betweenthesouthwesternspit(TrabucadorandLaBanya)andthemainland,thereisasemi-enclosedbaycalledAlfacsbay.Itisnearlyrectangularabout16 kmlongand4 kmwide,withanaveragedepthofabout4 m (6.5inthecenterofthebay).TheAlfacsbaycanbegeomorphologicallyclassifiedasabar-builtsemi-enclosedlagoonformedastheresultofthereworkingofanabandonedlobeofthefirstEbroRivermouth,directedtothesouthby waves(Pritchard,1952,Nienhuisetal.,2017). The bay receivesconsiderableamount of freshwater from rice crops(oftheorder of10m3·s‐−1during openedchannels)leading an almost stablestratificationthroughouttheyear.So, accordingtoconductivity,temperatureanddensityprofilersduringdifferentperiods,thesemi-enclosedlagoonisclassifiedhydrodynamicallyasasalt-wedgeestuary(Llebotetal.,2014;Cerralboetal.,2015a).Thebedismadeofsiltand
clay(largestpercentagesinthemiddleofthebay)(Palacinetal.,1991).However,nexttotheinnersideoftheTrabucadorbarrierbeachthereisasandyandshallowshelf,seeFigure.2,thatiscomposedoffinesandwith D 50 ≈ 0.15 − 0.3 mm(Falqués,1989).Thisshelfdeepensfrom0toroughly0.7 mwithinabout125 mwidthand,presumably,ithasbeenbuiltoutofsedimentsfromtheoutersideofthebeachduringoverwashevents. ThewindregimeintheEbroDeltaischaracterizedbytheprevalenceofNorth-Westerlywindsduringautumnandwinterandonshore/offshorewinds(dailysea-breeze)duringsummer(Cerralboetal.,2015b;Grifolletal.,2016). Usingone-yearobservedwinddata,Grifolletal.(2015)estimatethemonthlypercentageoftimeofNWwinds:duringwinterthispercentagewasabout80%decreasingduringsummerto20%.Seebreezeoccursdailyduringsummer andinalesserextendduringspring.North-Westerlywindsarise15 m·s‐−1andtypicalvaluesofsea-breezeareoftheorderof8 m·s‐−1.Accordingtothewindclimate,andbecausethisbeachissheltered,wavesareexpectedtobevery smallincomparisontotheopen-seasideoftheTrabucador. Thisbeachismicrotidal,withtidalrange~0.25 m,butsea-leveloscillationduetoseichesactivityisrelevant(Cerralboetal.,2014)withamplitudesofmorethan>0.5 mmeasuredwithinthebaywithcharacteristicperiodsof1 h and3 h(0thand1stseichemodes).ThewatercirculationinAlfacsbayatshorttemporalscales(ordersofhours)isdominatedbythelocalwindandtheoccasionalseicheactivity(Cerralboetal.,2015a).Atlargetemporalscalesthe Figure1Fig.1Left:theEbrodeltawiththeTrabucadorbeach(inyellow).Othergeographicallocationsmentionedinthetextarealsoshown.Right:TheTrabucadorbarrierbeachin2012.Theshallowshelfwiththetransversebarscanbeseenattheinnerside.Source:ICGC.(For interpretationofthereferencestocolourinthisfigurelegend,thereaderisreferredtothewebversionofthisarticle.) alt-text:Fig.1 Figure2Fig.2a)Detailedviewofthetransversebarsystemin2012.b)ViewoftheinnerTrabucadorbeachin2013wherethelargescaleundulationsarevisible.Source:ICGC.WhitearrowrepresentstheNorth. alt-text:Fig.2
waterresponseconsistsofacombinedeffectofagravitationalcirculationduetofreshwaterinputsandtheresidualeffectofwinds.Duringseicheactivitymaximumvelocitiesinthenodelocationreach0.5 m/s.Numericalmodelling excercises presented by Cerralbo et al. (2014) identify the antinode location near Trabucador beach shore, so that the expected velocities due to seiche activity are negligible. No measurements of currents are available at the TrabucadorinnerbeachwherebarsarelocatedbutthecurrentswereestimatedwithpassivetracersduringthefieldcampaignofSeptember1988(Falqués,1989).WithNorthwindof7.3 m/samaximumcurrentof0.30‐–0.35 m/stothe SWwasobservedoverthecrestofabarwithwaterdepthof0.05‐–0.1 m.Thecurrentwasclearlyfavoredbythebreakingofthesmallwavesoverthebar.Atthesametime,atthedeepsbetweenthebarswith0.3‐–0.5 mdepthcurrents of0.15‐–0.20 m/swereobserved.Themaximumcurrentintheoppositedirectionwasabout0.2 m/swith6.1 m/swindfromtheS-SE. 2.2.2.2Thealongshorerhythmicmorphologicalsystem ThissystemwasfirstdescribedbyFalqués(1989) on thebasis offieldobservations andtheanalysis ofaerial photos and maps.Periodicvisual inspectionssuggestthat the rhythmicsystem isverypersistent,althoughits prominencevariesfromveryobvioustonearlynonexistent.ThecharacteristicmorphologicalunitsareshowninFigure.2:barsandassociatedmegacuspsinFigure.2aandlargescaleshorelineundulationsinFigure.2b. 2.2.1.2.2.1Barsandmegacusps TheinnerTrabucadorbeachusuallyfeaturesmanylongfingertransversebars(LFTB).Theyareelongatedwithdifferentorientationsbutmostofthemtrendnearlyperpendiculartothecoastline(seeFigure.3)andtheyaretypicallyrotatedtothe leftwithrespecttoshorenormalviewingfromthebeach(anangle ϕ ≈ 10 − 40°),i.e.,theytrendapproximatelytotheW.Manyofthemreachthecoastlineandamegacuspdevelopsattheattachmentpoint(seeFigure.3c).Similarmegacuspshavebeen describedbyEvans(1939,1938).
ThelengthandalongshorespacingoftheLFTBareofthesameorderofmagnitude,theformertendingtobelargerthanthelatter.Thespacingroughlyrangesbetween5 mand100 mbutistypicallyabout20 m.Theirverticalreliefcanreach about0.4 m.Typically,theyareinter-tidal,becomingmostlyemergedatlowtide.Ascanbeseenintheaerialphotos,thewholebarsystemusuallyextendsoffshoreintotheAlfacsBayacrossasignificantportionoftheshallowshelf.Accordingtolocal fishermen,thebarshavebeentherefor,atleast,thelast70 years.Indeed,theyarealreadyvisibleinanaerialphotoof1946. Figure.3cshowsabathymetricsurveytakeninAugust1986(Falqués,1989)anddisplaying4mainbars.Theangleofthemainbarswiththeshorenormalis14°,24°,19°totheleftviewingfromthebeach.Thesebarshaveanasymmetriccross-bar profilesteeperattheNEsideandmilderattheSWside.Theiralongshorespacingisabout60‐–70 m.Atthetroughbetweentwomainbarsthereisaminorbarnotconnectedtotheshorelineandshowingaweakoppositeasymmetry.Thesesmallerbarsin betweenthelargerones(inlengthandinalongshorewavelength)arerelativelycommonallalongthebeach.Thebarsystemisdynamicandactive.AclearindicationwasobtainedinSeptember1988whenitwasobservedhowtheshapeofabarchanged andhowthecrestmigratedabout0. 4 mwithin20 hourshduringadailyfairweatherconditions.Themigrationwasconsistentwiththeasymmetryincross-barshape,i.e.,towardsthesteepestside.Itwasalsoobservedhowtheripplesatthebarcrestwere formed,changedtheshapeandsizeanddisappeared. Althoughtherearesomedifferences(e.g.,barorientationorlengthvsspacing)wethinkthesebarsshareimportantcharacteristicswiththosedescribedbyEvans(1938)astype2.Theyweredescribedas'‘largecuspswhichhavetheirapexes Figure3Fig.3a)Transversebarsatlowtideduringtheperiod1986‐–1988.b)Transversebarandtherefractivewavefocusingbyit.c)BathymetricmapofastretchoftheTrabucadorinnerbeachfromAugust1986showingaseriesof5barsand3megacusps.Thedepthsareincm. alt-text:Fig.3
continuingoutintothelakeasaridgeofsandonthelakebottom'’.Itwasclaimedtheyformonlowenergybeacheshavingaplentifulsupplyofsandandwithacross-shorebathymetricprofilewhichisabovetheequilibriumprofile,featuringashallow terrace. 2.2.2.2.2.2Largescaleshorelineundulations Accordingtovisual inspectionsoftheaerialphotos,inadditiontothesmall amplitude shoreline undulationsdeterminedbythemegacusps,therearelarger shoreline undulations(seeFigure.2b). The apexes of these undulations sometimes coincidewiththeattachmentofoneofthelargestbars.However,thereisnotaone-to-onerelationbetweenbarsandlargescaleundulationssince:i)therearealsosmallerbarswithmegacuspsattheembaymentsandii)sometimestheapexesofsome undulationsareapparentlynotassociatedtoanybar.Thewavelengthsoftheseundulationsareintherange100‐–250 m,i.e.,significantlylargerthanthetypicalalongshorespacingofthebars.Therefore,theycanbeidentifiedasKm-scaleshorelinesand wavesalthoughtheirwavelengthissmallerthanthatofthesimilarfeaturesdevelopingonopenoceanbeaches. 3.3Methods 3.1.3.1Imageprocessing ThissectiondescribesthemethodologyweusetoextracttheinformationofthebarpositionsfromtheimagesalongwiththealongshorewavelengthsusingaFourieranalysis.Aerialorthophotosfromthe InstitutCartogràfici GeològicdeCatalunya (ICGC)duringtheperiod1946‐–2014areanalyzedtoassessthetimeevolutionofthemorphologyandtoobtainitsalongshorewavelengths. TheavailableimagesarelistedinTable1.TheshorelineisclearlyvisibleinallofthembutidentifyingtheLFTBisnotsostraightforward.WehaveselectedthecentralstretchoftheTrabucadorBeach(2Km)wherethebarsare morevisible.Usually,thebarscanbedetectedfromelongatedpatcheswithacolourthatisinbetweenthecolourofthedrybeachandthecolourofthedeeperwater.Incaseofoldblackandwhiteimages,thedrybeachappearsin verylightgrey,thebarsabitdarkerandthetroughsdarkgrey.Incaseofcolourimages,thedrybeachappearsinlightyellow/brown,thebarsinthesamecolourscolorsbutdarker,andthedeeperwatersinblue.AscanbeseeninTable 1theimageshavedifferentresolution.Nevertheless,mostofthemareusefultolocateboththebarsandtheshorelineundulationsregardlessoftheirresolution,asshowninFigure.4a–b.However,theimagesfrom1983(Figure.4c), 1996,2000(Figure.4d)and2004arenotgoodenoughtocapturethelengthofthebarsduetothepoorqualityoftheorthophoto,butcanstillbeusedtoassessevaluatetheshorelineundulationsandthepositionofthebars. Table1OrthophotosavailableintheICGCwiththecorrespondingresolution.TheyearsthatarenotpresentinthetableareyearswheretheICGCdidnotperformaflightintheTrabucadorzone. alt-text:Table1 Year Resolution(cm/px) LFTBbars AlongshoreUnd. 1946 100 YES YES 1956 50 YES YES 1983 50 YES YES 1993 2500 NO YES 1994 50 YES YES 1996 2500 NO YES 2000 50 NO YES 2004 50 NO YES 2007 50 YES YES 2008 50 YES YES 2009 25 YES YES 2010 25 YES YES 2011 25 YES YES 2012 25 YES YES
2013 25 YES YES 2014 25 YES YES Themethodisbasedontrackingthesharpchangesincolourassociatedtothebarsandtotheshorelineinbetweenbarsbyamanualdigitalization.Inthisway,an ApparentShoreline/Bar signal(ASB)alongtheshoreisobtained. Thissignalis,therefore,influencedmainlybythesealevelandthelightduringtheflightoftheorthophoto.Thus,wedonotattempttoextracttheexactshorelinefromtheorthophotos,butasignalthatisrepresentativeenoughto capturethepositionofthemainbarsandtheundulationspresentintheshoreline.Inordertominimizetheinherenterrorsofamanualdigitalization,thesameprocedurewasrepeatedseveraltimes,obtainingdifferentsignalsfrom thesameorthophoto.Thesesensitivityanalysisonthemanualprofilingparametersshowednosignificantinfluencesonthefinalresults,sincethedifferencesbetweenwavelengthswerelowerthan5 metersm. Morespecifically,theASBsignalisobtainedfromtheimagesdetailedinTable1asfollows: 1. Manualprofiling:imagesareloadedtoaGISspecificsoftwareandaseriesofdotsismanuallyassignedtotheedgesofthebarcreststhatareattachedtotheshoreline(identifiedasasharpchangeincolour)andtotheshorelinestretchinbetweenbars(Figure.5a). Whilethereistypicallyasharpchangeincolouracrossthebarsthatallowsidentifyingrelativelywellbothedgesofabarcrest,thechangesalongthebarcrestsarequitegradual.Thus,thelengthofthebarsisnotwelldefinedfromtheimagesandisnotdiscussed here.Accordingtothehistoricalorthophotos,twoconcreteplatformswereconstructedalongtheinnerbeachandwerepresentintheshorelineprofilefrom1994to2000.InordertominimizethenoiseofthefinalFourieranalysis,theseplatformshavebeenremoved fromtheoriginalsignal. 2. Rotation:sincetheTrabucadorbarrierbeachisNorth-Eastoriented(55°)andtheFouriertransformneedsasingle-valuedsignal,theshoreline/barlineisrotatedtoapproximatethemeanshorelinealongthex-axis,usingaroadthatgoesalongthebarrierbeach(see Figure.5b).However,thebarsarestilltiltedtowardstheSouth-Westandasingle-valuedsignalisnotobtainedforallthebars.Thiscreatessomeproblemsinthenextstep. 3. Interpolation:theFourieranalysisrequiresanequispaceddigitalsignal.Thus,anewsignalisobtainedfromthepreviousstepbyinterpolatingtoanequispacedmesh.Alinearinterpolationwithaspacing ∆x = 0.5mturnedouttobesatisfactory,althoughweshould bearinmindthatitintroducesspuriouspeaksduetotheslightinclinationwithrespecttothenormaldirectionoftheTrabucadorbar,seeFigure.5c–d.ThesensitivitytodifferentinterpolationmethodsandgridspacinghasalsobeenexaminedinFigure.5dandthebest Figure4Fig.4Zoominoffourorthophotosduringdifferentyears.a)b)highqualityimages;c)d)lowqualityimages.Source:ICGC. alt-text:Fig.4
resultsareprovidedbythelinearinterpolation. Itisimportanttopointoutthatthepresentmethodologyincorporatessignificantnoisebelow10‐–15 mwavelengthduetothespuriouspeaksoftheASBsignaladdedduringtheinterpolationstep.Therefore,someexisting peaksaround10 mwillnotbecapturedaccuratelyusingourmethodology,asshowninFigure.5d.TheresultingASBsignalscanbeseeninFigure.6. Figure5Fig.5BrownpatchrepresentsthebeachcontourgivenbytheASBsignal.a)RawASBsignalfromthemanualprofiling.Theyellowlineindicatesthedirectionofthexaxisafterrotation;b)RotatedASBsignal;c)SensitivityanalysisoftheinterpolatedASBsignalasa functionoftheinterpolationspacing, ∆x .Thedifferentinterpolatedsignalsareshownindifferentcolours.;d)SensitivityanalysisoftheASBsignalinfunctionoftheinterpolationmethod.(Forinterpretationofthereferencestocolourinthisfigurelegend,thereaderisreferredtothe webversionofthisarticle.) alt-text:Fig.5
TheDiscreteFouriertransformofacontinuoussignal, f ( ξ ),inthespatialdomainisdefinedby: wherethecoefficientsare: andtheenergyofeachwavenumberis where λ 1 is the fundamental wavelength and ξ is the spatial coordinate. The interpolation of the original ASB signal of the images provides a discrete signal and must be analyzed through a Discrete Fourier Transform (DFT) algorithm.SeveralalgorithmscanbefoundintheliteraturetocomputetheDFTandmostofthemareclearlyaffectedbythesamplingfrequencyandotherinherentalgorithmparameters.Thepresentresearchusesthe FastFourier Transform (FFT)assumingthatvalidwavenumberswillbebelow1 m‐−1.TheresultsobtainedaftertheFFTcomputationareveryclearforhighwavelengths(Figure.7a)butdonotshowrobustpeaksinthelowwavelengthrange,see Figure.7b. Therefore, low wavelength scales are analyzed using Bartlett’'s Method. Bartlett (1948) estimates the power spectra by averaging the DFT results of N segments of the original signal. This is useful if the wavelength distributionisuniformalongthesignalbutmaysmoothdownsomepeaksthatarepresentataparticularlocationofthesignal.Therefore,Bartlett'sresultswillenhancetheexistenceofmainwavelengths,presentrepetitivelyatthe Figure6Fig.6Manualprofilingresultsoftheorthophotoimages. alt-text:Fig.6 (1) (2) (3)
Thedepthofclosure,Dc,isanimportantparameterforshorelinesandwaveformationoverwhichthereissomeuncertainty(Falquésetal.,2017).FromthemorphologyoftheinnerTrabucadorbeachencompassingfirstashallowterraceupto 0. 7 mdepthandthenasteepbreaktowardstheAlfacsbay(Figure.1),avalueDc ~ 0. 7 mseemsreasonable.Moreimportantly,theaerialphotosrevealthatthealongshorerhythmicbathymetryextendsuptowaterdepthofthisorderandthendisappears (Figure.1).Thus,wehereexploretwovalues,Dc = 0.6 mandDc = 0.75 m.Regardingtheshapeofthebathymetricperturbationalineardecayinbedlevelfrom1attheshorelineto0at Dc isconsidered.Theeffectofdifferentshapefunctionsonmodel shorelinesandwaveformationisdiscussedinIdieretal.(2017).RegardingthewaveconditionsweusetheSWANoutputsincaseofSWbreezeasareference.ThewaveheightisarobustoutputoftheSWANmodelsothatwefixHs = 0.25 mandwe exploreTpfrom2.5 sdownto1.25 s.Regardingthewaveangle,itisadvisabletoexploreawiderangeincludingtheoutputofSWAN(inordertoobtainathresholdanglefortheinstability)sothatweconsiderwaveanglesfromθ = 20°uptoθ = 70°(at D = 1.6 m).Finally,concerningthewavelengthoftheundulationsweexaminetherange100 m ≤ λ ≤ 500 m.Forthatrangeofwavelengths,theinstabilitycurves(growthrateσrasafunctionofλ)lookliketheonerepresentedinFigure.11bforHs = 0.25 m, Tp = 1.25 s,θ = 50°andDc = 0.6 m.Inthiscase,themaximumgrowthrateisσr=0.032d‐−1whichdefinesthedominantwavelength,λM = 180 m.Thismeansthatthecharacteristictimeformationofthoseundulationsisaboutσr‐−1 ≈ 31dassumingconstant waveconditions.AsshowninFigure.11c,thedominantwavelengthhasbeenexploredasafunctionofθforvariousDcandTp.Firstitisseenthattheformationofshorelineundulationsneedsaminimumwaveangleinaccordancewithhigh-anglewave instabilitytheory(Ashtonetal.,2001;FalquésandCalvete,2005).ForTp = 1.25 s,thisthresholdangleisabout40°forDc = 0.75 mand45°forDc = 0.6 m.ForTp = 2.5 sandDc = 0.6 mitis50°.Itisfoundthatthedominantwavelengthsrangebetween180 and500 mandtheyincreasewithDc.ForDc = 0.6 mandexcludingthosecorrespondingtotheanglesnearthethreshold,theyareintherange180‐–280 m.ForDc = 0.75 mtheyareintherange290‐–500 m.Therefore,therangeofvariabilityofthese wavelengthsasafunctionofthewaveparametersmatchesquitewelltherangeoflargealongshorelengthscalesrevealedbytheimageprocessingpresentedinSection4.1.Inparticular,therange150‐–250 mseemsquitewellcapturedbytheDc = 0.6 m modelresults.Thewavelengthsabout400 mcouldberoughlyassociatedtoDc = 0.75 m.Importantly,theexistenceofathresholdangleforshorelineinstabilityrulesoutthepossibilitythatNWwindeventscouldgeneratesuchlargescaleundulations. 5.5Discussionandpossibledrivingmechanisms 5.1.5.1Morphodynamiccharacteristicsoftheobservations. Previously,itisimportanttoidentifysomeofthecharacteristicsofthemorphologicalsystemthatcanprovidecluesintothedrivingprocesses: •Thesystemisnowadaysactive,itisnotrelict.ThisisveryobviousfromFigure.6.However,thefieldobservationsduringtheperiod1986topresentareveryscarceandtheaerialorthophotosareavailableonceayearatmost.Also,itissometimesdifficulttoascertain whetherthebarsarenotthereortheyaresimplynotvisible.Thus,withtheavailableinformationitisimpossibletoaccuratelyassessthetimevariabilityofthesystemandtrytolinkittoparticularmeteorologicalevents. Figure11Fig.11Modellinglargescaleshorelineundulationswith1D-morfoforstrongbreezefromtheSW.(a)cross-shoremeanbathymetricprofile,observed(diamonds)andanalyticalapproximationusedby1Dmorfo(thickline).(b)typicalinstabilitycurve,growthrateasafunction ofwavelength,showingadominantwavelengthabout180 m.TheparametersareHs = 0.25 m,Tp = 1.25 s,θ = 50°withDc = 0.6 m.(c)dominantwavelengthasafunctionofwaveangle,θ,forvariousDcandTpandforHs = 0.25 m.Belowθ = 40°thereisnoundulationgrowth. alt-text:Fig.11
•Therangeofwavelengths(Section4.1).Forthetransversebars/megacuspsitis15‐–65 mandtheoverallmostcommonisabout20 m.Forthelarge-scaleundulationsitisabout150‐–400 mbutthemosttypicalwavelengthsareabout150‐–250 m. •Refractivewavefocusingbythebars.AsshowninFigure.3b,thesmallwavesapproachingthebarsstronglyrefract,thewavecreststurntowardsthebarcrestandcrossthewavecrestscomingfromtheotherbarflank.Thiscreatesazonewithsmallbreakersoverthe barcrestandaquitestrongonshorecurrentoverthebar.Currentsofabout0.35 m/sandsedimentmobility(mainlybed-load)associatedtoithavebeenobservedbyFalqués(1989).Thisisaverystrikingandpersistentprocessonthisbeach.Thisprocesswasclearly describedbyNiederodaandTanner(1970)andwassuggestedtobeoneofthedrivingmechanismsofsometransversebars. •Obliquityandorientation.Accordingtotheaerialphotosthebarsarenotshorenormalbutoblique,i.e.,theyaretypicallyrotatedtotheleftviewingfromtheshoreananglerangingfrom10°to40°. •Alongshoreasymmetryofthecross-barbedprofile.Abathymetricmaptakenin1986(Figure.3c)showsanasymmetricprofileofallthebars,steeperattheNEflankofthebarsandmilderattheSWflank(exceptaminorone).Additionalbathymetricobservationsare notavailablesothatwedonotknowwhetherthisisalwaysthecaseornot.Fromthedistributionofdarkandlightcolorsintheaerialphotossomesymmetry/asymmetrycanbesometimesguessedanditseemsthattheasymmetryobservedinFigure.3cdoesnot alwaysoccur. 5.2.5.2Transversebars InSection4.2.2wehaveseenhowthemorfo55modeldescribestheself-organizedformationoftransversebarsoutofinitialbathymetricnoiseduringsome35 hourshofstrongSWbreeze(9 m/s).Ofcourse,35 hourshofsuch constantwindisunrealistic,butatleast6 hourshadayisquitecommonduringmanysummerdays.ThismeansthatthebarscouldformduringafewweeksofcommonsummerbreezeregimeatElTrabucador.Theobservedrangeof wavelengths,theobliqueorientation,thealongshoretranslationandeventhepresenceofdifferentbarsizesarecapturedbythemodel.Aquantitativecomparisononthealongshorecelerityisimpossiblewithoutknowingthewave conditionsfortheobservedtranslation.Accordingtothemorfo55equations,thephysicalprocessesthatareessentialfortheformationofthebarswouldbeasfollows.TheSWwinddriveswavesincomingobliquelytotheTrabucador beachfromtheW.ThebreakingofthesewavesdrivesalongshorecurrenttotheNE.Theincipientbarsactonthecurrentintwoways:a)theycauseameanderingofthecurrent,seawardatthecrestsandshorewardatthetroughs (bedfloweffect)andb)duetotherefractivefocusing,thereismorewaveenergyhencemorebreakingoverthecrestsandlesswaveenergyhencelessbreakingatthetroughs(bedsurfeffect).Thiscausesgradientsinwaveradiation stressesandinwave-setupthatalsoaltersthelongshorecurrent.Morfo55,then,computestheperturbedsedimentfluxanditturnsoutthatthegradientsinsedimentfluxcreatesedimentdepositionatthecrestsandsedimenterosion atthetroughstherebyapositivefeedbackoccurs.Thedetailsofthephysicalmechanismthatcandriveup-currentorientedtransversebarshadbeengenericallystudiedbyRibasetal.(2003),Garnieretal.(2006)andRibasetal.(2015). Regardingitsspecificapplicationtonature,Ribasetal.(2012)foundthatthemediumenergytransversefingerbarsobservedatNoordwijkbeach,theNetherlands,couldbeexplainedbythismechanism.However,thepresentresearch isthefirsttimethisfeedbackmechanismisfoundtodrivelongfingertransversebars(LFTB).AlthoughthephysicalmechanismissimilartoRibasetal.(2012)thereareanumberofdifferenceswiththepresentstudy(seeSection1).The mostessentialarei)thelevelofwaveenergy,ii)whethertheyareephemeralorpersistentandiii)whetherthereisashore-parallelbarornot.Moreover,thepresentstudyusesnonlinearmorphodynamictimeevolutionwhileRibasetal. (2012)uselinearstabilityanalysis.LFTBbarshavealsobeenobservedatElPuntalbeach,Santander(Cantabria,Spain)byPellónetal.(2014).However,therearetwoimportantdifferenceswiththeTrabucadorsystem:thosebarsare down-currentorientedandthetidalrangeissignificant,upto5 m.Moreover,amorphodynamicmodellingoftheirformationisnotavailablesothattheirformationmechanismisstillunknown. Thedrivingforceforthecurrentinthepresentsetupforthemorfo55modelisjustwavebreaking.OtherprocessescouldhelpindrivingacurrenttotheNEatElTrabucadorbeach.Firstofall,theSWwinditselfcoulddothis job.Pellónetal.(2014)comparedthewindandthewaveforcesonthewaterforatransversebarsysteminatidallow-energybeachandconcludedthattheywereofthesameorderofmagnitude.Cerralboetal.(2016)investigatethe watercirculationpatterninAlfacsBayusingafully3Doceannumericalmodel.TheirnumericalexperimentsshowedacurrentjetintheinnershoreoftheTrabucadorbeachassociatedtoNWwind.Also,lessfrequentNEwindsmaybe theoriginofasouthwestwardcurrentasalocalresponsetothewind.Asfarasweknow,thereisonlyonepublishedstudyonmorphodynamicinstabilitiesinducedbycurrentspartiallydrivenbywindinthesurfzone(Ribasetal.,2012). AlthoughtheenvironmentisquitedifferentfromElTrabucador(openoceanbarredbeach,mediumwaveenergyconditions)thatstudysuggeststhatwindeffectscouldreinforcetheself-organizationprocesspresentedinSection4.2.2. Certainly,moredetailedobservationsareneededalongwithspecificmodellingstudiesonmorphodynamicinstabilitiesgeneratedbypartiallywind-drivencurrent. TherefractivewavefocusingandbreakingbythebarswaspointedoutasthemainformationmechanismforLFTBbarsbyNiederodaandTanner(1970).Thiseffectisincludedinourmorfo55simulation(bedsurfeffect)butin combinationwiththelongshorecurrentdeflection(bedfloweffect)whichessentiallyoccursincaseofobliquewaveincidence.SincetheTrabucadorbarsareobliqueandfrequentlyasymmetrictheyverylikelyhaveformedunderthe actionofacurrent.IncontrastandaccordingtoNiederodaandTanner(1970),thebarswhichareprimarilyformedbytherefractivewavefocusingcanemergeevenincaseofnormalwaveincidencebeingshore-normalandcross-bar symmetric.Thus,weconcludethatthepresentbarsarenotformedonlybytherefractivewavefocusingbutbythecombinationofboththisprocessandthelongshorecurrentinteractionwiththebars. Atthispoint,wewonderwhichcouldbethepossibleroleofNWwindeventsonthebars.AmodelstudysimilartothatcarriedoutforSWwindisbeyondthescopeofthepresentworkandisleftforfutureresearch.Inthat case,thewaveincidenceisnearlyshore-normalandtherefractivewavefocusingcouldplayaroleininitiatingbarformationwithoutcross-barasymmetryandwithanearlyshore-normalorientation,accordingtoNiederodaandTanner (1970).Also,theNWwindcouldcontributetoreshapethebarsonceformedduringSWwindevents.Thus,wefinallyconcludethattheSWbreezewouldbetheprimarycauseofthetransversebarsalthoughtherecouldbesome
influenceofNWwindeventsthatstillremainsunknown. WeshouldfinallymentionanimportantoccasionalprocessthattakesplaceattheTrabucadorbarrierbeach.DuringseverestormsfromtheEorNEtheTrabucadorbarrierbeachissometimesfloodedandbreached(Graciaet al.,2013).Overwashfanscanthenappearandtheyperhapsmightactasaninitialperturbationthatwouldinitiatethefeedbackmechanismsbetweenmorphologyandhydrodynamics.Theself-organizedprocesseswouldlateroncreate thelongfingertransversebarsandinducethealongshoreperiodicity.Sincethebarsaremuchthinnerthantheoverwashfansitisunlikelythatthesefansweretheseedforthebars. 5.3.5.3Large-scaleshorelineundulations InSection4.2.3wehaveseenthatifweignorethetransversebarsandthedetailsofthesurfzonemorphologyapositivefeedbackbetweenhydrodynamicsandthesmoothedmorphologycanstilloccuratwavelengthsoneorder ofmagnitudelargerthantransversebarsiftheSWwaveincidenceangleislargeenough.Thisfeedbackisbasedonthealongshoregradientsinthecross-shoreaveragedbedelevationwhenthealongshoregradientsoftotallongshore transportcausedepositionattheshoalsoftheaveragedbathymetryanderosionatthedepressions.Thesegradientsoccurbecauseofthealongshorevariationsinwaverefractionandshoalingintheshoalingzone,whichmustbe coupledtothesurfzone.Dependingonsomeparameters,the1D-morfocomputationsgivearangeofwavelengthsλ ~ 180‐–500 mwhichmatchesquitewelltherangeoftheobservedlargescaleshorelineundulations.Wethusconclude thathigh-anglewaveinstabilityforSWbreezescouldprovideaplausibleexplanationforthelargescaleshorelineundulations. 6.6Concludingremarks TheanalysisoftheaerialorthophotosanddirectobservationsoftheTrabucadorbeachhaveshownthattheinnersidecommonlyfeaturesanintricateandcomplexmorphologythatisalongshorerhythmicatanumberof differentlengthscales.Themostapparentmorphologicalunitisthelongfingertransversebars(LFTB).Thesebarstypicallyattachtotheshorelinebyamegacuspandboththemeanandthemostfrequentalongshorespacingsare about L = 20m.Thespectralanalysisofthedigitalizedsignalsshowsmanypeakswithinter-annualvariability.However,thereisalwaysasignificantconcentrationofpeaksatthe15‐–25 mband.Thispreferentlenghscaleisalso confirmedbytheinsituobservations.Therefore,weconcludethatthedominantalongshorespacingoftheLFTBistypicallyabout20 m.Thespectralanalysisalsoshows(lessfrequent)peaksintheband L = 30 − 65m.Boththe orthophotosandthedirectobservationsshowthattheycorrespondtothespacingbetweenthelargestbarshavingminorbarsinbetween. Anothermorphologicalunitisthelargescaleshorelineundulations.TheirsystematicorrhythmicbehaviorisrevealedbythespectralanalysisoftheASBsignalsandbecomesveryclearinsomeoftheorthophotos.These undulationsseemtobeadifferentelementnotdirectlyrelatedtotheLFTB.Thecorrespondingalongshorewavelengthsofthoseundulationsrangebetweenabout150 mupto700 mbutthehighestconcentrationofspectralpeaks associatedtothemisaround200 m. Althoughthemorphologyishighlydynamic,itsmaincharacteristicsarepersistentintimeatleastforthelast70 years.Itisimportanttobearinmindthatthelackofsystematicexperimentalinformationhasnotyetalloweda detailedquantitativemodellingofthedevelopmentofthemorphologyduringthepertinentobservedmeteorologicalandhydrodynamicconditions.Despiteofit,wehaveusedthe2DHmorfo55andtheone-linelinearstability1D-morfo modelstoshowthatboththetransversebars/megacuspsandthelargescaleshorelineundulationscouldbeoriginatedbythewave-drivenalongshorecurrentduringthestrongSWbreezes.Thus,ourmodelresultssuggestthatthe primarycauseoftheobservedmorphologicalsystemisthefeedbacksbetweeni)wavesandcurrentsandii)morphology,throughsedimenttransportduringSWwindevents.Althoughtheinfluenceofotherweatherconditions,in particularNWorNEwind,couldbeimportantwethinkitisnotessentialforthegenesisofthemorphologyandisleftforfutureresearch.NotwithstandingthatourhypothesisthatthemaindriverisSWwindeventsseemshighly plausible,manyaspectsdeservefurtherattention.Fromthemodellingside,theroleofNWwindeventsandthewindforcingonthecurrentshouldbeinvestigated.Fromtheobservationalside,thestateofthebarsystemshouldbe monitoredatleastduringoneyeartoseewhetherbarsaremoreprominentinsummerorinwinter,andtoseekforpossiblecorrelationwithmeteorologicalevents,inparticularwithprevailingSWwinds.Also,directmeasurementof wavesandcurrentsatthesitewouldbehighlydesirableasapreviousstepofamorespecificmorphodynamicmodelling.Otherfutureworkswillbeaddressedtoevaluatethesignificanceoftheresultsatmethodologicalstepssuchas thecoastlineprofilingortheFourieranalysis. Uncitedreferences Calveteetal.,2001 Hallermeier,1978 Acknowledgements ThisresearchispartoftheSpanishGovernmentprojectCTM2015-63166225-C2-1-P(MINECO/FEDER).TheauthorswanttogratefullyacknowledgeDr.J.Guillénforprovidingbathymetricdataonthecrossshore beach profile at El Trabucador, Dr. R. Garnier for his support for the morfo55 runs and Dr. Cerralbo for providing numerical mesh information. The first author’'s funding comes from the research group
LIM/UPC. References AshtonA.,MurrayA.B.andArnaultO.,Formationofcoastlinefeaturesbylarge-scaleinstabilitiesinducedbyhigh-anglewaves, Nature 414,2001,296–300. BartlettM.S.,Smoothingperiodogramsfromtimeserieswithcontinuousspectra, Nature 161,1948,686–687. BooijN.,RisR.C.andHolthuijsenL.H.,Athird-generationwavemodelforcoastalregions:1.Modeldescriptionandvalidation, J.Geophys.Res.Oceans 104,1999,7649–7666,https://doi.org/10.1029/98JC02622. BrunerK.R.andSmosnaR.A.,Themovementandstabilizationofbeachsandontransversebars,AssateagueIsland, Virginia.J.Coast.Res. 5,1989,593–601. CalveteD.,FalquésA.,deSwartH.E.andWalgreenM.,Modellingtheformationofshoreface-connectedsandridgesonstorm-dominatedinnershelves, J.FluidMech. 441,2001,169–193, https://doi.org/10.1017/S0022112001004815. CalveteD.,DoddN.,FalquésA.andvanLeeuwenS.M.,MorphologicalDevelopmentofRipChannelSystems:NormalandNnearNormalWaveIncidence, J.Geophys.Res. 110,2005. CerralboP.,GrifollM.,Valle-LevinsonA.andEspinoM.,Tidaltransformationandresonanceinashort,microtidalMediterraneanestuary(AlfacsBayinEbredelta), Estuar.Coast.ShelfSci. 145,2014,57–68. CerralboP.,GrifollM.andEspinoM.,Hydrodynamicresponseinamicrotidalandshallowbayunderenergeticwindandseicheepisodes, J.Mar.Syst. 149,2015a,1–13,https://doi.org/10.1016/J.JMARSYS.2015.04.003. CerralboP.,GrifollM.,MoréJ.,BravoM.,AfifA.S.andEspinoM.,Windvariabilityinacoastalarea(AlfacsBay,EbroRiverdelta), Adv.Sci.Res. 12,2015b,11–21. CerralboP.,EspinoM.andGrifollM.,Modelingcirculationpatternsinducedbyspatialcross-shorewindvariabilityinasmall-sizecoastalembayment, OceanModel.OceanModel 104,2016,84–98, https://doi.org/10.1016/J.OCEMOD.2016.05.011. CocoG.andMurrayA.B.,Patternsinthesand:Ffromforcingtemplatestoself-organization, Geomorphology 91,2007. CocoG.,HuntleyD.A.andO’'HareT.J.,Investigationofaself-organizationmodelforbeachcuspformationanddevelopment, J.Geophys.Res. 105,2000,21991–22002. DeigaardR.,DrønenN.,FredsoeJ.,JensenJ.H.andJørgensenM.P.,Amorphologicalstabilityanalysisforalongstraightbarredcoast, Coast.Eng. 36,1999,171–195. DoddN.,StokerA.,CalveteD.andSriariyawatA.,OnBeachCuspFormation, J.FluidMech. 597,2008,145–169. EvansO.F.,Theclassificationandoriginofbeachcusps, J.Geol. 46,1938,615–627. EvansO.F.,Masstransportationofsedimentsonsubaqueousterraces, J.Geol. 47,1939,325–334. FalquésA.,FormacióndeTopografíaarítmicaenelDeltadelEbro, Rev.Geofísica 45,1989,143–156. FalquésA.andCalveteD.,Largescaledynamicsofsandycoastlines.Diffusivityandinstability, J.Geophys.Res. 110,2005. FalquésA.,CocoG.andHuntleyD.A.,Amechanismforthegenerationofwave-drivenrhythmicpatternsinthesurfzone, J.Geophys.Res. 105,2000,24071–24088. FalquésA.,RibasF.,IdierD.andArriagaJ.,Formationmechanismsforself-organizedkm-scaleshorelinesandwaves, J.Geophys.Res.EarthSurf. 2017,122. GarnierR.,CalveteD.,FalquésA.andCaballeriaM.,Generationandnonlinearevolutionofshore-oblique/transversesandbars, J.FluidMech. 567,2006,327–360. GarnierR.,CalveteD.,FalquésA.andDoddN.,Modellingtheformationandthelong-termbehaviorofripchannelsystemsfromthedeformationofalongshorebar, J.Geophys.Res. 113,2008. GelfenbaumG.andBrooksG.R.,Themorphologyandmigrationoftransversebarsoffthewest-centralFloridacoast, Mar.Geol. 200,2003,273–289. GraciaV.,GarciaM.,GrifollM.andSánchez-ArcillaA.,Breachingofabarrierunderextremeevents.Theroleofmorphodynamicsimulations, J.Coast.Res. 65,2013,951–956.
GrifollM.,AretxabaletaA.L.andEspinoM.,Shelfresponsetointenseoffshorewind, J.Geophys.Res.Oceans 120,2015. GrifollM.,NavarroJ.,PallaresE.,RàfolsL.,EspinoM.andPalomaresA.,Ocean-atmosphere-wavecharacterisationofawindjet(Ebroshelf,NWMediterraneanSea), NonlinearProcess.Geophys. 23,2016,143–158. Guillén,J.,Acosta,J.,Chiocci,F.L.,Palanques,A.,2017.AtlasofBedformsintheWesternMediterranean.SpringerInternationalPublishing,Cham.doi:https://doi.org/10.1007/978-3-319-33940-5. Hallermeier,R.J.,1978.Usesforacalculatedlimitdepthtobeacherosion,in:CoastalEng.1978.Am.Soc.ofCiv.Eng.,pp.1493–1512. HolmanR.A.andBowenA.J.,Bars,bumps,andholes:modelsforthegenerationofcomplexbeachtopography, J.Geophys.Res. 87,1982,457–468. IdierD.andFalquésA.,Howkilometricsandyshorelineundulationscorrelatewithwaveandmorphologycharacteristics:preliminaryanalysisonthe{A}tlanticcoastof{A}frica, Adv.Geosci. 39,2014,55–60. IdierD.,FalquésA.,RohmerJ.andArriagaJ.,Self-organizedkilometre-scaleshorelinesandwavegeneration:sensitivitytomodelandphysicalparameters, J.Geophys.Res. 122,2017. JiménezJ.,Sánchez-ArcillaA.,ValdemoroH.I.,GraciaV.andNietoF.,ProcessesreshapingtheEbrodelta, Mar.Geol. 144,1997,59–79,https://doi.org/10.1016/S0025-3227(97)00076-5. KaergaardK.andFredsoeJ.,Numericalmodelingofshorelineundulationspart1:Cconstantwaveclimate, Coast.Eng. 75,2013,64–76. Khabidov,A.,2001.Transversebarsformationonatidelessbeach,in:ofCivilEngineers.Lund,S.(Ed.),ProceedingsCoastalDynamics2001.pPp.666–672. KomarP.D.,Beach{P}rocessesand{S}edimentation,Second.ed.,1998,PrenticeHall;EnglewoodCliffs,N.J. KonickiK.M.andHolmanR.A.,Thestatisticsandkinematicsoftransversebarsonanopencoast, Mar.Geol. 169,2000,69–101. LevoyF.,AnthonyE.J.,MonfortO.,RobinN.andBretelP.,Formationandmigrationoftransversebarsalongatidalsandycoastdeducedfrommulti-temporalLidardatasets, Mar.Geol. 342,2013,39–52. LinW.,SanfordL.P.andSuttlesS.E.,WavemeasurementandmodelinginChesapeakeBay, Cont.ShelfRes. 18,2002,2673–2686. LlebotC.,RuedaF.J.,SoléJ.,ArtigasM.L.andEstradaM.,Hydrodynamicstatesinawind-drivenmicrotidalestuary(AlfacsBay), J.SeaRes. 85,2014,263–276,https://doi.org/10.1016/J.SEARES.2013.05.010. MoeiniM.H.andEtemad-ShahidiA.,ApplicationoftwonumericalmodelsforwavehindcastinginLakeErie, Appl.OceanRes. 29,2007,137–145. NiederodaA.W.andTannerW.F.,Preliminarystudyontransversebars, Mar.Geol. 9,1970,41–62. NienhuisJ.H.J.H.,AshtonA.D.A.D.,KettnerA.J.A.J.andGiosanL.,Large-scalecoastalandfluvialmodelsconstrainthelateHoloceneevolutionoftheEbroDelta, EarthSurf.Dyn. 5,2017,585–603. PalacinC.,MartinD.andGiliJ.M.,FeaturesofspatialdistributionofbenthicinfaunainaMediterraneanshallow-waterbay, Mar.Biol. 110,1991,315–321. PellónE.,GarnierR.andMedinaR.,IntertidalfingerbarsatElPuntalspit,BayofSantander,Spain:observationandforcinganalysis, EarthSurf.Dyn. 2,2014,349–361. Pritchard,A.,1952.EstuarineHydrography.AdvancesinGeophysics,vol.1.AcademicPress,Inc.,NewYork,USA. RibasF.andKroonA.,Characteristicsanddynamicsofsurfzonetransversefingerbars, J.Geophys.Res. 112,2007. RibasF.,FalquésA.andMontotoA.,Nearshoreobliquesandbars, J.Geophys.Res. 108,2003. RibasF.,deSwartH.E.,CalveteD.andFalquésA.,Modelingandanalyzingobservedtransversesandbarsinthesurfzone, J.Geophys.Res. 117,2012. RibasF.,FalquésA.,deSwartH.E.,DoddN.,GarnierR.andCalveteD.,Understandingcoastalmorphodynamicpatternsfromdepth-averagedsedimentconcentration, Rev.Geophys. 53,2015. ShepardF.P.,Revisednomenclaturefordepositionalcoastalfeatures, Bull.Am.Assoc.Pet.Geol. 36,1952,1902–1912. ShortA.D.,HandbookofbBeachandsShorefacemMorphodynamics,1999,JohnWiley. vandenBergN.,FalquésA.andRibasF.,Modellinglargescaleshorelinesandwavesunderobliquewaveincidence, J.Geophys.Res. 117,2012.