Temporal and spatial variation in the Nazaré Canyon (Western Iberian margin): Inter-annual and canyon heterogeneity effects on meiofauna biomass and diversity

The Nazaré Canyon on the Portuguese Margin (NE Atlantic) was sampled during spring-summer for three consecutive years (2005–2007), permitting the first inter-annual study of the meiofaunal communities at the Iberian Margin at two abyssal depths (~3500 m and ~4400 m). Using new and already published data, the meiofauna standing stocks (abundance and biomass) and nematode structural and functional diversity were investigated in relation to the sediment biogeochemistry (e.g. organic carbon, nitrogen, chlorophyll a, phaeopigments) and grain size. A conspicuous increase in sand content from 2005 to 2006 and decrease of phytodetritus at both sites, suggested the occurrence of one or more physical disturbance events. Nematode standing stocks and trophic diversity decreased after these events, seemingly followed by a recovery/recolonisation period in 2007, which was strongly correlated with an increase in the quantity and bioavailability of phytodetrital organic matter supplied. Changes in meiofauna assemblages, however, also differed between stations, likely because of the contrasting hydrodynamic and food supply conditions. Higher meiofauna and nematode abundances, biomass and trophic complexity were found at the shallowest canyon station, where the quantity, quality and bioavailability of food material were higher than at the deeper site. The present results suggest that even though inter-annual variations in the sedimentary environment can regulate the meiofauna in the abyssal Nazaré Canyon, heterogeneity between sampling locations in the canyon were more pronounced.

Biomass changes and trophic amplification of plankton in a warmer ocean

Ocean warming can modify the ecophysiology and distribution of marine organisms, and relationships between species, with nonlinear interactions between ecosystem components potentially resulting in trophic amplification. Trophic amplification (or attenuation) describe the propagation of a hydroclimatic signal up the food web, causing magnification (or depression) of biomass values along one or more trophic pathways. We have employed 3-D coupled physical-biogeochemical models to explore ecosystem responses to climate change with a focus on trophic amplification. The response of phytoplankton and zooplankton to global climate-change projections, carried out with the IPSL Earth System Model by the end of the century, is analysed at global and regional basis, including European seas (NE Atlantic, Barents Sea, Baltic Sea, Black Sea, Bay of Biscay, Adriatic Sea, Aegean Sea) and the Eastern Boundary Upwelling System (Benguela). Results indicate that globally and in Atlantic Margin and North Sea, increased ocean stratification causes primary production and zooplankton biomass to decrease in response to a warming climate, whilst in the Barents, Baltic and Black Seas, primary production and zooplankton biomass increase. Projected warming characterized by an increase in sea surface temperature of 2.29 ± 0.05 °C leads to a reduction in zooplankton and phytoplankton biomasses of 11% and 6%, respectively. This suggests negative amplification of climate driven modifications of trophic level biomass through bottom-up control, leading to a reduced capacity of oceans to regulate climate through the biological carbon pump. Simulations suggest negative amplification is the dominant response across 47% of the ocean surface and prevails in the tropical oceans; whilst positive trophic amplification prevails in the Arctic and Antarctic oceans. Trophic attenuation is projected in temperate seas. Uncertainties in ocean plankton projections, associated to the use of single global and regional models, imply the need for caution when extending these considerations into higher trophic levels.

Role of zooplankton dynamics for Southern Ocean phytoplankton biomass and global biogeochemical cycles

Global ocean biogeochemistry models currently employed in climate change projections use highly simplified representations of pelagic food webs. These food webs do not necessarily include critical pathways by which ecosystems interact with ocean biogeochemistry and climate. Here we present a global biogeochemical model which incorporates ecosystem dynamics based on the representation of ten plankton functional types (PFTs); six types of phytoplankton, three types of zooplankton, and heterotrophic bacteria. We improved the representation of zooplankton dynamics in our model through (a) the explicit inclusion of large, slow-growing zooplankton, and (b) the introduction of trophic cascades among the three zooplankton types. We use the model to quantitatively assess the relative roles of iron vs. grazing in determining phytoplankton biomass in the Southern Ocean High Nutrient Low Chlorophyll (HNLC) region during summer. When model simulations do not represent crustacean macrozooplankton grazing, they systematically overestimate Southern Ocean chlorophyll biomass during the summer, even when there was no iron deposition from dust. When model simulations included the developments of the zooplankton component, the simulation of phytoplankton biomass improved and the high chlorophyll summer bias in the Southern Ocean HNLC region largely disappeared. Our model results suggest that the observed low phytoplankton biomass in the Southern Ocean during summer is primarily explained by the dynamics of the Southern Ocean zooplankton community rather than iron limitation. This result has implications for the representation of global biogeochemical cycles in models as zooplankton faecal pellets sink rapidly and partly control the carbon export to the intermediate and deep ocean.

Difference equation approach to two-thermocouple sensor characterisation in constant velocity flow environments

Thermocouples are one of the most popular devices for temperature measurement due to their robustness, ease of manufacture and installation, and low cost. However, when used in certain harsh environments, for example, in combustion systems and engine exhausts, large wire diameters are required, and consequently the measurement bandwidth is reduced. This article discusses a software compensation technique to address the loss of high frequency fluctuations based on measurements from two thermocouples. In particular, a difference equation sDEd approach is proposed and compared with existing methods both in simulation and on experimental test rig data with constant flow velocity. It is found that the DE algorithm, combined with the use of generalized total least squares for parameter identification, provides better performance in terms of time constant estimation without any a priori assumption on the time constant ratios of the thermocouples.

Exploiting a priori time constant ratio information in difference equation two-thermocouple sensor characterisation

The characterization of thermocouple sensors for temperature measurement in varying-flow environments is a challenging problem. Recently, the authors introduced novel difference-equation-based algorithms that allow in situ characterization of temperature measurement probes consisting of two-thermocouple sensors with differing time constants. In particular, a linear least squares (LS) lambda formulation of the characterization problem, which yields unbiased estimates when identified using generalized total LS, was introduced. These algorithms assume that time constants do not change during operation and are, therefore, appropriate for temperature measurement in homogenous constant-velocity liquid or gas flows. This paper develops an alternative ß-formulation of the characterization problem that has the major advantage of allowing exploitation of a priori knowledge of the ratio of the sensor time constants, thereby facilitating the implementation of computationally efficient algorithms that are less sensitive to measurement noise. A number of variants of the ß-formulation are developed, and appropriate unbiased estimators are identified. Monte Carlo simulation results are used to support the analysis.

Solving the advection-dispersion equation using the discrete puff particle method

A flexible, mass-conservative numerical technique for solving the advection-dispersion equation for miscible contaminant transport is presented. The method combines features of puff transport models from air pollution studies with features from the random walk particle method used in water resources studies, providing a deterministic time-marching algorithm which is independent of the grid Peclet number and scales from one to higher dimensions simply. The concentration field is discretised into a number of particles, each of which is treated as a point release which advects and disperses over the time interval. The dispersed puff is itself discretised into a spatial distribution of particles whose masses can be pre-calculated. Concentration within the simulation domain is then calculated from the mass distribution as an average over some small volume. Comparison with analytical solutions for a one-dimensional fixed-duration concentration pulse and for two-dimensional transport in an axisymmetric flow field indicate that the algorithm performs well. For a given level of accuracy the new method has lower computation times than the random walk particle method.

A Van der Pol-Mathieu equation for the dynamics of dust grain charge in dusty plasmas

The chaotic profile of dust grain dynamics associated with dust-acoustic oscillations in a dusty plasma is considered. The collective behaviour of the dust plasma component is described via a multi-fluid model, comprising Boltzmann distributed electrons and ions, as well as an equation of continuity possessing a source term for the dust grains, the dust momentum and Poisson's equations. A Van der Pol–Mathieu-type nonlinear ordinary differential equation for the dust grain density dynamics is derived. The dynamical system is cast into an autonomous form by employing an averaging method. Critical stability boundaries for a particular trivial solution of the governing equation with varying parameters are specified. The equation is analysed to determine the resonance region, and finally numerically solved by using a fourth-order Runge–Kutta method. The presence of chaotic limit cycles is pointed out.