Shallow-water longshore drift-fed submarine fan deposition (Moisie River Delta, Eastern Canada)

Submarine canyons and associated submarine fans are in some cases located at the end of a littoral cell where they act as conduits for the transfer of eroded terrigenous sediments to the marine environment. Such fans are generally found in deep-water settings at >500 m water depth. Offshore the Moisie River Delta (NW Gulf of St. Lawrence, Eastern Canada), high-resolution multibeam bathymetry and seismic data led to the discovery of an unusually shallow submarine fan (≤60 m) located at the end of a littoral cell. Sediment is transported westward on the shallow coastal shelf, as demonstrated by the downcurrent displacement of oblique nearshore sandbars where the shelf narrows to less than 1 km. The steep slope near the end of the littoral cell is incised by a channel that feeds a submarine fan composed of smaller channels and depositional lobes. According to existing Holocene evolution models for the region, the fan formed within the last 5,000 years. Its evolution is largely due to the transport of sediment by longshore drift. Multibeam echosounder and seismic data also reveal that the gravity-driven accretion of the submarine fan is characterized mainly by two processes, i.e., frequent small-scale, downslope migration of sandwaves on the slope, and more episodic slumping/turbidity-current activity in the deeper part of the fan. This study documents that, besides their common deep-water location, smaller-scale submarine fans can occur also in very shallow water, implying that they could be more frequent than previously thought both in modern environments and in the rock record.     

Spatiotemporal aspects of silica buffering in restored tidal marshes

Losses of pelagic diatom production resulting from silica limitation have not only been blamed for toxic algal blooms, but for the reduction in ability of coastal food webs to support higher trophic levels. Recent research has shown the importance of advective seepage water fluxes of dissolved silica (DSi) from freshwater marshes to pelagic waters during moments of riverine Si-limitation. In this study, we investigated the potential impact of recently installed new tidal areas along the Schelde estuary, located in former polder areas and characterized by so-called controlled reduced tidal regimes (CRT). Nine mass-balance studies were conducted in a newly constructed CRT in the freshwater Schelde estuary. During complete tidal cycles both DSi and amorphous silica (ASi) concentrations were monitored at the entrance culverts and in different habitats in the marsh. A swift DSi-delivery capacity was observed despite the shifted spatiotemporal frame of exchange processes compared to reference marshes. As silica-accumulating vegetation is not yet present, and difference with reference marshes’ deliveries is surprisingly small, we indicate diatomaceous debris and phytoliths to be the main silica source. Although further research is necessary on the driving forces of the different processes involved, restoration of former agricultural areas under CRT-regime provide the potential to buffer silica in the estuary.     

Spurious diapycnal mixing in terrain-following coordinate models: The problem and a solution

In this paper, we identify a crucial numerical problem in sigma coordinate models, leading to unacceptable spurious diapycnal mixing. This error is a by-product of recent advances in numerical methods, namely the implementation of high-order diffusive advection schemes. In the case of ROMS, spurious mixing is produced by its third-order upwind advection scheme, but our analysis suggests that all diffusive advection schemes would behave similarly in all sigma models. We show that the common idea that spurious mixing decreases with resolution is generally false. In a coarse-resolution regime, spurious mixing increases as resolution is refined, and may reach its peak value when eddy-driven lateral mixing becomes explicitly resolved. At finer resolution, diffusivities are expected to decrease but with values that only become acceptable at resolutions finer than the kilometer. The solution to this problem requires a specifically designed advection scheme. We propose and validate the RSUP3 scheme, where diffusion is split from advection and is represented by a rotated biharmonic diffusion scheme with flow-dependent hyperdiffusivity satisfying the Peclet constraint. The rotated diffusion operator is designed for numerical stability, which includes improvements of linear stability limits and a clipping method adapted to the sigma-coordinate. Realistic model experiments in a southwest Pacific configuration show that RSUP3 is able to preserve low dispersion and diffusion capabilities of the original third-order upwind scheme, while preserving water mass characteristics. There are residual errors from the rotated diffusion operator, but they remain acceptable. The use of a constant diffusivity rather than the Peclet hyperdiffusivity tends to increase these residual errors which become unacceptable with Laplacian diffusion. Finally, we have left some options open concerning the use of time filters as an alternative to spatial diffusion. A temporal discretization approach to the present problem (including implicit discretization) will be reported in a following paper.     

Characterisation of Soot Emitted by Domestic Heating,Aircraft and Cars Using Diesel or Biodiesel

The characterisation of aggregates, like soot, firstly requires the determination of the size distribution of the primary particles. The primary particle size of combustion generated aggregates depends upon the combustion environment and the formation conditions, such as temperature, pressure and fuel-to-air ratio, among others. Since the combustion characteristics are different in the different types of burners, the characterisation of primary particles may offer the possibility to distinguish soot from different sources. In this paper, we present the signature of the primary particles and the aggregates of soot emitted by cars using diesel or biodiesel, by domestic heating, and by aircraft exhausts, which can be considered as the major sources as derived from measurements on transmission electron micrographs. The size distributions of all aggregates types with different aerodynamic diameter were log-normal and quasi-monodisperse. The size distribution of the primary particles for soot emitted by different sources showed minor differences. However, a comparison between the diameter of the primary particles and those obtained using a standard method for carbon black revealed discrepancies. The median diameter of the primary particles was combined with the median number of primary particles in an aggregate to calculate the relative particle surface area available for adsorption. In a similar way, the relative specific surface area was determined. The surface area was measured using the Brunauer-Emmett-Teller (B.E.T.) nitrogen adsorption method and the relative surface area available for adsorption was calculated.     

The chemical changes of DOM from black waters to coastal marine waters by HPLC combined with ultrahigh resolution mass spectrometry

How dissolved organic matter (DOM) undergoes chemical changes during its transit from river to ocean remains a challenge due to its complex structure. In this study, DOM along a river transect from black waters to marine waters is characterized using an offline combination of reversed-phase high performance liquid chromatography (RP-HPLC) coupled to electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry (ESI-FTICR-MS), as well as tandem ESI-FTICR-MS. In addition, a water extract from degraded wood that mainly consists of lignins is used for comparison to the DOM from this transect. The HPLC chromatograms of all DOM samples and the wood extract show two major well-separated components; one is hydrophilic and the other is hydrophobic, based on their elution order from the C₁₈ column. From the FTICR-MS analysis of the HPLC fractions, the hydrophilic components mainly contain low molecular weight compounds (less than 400 Da), while the hydrophobic fractions contain the vast majority of compounds of the bulk C₁₈ extracted DOM. The wood extract and the DOM samples from the transect of black waters to coastal marine waters show strikingly similar HPLC chromatograms, and the FTICR-MS analysis further indicates that a large fraction of molecular formulas from these samples are the same, existing as lignin-like compounds. Tandem mass spectrometry experiments show that several representative molecules from the lignin-like compounds have similar functional group losses and fragmentation patterns, consistent with modified lignin structural entities in the wood extract and these DOM samples. Taken together, these data suggest that lignin-derived compounds may survive the transit from the river to the coastal ocean and can accumulate there because of their refractory nature.     

Contaminant concentrations and toxicity of sea-surface microlayer near Los Angeles, California

Sea-surface microlayer samples were collected from six nearshore areas receiving different amounts of anthropogenic inputs. The samples were analyzed for selected trace metals, chlorinated hydrocarbons, and polycyclic aromatic hydrocarbons. The relative toxicities of the samples were determined with fish embryo bioassays.Contaminant concentrations generally increased from offshore to the inshore stations. Contaminant concentrations were several orders of magnitude higher in microlayer samples from the highly industrialized Los Angeles and Long Beach harbors compared to samples from a site 15km offshore. Microlayer samples from the inshore stations were significantly more toxic, and induced significantly more developmental abnormalities and chromosome aberrations, than samples from the offshore stations.     

Wind-driven interannual variability over the northeast Pacific Ocean

Interannual variability of the sea surface height (SSH) over the northeast Pacific Ocean is hindcast with a reduced-gravity, quasi-geostrophic model that includes linear damping. The model is forced with monthly Ekman pumping fields derived from the NCEP reanalysis wind stresses. The numerical solution is compared with SSH observations derived from satellite altimeter data and gridded at a lateral resolution of 1 degree. Provided that the reduced gravity parameter is chosen appropriately, the results demonstrate that the model has significant hindcast skill over interior regions of the basin, away from continental boundaries. A damping time scale of 2 to 3 years is close to optimal, although the hindcast skill is not strongly dependent on this parameter.A simplification of the quasi-geostrophic model is considered in which Rossby waves are eliminated, yielding a Markov model driven by local Ekman pumping. The results approximately reproduce the hindcast skill of the more complete quasi-geostrophic model and indicate that the interannual SSH variability is dominated by the local response to wind forcing. There is a close correspondence the two leading empirical orthogonal modes of the local model and those of the observed SSH anomalies. The latter account for over half of the variance of the interannual signal over the region.     

The Regional Ocean Modeling System (ROMS) 4-dimensional variational data assimilation systems: Part II - Performance and application to the California Current System

The Regional Ocean Modeling System (ROMS) 4-dimensional variational (4D-Var) data assimilation systems have been systematically applied to the mesoscale circulation environment of the California Current to demonstrate the performance and practical utility of the various components of ROMS 4D-Var. In particular, we present a comparison of three approaches to 4D-Var, namely: the primal formulation of the incremental strong constraint approach; the dual formulation “physical-space statistical analysis system”; and the dual formulation indirect representer approach. In agreement with theoretical considerations all three approaches converge to the same ocean circulation estimate when using the same observations and prior information. However, the rate of convergence of the dual formulation was found to be inferior to that of the primal formulation. Other aspects of the 4D-Var performance that relate to the use of multiple outer-loops, preconditioning, and the weak constraint are also explored. A systematic evaluation of the impact of the various components of the 4D-Var control vector (i.e. the initial conditions, surface forcing and open boundary conditions) is also presented. It is shown that correcting for uncertainties in the model initial conditions exerts the largest influence on the ability of the model to fit the available observations. Various important diagnostics of 4D-Var are also examined, including estimates of the posterior error, the information content of the observation array, and innovation-based consistency checks on the prior error assumptions. Using these diagnostic tools, we find that more than 90% of the observations assimilated into the model provide redundant information. This is a symptom of the large percentage of satellite data that are used and to some extent the nature of the data processing employed. This is the second in a series of three papers describing the ROMS 4D-Var systems.