Data assimilation for a coastal area morphodynamic model: Morecambe Bay

This paper investigates the use of data assimilation in coastal area morphodynamic modelling using Morecambe Bay as a study site. A simple model of the bay has been enhanced with a data assimilation scheme to better predict large-scale changes in bathymetry observed in the bay over a 3-year period. The 2DH decoupled morphodynamic model developed for the work is described, as is the optimal interpolation scheme used to assimilate waterline observations into the model run. Each waterline was acquired from a SAR satellite image and is essentially a contour of the bathymetry at some level within the inter-tidal zone of the bay. For model parameters calibrated against validation observations, model performance is good, even without data assimilation. However the use of data assimilation successfully compensates for a particular failing of the model, and helps to keep the model bathymetry on track. It also improves the ability of the model to predict future bathymetry. Although the benefits of data assimilation are demonstrated using waterline observations, any observations of morphology could potentially be used. These results suggest that data assimilation should be considered for use in future coastal area morphodynamic models.     

Three-Dimensional Velocity Field and Cross-Frontal Water Exchange in the Kuroshio Extension

A spread of warm water from the first crest of the Kuroshio Extension is periodically enhanced by northward warm water intrusions from the main current. The water type in the spread area was previously found to be the same as that in the Kuroshio front at depth. In looking for the possible mechanism responsible for the northward warm water intrusions, a dynamic analysis in the Kuroshio front was carried out by using CTD, ADCP, AVHRR and ARGOS buoy data, obtained in 1996 by the R.V. Hakuho Maru. Downstream, cross-stream and vertical velocities in the Kuroshio Extension were found by using a "stream coordinate system". The velocity field in the Kuroshio front at the first crest showed a double structure with two surface velocity maxima. In the inner part of the front, relatively high cross-stream (northward) and vertical (upward) velocities were found. Thus, this study suggests that while water particles flow downstream along the first stationary meander of the Kuroshio Extension, they also experience lateral and vertical movements which allow the deeper water from an upstream location to rise to the surface layer, and in certain locations to deflect northward. By assuming isopycnal movement and conservation of potential vorticity, it was found that in those locations where anticyclonic curvature of the meander increases, warm water is more likely to deflect northward. High ageostrophic components observed in the first 300 m of the water column are probably related to the relatively high cross-stream and vertical velocities in the inner part of the front.     

A preliminary study of carbon system in the East China Sea

In the central part of the East China Sea, the activity of CO₂ in the surface water and total carbonate, pH and alkalinity in the water column were determined in winter and autumn of 1993. The activity of CO₂ in the continental shelf water was about 50 ppm lower than that of surface air. This decrease corresponds to the absorption of about 40 gC/m²/yr of atmospheric CO₂ in the coastal zone or 1 GtC/yr in the global continental shelf, if this rate is applicable to entire coastal seas. The normalized total carbonate contents were higher in the water near the coast and near the bottom. This increase toward the bottom may be due to the organic matter deposited on the bottom. This conclusion is supported by the distribution of pH. The normalized alkalinity distribution also showed higher values in the near-coast water, but in the surface water, indicating the supply of bicarbonate from river water. The residence time of the East China Sea water, including the Yellow Sea water, has been calculated to be about 0.8 yr from the excess alkalinity and the alkalinity input. Using this residence time and the excess carbonate, we can estimate that the amount of dissolved carbonate transported from the coastal zone to the oceanic basin is about 70 gC/m²/yr or 2 GtC/yr/area-of-global-continental-shelf. This also means that the rivers transport carbon to the oceans at a rate of 30 gC/m²/yr of the coastal sea or 0.8 GtC/yr/ area-of-global shelf, the carbon consisting of dissolved inorganic carbonate and terrestrial organic carbon decomposed on the continental shelf.     

砷酸根的液—固界面交换吸附研究——pH和Ca(Ⅱ)的影响

研究ｐＨ和Ｃａ（Ⅱ）对砷酸根在海水与ＮａＣｌ介质中对水合氧化物和粘土矿物上交换吸附的影响：（１）在酸性条件下和海水及ＮａＣｌ介质中，砷酸根在针铁矿上的交换率Ｅ（％）—ｐＨ曲线呈平顶峰型、在δ－ＭｎＯ呈峰型、在粘土上分别呈双峰型和肩峰型。这是砷酸根存在多种形式与固体表面发生离子／配位子交换以及介质中阴离子竞争交换吸附共同作用的结果。（２）在碱性条件下的海水介质中，砷酸根与水合氧化物和粘土矿物的Ｅ（％）—ｐＨ曲线都呈Ｖ型，这是因为生成了具有一定结构的液—固界面三元络合物之故。（３）在ＮａＣｌ介质中，Ｃａ加入起促进作用，再次证明液—固界面三元络合物生成机理与Ｖ型曲线的内在关系     

A Reconstruction of Observed Profiles in the Sea East of Japan Using Vertical Coupled Temperature-Salinity EOF Modes

It is important to estimate hard-to-observe parameters in the ocean interior from easy-to-observe parameters. This study therefore demostrates a reconstruction of observed temperature and salinity profiles of the sea east of Japan (30°≈40°N, 140°≈150°E). The reconstruction was done by estimating suboptimal state from several values of the observed profiles and/or sea surface dynamic height (SDH) calculated from the profiles. The estimation used a variational method with vertical coupled temperature-salinity empirical orthogonal function (EOF) modes. Profiles of temperature and salinity in the subtropical region are effectively reconstructed from in situ temperature profile data, or sea surface temperature (SST) and SDH. For example, the analyzed temperature field from SST and SDH has an accuracy to within 1°C in the subtropical region. Salinity in the sea north of Kuroshio, however, is difficult to estimate because of its complex variability which is less correlated with temperature than in the subtropical region. Sea surface salinity is useful to estimate the subsurface structure. We also show the possibility that the estimation is improved by considering nonlinearity in the equation calculating SDH from temperature and salinity analysis values in order to examine the misfit between analysis and observation. Analysis using TOPEX/POSEIDON altimetry data instead of SDH was also performed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Dynamics and Variability of Terra Nova Bay Polynya

Abstract. We present a process study on the dynamics and variability of the Terra Nova Bay polynya in the western sector of the Ross Sea. The air-sea heat exchange is known to be particularly large in polynya during the winter, when differences between air and sea temperatures are large. We apply a 1-D model (Pease, 1987; Van Woert, 1999a, 1999b), which is modified in the latent heat parameterisation in order to account for time-dependent relative humidity and cloud coverage. Furthermore, the Ice Collection Depth is correlated linearly with a variable wind speed. The model is forced with two different meteorological data sets: the operational analysis of the European Center for Medium Range Weather Forecasts atmospheric data set and the meteorological parameters measured by an Automatic Weather Station located on the coast of Terra Nova Bay. The results are compared in terms of polynya extension, ice, and High Salinity Shelf Water production. According to the two different wind velocities, the results obtained from the different data sets clearly differ. Qualitatively, however, the results are in good agreement.     

Seakeeping assessment of fishing vessels in conceptual design stage

The main idea of this paper is to identify functional relations between seakeeping characteristics and hull form parameters of Mediterranean fishing vessels. Multiple regression analysis is used for quantitative assessment through a computer software that is based on the SQL Server Database. The seakeeping attributes under investigation are the transfer functions of heave and pitch motions and of absolute vertical acceleration at stern, while the ship parameters influencing motion dynamics have been classified into two groups: displacement (Δ) and main dimensions (L, B, T), coefficients that define the details of the hull form (C_WP, C_VP, LCB, LCF, etc.).Four multiple regression models having different parameter combinations are here investigated and discussed, giving way to the so-called ‘Simple Model’, ‘Intermediate Model’, ‘Enhanced 1 Model’ and ‘Enhanced 2 Model’. The obtained results are more than satisfactory for seakeeping predictions during the conceptual design stage.     

Exchange flow of oil and sea-water in a ruptured submarine pipeline

The rupture of a submarine oil pipeline starts various mechanisms leading to an oil spill. Among these mechanisms the leakage of oil driven by the difference in specific gravities of oil and sea-water is difficult to estimate. A simple mathematical model has been developed and laboratory experiments have been carried out to obtain an insight into the density-driven exchange flow and to determine the leak rate. The mathematical model is predictive and takes account of the effects of friction, inclination of the pipeline, and inertia of the fluid. The experiments were done in a horizontal model pipeline. Theoretical and experimental results are in satisfactory agreement.     

Preconditioned Optimizing Utility for Large-dimensional analyses (POpULar)

I present the derivation of the Preconditioned Optimizing Utility for Large-dimensional analyses (POpULar), which is developed for adopting a non-diagonal background error covariance matrix in nonlinear variational analyses (i.e., analyses employing a non-quadratic cost function). POpULar is based on the idea of a linear preconditioned conjugate gradient method widely adopted in ocean data assimilation systems. POpULar uses the background error covariance matrix as a preconditioner without any decomposition of the matrix. This preconditioning accelerates the convergence. Moreover, the inverse of the matrix is not required. POpULar therefore allows us easily to handle the correlations among deviations of control variables (i.e., the variables which will be analyzed) from their background in nonlinear problems. In order to demonstrate the usefulness of POpULar, we illustrate two effects which are often neglected in studies of ocean data assimilation before. One is the effect of correlations among the deviations of control variables in an adjoint analysis. The other is the nonlinear effect of sea surface dynamic height calculation required when sea surface height observation is employed in a three-dimensional ocean analysis. As the results, these effects are not so small to neglect.     

Consequences of winter upwelling events on biogeochemical and phytoplankton patterns in a western Galician ria (NW Iberian peninsula)

The consequences of two upwelling events in mid- (MW) and late (LW) winter on biogeochemical and phytoplankton patterns were studied in the Pontevedra Ria and compared with the patterns measured under typical winter conditions and under a summer upwelling event. Thermohaline patterns measured during the mid-winter upwelling event (MW-up) revealed the intrusion of saltier seawater (35.9) into the ria associated with the Iberian Poleward Current (IPC). During the late-winter upwelling event (LW-up), the seawater which had welled up into the ria showed characteristics of the Eastern North Atlantic Central Water mass (ENACW). In both cases the measured water residence time (4 days during MW-up and 10 days during LW-up) was related to both meteorological and fluvial forcing. This residence time contrasts with that of summer upwelling (7 days) and with that estimated under unfavorable upwelling atmospheric conditions (2–4 weeks). During MW-up, the ria became poor in nutrients due to continental freshwater dilution, associated with the shorter residence time of the water, and the intrusion of IPC, which is a water body poor in nutrient salts: 2.9 μM of nitrate, 0.1 μM of phosphate and 1.5 μM of silicate. During this event, the ria exported 3.4 mol_DIN s⁻¹, compared with 6.9 mol_DIN s⁻¹ in non-upwelling conditions. Phytoplankton showed a uniform distribution throughout the ria, as during unfavorable upwelling conditions, and was characterized by the dominance of diatoms, mainly Nitzschia longissima and Skeletonema costatum. During LW-up, a nutrient depletion in the photic layer also occurred, but as a result of a phytoplankton spring bloom developing at this time. The ria was a nutrient trap where 4.1 mol_DIN s⁻¹ were processed by photosynthesis. This budget is three times higher than the one under non-upwelling conditions. In contrast with the MW-up, which had no effect on primary production, during LW-up the ria became more productive, although not as productive as during a summer upwelling event (9.9 mol_DIN s⁻¹). The taxonomic composition of the phytoplankton community did not change noticeably during LW-up and the summer upwelling, with the same species present and changing only in relative proportions. Diatoms were always the dominant microphytoplankton community, with Pseudonitzschia pungens, Thalassionema nitzschioides and several species of Chaetoceros as characteristic taxons.