Tidal and subtidal sea level variability at the northern shelf of the Brazilian Northeast Region

A characterization of the sea level variability at tidal and sub-tidal frequencies at the northern shore of the Brazilian Northeast shelf for the period 2009–2011 is presented. The sea level data used was obtained from the Permanent Geodetic Tide Network from the Brazilian Institute of Geography and Statistics for the Fortaleza gauge station. Local wind data was also used to assess its effects on the low-frequency sea level variability. The variability of the sea level was investigated by classical harmonic analysis and by morphology assessment over the tidal signal. The low frequencies were obtained by low-pass filtering. The tidal range oscillated with the highest value of 3.3 m during the equinox and the lowest value of 0.7 m during the solstice. Differences between the spring and neap tides were as high as 1 m. A total of 59 tidal constituents were obtained from harmonic analysis, and the regional tide was classified as semi-diurnal pure with a form number of 0.11. An assessment of the monthly variability of the main tidal constituents (M2, S2, N2, O1, and K1) indicated that the main semi-diurnal solar S2 presented the highest variability, ranging from 0.21 to 0.41 m; it was the main element altering the form number through the years. The low frequency sea-level variability is negligible, although there is a persistent signal with an energy peak in the 10–15 day period, and it cannot be explained by the effects of local winds.

The tidal energy spans on different frequency bands (Munk and Cartwright 1966), although most of the energy is concentrated in the semi-diurnal (12.42 h) and diurnal (24.83 h) frequencies, called tidal frequencies (TF) to distinguish them from those that occur at higher frequencies, called supratidal frequencies (SupTF), and those that occur at lower frequencies, called sub-tidal frequencies (SubTF).Most of the coastal sea level variability at SubTF is related to meteorological events, trapped shelf waves, and river discharge (Miranda et al. 2002, Truccolo et al. 2006).
In spite of its importance for coastal management, studies on the sea level variability along the Brazilian Northeast's shores are scarce.In fact, there are only three published studies regarding tides and sea level variability in this region.Lessa et al. (2001) reported the tide characteristics in the Todos os Santos Bay, northern coast of Bahia State; Lessa and Cirano (2006) assessed the hydrodynamics of the Canal Sueste, south of Bahia State, and also provided a characterization of the tides; Schettini et al. (2011)  in tidal records from harbors in restricted seas (e.g., Pugh 1987).

phYsiCal seTTiNG
The northern continental shelf of the Brazilian Northeast (NCS-NE) region encompasses the states of Rio Grande do Norte (RN), Ceará (CE), Piauí (Pi), and Maranhão (MA). it can be separated into two main sectors, the Maranhão shelf and the shelf of the other states.The Maranhão shelf is wider with a general E-W orientation and has the "reentrâncias Maranhenses" which is the landward coastline indentation that forms the São Marcos and São José do Ribamar bays.The other sector shows a variable orientation, being E-W in the RN, SE-NE in the east portion of CE, and E-W again in the western portion of CE and Pi (Figure 1a).The relief is mild in general, with localized reef formations close to the shore.The shelf width is of the order of 60 km, and the shelf break occurs at nearly 70 m with an average slope of 2 m/km (Freire 1985).
The regional wind regime is strongly modulated by the zonal displacement of the Inter-Tropical Convergence Zone (ITCZ).During the period from August to November, the iTCZ is at a northern position, with the occurrence of strong southeasterly trade winds in the NCS-NE and lower or even a null precipitation rate.During March to May the iTCZ is over the NCS-NE, resulting in weaker northeast trade winds and higher precipitation rates as well (Uvo and Nobre 1989, Castro and Miranda 1998, Campos 2003).The intensification of the southeasterly winds during the dry period intensifies the thermal gradient between land and sea, favoring the occurrence of land/sea breezes (Ferreira 1996).The iTCZ, and as a consequence the wind pattern, is influenced by large-scale oceanographic phenomena such as the Atlantic Dipole and the El Niño-Southern Oscillation (ENSO) cycle (e.g., Cane and Zebiak 1985, Uvo and Berndtsson 1996).The El Niño events -the ENSO warm phase -is associated with dry conditions due to displacement of the ITCZ further north and a high pressure area in the Northeast region (e.g., Cavalcanti 1996), which also alters the main atmospheric circulations.
The tides in the Tropical Atlantic Ocean are characterized by two main amphidromic systems, one at the north, in the Gulf of Mexico, and the other at the south, in the central portion of the Atlantic basin (Trujillo and Thurman 2011).Figure 2 shows the co-tidal charts of the main semidiurnal constituent, M 2 , for the tropical Atlantic Ocean provided by Cartwright et al. (1991).Table I shows the main harmonic tidal constituents for Natal (RN), Fortaleza (CE), and Luis Correia (Pi).These harmonic constituents were obtained from data sets of 30 days (FEMAR 1999), and show high similarity in terms of relative weight of the semidiurnal ones, and a trend of increase in the range eastward as indicated by the co-tidal chart (Figure   2).The tide in Fortaleza is pure semi-diurnal, with mean neap and spring tide ranges of 1.5 and 2.8 m, respectively (Schettini et al. 2011).

MaTeRials aND MeThoDs
Sea level data was obtained from the tide gauge station at Fortaleza, CE.This tide gauge is part of the Permanent Tidal Monitoring Network for Geodetics of the Brazilian Institute of Geography and Statistics (RMPG-iBGE, acronyms from Portuguese). it is located in the Fortaleza harbor at 03° 42.9' S and 38° 28.6' W (Figure 1b). in contrast to most of the tide gauges along the Brazilian shore, which are located in estuaries, this one is located in a semi-sheltered area that still presents open shore conditions.
The Fortaleza tide gauge station was commissioned in September 2007, with a conventional analogical tidal gauge, providing hourly sea level records.A digital tidal gauge was added in April 2008, providing sea level records at 5 minute intervals.A geodetic GPS was added in August 2008.The station also has a standard RMPG tidal staff, and the tide gauges are assembled in two resting wells.The station is leveled to the national reference level of iBGE.Further details can be accessed at http://www.ibge.gov.br/home/geociencias/geodesia/default.shtm.In the present study, we analyzed the sea level data recorded from January 2009 to December 2011.
The temporal sea level series were visually inspected to identify spurious signals, missing values, and step changes in the mean sea level.The total amount of errors accounted for less than 0.1% of the entire record.The missing data periods ranged between a few hours to a day, and were filled by the harmonic signal obtained by the harmonic prediction for the missing period using the Matlab (Mathworks Inc.) T-Tide routines (Pawlowicz et al. 2002, available at http://www2.ocgy.ubc.ca/~rich/#T_Tide).Sea levels shifts, positive or negative, were corrected by the adjustment of the averaged sea level referenced to the long-term sea level.Figure 3 (a1, b1 and c1) shows the time series of the observed sea level at Porto de Fortaleza for the years 2009, 2010 and 2011 after corrections.
The sea level data, recorded at 5-minute intervals, were reduced to hourly intervals and leveled to the mean level.The sea level data was initially analyzed for its morphometric characteristics by extracting the extreme high (hW) and low (LW) sea level values to obtain a time series of tidal height.The tidal height was obtained by the difference of consecutive LW and hW, and the reference time for a given set of values was the mean time between the respective LW and hW.
The harmonic analysis was performed for each year (2009, 2010, and 2011), and also for each month using a window of 29 days starting always at the first day of each month.For the month of February, data for March 1 st was added to complete the window.
The tidal form number (FN), introduced by Courtier (1938( , apud Defant 1961)), is calculated as the ratio of the sum of the amplitudes, , of the main diurnal harmonic constituents (O 1 -main diurnal lunar, and K 1 -main diurnal solar) to the sum of the main semi-diurnal harmonic constituents (M 2 -main semi-diurnal lunar, and S 2 -main semidiurnal solar), aO aK FN aM aS FN is used to classify the tide dominance.When FN < 0.25, the main tidal regime is semidiurnal pure; when 0.25 < FN < 1.5, the tidal regime is mixed with semi-diurnal dominance; when 1.5 < FN < 3, the tidal regime is mixed with diurnal dominance; and when FN > 3 the tidal regime is diurnal pure (Defant 1961).
The non-astronomical sea level signal (SubTF) was obtained through the filtration of the original dataset with a temporal window (Walters and heston 1982), applying a cut frequency smaller than the Nyquist frequency to avoid signal aliasing (Jenkins andWatts 1968, Newland 1975).The cutting frequency used was 48 hours.The signal was also detrended by a linear relationship with time.
Wind speed and direction data were obtained from a meteorological station located in the campus of the Federal University of Ceará (UFC),  Cienc (2016) 88 (3) TiDES AT ThE NORThERN ShELF OF BRAZiL 1377 positioned at 03° 43' S and 38° 33' W (Fig. 1b).The data is recorded three times per day, at 9, 15, and 21 hours.The location and acquisition rate are not ideal; however, it was the only source of data available.The station is located a few kilometers inshore, and the local topography, the presence of natural and artificial obstacles, variations in the heat fluxes among distinct surfaces, etc., can decrease the wind speed compared to the wind observed offshore (Silva 2003).however, it is still useful since we are looking for some phase agreement between the wind and sea level.

An Acad Bras
The wind was decomposed in the along and cross-shore components after a rotation by 53°, which is the general orientation of the coastline.The wind components were interpolated by applying the cubic spline method (e.g., Emery and Thomson 1997) to obtain an hourly signal.The spline method segments the time series in smaller parts and applies small order polynomials, and the sum of these polynomials build the interpolated signal over the domain (Mazzini and Schettini 2009).The wind signal was also filtered to SubTF following the same procedure described for the sea level.
The relationship between the wind and sea level was assessed in terms of the yield shear stress over the sea interface.The wind stress, , was obtained by where is the density of air (~1.22 kgm −3 ), is the drag coefficient (1.153; Stech and Lorenzetti 1992), and is the wind speed.
The relation between the SubTF wind stress and sea level was investigated by applying the wavelets analysis which allows the identification of not only the frequency energy distribution for a given period, but also the visualization of its variability with time.Like the Fourier analysis, the wavelet gives the wavelet power spectra, cross-wavelet spectra and wavelet coherence and phase spectra, which were calculated following the procedure described in Torrence and Compo (1998), Torrence and Webster (1999) and Grinsted et al. (2004).The wavelets analysis was adapted from the Matlab (Mathworks Inc.) routines available at http://www.pol.ac.uk/home/research/waveletcoherence/.It uses the Morlet wavelet, which uses a continuous function modulated by a Gaussian wave.

ResulTs
Figure 3 (a2, b2 and c2) shows the time series of the tidal height for the years 2009, 2010, and 2011, with indications (black dots) of the extreme values of the tidal height for the spring and neap tides.The average tidal height was 1.95 m for the three years.The mean spring tide peak height was 2.79 m, ranging from 2.23 to 3.34 m.The mean neap tide peak height was 1.12 m, ranging from 0.72 to 1.72 m.The largest spring tide heights occurred during the equinox period (March and September), which had the smallest neap tide heights, while the smallest spring tide heights occurred during the solstice period (June and December), which had the largest neap tide heights.
Figure 4a shows the relationship between the observed and predicted tidal signals for the year 2009, at hourly intervals.The agreement is very high with the r 2 nearly at unity.Therefore, there is some scattering from the 1:1 ratio toward the largest heights.For the observed height of 2.5 m, the respective predicted heights range from 2.0 to 3.0 m, although for a limited number of occurrences.Figure 4b shows the tidal height occurrence distribution from the observed and predicted signals for the year 2009 and Figure 4c shows the neap and spring tide occurrence distribution from the observed and predicted signals.There is a very good agreement for the tide height and spring tide height occurrences, and therefore, the agreement is not as good for the neap tide heights, where the FELiPE F. FROTA, ELiANE C. TRUCCOLO and CARLOS A.F. SChETTiNi mode of the observed heights is smaller than the mode of the predicted heights.The tidal height distribution is asymmetric, with a dominance of heights smaller than 2 m.The main neap tide height mode was about 1.2 m, and the main spring tide height mode was about 2.4 m.
The results of the tidal harmonic analysis performed for the datasets of each year resulted in 59 harmonic constituents.The main diurnal and semi-diurnal harmonic constituents are shown in Table II.There is an excellent agreement among the different years and also with the published harmonic constituents by FEMAR, which is expected because of the constancy of the tidal generating forces and periodicity (Pugh 1987).The averaged maximum potential tidal amplitude was 1.82 m, which can be stated as the sum of the amplitudes of the main constituents shown in Table III.The M 2 accounted for 52.5% of the maximum potential amplitude, while the S 2 and the N 2 , the larger lunar elliptic, have about the same weight, 17.1 and 11.0%, respectively.The O 1 and K 1 accounted for less than 4% each.The averaged explained variance of the sea level due to the harmonic constituents was 99.7%.Similar results, in terms of the most important tidal constituents, were reported for the Amazon shelf by Beardsley et al. (1995).Table III shows these parameters derived from the tidal constituents.
it is noticeable in Figure 3 that there is a strong anomalistic modulation caused by the N 2 , which implies synodical inequalities.This means that there is a height difference between consecutive spring or neap tides.Figure 5 shows the temporal variability of the height difference between consecutive spring and neap tides, depicted as dots in Figure 3 (a2, b2 and c2).There is a polarity inversion between the maximum spring tide height differences and the neap tide ones, or, when the maximum difference occurs in the spring tides, the smallest is observed in the neap tides, and viceversa.The inversions occur with a period larger than 6 months.During the period analyzed, the maximum height differences of the spring tide occurred in January and July of 2009, February and September of 2010, and April and November of 2011.This effect is due to the nodal regression of the Moon, which has a period of 18.61 years (Pugh 1987).Taking the mean spring tide height from the morphometry (2.8 m) as reference, the spring tide height variability can be as high as 1 m, or 36% of the mean value.
The time series of along and cross-shore wind velocity components and the sea level at SubTF for 2009 are shown in Figure 6 (a, b and c).The maximum along, , and cross, , -shore wind velocity were 4.3 and 9.6 ms −1 , respectively, and the mean values were nearly null and 2.2 ms −1 , respectively.There is an annual modulation of the wind velocity.
The lowest velocity of and were recorded between April and June.The highest wind velocities were recorded between September and November.The sea level variability at SubTF was very low, with a maximum of 0.12 m.Therefore, there are still oscillations apparently not related to the wind.The wavelet power spectra of the sea level and the wind shear stress, and the wavelet coherence and phase spectrum between them are shown in Figure 6 (c, d and e).The x-axis represents time and the y-axis the period in 'days'.The side areas with transparency indicate the cone of influence (the inner part) where the analysis is valid.The color pattern indicates the energy, increasing to the red.A black line separates the regions where the levels of significance are higher than 5%.The level of significance of the coherence spectrum was estimated using the Monte Carlo method.The arrows in the coherence spectrum indicate the phase between the series.Arrows pointing to the right denote that the series are in phase, while arrows pointing to the left denote that the series are 180° out of phase.The coherence spectrum is the square of the smoothed cross spectrum, normalized with the smoothed wavelet spectrum.This procedure provides the interval of coherence between 0 and 1, where a value of one indicates maximum coherence (Torrence and Webster 1999).The application of the coherence and phase spectrum allows the visualization of the regions of higher covariance, in contrast to the cross-wavelet, which reveals common regions of higher energy between series (Truccolo 2009).
The sea level power spectrum shows a region of higher energy within the periods between 10 and 15 days for most of the time, and other scattered to lower time-periods during all the year.The wind stress power spectrum shows higher energy for the periods below 10 days, except for the period of April/May when the energy was the lowest, which is related to the lowest wind velocity (Figure 6c, d and e).There was an increase of energy at higher periods for the period of 10 to 20 days for the period between August and October.The only region of the coherence spectrum coinciding with higher energy in the sea level power spectrum was for March-April with periods of 10-15 days.Therefore, there is no significant energy in the wind power spectrum.

DisCussioN
The world's shores can be classified according to the tidal range, , as microtidal, mesotidal and macrotidal, when < 2.0 m, 2.0 < < 4.0 m, and > 4.0 m, respectively (Davies 1964), or, in a similar fashion, as microtidal, lower mesotidal, upper mesotidal, lower macrotidal and upper macrotidal, when <  1.0 m, 1.0 < < 2.0 m, 2.0 < < 3.5 m, 3.5 < < 5.5 m, and > 5.5 m, respectively (hayes 1979).These classifications are useful since they qualitatively weight the role of tides, and consequently the tidal energy controlling the coastal processes, providing a basis for coastal management.This is particularly relevant especially when one considers the scenarios of sea-level rise (e.g., Pethick 2001).
The basic approach to classify a given stretch of shore is to obtain a time series of sea level for a period of at least 30 days, from which the main harmonic tidal constituents can be determined, and from them, the basic descriptive parameters of the tide (Table III).These parameters provide a robust characterization of the tide in average terms, in error bounds of about 10% compared with observed values in the present case.Therefore, there is an intrinsic tidal variability that must be accounted to better describe the phenomenon.In the present study, we analyzed a time series of three years of sea level recorded in nearly open sea, which represent the shelf tides very well.Almost all the sea level variance was explained by tidal harmonic constituents, which were compared with tidal morphologic parameters.The tidal harmonic analysis is robust in obtaining coefficients from a time series, although it will never reproduce the original time series perfectly.
With the application of the traditional tidal descriptive parameters (Table iii), we obtained a mean tidal height of 2.05 m, which is about 7% larger than the value obtained from the tidal morphology.The calculated mean spring and neap tide height were about 10% smaller and larger, respectively, from those obtained from morphology.In spite of the relatively small difference between the calculated and observed values, this difference is, fortuitously in the present case, enough to change the regional shore classification.Even this is not necessarily an important issue regarding physical coastal processes the misclassification may suggest to someone that the tides are more important than it really is.
The variability around the mean spring and neap tide values is usually not considered.In the present case, the neap tides can be 52% higher or 36% lower than the mean, and the spring tides can be about 20% higher or lower than the mean value.As shown in Figure 7, the tidal height difference between consecutive spring tides can be as high as 1 m, and 0.7 between consecutive neap tides.The main cause of such variability is related to the importance of the N 2 , which is nearly equivalent to the S 2 .The N 2 represents the effects of the elliptic orbit of the Moon, with a period of 27.55 days.Thus, in the period of a lunar month, the distance between the Earth and Moon varies, and when the Moon is closest to the Earth (apogee), the tidal generating forces will be stronger, while in the opposite case (perigee), they will be weaker (Pugh 1987).
The diurnal inequalities can be important as well.Diurnal inequalities are due to the interactions between diurnal and semi-diurnal constituents, which are weighted by the tidal form number.The tidal constituents are not constant through the year, and consequently the form number changes too. Figure 7a shows the annual variability of the main diurnal (O 1 , K 1 ) and semi-diurnal (M 2 , S 2 and N 2 ) harmonic tidal constituents.The O 1 , K 1 , M 2 , and N 2 show a relatively small variability around their mean values, while the S 2 varies much more, between 0.2 and 0.4.The variability ratio, given by the ratio of the standard deviation and the mean, was 0.11, and due to the annual variation of the S 2 amplitude, it ranges from 0.08 in March to 0.13 in June (Figure 7b).
A practical aspect concerning the understanding of the tidal variability is related to the sampling strategies assessing the hydrodynamics and scalar transport processes in estuaries (e.g., Kjerfve 1990), especially when the continent-ocean biogeochemical fluxes are being estimated (e.g., Pereira-Filho et al. 2001, Schettini andMiranda 2010).The sampling strategy is usually based on sampling sea level, current velocity and direction, salinity, and several other scalars in a cross-section during at least one complete tidal cycle to obtain the instantaneous volume transport and the scalar mass transport.This approach furnishes the intratidal variability (e.g., Schettini et al. 2006), and the time-averaged values give the residual transport by elimination of the oscillatory character of the tidal variability.The latter is a very important result as it can highlight whether a system is a source or a trap for a specific scalar (Shubel andCarter 1984, Schettini et al. 2013), which can be extrapolated to obtain long-term trends.Therefore, the robustness of this result is strongly dependent on the sampling strategy.
Along shores subjected to pure semi-diurnal tides (FN < 0.4), it has been demonstrated that sampling during a single tidal cycle (~13 hours) is sufficient to obtain robust residual fluxes, since the semi-diurnal inequalities are small (Kjerfve 1990).This strategy has many logistic advantages, since an experiment can be performed during daylight (it is safer), can be run by a single team, and it is cheaper than a 25-hour or longer experiments.Therefore, especially in arid and semi-arid shores where the estuaries are under a low-inflow regime (e.g., Largier 2010) and the residual fluxes are due mainly to the tidal regime and are expected to be small, even small tidal inequalities can introduce a bias in the residual flux estimation.Studies aiming to obtain fluxes based on 13-hour surveys would achieve better results if sampling is conducted in May and October when the FN number and tidal inequalities are smaller (Figure 4).
Another important aspect associated with tidal experiments to estimate fluxes is that they must be performed at hydrological and tidal modal conditions.This means that experiments should be performed during the lower and higher periods of river inflow to the estuary, which depends on the annual hydrological regime (e.g., Pereira et al. 2010), and at spring and neap tides, which present current velocity extremes (Byun et al. 2013).Regarding the tides, usually the most extreme conditions are sought to visualize the wider spectrum of tidal effects, or, sampling is performed during the larger spring and smallest neap tide, consecutively (the order doesn't matter).These extremes occur during the equinoxes, at March and September, and do not coincide with the minimum FN periods.Additionally, the synodical inequalities should also be considered, since the range difference between consecutive spring/neap tides can be important, and will change on a yearly basis in a cycle of about 20 years.
Non-astronomical effects on the sea level are irrelevant at the study area. in spite of this, it is notable that there still is a SubTF that produces oscillations of 0.1 m in amplitude, and it is even more interesting that these oscillations have spectral energy concentrated in the 10-15 day period.The explanation of this variability is beyond the scope of the present paper, although some hypothesis can be put forth.The coherence spectrum between the local wind stress and sea level indicates that these oscillations are not due to the local response through Ekman mechanism, however, this result doesn't discard the idea that remote winds could have produced them (e.g., Truccolo et al. 2006).According to Shetye and Vijith (2013), there are two known mechanisms for generating sea level oscillations at the continental shelf: geostrophic adjustment of a remotely forced current and oscillations forced by local winds that could be ageostrophic.Since the continental shelf is relatively narrow along this stretch, it is possible that some sea level oscillations could be related to the continental shelf waves (due to remotely wind effect) and the variability of the North Brazil current-NBC (e.g., Stramma and Schott 1999).The NBC flowing westward at the continental margin produces rings (Philander and Pacanowski 1986), and the passage of these rings could affect somewhat the shelf sea level.At the shore of the Amazon, the synodical constituent MSf is relatively important (Gallo and Vinzon 2005).The SubTF oscillations recorded in Fortaleza (~ 1,000 km eastward) may be associated with both effects, although with the energy spanned from the harmonic frequency due to non-linear interactions with the NBC and shelf morphology.

CoNClusioNs
The tidal regime along the northern shelf of the Brazilian Northeast region can be classified as an intermediate state between micro and mesotidal (e.g., Davies 1964) or upper mesotidal regimes (e.g., hayes 1979).The tidal height ranges from 0.7 (the minimum neap tide height) to 3.3 m (the maximum spring tide height). in the present study, the use of a relatively simple morphology assessment over the observed and/or predicted (based on a 365-day analysis) data provided a useful insight into the harmonic tidal variability.
The morphology assessment from both, observations and harmonic prediction, furnish similar results in terms of the frequency distribution of height, although some scattering at higher tidal heights is present.The histograms of the neap and spring tide also do not correlate very well.
The descriptive tidal parameters from the harmonic constituents work well for the mean tidal range; however, there are overestimations/ underestimations of about 10% for the mean spring/ neap tidal height.
The S 2 annual variation produces a reflex in the form number, with the maximum semi-diurnal inequalities occurring during May and October.The N 2 synodical inequalities can be as large as 1 m between consecutive spring tides.
The non-harmonic sea level oscillations are negligible.Nevertheless, there is a persistent oscillation of about 0.1 m with energy concentrated in the 10-15 day period.
investigated the sea level variability in terms of the joint probability of tides and waves at the Fortaleza coast.This study aims to provide a comprehensive characterization of the sea level variability at the northern shore of the Brazilian Northeast through analysis of the sea level data recorded by the Permanent Tidal Monitoring Network for Geodetics of the Brazilian institute of Geography and Statistics (iBGE) at Fortaleza.The tide gauge station is located at Porto do Mucuripe, and reflects open sea conditions with very low shallow sea distortions usually observed

Figure 1 -
Figure 1 -(a) The northern Brazilian Northeast Continental Shelf and (b) details of the Fortaleza shore showing the tide gauge and meteorological station.

Figure 2 -
Figure 2 -Co-Tidal chart for the M 2 constituent in the tropical Atlantic Ocean.(a) Amplitude in centimeters and (b) phase in degrees to Greenwich.Adapted from Beardsley et al. 1995.

Figure 3 -
Figure 3 -Observed hourly sea level time series (a1, b1 and c1) and respective tidal height (a2, b2 and c2) at Fortaleza for the years 2009 (a), 2010 (b) and 2011 (c).The black dots at tidal height time series indicate the extreme higher and lower spring and neap tide heights.

Figure 4 -
Figure 4 -(a) Relation between the observed and predicted tidal heights.(b) Frequency distribution of observed (full line) and predicted (dashed line) tidal height.(c) Frequency distribution of observed (full line) and predicted (dashed line) of neap and spring tidal height.

Figure 5 -
Figure 5 -Temporal variability between consecutive neap (full line) and spring (dashed line) tides for the years 2009, 2010 and 2011.

Figure 6 -
Figure 6 -Time series of sub-tidal signals of (a) alongshore and (b) cross-shore wind velocity components and (c) sea level.Wavelet power spectra (PS) of (d) sea level (SL) and (e) alongshore wind stress (LWS), and (f) wavelet coherence and phase spectrum (CS) between wind stress and sea level.

Figure 7 -
Figure 7 -(a) Monthly variation of the amplitude of the main diurnal (O 1 , K 1 ) and semidiurnal (M 2 , S 2 and N 2 ) tidal constituents averaged for the years 2009, 2010 and 2011.(b) Monthly variation of the tidal form number anomaly (full line: 2009; dashed line: 2010; dotted line: 2011).

aCKNoWleDGMeNTs
We are thankful to the instituto Brasileiro de Geografia e Estatística (iBGE) for making the sea level data freely available; to the Universidade Federal do Ceará (UFC) for the meteorological data.FFF has a Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) scholarship; ECT has a CNPq and Fundação de Amparo à Ciência e Tecnologia do Estado de Pernambuco (FACEPE) DCR scholarship; CAFS has a CNPq scholarship PQ 308665/2013-9.

TaBle iii Descriptive tidal parameters derived from harmonic tidal constituents (Values-h) compared with the values derived from morphology (Values-M). FN: form number; Mh: mean tidal height; Msh: mean spring tide height; MNh: mean neap tide height.
FELiPE F. FROTA, ELiANE C. TRUCCOLO and CARLOS A.F. SChETTiNi