Numerical modeling of tidal currents, sediment transport and morphological evolution in Hangzhou Bay, China

A 2D depth-averaged numerical model is set up to simulate the macro-scale hydrodynamic characteristics, sediment transport patterns and morphological evolution in Hangzhou Bay, a large macro-tidal estuary on the eastern coast of China. By incorporating the shallow water equations, the suspended sediment transport equation and the mass-balance equation for sediment; short-term hydrodynamics, sediment transport and long-term morphological evolution for Hangzhou Bay are simulated and the underlying physical mechanisms are analyzed. The model reproduces the spatial distribution patterns of suspended sediment concentration (SSC) in Hangzhou Bay, characterized by three high SSC zones and two low SSC zones. It also correctly simulates the residual flow, the residual sediment transport and the sediment accumulation patterns in Hangzhou Bay. The model results are in agreement with previous studies based on field measurements. The residual flow and the residual sediment transport are landwards directed in the northern part of the bay and seawards directed in the southern part. Sediment accumulation takes place in most areas of the bay. Harmonic analysis revealed that the tide is flood-dominant in the northern part of the bay and ebb-dominant in the southern part of the bay. The strength of the flood-dominance increases landwards along the northern Hangzhou Bay. In turn sediment transport in Hangzhou Bay is controlled by this tidal asymmetry pattern. In addition, the direction of tidal propagation in the East China Sea, the presence of the archipelago in the southeast and the funnel-shaped geometry of the bay, play important roles for the patterns of sediment transport and sediment accumulation respectively.     

三门湾沿海声层析潮流观测实验

2009年9月6日至9日在三门湾进行了沿海声层析(Coastal Acoustic Tomography,CAT)潮流观测实验.实验由7台沿海声层析仪组网进行,并分别由渔船定点抛锚于7个站位.实验期间,还进行了定点ADCP(Acoustic Doppler Current Profiler)观测.通过建立逆模式对声传播时间差进行解析,引入权重因子,用L-curve法确定阻尼因子的最佳值,继而根据阻尼最小二乘法得到流速的最佳解.根据逆模式得到的流速分布可知该海区的潮流以半日潮(M2)为主,M2潮流椭圆呈东南-西北走向,潮流基本都是顺着水道,即涨潮为西北流向,退潮为东南流向.西北向与东南向最大流速分别为1.03 m·s-1和1.09 m·s-1.实验期间该区域的余流是从湾外流入湾内,平均流速约为0.05 m·s-1.CAT与定点ADCP流速的东分量和北分量的均方差均小于0.18 m·s-1.这样大面积的潮流和余流水平分布的同步观测,用传统观测手段很难实现.通过以上结果可以得出,沿海声层析技术可以作为一种新的测流方法对强潮海区进行大面积潮流观测,可在我国沿海的海洋环境监测等方面发挥重要作用.     

Hemolytic activity and fatty acids composition in the ichthyotoxic dinoflagellate Cochlodinium polykrikoides isolated from Bahía de La Paz, Gulf of California

The hemolytic activity of the dinoflagellate Cochlodinium polykrikoides from Bahía de La Paz, Gulf of California was investigated as part of the ichthyotoxic mechanism of this microalga. Two different kinds of erythrocytes, fish and human, were tested for the hemolytic assay. Since fatty acids have been associated with hemolytic activity in C. polykrikoides, the composition of fatty acids of this dinoflagellate was also analyzed. The concentration of C. polykrikoides causing 50% hemolysis (HE₅₀) was 4.88 and 5.27 × 10⁶ cells L⁻¹, for fish and human erythrocytes, respectively. According to the standard curve of saponin, an equivalence between the hemolytic activity of saponin and the dinoflagellate concentration was found with 1 μg saponin mL⁻¹ equivalent to 1 × 10⁶ cells L⁻¹ of C. polykrikoides. The polyunsaturated fatty acids: hexadecaenoic (16:0), docosahexaenoic (22:6 n3) and octadecapentaenoic (18:5 n3) were found in an abundance of 62% of total fatty acids.     

Active segment of the 12 November 2003 Mw 5.6 earthquake at Salsipuedes oceanic basin,Gulf of California. Mexico

We analyzed records of eight seismic stations of the autonomous broadband seismograph network of a joint project between Utrecht University (the Netherlands), California Institute of Technology, and Centro de Investigación Científica y de Estudios Superiores de Ensenada (CICESE). These stations recorded the Mw 5.6 earthquake that occurred on 12 November 2003 at Salsipuedes basin in the middle of the Gulf of California 2 km west of the island Angel de la Guarda. This event was located at 29.16º N and 113.37º W, 30 km northeast of Bahia de los Angeles. A foreshock and hundreds of aftershocks were recorded in the 48 hours after its origin time. With the location of 29 earthquakes we identified the active segment, perpendicular to the main transform fault NW–SE of Canal de Ballenas, representing the transtensional boundary between the Pacific and North American plates. The direction of the active fault described is consistent with the normal fault mechanism reported by the National Earthquake Information Center (strike=39º, dip=34º, slip=–44º).From the duration magnitude of 456 aftershocks, we calculated a b-value of 1.14±0.28; furthermore, we calculated a seismic moment of (3.5 ±3.3) X10¹⁷Nm, a source radius of 3.7 ± 2.63 km, and a static stress drop of 3.94 ± 1.15 MPa (39.4 ± 11.5 bar.).     

Benthic foraminifera show some resilience to ocean acidification in the northern Gulf of California,Mexico

Extensive CO₂ vents have been discovered in the Wagner Basin, northern Gulf of California, where they create large areas with lowered seawater pH. Such areas are suitable for investigations of long-term biological effects of ocean acidification and effects of CO₂ leakage from subsea carbon capture storage. Here, we show responses of benthic foraminifera to seawater pH gradients at 74–207 m water depth. Living (rose Bengal stained) benthic foraminifera included Nonionella basispinata, Epistominella bradyana and Bulimina marginata. Studies on foraminifera at CO₂ vents in the Mediterranean and off Papua New Guinea have shown dramatic long-term effects of acidified seawater. We found living calcareous benthic foraminifera in low pH conditions in the northern Gulf of California, although there was an impoverished species assemblage and evidence of post-mortem test dissolution.     

Near-inertial currents in the DeSoto Canyon region

Near-inertial currents in the DeSoto Canyon region are described using current and wind observations taken between April 1997 and March 1998 for the “DeSoto Canyon Eddy Intrusion Study”. Distinct energy peaks are present at near-inertial frequencies for the clockwise spectrum and there is little energy at the same frequencies for the counterclockwise current spectrum. In this region, amplitudes of the near-inertial currents can be as high as 40 cm s⁻¹. These currents are surface-intensified and display an increase in amplitude from the shelf break to offshore. Between November 1997 and March 1998, they were effectively generated by shifting winds accompanying passages of cold fronts. For this time period, near-inertial currents are reasonably well-simulated by a mixed-layer model forced by observed winds. During summer 1997, however, enhanced near-inertial motions often resulted from resonance between winds and existing currents.     

Modeling the tidal channel morphodynamics in a macro-tidal embayment,Hangzhou Bay,China

The Hangzhou Bay is a macro-tidal bay located to the south of the Changjiang estuary in China. Along its northern shore, a large-scale tidal channel system has developed, which includes a main northern tidal channel, with a length of more than 50 km and a width up to 10 km, and a secondary southern tidal channel. A process-based morphodynamic model, incorporating the cohesive sediment transport module of Delft3D, is used to analyze the physical processes and mechanisms underlying the formation and evolution of this tidal channel system. The results show that spatial gradients of flood dominance, caused by boundary enhancement via current convergences, is responsible for the formation of the channel system, due to a combination of the various factors such as funnel-shaped geometry hindering associated with the presence of islands, and flow deviation by the southern tidal flat and so on. The model results agree well with the real morphological features. This study also indicates that the reclamation of the southern tidal flat imposes a profound influence on the morphological evolution of the tidal channel system in the Hangzhou Bay. It is feasible to use the model to simulate long-term estuarine morphological changes with cohesive sediment settings.     

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.     

Geology, ore deposits and hydrothermal venting in Bahía Concepción, Baja California Sur, Mexico

Abstract Bahía Concepción is located in the eastern coast of the Baja California peninsula and it is shaped by northwestern–southeastern normal faults. These are associated with a 12–6 Ma rifting episode, although some have been reactivated since the Pliocene. The most abundant rocks correspond to the arc related Comondú Group, Oligocene to Miocene, which forms a mainly calc‐alkaline volcanic and volcaniclastic sequence. There are less extensive outcrops of sedimentary rocks, lava flows, domes and pyroclastic rocks of Pliocene to Quaternary ages. The Neogene volcanism in the area indicates a shift from a subduction regime to an intraplate volcanism related to continental extension and the opening of an oceanic basin. The Bahía Concepción area contains numerous Mn ore deposits, being the biggest at El Gavilán and Guadalupe. The Mn deposits occur as veins, breccias and stockworks, and are composed by Mn oxides (pyrolusite, coronadite, romanechite), dolomite, quartz and barite. The deposits are hosted in volcanic rocks of the Comondú Group and, locally, in Pliocene sedimentary rocks. Thus, the Mn deposits formed between the Middle Miocene and the Pliocene. The mineralized structures are associated with Miocene northwestern–southeastern fault systems, which are analogous to those associated with the Cu‐Co‐Zn‐Mn deposits of El Boleo. The Bahía Concepción area also bears subaerial and submarine hot springs, which are associated with the same fault systems and host rocks. The submarine and subaerial geothermal manifestations south of the bay are possibly related with recent volcanism. The geothermal manifestations within the bay are intertidal hot springs and shallow submarine diffuse venting areas. Around the submarine vents (5–15 m deep, 87°C), Fe‐oxyhydroxide crusts with pyrite and cinnabar precipitate. In the intertidal vents (62°C), aggregates of opal, calcite, barite and Ba‐rich Mn oxides occur covered by silica‐carbonate stromatolitic sinters. Some 10–30 cm thick crustiform veins formed by chalcedony, calcite and barite are also found close to the vents. The hydrothermal fluids exhibit mixed isotopic compositions between δ¹⁸O‐enriched meteoric and local marine water. The precipitation of Ba‐rich Mn oxides around the vent sites could be an active analog for the processes that produced Miocene to Pliocene hydrothermal Mn‐deposits.     

Intertidal mudflat properties, currents and sediment erosion in the partially mixed Tamar Estuary, UK

Results are presented from a 1-year campaign to measure the seasonal variability of some key physical and biological properties of intertidal mudflats over a section of the central Tamar Estuary and to relate these to the physical environment. Seasonal variations in physical mudflat properties, such as grain size, density and moisture content were relatively small. With the exception of the particulate organic carbon content in the upper 0.002 m of surface sediment, biological variations were large. Redox potential exhibited considerable seasonal variation and showed that the sediments were less reduced in winter and more reduced in summer. Chlorophyll a and extracellular polymeric substances (EPS) within the surface 0.002 m of sediment (due to the presence of benthic diatoms) were strongly correlated and exhibited a pronounced seasonal pattern, with smallest values during winter and greatest values during late summer and early autumn. EPS had a dominating influence on sediment erosion, as determined from annular flume measurements. Velocity measurements and velocity modelling indicated that during the flood, and for much of the time during benthic diatom bloom conditions of high chlorophyll a and EPS sediment contents, the stresses exerted by tidal currents were too small to cause significant suspension of sediments over much of the middle and upper mudflats. Suspended fine sediment in the turbidity maximum zone was transported down-estuary and deposited in the main channel at low-water (LW) slack. Some of this sediment, in the form of relatively large aggregates, was subsequently transported onto the mudflats during the flooding tide, where slack currents and fast settling velocities may have enhanced sediment deposition there.Responsible Editor: Hans Burchard     

Observations of sediment resuspension and settling off the mouth of Jiaozhou Bay,Yellow Sea

We deployed bottom-mounted quadrapod equipped with acoustic Doppler current profiler (ADCP), acoustic Doppler velocimeter (ADV), and optical backscatter sensor (OBS) over two semidiurnal tidal cycles along the western coast of the Yellow Sea, China. In combination with shipboard profiling of CTD and LISST-100, we resolved the temporal and spatial distributions of tidal currents, turbulent kinetic energy (TKE), suspended sediment concentration (SSC) and particle size distributions. During the observations, tidal-induced bottom shear stress was the main stirring factor. However, weak tidal flow during the ebb phase was accompanied by two large SSC and median size events. The interactions of seiche-induced oscillations with weak ebb flow induced multiple flow reversals and provided a source of turbulence production, which stripped up the benthic fluff layers (only several millimeters) around the Jiaozhou Bay mouth. Several different methods for inferring mean suspended sediment settling velocity agreed well under peak currents, including estimates using LISST-based Stokes’ settling law, and ADCP-based Rouse profiles, ADV-based inertial-dissipation balance and Reynolds flux. Suspended particles in the study site can be roughly classified into two types according to settling behavior: a smaller, denser class consistent with silt and clay and a larger, less dense class consistent with loosely aggregated flocs. In the present work, we prove that acoustic approaches are robust in simultaneously and non-intrusively estimating hydrodynamics, SSC and settling velocities, which is especially applicable for studying sediment dynamics in tidal environments with moderate concentration levels.     

Tidal currents,energy flux and bottom boundary layer thickness in the Clyde Sea and North Channel of the Irish Sea

A high-resolution three-dimensional model of the Clyde Sea and the adjacent North Channel of the Irish Sea is used to compute the major diurnal and semidiurnal tides in the region, the associated energy fluxes and thickness of the bottom boundary layer. Initially, the accuracy of the model is assessed by performing a detailed comparison of computed tidal elevations and currents in the region, against an extensive database that exists for the M₂, S₂, N₂, K₁ and O₁ tides. Subsequently, the model is used to compute the tidal energy flux vectors in the region. These show that the major energy flux is confined to the North Channel region, with little energy flux into the Clyde Sea. Comparison with the observed energy flux in the North Channel shows that its across-channel distribution and its magnitude are particularly sensitive to the phase difference between elevation and current. Consequently, small changes in the computed values of these parameters due to slight changes of the order of the uncertainty in the open-boundary values to the model, can significantly influence the computed energy flux. The thickness of the bottom boundary layer in the region is computed using a number of formulations. Depending upon the definition adopted, the empirical coefficient C used to determine its thickness varies over the range 0.1 to 0.3, in good agreement with values found in the literature. In the North Channel, the boundary layer thickness occupies the whole water depth, and hence tidal turbulence produced at the sea bed keeps the region well mixed. In the Clyde Sea, the boundary layer thickness is a small fraction of the depth, and hence the region stratifies.Responsible Editor: Phil Dyke     

A Numerical Simulation of the Seasonal Cycle of Temperature, Salinity and Velocity Fields in Tokyo Bay

To understand when oxygen-depleted waters occur, how they develop and when they dissipate in inner Tokyo Bay, realistic simulations were attempted with fine spatial and temporal resolution by applying realistic time dependent external forcing. A 3D hydrodynamic model was driven by time-dependent external forcing factors/parameters such as solar radiation, air temperature, relative humidity, wind velocity, and fluvial discharge, under the open boundary conditions of 1995. A simulated time series of salinity and temperature agreed fairly well with observed data, except in summer. The model failed to reproduce the development of the surface mixed layer in summer. Several sensitivity analyses on the external forcing parameters such as wind velocity and vertical diffusivity were conducted to reproduce the mixed layer. However, changing these parameter values did not improve the model results.     

Tidal and residual circulation in the St. Andrew Bay system,Florida

Two 24-h surveys were conducted in St. Andrew Bay, Florida, during spring and neap tides to describe the tidal and non-tidal circulation patterns and to determine the factors that affect these patterns. In particular, the effect of tidal forcing in modulating such circulation patterns was explored. Observed velocities were fitted to diurnal and semidiurnal harmonics separating tidal motions from sub-tidal motions. Residual flows were compared with an analytic model that allowed variations in the relative contributions from Coriolis acceleration and friction using the Ekman number. A solution with an Ekman number of 0.04 resembled the observations best and indicated that the hydrodynamics were governed by pressure gradient, Coriolis and friction. Locally, advective accelerations became important around headlands in sub-estuaries in the system. The consistency of the mean pattern from October to March suggests that tides play a minor role in modulating the exchange flow. Deviations from the long-term mean are mainly caused by wind-driven coastal setup and setdown.     

Impact of tidal power dams upon tides and environmental conditions in the Sea of Okhotsk

The transformation of natural tidal sea-level and currents is studied resulting from large-scale tidal power plant (TPP) dams in bays of the Sea of Okhotsk (SO). Some effects due to this transformation are estimated based on predictive modelling and a number of expected changes in amplitudes and phases, and spectral composition of tidal oscillations are described. Changes of morphometric properties of basins change the character of tidal motions even on significant distance from a dam. That is why, it is impossible to estimate this impact as usual boundary-value problems. The problem is solved based on “impedance” conditions on the open boundary of the model area, allowing to take into account the radiation of the additional perturbations induced by both waves reflected from the dam and nonlinear effects inside the area. In general, the transformation effects are proportional to the dam size and depend essentially on the dam location, the creation of which can change dissipative and resonance properties of the bays. The changes in tidal energetics of SO due to the dam construction are also considered to show noticeable reconstruction of horizontal energy fluxes and changes in the energy dissipation. Possible environmental consequences are related mainly to the transformation of tidal currents.     

Dynamics of a Prorocentrum minimum bloom along the northern coast of Sinaloa,Mexico

To investigate the relative importance of mesoscale physical events, such as upwellings and physical and chemical variables during an algae bloom of Prorocentrum minimum, 25 sampling sites were established offshore of the Navachiste Lagoon Complex on the east side of the Gulf of California. Samples were analyzed for phytoplankton concentration, water chemistry, and temperature during November 1999, January, March, April, May, and August 2000. Satellite imagery of sea surface temperature (SST) for April 2000 was processed to obtain a synoptic view of the area during the extraordinary bloom of P. minimum in the open waters of the Gulf of California. The bloom was associated with change of oceanographic conditions from moderate winds to calm period, temperature increase and high nitrate (NO₃⁻-N) and ammonia (NH₄⁺-N) content in the offshore waters. Depletion of these nutrients during the bloom suggests that this species uses both types of nitrogen substrates. Cysts in the northernmost sampling stations in January and March indicate that upwelling water, rich in nitrates, also carried a seed stock population of P. minimum. SST patterns in the satellite imagery suggest wind-forcing as the responsible mechanisms triggering the algal bloom offshore of the Navachiste Lagoon Complex.     

Wind effects on the lateral structure of density-driven circulation in Chesapeake Bay

The response of the density-driven circulation in the Chesapeake Bay to wind forcing was studied with numerical experiments. A model of the bay with realistic bathymetry was first applied to produce the density-driven flow under average river discharge and tidal forcing. Subsequently, four spatially uniform wind fields (northeasterly, northwesterly, southwesterly, and southeasterly) were imposed to examine the resulting cross-estuary structure of salinity and flow fields. In general, northeasterly and northwesterly winds intensified the density-driven circulation in the upper and middle reaches of the bay, whereas southeasterly and southwesterly winds weakened it. The response was different in the lower bay, where downwind flow from the upper and middle reaches of the bay competed with onshore/offshore coastal flows. Wind remote effects were dominant, over local effects, on volume transports through the bay entrance. However, local effects were more influential in establishing the sea-level slopes that drove subtidal flows and salinity fields in most of the bay. The effect of vertical stratification on wind-induced flows was also investigated by switching it off. The absence of stratification allowed development of Ekman layers that reached depths of the same order as the water depth. Consequently, bathymetric effects became influential on the homogeneous flow structure causing the wind-induced flow inside the bay to show a marked transverse structure: downwind over the shallow areas and upwind in the channels. In the presence of stratification, Ekman layers became shallower and the wind-induced currents showed weaker transverse structure than those that developed in the absence of stratification. In essence, the wind-driven flows were horizontally sheared under weak stratification and vertically sheared under stratified conditions.     

The impact of hurricanes on sedimenting particulate matter in the semi-arid Bahía de La Paz,Gulf of California

From 2002 through 2004, time-series sediment trap samples were collected from a depth of 410 m in Cuenca Alfonso, Bahía de La Paz, on the SW coast of the Gulf of California. The instrument recorded the impact of the local passage of hurricanes “Ignacio” (24–26 August) and “Marty” (21–23 September) in 2003. These two events accounted for 82% of the total rainfall measured in 2003, equivalent to the annual average precipitation in years without hurricanes. Mean total mass fluxes (TMFs) of 2.88 and 3.58 g m⁻² d⁻¹ were measured during the week of each hurricane as well as the following week. This may have been enough to produce a lamina in the underlying sediment with characteristics peculiar to such events. The terrigenous component was particularly abundant, with notably higher concentrations of Fe, Sc, Co and Cs and REEs. In contrast, TMFs throughout 2002–2004 (excluding the hurricane periods) averaged only 0.73 g m⁻² d⁻¹ and had a larger marine biogenic component. The extraordinary elemental fluxes during the 29 days of hurricane-influenced sedimentation represented a great proportion of the totals over an entire “normal” year: Co (67.8%) >Sc (62.6) >Fe (59.6) >Cs (53.4)>Lu (51.5)>La (51.3)>Yb (51.0)>Ce (49.5) >Tb (48.4) >Sm (44.7)>Cr (36.5) >Ca (31.0)>Eu (25.4%). The terrigenous fraction was calculated using (a) TMF minus the sum of CaCO₃, biogenic silica and organic matter and (b) the ratio of Sc in the trap samples to the average in the Earth's crust. The latter was consistently smaller, but the two methods offered similar results following hurricanes (78% vs. 63%, respectively). For normal sedimentation, however, the difference method yielded values twice as large as the Sc method (58% vs. 30%) This suggests that the mineralogy of the terrigenous fraction may also vary, with unsorted dessert soil being carried to sea by the powerful flash floods associated with hurricanes. Eolian supply of particles, particularly Sc-free quartz grains, possibly from beyond the limited fluvial drainage basin, apparently dominates normal sedimentation.     

Simulating the migration of drifters deployed in the Bay of Biscay, during the Prestige crisis

A main conclusion following the oil spill from the Prestige tanker was that improvements in ocean circulation models were necessary; this was in order to predict, more accurately, the trajectories followed by the oil slicks and hence assist in fight against oil pollution operations. In this contribution, the results of the validation of a semi-empirical ocean circulation model, parameterised for the Bay of Biscay and forced with operational oceano-meteorological remote sensing observations, are shown. The model results have been validated with observations from drifting buoys, deployed in the Bay of Biscay during the crisis. The results show that the model explains a relatively large percentage of the current variability. The comparisons between the real and the estimated drifter trajectories indicate that for 3, 5 and 7 day-long trajectories, the drifter position is estimated with errors of approximately 23, 35 and 46km, respectively. The model reproduces relatively well the trajectory followed by the drifter with the shortest period (23 days).