Investigation and sensitivity analysis of a mechanistic phytoplankton model implemented in a new modular numerical tool (Eco3M) dedicated to biogeochemical modelling

A new class of phytoplankton models with a mechanistic basis has been presented in a companion paper (Baklouti, M., Diaz, F., Pinazo, C., Faure, V., Queguiner, B., 2006. Investigation of mechanistic formulations depicting phytoplankton dynamics for models of marine pelagic ecosystems. Progress in Oceanography). It is the default class of models implemented in our new numerical tool Eco3M, which is dedicated to Ecological, Mechanistic and Modular Modelling. A brief overview of its main features is given in Section 2 of the present paper. In the next sections, a particular phytoplankton model among the aforementioned class has been tested with special emphasis on the mechanistic photosynthesis component relating the photosynthetic rate to the proportion of open photosystems II. The present study encompasses several essential steps that are inherent to any modelling, including model reduction, model sensitivity analysis and comparison of model outputs with experiments. The global sensitivity analysis of the plankton model for one-at-a-time parameter perturbations revealed a restricted set of parameters having major influence on the model outputs. Sensitivity tests involving simultaneous parameter perturbations within the range actually encountered in the literature provided a confidence interval for the outputs. Chemostat experiments performed on nitrate-limited diatoms grown under low (LL) and high-light (HL) conditions have been used for comparison with model outputs. The good fit between measured data and model outputs using the same parameter values in both the LL and HL cases demonstrates the ability of our model to represent the main features of phytoplankton dynamics including photoacclimation. Finally, Eco3M is ultimately intended to include explicit bacterial and zooplankton compartments, as well as to be coupled with ocean circulation models, but the intrinsic behavior of the phytoplankton model has been investigated first, independently of physical forcing.     

GPS测高技术在无验潮水深测量中的应用

应用双频GPS动态后处理高精度定位技术，建立了一套完整的GPS无验潮海洋深度测量作业模式，通过海上试验与传统作业模式作了数值分析比较，结果表明，该作业模式不仅无需验潮，而且能够有效消除传统作业模式中船只动态吃水和涌浪等因素对测量成果的影响，显著提高水深测量成果的精度。     

A Sensitivity Analysis of the Waterline Method of Constructing a Digital Elevation Model for Intertidal Areas in ERS SAR scene of Eastern England

A sensitivity analysis of the waterline method of constructing a Digital Elevation Model (DEM) of an intertidal zone using remote sensing and hydrodynamic modelling is described. Variation in vertical height accuracy as a function of beach slope is investigated using a set of nine ERS Synthetic Aperture Radar (SAR) images of the Humber/Wash area on the English east coast acquired between 1992 and 1994. Waterlines from these images are heighted using a hydrodynamic tide-surge model and interpolated using block kriging. On 1:500 slope beaches, an average block height estimation standard deviation of 18–22 cm is achieved. This rises to 27 cm on 1:100 slope beaches, and 32 cm on 1:30 slope beaches. The average heighting error at different slopes is decomposed into components due to waterline heighting error, inadequate sensor resolution and interpolation inaccuracy. It is shown that, at 1:500 slope, waterline heighting error and interpolation inaccuracy are the main error sources, whilst at 1:30 slope, errors due to inadequate sensor resolution become dominant. The ability of the technique to generate intertidal DEMs for almost the entire coastal zone in a complete ERS SAR scene covering 100×100 km is demonstrated.     

走航测深数据的自动采集

本文介绍应用IBM微机实现走航实时水深及定位数据的自动采集系统。给出了联机系统方框图及采集、传送数据的程序流程图。并对IBM-PC/XT微机与测深仪之间的接口、数码转换和数据处理作了简要的叙述。     

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.     

Beneath the surface: Characteristics of oceanic ecosystems under weak mixing conditions – A theoretical investigation

This study considers an important biome in aquatic environments, the subsurface ecosystem that evolves under low mixing conditions, from a theoretical point of view. Employing a conceptual model that involves phytoplankton, a limiting nutrient and sinking detritus, we use a set of key characteristics (thickness, depth, biomass amplitude/productivity) to qualitatively and quantitatively describe subsurface biomass maximum layers (SBMLs) of phytoplankton. These SBMLs are defined by the existence of two community compensation depths in the water column, which confine the layer of net community production; their depth coincides with the upper nutricline. Analysing the results of a large ensemble of simulations with a one-dimensional numerical model, we explore the parameter dependencies to obtain fundamental steady-state relationships that connect primary production, mortality and grazing, remineralization, vertical diffusion and detrital sinking. As a main result, we find that we can distinguish between factors that determine the vertically integrated primary production and others that affect only depth and shape (thickness and biomass amplitude) of this subsurface production layer. A simple relationship is derived analytically, which can be used to estimate the steady-state primary productivity in the subsurface oligotrophic ocean. The fundamental nature of the results provides further insight into the dynamics of these “hidden” ecosystems and their role in marine nutrient cycling.     

Elucidating the dynamics and mixing agents of a shallow fjord through age tracer modelling

The mixing agents and their role in the dynamics of a shallow fjord are elucidated through an Eulerian implementation of artificial tracers in a three-dimensional hydrodynamic model. The time scales of vertical mixing in this shallow estuary are short, and the artificial tracers are utilized in order to reveal information not detectable in the temperature or salinity fields. The fjord's response to external forcing is investigated through a series of model experiments in which we quantify vertical mixing, transport time scales of fresh water runoff and estuarine circulation in relation to external forcing.Using age tracers released at surface and bottom, we quantify the time scales of downward mixing of surface water and upward mixing of bottom water. Wind is shown to be the major agent for vertical mixing at nearly all depth levels in the fjord, whereas the tide or external sea level forcing is a minor agent and only occasionally more important just close to the bottom. The time scale of vertical mixing of surface water to the bottom or ventilation time scale of bottom water is estimated to be in the range 0.7 h to 9.0 days, with an average age of 2.7 days for the year 2004.The fjord receives fresh water from two streams entering the innermost part of the fjord, and the distribution and age of this water are studied using both ageing and conservative tracers. The salinity variations outside this fjord are large, and in contrast to the salinity, the artificial tracers provide a straight forward analysis of river water content. The ageing tracer is used to estimate transport time scales of river water (i.e. the time elapsed since the water left the river mouth). In May 2004, the typical age of river water leaving the fjord mouth is 5 days. As the major vertical mixing agent is wind, it controls the estuarine circulation and export of river water. When the wind stress is set to zero, the vertical mixing is reduced and the vertical salinity stratification is increased, and the river water can be effectively exported out of the fjord.We also analyse the river tracer fields and salinity field in relation to along estuary winds in order to detect signs of wind-induced straining of the along estuary density gradient. We find that events of down estuary winds are primarily associated with a reduced along estuary salinity gradient due to increased surface salinity in the innermost part of the fjord, and with an overall decrease in vertical stratification and river water content at the surface. Thus, our results show no apparent signs of wind-induced straining in this shallow fjord but instead they indicate increased levels of vertical mixing or upwelling during down estuary wind events.     

Longshore drift estimation using fluorescent tracers: New insights from an experiment at Comporta Beach, Portugal

A fluorescent sand-tracer experiment was performed at Comporta Beach (Portugal) with the aim of acquiring longshore sediment transport data on a reflective beach, the optimization of field and laboratory tracer procedures and the improvement of the conceptual model used to support tracer data interpretation.

The field experiment was performed on a mesotidal reflective beach face in low energetic conditions (significant wave height between 0.4 and 0.5 m). Two different colour tracers (orange and blue) were injected at low tide and sampled in the two subsequent low tides using a high resolution 3D grid extending 450 m alongshore and 30 m cross-shore. Marked sand was detected using an automatic digital image processing system developed in the scope of the present experiment.

Results for the two colour tracers show a remarkable coherence, with high recovery rates attesting data validity. Sand tracer displayed a high advection velocity, but with distinct vertical distribution patterns in the two tides: in the first tide there was a clear decrease in tracer advection velocity with depth while in the second tide, the tracer exhibited an almost uniform vertical velocity distribution. This differing behaviour suggests that, in the first tide, the tracer had not reached equilibrium within the transport system, pointing to a considerable time lag between injection and complete mixing. This issue has important implications for the interpretation of tracer data, indicating that short term tracer experiments tend to overestimate transport rates. In this work, therefore, longshore estimates were based on tracer results obtained during the second tide.

The estimated total longshore transport rate at Comporta Beach was 2 × 10^− 3 m³/s, more than four times larger than predicted using standard empirical longshore formulas. This discrepancy, which results from the unusually large active moving layer observed during the experiment, confirms the idea that most common longshore transport equations under-estimate total sediment transport in plunging/surging waves.     

3D gravity modelling for Anyongbok Seamount in the East Sea

A seamount chain with an approximately WNW trend is observed in the northeastern Ulleung Basin. It has been argued that these seamounts, including two islands called Ulleung and Dok islands, were formed by a hotspot process or by ridge related volcanism. Many geological and geophysical studies have been done for all the seamounts and islands in the chain except Anyongbok Seamount, which is close to the proposed spreading ridge. We first report morphological characteristics, sediment distribution patterns, and the crustal thickness of Anyongbok Seamount using multibeam bathymetry data, seismic reflection profiles, and 3D gravity modeling. The morphology of Anyongbok Seamount shows a cone shaped feature and is characterized by the development of many flank cones and flank rift zones. The estimated surface volume is about 60 km³, and implies that the seamount is smaller than the other seamounts in the chain. No sediments have been observed on the seamount except the lower slope, which is covered by more than 1,000 m of strata. The crustal structure obtained from a 3D gravity modeling (GFR = 3.11, SD 3.82 = mGal) suggests that the seamount was formed around the boundary of the Ulleung Plateau and the Ulleung Basin, and the estimated crustal thickness is about 20 km, which is a little thicker than other nearby seamounts distributed along the northeastern boundary of the Ulleung Basin. This significant crustal thickness also implies that Anyongbok Seamount might not be related to ridge volcanism.