Tide model comparison over the Southwestern Atlantic Shelf

Sea surface height (SSH) as measured by satellites has become a powerful tool for oceanographic and climate related studies. Whereas in the open ocean good accuracy has been achieved, more energetic dynamics and a number of calibration problems have limited applications over continental shelves and near the coast. Tidal ranges in the Southwestern Atlantic (SWA) continental shelf are among the highest in the world ocean, reaching up to 12 m at specific locations. This fact highlights the relevance of the accuracy of the tidal correction that must be applied to the satellite data to be useful in the region. In this work, amplitudes and phases of tidal constituents are extracted from five global tide models and three regional models and compared to the corresponding harmonics estimated from coastal tide gauges (TGs) and satellite altimetry data. The Root Sum Square (RSS) of the misfit of the common set of the five tidal constituents solved by the models (M₂, N₂, S₂, K₁ and O₁) is higher than 18 cm close to the coast for two of the regional models and higher than 24.5 cm for the rest of the models considered. Both values are too high to provide an accurate estimation of geostrophic non-tidal currents from satellite altimetry in the coastal region. On the other hand, the global model with the highest spatial resolution has a RSS lower than 4.5 cm over the continental shelf even when the non-linear M₄ overtide is considered. Comparison with in-situ current measurements suggests that this model can be used to de-tide altimetry data to compute large-scale patterns of SSH and associated geostrophic velocities. It is suggested that a local tide model with very high resolution that assimilates in-situ and satellite data should meet the precision needed to estimate geostrophic velocities at a higher resolution both close to the coast and over the Patagonian shelf.     

Tephrostratigraphy of the last 170 ka in sedimentary successions from the Adriatic Sea

In this study are discussed new SEM-EDS analyses performed on glass shards from five cores collected in the Central Adriatic Sea and two cores recovered from the South Adriatic Sea. A total of 26 tephra layers have been characterized and compared with the geochemical features of terrestrial deposits and other tephra archives in the area (South Adriatic Sea and Lago Grande di Monticchio, Vulture volcano). The compositions are compatible with either a Campanian or a Roman provenance. The cores, located on the Central Adriatic inner and outer shelf, recorded tephra referred to explosive events described in the literature: AP3 (sub-Plinian activity of the Somma-Vesuvius, 2710 ± 60 ¹⁴C years BP); Avellino eruption (Somma–Vesuvius, 3548 ± 129 ¹⁴C years BP); Agnano Monte Spina (Phlegrean Fields, 4100 ± 400 years BP); Mercato eruption (Somma–Vesuvius, 8010 ± 35 ¹⁴C years BP; Agnano Pomici Principali eruption (Phlegrean Fields, 10,320 ± 50 ¹⁴C years BP); Neapolitan Yellow Tuff (Phlegrean Fields, 12,100 ± 170 ¹⁴C years BP). Some of these layers were also observed in the South Adriatic core IN68-9 in addition to younger (AP2, sub-Plinian eruption, Somma–Vesuvius, 3225 ± 140 ¹⁴C years BP), and older layers (Pomici di Base eruption, Somma–Vesuvius, 18,300 ± 150 ¹⁴C years BP). Significant is the tephra record of core RF95-7 that, for the first time in the Adriatic Sea, reports the occurrence of tephra layers older than 60 ka: the well known Mediterranean tephra layers X2 (ca. 70 ka), W1 (ca. 140 ka) and V2 (Roman origin, ca. 170 ka) as well as other tephra layers attributed, on the basis of geochemistry and biostratigraphy, to explosive eruptions occurred at Vico (138 ± 2 and 151 ± 3 ka BP) and Ischia (147–140 ka BP).     

Ingestion of marine debris by loggerhead sea turtles, Caretta caretta, in the Adriatic Sea

We examined the occurrence of marine debris in the gastrointestinal tract of 54 loggerhead sea turtles (Caretta caretta) found stranded or incidentally captured dead by fisheries in the Adriatic Sea, with a curved carapace length of 25.0-79.2 cm. Marine debris was present in 35.2% of turtles and included soft plastic, ropes, Styrofoam and monofilament lines found in 68.4%, 42.1%, 15.8% and 5.3% of loggerheads that have ingested debris, respectively. The dry mass of debris per turtle was low, ranging from <0.01 to 0.71 g, and the ingestion was not significantly affected by sex or body size (all p > 0.05). Marine debris averaged 2.2 ± 8.0% of dry mass of gut content, with a maximum of 35% found in a juvenile turtle that most likely died due to debris ingestion. Considering the relatively high occurrence of debris intake and possible sub-lethal effects of even small quantities of marine debris, this can be an additional factor of concern for loggerheads in the Adriatic Sea.     

Variability of baroclinic tidal currents on the Mackenzie Shelf,the Southeastern Beaufort Sea

Semidiurnal tidal currents on the outer shelf of the Mackenzie Shelf in the Beaufort Sea were found to be strongly influenced by the locally generated baroclinic tide. Two primary factors are involved in this process: (1) the sharp shelf break along the northeastern Mackenzie Shelf, promoting the generation of vigorous internal tidal waves; and (2) the proximity to critical latitudes for M₂ and N₂ motions locking these waves and preventing them from leaving the source region. As a result, internal tides are resonantly trapped between the shelf and critical latitudes. The physical properties and temporal variations of tidal motions were examined using current meter measurements obtained from 1987–1988 at four sites (SS1, SS2, SS3, and SS4) offshore of the shelf break at depths of ∼200 m. Each mooring had Aanderaa RCM4s positioned at ∼35 m below the surface and ∼50 m above the bottom. Complex demodulation was used to compute the envelopes (amplitude modulation) of these components. A striking difference in the variability of clockwise (CW) and counterclockwise (CCW) tidal currents was found. The CW tides are highly variable, have greater amplitude, exhibit a burst-like character associated with wind events and contain about 80% of the total energy of the semidiurnal tidal currents. In contrast, the CCW components have a more regular temporal regime with distinct monthly, fortnightly and 10-day modulation at astronomical periodicities associated with frequency differences M₂–N₂ (0.03629 cpd), S₂–M₂ (0.06773 cpd), and S₂–N₂ (0.10402 cpd). Significant horizontal correlation of the CW current envelopes was found only between stations near the northeast Mackenzie Shelf, indicating this to be the main area of baroclinic internal wave generation.     

The oceanic tides in the South Atlantic Ocean

The finite element ocean tide model of Le Provost and Vincent (1986) has been applied to the simulation of the M₂ and K₁ components over the South Atlantic Ocean. The discretisation of the domain, of the order of 200 km over the deep ocean, is refined down to 15 km along the coasts, such refinement enables wave propagation and damping over the continental shelves to be correctly solved. The marine boundary conditions, from Dakar to Natal, through the Drake passage and from South Africa to Antarctica, are deduced from in situ data and from Schwiderski’s solution and then optimised following a procedure previously developed by the authors. The solutions presented are in very good agreement with in situ data: the root mean square deviations from a standard subset of 13 pelagic stations are 1.4 cm for M₂ and 0.45 cm for K₁, which is significantly better overall than solutions published to date in the literature. Zooms of the M₂ solution are presented for the Falkland Archipelago, the Weddell Sea and the Patagonian Shelf. The first zoom allows detailing of the tidal structure around the Falklands and its interpretation in terms of a stationary trapped Kelvin wave system. The second zoom, over the Weddell Sea, reveals for the first time what must be the tidal signal under the permanent ice shelf and gives a solution over that sea which is generally in agreement with observations. The third zoom is over the complex Patagonian Shelf. This zoom illustrates the ability of the model to simulate the tides, even over this area, with a surprising level of realism, following purely hydrodynamic modelling procedures, within a global ocean tide model. Maps of maximum associated tidal currents are also given, as a first illustration of a by-product of these simulations.     

Tides and mixing in the northwestern East China Sea Part I: Rotating and reversing tidal flows

Bottom-mounted ADV and ADCP instruments in combination with CTD profiling measurements taken along the Chinese coast of the East China Sea were used to study the vertical structure of temperature, salinity, and velocity in reversing tidal currents on a shallow inner shelf and in rotating tidal flows over a deeper sloping bottom of the outer shelf. These two regimes of barotropic tide affect small-scale dynamics in the lower part of the water column differently. The reversing flow was superimposed by seiches of ∼2.3 h period generated in semienclosed Jiaozhou Bay located nearby. As the tidal vector rotates over the sloping bottom, the height of the near-bottom logarithmic layer is subjected to tidal-induced variations. A maximum of horizontal velocity U_max appears at the upper boundary of the log layer during the first half of the current vector rotation from the minor to the major axis of tidal ellipse. In rotating tidal flow, vertical shear generated at the seafloor, propagated slowly to the water interior up to the height of U_max, with a phase speed of ∼5 m/h. The time-shifted shear inside the water column, relative to the shear at the bottom, was associated with periodically changing increases and decreases of the tidal velocity above the log layer toward the sea surface. In reversing flows, the shear generated near the bottom and the shear at the upper levels were almost in phase.     

The tidal and wind induced hydrodynamics of the composite system Adriatic Sea/Lagoon of Venice

A shallow water hydrostatic 2D hydrodynamic numerical model, based on the boundary conforming coordinate system, was used to simulate aspects of both general and small scale oceanic features occurring in the composite system constituted by the Adriatic Sea and the Lagoon of Venice (Italy), under the influence of tide and realistic atmospheric forcing. Due to a specific technique for the treatment of movable lateral boundaries, the model is able to simulate efficiently dry up and flooding processes within the lagoon. Firstly, a model calibration was performed by comparing the results of the model, forced using tides and ECMWF atmospheric pressure and wind fields, with observations collected for a set of 33 mareographic stations uniformly distributed in the Adriatic Sea and in the Lagoon of Venice. A second numerical experiment was then carried out by considering only the tidal forcing. Through a comparison between the results obtained in the two experiments it was possible to assess the reliability of the estimated parameter through the composite forcing. Model results were then verified by comparing simulated amplitude and phase of each tidal constituent as well as tidal velocities simulated at the inlets of the lagoon and in the Northern Adriatic Sea with the corresponding observed values. The model accurately reproduces the observed harmonics: mean amplitude differences rarely exceed 1 cm, while phase errors are commonly confined below 15°. Semidiurnal and diurnal currents were correctly reproduced in the northern basin and a good agreement was obtained with measurements carried out at the lagoon inlets. On this basis, the outcomes of the hydrodynamic model were analyzed in order to investigate: (i) small-scale coastal circulation features observed at the interface between the adjoining basins, which consist often of vortical dipoles connected with the tidal flow of Adriatic water entering and leaving the Lagoon of Venice and with along-shore current fields connected with specific wind patterns; (ii) residual oscillations, which are often connected to meteorological forcing over the basin. In particular, it emerges that small-scale vortical features generated near the lagoon inlet can be efficiently transported toward the open sea, thus contributing to the water exchange between the two marine regions, and a realistic representation of observed residual oscillations in the area would require a very detailed knowledge of atmospheric as well as remote oceanic forcing.     

An inter-comparison of tidal solutions computed with a range of unstructured grid models of the Irish and Celtic Sea Regions

Three finite element codes, namely TELEMAC, ADCIRC and QUODDY, are used to compute the spatial distributions of the M₂, M₄ and M₆ components of the tide in the sea region off the west coast of Britain. This region is chosen because there is an accurate topographic dataset in the area and detailed open boundary M₂ tidal forcing for driving the model. In addition, accurate solutions (based upon comparisons with extensive observations) using uniform grid finite difference models forced with these open boundary data exist for comparison purposes. By using boundary forcing, bottom topography and bottom drag coefficients identical to those used in an earlier finite difference model, there is no danger of comparing finite element solutions for “untuned unoptimised solutions” with those from a “tuned optimised solution”. In addition, by placing the open boundary in all finite element calculations at the same location as that used in a previous finite difference model and using the same M₂ tidal boundary forcing and water depths, a like with like comparison of solutions derived with the various finite element models was possible. In addition, this open boundary was well removed from the shallow water region, namely the eastern Irish Sea where the higher harmonics were generated. Since these are not included in the open boundary, forcing their generation was determined by physical processes within the models. Consequently, an inter-comparison of these higher harmonics generated by the various finite element codes gives some indication of the degree of variability in the solution particularly in coastal regions from one finite element model to another. Initial calculations using high-resolution near-shore topography in the eastern Irish Sea and including “wetting and drying” showed that M₂ tidal amplitudes and phases in the region computed with TELEMAC were in good agreement with observations. The ADCIRC code gave amplitudes about 30 cm lower and phases about 8° higher. For the M₄ tide, in the eastern Irish Sea amplitudes computed with TELEMAC were about 4 cm higher than ADCIRC on average, with phase differences of order 5°. For the M₆ component, amplitudes and phases showed significant small-scale variability in the eastern Irish Sea, and no clear bias between the models could be found. Although setting a minimum water depth of 5 m in the near-shore region, hence removing wetting and drying, reduced the small-scale variability in the models, the differences in M₂ and M₄ tide between models remained. For M₆, a significant reduction in variability occurred in the eastern Irish Sea when a minimum 5-m water depth was specified. In this case, TELEMAC gave amplitudes that were 1 cm higher and phases 30° lower than ADCIRC on average. For QUODDY in the eastern Irish Sea, average M₂ tidal amplitudes were about 10 cm higher and phase 8° higher than those computed with TELEMAC. For M₄, amplitudes were approximately 2 cm higher with phases of order 15° higher in the northern part of the region and 15° lower in the southern part. For M₆ in the north of the region, amplitudes were 2 cm higher and about 2 cm lower in the south. Very rapid M₆ tidal-phase changes occurred in the near-shore regions. The lessons learned from this model inter-comparison study are summarised in the final section of the paper. In addition, the problems of performing a detailed model–model inter-comparison are discussed, as are the enormous difficulties of conducting a true model skill assessment that would require detailed measurements of tidal boundary forcing, near-shore topography and precise knowledge of bed types and bed forms. Such data are at present not available.     

A numerical study of the barotropic tides and tidal energy distribution in the Indonesian seas with the assimilated finite volume coastal ocean model

The tides and tidal energetics in the Indonesian seas are simulated using a three-dimensional finite volume coastal ocean model. The high-resolution coastline-fitted model is configured to better resolve the hydrodynamic processes around the numerous barrier islands. A large model domain is adopted to minimize the uncertainty adjacent to open boundaries. The model results with elevation assimilation based on a simple nudge scheme faithfully reproduced the general features of the barotropic tides in the Indonesian Seas. The mean root-mean-square errors between the observed and simulated tidal constants are 2.3, 1.1, 2.4, and 1.5 cm for M₂, S₂, K₁, and O₁, respectively. Analysis of the model solutions indicates that the semidiurnal tides in the Indonesian Seas are primarily dominated by the Indian Ocean, whereas the diurnal tides in this region are mainly influenced by the Pacific Ocean, which is consistent with previous studies. Examinations of tidal energy transport reveal that the tidal energy for both of the simulated tidal constituents are transported from the Indian Ocean into the IS mainly through the Lombok Strait and the Timor Sea, whereas only M₂ energy enters the Banda Sea and continues northward. The tidal energy dissipates the most in the passages on both sides of Timor Island, with the maximum M₂ and K₁ tidal energy transport reaching about 750 and 650 kW m^–1, respectively. The total energy losses of the four dominant constituents in the IS are nearly 338 GW, with the M₂ constituent dissipating 240.8 GW. It is also shown that the bottom dissipation rate for the M₂ tide is about 1–2 order of magnitudes larger than that of the other three tidal components in the Indonesian seas.     

Biomonitoring of organochlorine compounds and trace metals along the Eastern Adriatic coast (Croatia) using Mytilus galloprovincialis

The paper presents the results of monitoring trace organic contaminants (17 PCBs and 7 OCPs) and 6 trace metals (Cd, Cr, Cu, Hg, Pb and Zn) in the coastal waters of the eastern Adriatic in 2006 using Mytilus galloprovincialis as indicator species. OCPs were found in all samples (0.07-14.3 ng g⁻¹ d.wt.), while in positive samples concentrations of PCBs ranged between 0.16 and 20.5 ng g⁻¹ d.wt. The α-HCH/γ-HCH and DDE/DDT ratios indicated recent input of γ-HCH and DDT. Spatial distributions of HCB, PCBs, DDTs and Hg were significantly correlated with population density. Although elevated concentrations of all trace metals and organic contaminants were recorded in the densely populated and industrialized areas, overall data point to relatively low level of contamination of the eastern Adriatic coast. Evaluation of the risks to human health associated with consumption of the mussels suggested that there is no health risk for moderate shellfish consumers.     

Ostreopsis cf. ovata bloom in the northern Adriatic Sea during summer 2009: ecology, molecular characterization and toxin profile

Intense blooms of the benthic dinoflagellate Ostreopsis cf. ovata have occurred in the northern Adriatic Sea since 2006. These blooms are associated with noxious effects on human health and with the mortality of benthic organisms because of the production of palytoxin-like compounds. The O. cf. ovata bloom and its relationships with nutrient concentrations at two stations on the Conero Riviera (northern Adriatic Sea) were investigated in the summer of 2009. O. cf. ovata developed from August to November, with the highest abundances in September (1.3 × 10⁶ cells g⁻¹ fw corresponding to 63.8 × 10³ cells cm⁻²). The presence of the single O. cf. ovata genotype was confirmed by a PCR assay. Bloom developed when the seawater temperature was decreasing. Nutrient concentrations did not seem to affect bloom dynamics. Toxin analysis performed by high resolution liquid chromatography-mass spectrometry revealed a high total toxin content (up to 75 pg cell⁻¹), including putative palytoxin and all the ovatoxins known so far.     

River water and nutrient discharges in the Northern Adriatic Sea: Current importance and long term changes

Runoff and nutrient transport by rivers were analysed in the Northern Adriatic continental shelf, in order to evaluate their interannual and multidecal variability, as well as their current contribution to determine freshwater and nutrient budgets in this marine region. During the years 2004-2007, the runoff in the basin (34.1-64.6 km³ yr⁻¹) was highly imbalanced, being 84% of freshwater discharged along the western coast, because of the contributions of Po, Adige and Brenta rivers. In the northern and eastern sections of the coast, freshwater discharge by rivers was less important (10 and 6%, respectively), but not negligible in determining the oceanographic properties at sub-regional scales. The oscillations of the transport of biogenic elements (124-262×10³ t N yr⁻¹ for TN, 72-136×10³ t N yr⁻¹ for DIN, 4.5-11.1×10³t P yr⁻¹ for TP, 2.2-3.5×10³ t P yr⁻¹ for PO₄ and 104-196×10³ t Si yr⁻¹ for SiO₂) were strictly dependant to the differences in the annual runoff. A strong excess of N load in comparison to P load characterised all rivers, both in inorganic nutrient (DIN/PO4=37-418) and total (TN/TP=48-208) pools, particularly in the northern and eastern areas of the basin.The annual runoff showed significant oscillations for Po on multidecadal time scale, whereas a general decrease (−33%) was observed for the other N Adriatic rivers as the recent discharges were compared to those before the 1980s. During the dry years 2005-2007, a strong reduction of river water flows and nutrient loads was experienced by the N Adriatic ecosystem with respect to years characterised by medium-high regimes. An increased frequency of similar drought periods, due to ongoing climate changes or to a larger human usage of continental waters, would be easily able to significantly change the biogeochemistry of this basin.     

On the effect of the sampling frequency of sea level measurements on return period estimate of extremes—Southern European examples

Estimates of extreme sea levels and return periods have been based mainly on hourly sampling rates. Technological development has enabled the sampling rates to increase and sampling rates of 5–10 min are becoming increasingly common. In this paper we explore the relationship between extreme sea levels and estimated return periods based on hourly and shorter sampling periods in three tide-gauges one at the Atlantic coasts of Spain (Coruña), one in the western Mediterranean (Malaga) and one in the N. Adriatic (Trieste). Significant differences of several centimetres are found in the hourly and 5 min extremes. These reflect in significant underestimation of the 50-year return levels which in Trieste reach 38 cm. A theoretical relationship between the high and the low sampling rate of extremes is also tested. Thus updated 50-year return levels for the Mediterranean and the coasts of the Iberian peninsula are produced assuming that the differences identified in the various stations generalise to other tide-gauge (hourly) records for which hourly values have been analysed earlier.     

Fortnightly oscillations observed in the Adriatic Sea

This paper documents the occurrence of strong fortnightly oscillations in the Adriatic Sea frequently observed in the current measurements. For that purpose, we analyzed half-decadal sea level series and long-term currents collected within different parts and layers of the Adriatic Sea. Harmonic analysis and band-pass filter with cutoff periods at 10 and 20 days have been applied to the series. The sea level Mf tide is found to be the only significant tidal constituent over periods between the diurnal and semi-annual ones, having no significant phase variations in the Adriatic Sea. Moreover, the currents on periods between 10 and 20 days could only partially be explained by atmospheric forcing. The fact that they are in-phase with the Mf tide gives a hint on a possible mechanism responsible for the generation of these oscillations. Additionally, these currents are usually strongly amplified during a weakly stratified season, indicating a baroclinic mechanism responsible for the multiplication of the fortnightly currents.     

Inter-annual and seasonal variations in the air–sea CO2 balance in the central Baltic Sea and the Kattegat

We estimated the net annual air–sea exchange of carbon dioxide (CO₂) using monitoring data from the East Gotland Sea, Bornholm Sea, and Kattegat for the 1993–2009 period. Wind speed and the sea surface partial pressure of CO₂ (pCO₂^w), calculated from pH, total alkalinity, temperature, and salinity, were used for the flux calculations. We demonstrate that regions in the central Baltic Sea and the Kattegat alternate between being sinks (−) and sources (+) of CO₂ within the −4.2 to +5.2 mol m⁻² yr⁻¹ range. On average, for the 1994–2008 period, the East Gotland Sea was a source of CO₂ (1.64 mol m⁻² yr⁻¹), the Bornholm Sea was a source (2.34 mol m⁻² yr⁻¹), and the Kattegat was a sink (−1.16 mol m⁻² yr⁻¹). Large inter-annual and regional variations in the air–sea balance were observed. We used two parameterizations for the gas transfer velocity (k) and the choice varied the air–sea exchange by a factor of two. Inter-annual variations in pCO₂^w between summers were controlled by the maximum concentration of phosphate in winter. Inter-annual variations in the CO₂ flux and gas transfer velocity were larger between winters than between summers. This indicates that the inter-annual variability in the total flux was controlled by winter conditions. The large differences between the central Baltic Sea and Kattegat were considered to depend partly on the differences in the mixed layer depth.     

Air-sea CO2 fluxes in the southern Yellow Sea: An examination of the continental shelf pump hypothesis

The southern Yellow Sea (SYS), located to the north of the East China Sea (ECS), was considered part of the ECS when Tsunogai et al. (1999) proposed the “continental shelf pump” (CSP) hypothesis. However, the original CSP carbon dioxide (CO₂) uptake flux (2.9 mol C m⁻² yr⁻¹) appears to have been overestimated, primarily due to the differences between the SYS and the ECS in terms of their CO₂ system. In this paper, we estimated air-sea CO₂ fluxes in the SYS using the surface water partial pressure of CO₂ (pCO₂) measured in winter, spring, and summer, as well as that estimated in fall via the relationship of pCO₂ with salinity, temperature, and chlorophyll a. The results indicate that overall, the entire investigated area was a net source of atmospheric CO₂ during summer, winter, and fall, whereas it was a net sink during spring. Spatially, the nearshore area was almost a permanent CO₂ source, while the central SYS shifted from being a CO₂ sink in spring to a source in the other seasons of the year. Overall, the SYS is a net source of atmospheric CO₂ on an annual scale, releasing ∼7.38 Tg C (1 Tg=10¹² g) to the atmosphere annually. Thus, the updated CO₂ uptake flux in the combined SYS and ECS is reduced to ∼0.86 mol C m⁻² yr⁻¹. If this value is extrapolated globally following Tsunogai et al. (1999), the global continental shelf would be a sink of ∼0.29 Pg C yr⁻¹, instead of 1 Pg C yr⁻¹ (1 Pg=10¹⁵ g).The SYS as a net annual source of atmospheric CO₂ is in sharp contrast to most mid- and high-latitude continental shelves, which are CO₂ sinks. We argue that unlike the ECS and the North Sea where carbon on the shelf could be exported to the open ocean, the SYS lacks the physical conditions required by the CSP to transport carbon off the shelf effectively. The global validity of the CSP theory is thus questionable.     

Modeling studies of tidal dynamics and the associated responses to coastline changes in the Bohai Sea,China

Over the past 30 years, reclamation projects and related changes have impacted the hydrodynamics and sediment transport in the Bohai Sea. Three-dimensional tidal current models of the Bohai Sea and the Yellow Sea were constructed using the MIKE 3 model. We used a refined grid to simulate and analyze the effects of changes in coastline, depth, topography, reclamation, the Yellow River estuary, and coastal erosion on tidal systems, tide levels, tidal currents, residual currents, and tidal fluxes. The simulation results show that the relative change in the amplitude of the half-day tide is greater than that of the full-day tide. The changes in the tidal amplitudes of M₂, S₂, K₁, and O₁ caused by coastline changes accounted for 27.76–99.07% of the overall change in amplitude from 1987 to 2016, and water depth changes accounted for 0.93–72.24% of the overall change. The dominant factor driving coastline changes is reclamation, accounting for 99.55–99.91% of the amplitude changes in tidal waves, followed by coastal erosion, accounting for 0.05–0.40% of the tidal wave amplitude changes. The contribution of changes in the Yellow River estuary to tidal wave amplitude changes is small, accounting for 0.01–0.12% of the amplitude change factor. The change in the highest tide level (HTL) is mainly related to the amplitude change, and the correlation with the phase change is small. The dominant factor responsible for the change in the HTL is the tide amplitude change in M2, followed by S2, whereas the influence of the K1 and O1 tides on the change in the HTL is small. Reclamation resulted in a decrease in the vertical average maximum flow velocity (V_VAM) in the Bohai Sea. Shallower water depths have led to an increase in the V_VAM; deeper water depths have led to a decrease in the maximum flow velocity. The absolute value of the maximum flow velocity gradually decreases from the surface to the bottom, but the relative change value is basically constant. The changes in the tidal dynamics of the Bohai Sea are proportional to the degree of change in the coastline. The maximum and minimum changes in the tidal flux appear in Laizhou Bay (P-LZB) and Liaodong Bay (P-LDB), respectively. The changes in the tidal flux are related to the change in the area of the bay. Due to the reduced tidal flux, the water exchange capacity of the Bohai Sea has decreased, impacting the ecological environment of the Bohai Sea. Strictly controlling the scale of reclamation are important measures for reducing the decline in the water exchange capacity of the Bohai Sea and the deterioration of its ecological environment.     

Vertical Crustal Motion along the Mediterranean and Black Sea Coast Derived from Ocean Altimetry and Tide Gauge Data

Tide gauge (TG) data along the northern Mediterranean and Black Sea coasts are compared to the sea-surface height (SSH) anomaly obtained from ocean altimetry (TOPEX/Poseidon and ERS-1/2) for a period of nine years (1993–2001). The TG measures the SSH relative to the ground whereas the altimetry does so with respect to the geocentric reference frame; therefore their difference would be in principle a vertical ground motion of the TG sites, though there are different error sources for this estimate as is discussed in the paper. In this study we estimate such vertical ground motion, for each TG site, from the slope of the SSH time series of the (non-seasonal) difference between the TG record and the altimetry measurement at a point closest to the TG. Where possible, these estimates are further compared with those derived from nearby continuous Global Positioning System (GPS) data series. These results on vertical ground motion along the Mediterranean and Black Sea coasts provide useful source data for studying, contrasting, and constraining tectonic models of the region. For example, in the eastern coast of the Adriatic Sea and in the western coast of Greece, a general subsidence is observed which may be related to the Adriatic lithosphere subducting beneath the Eurasian plate along the Dinarides fault.     

Spatial distribution of polyfluoroalkyl compounds in dab (Limanda limanda) bile fluids from Iceland and the North Sea

The spatial distribution of polyfluoroalkyl compounds (PFCs) was investigated in dab (Limanda limanda) bile fluids collected from Iceland and the North Sea. Concentrations of various PFCs, including perfluorinated sulfonates (C₄-C₆, C₈ PFSAs), perfluorinated carboxylic acids (C₉-C₁₄ PFCAs) and n-methyl perfluorooctane sulfonamidoethanol (MeFOSE), were quantified. Perfluorooctane sulfonate (PFOS) was the predominant compound with highest concentrations along the Danish and German coast (mean 9.36 ng/g wet weight (ww)). Significantly lower PFOS concentrations were found at the other sampling stations in the North Sea and Iceland (p < 0.01, t-test). Conversely, the spatial distribution of the PFCAs in Iceland and the North Sea was more uniform. The most abundant PFCA was perfluorononanoic acid (PFNA), while the mean concentration decreased with increasing chain length from 4.7 ng/g ww for PFNA to 0.04 ng/g ww for perfluorotetradecanoic acid (PFTeDA). Overall, the different spatial distribution of PFCs indicates different origin of sources and different transportation mechanism.