Validating and assessing the sensitivity of the climate model with an ocean general circulation model developed at the Institute of Atmospheric Physics,Russian Academy of Sciences

A new version of the Institute of Atmospheric Physics, Russian Academy of Sciences (IAP RAS), climate model (CM) has been developed using an ocean general circulation model instead of the statistical-dynamical ocean model applied in the previous version. The spatial resolution of the new ocean model is 3° in latitude and 5° in longitude, with 25 unevenly spaced vertical levels. In the previous version of the oceanic model, as in the atmospheric model, the horizontal resolution was 4.5° in latitude and 6° in longitude, with four vertical levels (the upper quasi-homogeneous layer, seasonal thermocline, abyssal ocean, and bottom friction layer). There is no correction for the heat and momentum fluxes between the atmosphere and ocean in the new version of the IAP RAS CM. Numerical experiments with the IAP RAS CM have been performed under current initial and boundary conditions, as well as with an increasing concentration of atmospheric carbon dioxide. The main simulated atmospheric and oceanic fields agree quite well with observational data. The new version’s equilibrium temperature sensitivity to atmospheric CO₂ doubling was found to be 2.9 K. This value lies in the mid-range of estimates (2–4.5 K) obtained from simulations with state-of-the-art models of different complexities.     

Optimal interpolation of the data of remote sensing of the sea surface

We analyze the statistical structure of the data of remote sensing of the fields of surface temperature and the concentration of chlorophyll in the northwest part of the Black Sea with an aim to apply the method of optimal interpolation to the restoration of the maps of these fields according to the restricted amounts of data. We determine the scales of smoothing for the selection of the random components of the fields and estimate their difference from the statistical model uniform and isotropic over the space correlation coefficients. The functions of spectral densities are constructed along various sections and the best intervals of discretization of the fields required to perform the optimal interpolation with given accuracy are determined. It is shown that the construction of the maps of fields under the conditions of controlled accuracy of reconstruction of their values remains possible even if the amounts of the input satellite data are significantly decreased. We consider a procedure of evaluation of the best possible amounts of data required for the use of the method of optimal interpolation. Some examples of application of this procedure to the data of remote sensing of the fields of surface temperature and the concentration of chlorophyll are presented.     

Impact of dynamic feedbacks between sedimentation,sea-level rise,and biomass production on near-surface marsh stratigraphy and carbon accumulation

Salt marshes accrete both organic and inorganic sediments. Here we present analytical and numerical models of salt marsh sedimentation that, in addition to capturing inorganic processes, explicitly account for above- and belowground organic processes including root growth and decay of organic carbon. The analytical model is used to examine the bias introduced by organic processes into proxy records of sedimentation, namely ¹³⁷Cs and ²¹⁰Pb. We find that accretion rates estimated using ²¹⁰Pb will be less than accretion rates estimated using the ¹³⁷Cs peak in steadily accreting marshes if (1) carbon decay is significant and (2) data for ²¹⁰Pb extend below the ¹³⁷Cs peak. The numerical model expands upon the analytical model by including belowground processes such as compaction and root growth, and by explicitly tracking the evolution of aboveground biomass and its effect on sedimentation rates. Using the numerical model we explore how marsh stratigraphy responds to sediment supply and the rate of sea-level rise. It is calibrated and tested using an extensive data set of both marsh stratigraphy and measurements of vegetation dynamics in a Spartina alterniflora marsh in South Carolina, USA. We find that carbon accumulation in marshes is nonlinearly related to both the supply of inorganic sediment and the rate of sea-level rise; carbon accumulation increases with sea-level rise until sea-level rise reaches a critical rate that drowns the marsh vegetation and halts carbon accumulation. The model predicts that changes in carbon storage resulting from changing sediment supply or sea-level rise are strongly dependent on the background sediment supply: if inorganic sediment supply is reduced in an already sediment poor marsh the storage of organic carbon will increase to a far greater extent than in a sediment-rich marsh, provided that the rate of sea-level rise does not exceed a threshold. These results imply that altering sediment supply to estuaries (e.g., by damming upstream rivers or altering littoral sediment transport) could lead to significant changes in the carbon budgets of coastal salt marshes.     

The fate of terrestrial organic matter in two Scottish sea lochs

Sea lochs are zones of rapid organic matter (OM) turnover. Most of this OM is of allochthonous origin, being introduced into the lochs via freshwater input. In this study the behaviour of terrestrially derived OM was elucidated using a combination of parameters which indicate OM diagenesis in the near surface sediments from two Scottish sea lochs, Loch Creran and Loch Etive. Alkaline CuO oxidation was used to determine lignin phenols which serve as biomarkers for terrestrial OM in sediments. Stable carbon isotope, total carbon and nitrogen and total OM (including the labile and refractory fractions) compositions were also determined.     

Isotopic,trace element and nutrient characterization of coastal waters from Ubatuba inner shelf area,south-eastern Brazil

Stable isotopes, tritium, radium isotopes, radon, trace elements and nutrients data were collected during two sampling campaigns in the Ubatuba coastal area (south-eastern Brazil) with the aim of investigating submarine groundwater discharge (SGD) in the region. The isotopic composition (δD, δ¹⁸O, ³H) of submarine waters was characterised by significant variability and heavy isotope enrichment. The stable isotopes and tritium data showed good separation of groundwater and seawater groups. The contribution of groundwater in submarine waters varied from a few % to 17%. Spatial distribution of ²²²Rn activity concentration in surface seawater revealed changes between 50 and 200 Bq m⁻³ which were in opposite relationship with observed salinities. Time series measurements of ²²²Rn activity concentration in Flamengo Bay (from 1 to 5 kBq m⁻³), obtained by in situ underwater gamma-spectrometry showed a negative correlation between the ²²²Rn activity concentration and tide/salinity. This may be caused by sea level changes as tide effects induce variations of hydraulic gradients, which increase ²²²Rn concentration during lower sea level, and opposite, during high tides where the ²²²Rn activity concentration is smaller. The estimated SGD fluxes varied during 22–26 November between 8 and 40 cm d⁻¹, with an average value of 21 cm d⁻¹ (the unit is cm³/cm² per day). The radium isotopes and nutrient data showed scattered distributions with offshore distance and salinity, which implies that in a complex coast with many small bays and islands, the area has been influenced by local currents and groundwater–seawater mixing. SGD in the Ubatuba area is fed by coastal contaminated groundwater and re-circulated seawater (with small admixtures of groundwater), which claims for potential environmental concern with implications on the management of freshwater resources in the region.     

Water level variations in the Neuse and Pamlico Estuaries,North Carolina due to local and remote forcing

Water level time series records from the Neuse and Pamlico River Estuaries were statistically compared to local and distant wind field data, water level records within the Pamlico Sound and also coastal ocean sites to determine the relative contribution of each time series to water levels in the Neuse and Pamlico Estuaries. The objectives of this study were to examine these time series data using various statistical methods (i.e. autoregressive, empirical orthogonal function analysis (EOF), exploratory data analysis (EDA)) to determine short- and long-time-scale variability, and to develop predictive statistical models that can be used to estimate past water level fluctuations in both the Neuse Estuary (NE) and Pamlico Estuary (PE). Short- and long-time-scale similarities were observed in all time series of estuarine, Pamlico Sound and subtidal coastal ocean water level and wind component data, due to events (nor'easters, fronts and tropical systems) and seasonality. Empirical orthogonal function analyses revealed a strong coastal ocean and wind field contribution to water level in the NE and PE. Approximately 95% of the variation was captured in the first two EOF components for water level data from the NE, sound and coastal ocean, and 70% for the PE, sound and coastal ocean. Spectral density plots revealed strong diurnal signals in both wind and water level data, and a strong cross correlation and coherency between the NE water level and the North/South wind component. There was good agreement between data and predictions using autoregressive statistical models for the NE (R² = 0.92) and PE (R² = 0.76). These methods also revealed significant autoregressive lags for the NE (days 1 and 3) and for the PE (days 1, 2 and 3). Significant departures from predictions are attributed to local meteorological and hydrological events. The autoregressive techniques showed significant predictive improvement over ordinary least squares methods. The results are considered within the context of providing long time-scale hindcast data for the two estuaries, and the importance of these data for multidisciplinary researchers and managers.     

The onset of fish colonization in a coastal defence structure (Chioggia,Northern Adriatic Sea)

Coastal defence structures constitute the most extensive hard substrates of the Northwestern Adriatic Sea and are known to sustain rich benthic and nektonic communities. To appreciate the pattern of colonization, we studied the fish assemblage of a recently deployed breakwater. We compared observations from two years, the different sides (landward and seaward) of the barrier, and the two fringes, characterized by timing of work completion. The results indicate that colonization, still in process, follows different patterns among species. Benthic and necto-benthic species presented a striking increase in abundance and richness in the second year of colonization, while more mobile species did not evince any variation between years. Differences in mobility among species suggest that the latter group may have reached the breakwater from nearby artificial substrates, whereas the former colonized the new structure as recruits. In addition, fish assemblages differed between the two sides, likely due to variation in the environmental characteristics, and according to depth, reflecting species preferences.     

Lagrangian measurements and modelling of fluid advection in the inner surf and swash zones

New laboratory and field data are presented on fluid advection into the swash zone. The data illustrate the region of the inner surf zone from which sediment can be directly advected into the swash zone during a single uprush, which is termed the advection length. Experiments were conducted by particle tracking in a Lagrangian reference frame, and were performed for monochromatic breaking waves, solitary bores, non-breaking solitary waves and field conditions. The advection length is normalised by the run-up length to give an advection ratio, A, and different advection ratios are identified on the basis of the experimental data. The data show that fluid enters the swash zone from a region of the inner surf zone that can extend a distance seaward of the bore collapse location that is approximately equal to half of the run-up length. This region is about eight times wider than the region predicted by the classical swash solution of Shen and Meyer [Shen, M.C., Meyer, R.E., 1963. Climb of a bore on a beach. Part 3. Runup. Journal of Fluid Mechanics 16, 113–125], as illustrated by Pritchard and Hogg [Pritchard, D., Hogg, A.J., 2005. On the transport of suspended sediment by a swash event on a plane beach. Coastal Engineering 52, 1–23]. Measured advection ratios for periodic waves show no significant trend with Iribarren number, consistent with self-similarity in typical swash flows. The data are compared to recent characteristic solutions of the non-linear shallow water wave (NLSW) equations and both finite difference and finite volume solutions of the NLSW equations.     

Large-scale dune erosion tests to study the influence of wave periods

Large-scale physical model tests were performed to quantify the effects of the wave period on dune erosion. Attention was focussed on 2D cross-shore effects in a situation with sandy dunes and extreme water levels and wave conditions. Besides profile measurements, detailed measurements in time and space of water pressure, flow velocities and sediment concentrations were performed in the near near-shore area. It was concluded that a longer wave period leads to a larger dune erosion volume and to a larger landward retreat of the dune face. Tests with double-peaked wave spectra showed that the influence of the spectral shape on dune erosion was best represented by the T_m − 1,0 spectral mean wave period, better than the peak wave period, T_p. The effect of the wave period on dune erosion was implemented in a dune erosion prediction method that estimates erosion volumes during normative storm conditions for the Dutch coast. More details of the measurements and additional analyses of physical processes are described in an accompanying paper by Van Thiel de Vries et al. [Van Thiel de Vries, J.S.M., van Gent, M.R.A., Reniers, A.J.H.M. and Walstra, D.J.R., submitted for publication. Analysis of dune erosion processes in large scale flume experiments, In this volume of Coastal Engineering.].     

Morphodynamic response of a two-barred beach to a shoreface nourishment

The response to a shoreface nourishment of the two-bar system at Noordwijk (the Netherlands) is analyzed based on a daily data set of time-exposure video images collected during about 6 years, complemented with topographic and bathymetric surveys. The 1.7 Mm³ nourishment, implemented as a 3 km alongshore bump seaward of the outer bar, migrated more than 300 m onshore in 4 years before losing its integrity. Furthermore, the nourishment interrupted the autonomous seaward migration of both bars for the entire duration of the study period and, allaying earlier fears, did not intensify the three-dimensional patterns in the bars, such as the crescentic plan-shape and rip channels. The nourishment did result in clear head effects on both flanks, with the bar becoming discontinuous and the flank section decaying or becoming attached to an offshore-located bar, while the section of bar landward of the nourishment became attached to a landward-located bar. This sequence of morphologies is known as bar switching. Each switching episode took almost one year to complete and can therefore not be ascribed to individual wave events. We suspect that shoreface nourishments enhance the possibility of bar switching by creating alongshore variability in the position and depth of the outer bar and in its cross-shore migration rate and direction. The Noordwijk nourishment did not influence the shoreline position as its trend did not undergo distinctive variations after 1998. Differences in the response of the Noordwijk sandbars to the shoreface nourishment compared with other Dutch nourishments are attributed to the location and size (volume per unit length) of the nourishment with respect to the sandbars and to the median grain size of the nourished material.