The biogeochemical structuring role of horizontal stirring: Lagrangian perspectives on iron delivery downstream of the Kerguelen Plateau

Field campaigns are instrumental in providing ground truth for understanding and modeling global ocean biogeochemical budgets. A survey however can only inspect a fraction of the global oceans, typically a region hundreds of kilometers wide for a temporal window of the order of (at most) several weeks. This spatiotemporal domain is also the one in which the mesoscale activity induces through horizontal stirring a strong variability in the biogeochemical tracers, with ephemeral, local contrasts which can easily mask the regional and seasonal gradients. Therefore, whenever local in situ measures are used to infer larger-scale budgets, one faces the challenge of identifying the mesoscale structuring effect, if not simply to filter it out. In the case of the KEOPS2 investigation of biogeochemical responses to natural iron fertilization, this problem was tackled by designing an adaptive sampling strategy based on regionally optimized multisatellite products analyzed in real time by specifically designed Lagrangian diagnostics. This strategy identified the different mesoscale and stirring structures present in the region and tracked the dynamical frontiers among them. It also enabled back trajectories for the ship-sampled stations to be estimated, providing important insights into the timing and pathways of iron supply, which were explored further using a model based on first-order iron removal. This context was essential for the interpretation of the field results. The mesoscale circulation-based strategy was also validated post-cruise by comparing the Lagrangian maps derived from satellites with the patterns of more than one hundred drifters, including some adaptively released during KEOPS2 and a subsequent research voyage. The KEOPS2 strategy was adapted to the specific biogeochemical characteristics of the region, but its principles are general and will be useful for future in situ biogeochemical surveys.

Surface Salinity in the North Atlantic Subtropical Gyre During the STRASSE/SPURS Summer 2012 Cruise

We investigated a 100 × 100 km high-salinity region of the North Atlantic subtropical gyre during the Sub-Tropical Atlantic Surface Salinity Experiment/Salinity Processes in the Upper-ocean Regional Study (STRASSE/SPURS) cruise from August 21, 2012, to September 9, 2012. Results showed great variability in sea surface salinity (SSS; over 0.3 psu) in the mesoscale, over 7 cm of total evaporation, and little diapycnal mixing below 36 m depth, the deepest mixed layers encountered. Strong currents in the southwestern part of the domain, and the penetration of freshwater, suggest that advection contributed greatly to salinity evolution. However, it was further observed that a smaller cyclonic structure tucked between the high SSS band and the strongest currents contributed to the transport of high SSS water along a narrow front. Cross-frontal transport by mixing is also a possible cause of summertime reduction of SSS. The observed structure was also responsible for significant southward salt transport over more than 200 km.

Lateral diffusivity coefficients from the dynamics of a SF6 patch in a coastal environment

The dispersion of a patch of the tracer sulfur hexafluoride (SF6) is used to assess the lateral diffusivity in the coastal waters of the western part of the Gulf of Lion (GoL), northwestern Mediterranean Sea, during the Latex10 experiment (September 2010). Immediately after the release, the spreading of the patch is associated with a strong decrease of the SF6 concentrations due to the gas exchange from the ocean to the atmosphere. This has been accurately quantified, evidencing the impact of the strong wind conditions during the first days of this campaign. Few days after the release, as the atmospheric loss of SF6 decreased, lateral diffusivity coefficient at spatial scales of 10 km has been computed using two approaches. First, the evolution of the patch with time was combined with a diffusion-strain model to obtain estimates of the strain rate (γ = 2.5 10- 6 s- 1) and of the lateral diffusivity coefficient (Kh = 23.2 m2s− 1). Second, a steady state model was applied, showing Kh values similar to the previous method after a period of adjustment between 2 and 4.5 days. This implies that after such period, our computation of Kh becomes insensitive to the inclusion of further straining of the patch. Analysis of sea surface temperature satellite imagery shows the presence of a strong front in the study area. The front clearly affected the dynamics within the region and thus the temporal evolution of the patch. Our results are consistent with previous studies in open ocean and demonstrate the success and feasibility of those methods also under small-scale, rapidly-evolving dynamics typical of coastal environments.

20 Years of the Atlantic Meridional Transect — AMT

The AMT (www.amt-uk.org) is a multidisciplinary programme which undertakes biological, chemical, and physical oceanographic research during an annual voyage between the UK and a destination in the South Atlantic such as the Falkland Islands,South Africa, or Chile. This transect of >12,000 km crosses a range of ecosystems from subpolar to tropical, from euphotic shelf seas and upwelling systems, to oligotrophic mid-ocean gyres. The year 2015 has seen two milestones in the history of the AMT: the achievement of 20 years of this unique ocean going programme and the departure of the 25th cruise on the 15th of September. Both of these events were celebrated in June this year with an open science conference hosted by the Plymouth Marine Laboratory (PML) and will be further documented in a special issue of Progress in Oceanography which is planned for publication in 2016. Since 1995, the 25 research cruises have involved 242 sea-going scientists from 66 institutes representing 22 countries.

Initial Examination of AltiKa's Individual Echoes

The AltiKa altimeter records the reflection of Ka-band radar pulses from the Earth’s surface, with the commonly used waveform product involving the summation of 96 returns to provide average echoes at 40 Hz. Occasionally there are one-second recordings of the complex individual echoes (IEs), which facilitate the evaluation of on-board processing and offer the potential for new processing strategies. Our investigation of these IEs over the ocean confirms the on-board operations, whilst noting that data quantization limits the accuracy in the thermal noise region. By constructing average waveforms from 32 IEs at a time, and applying an innovative subwaveform retracker, we demonstrate that accurate height and wave height information can be retrieved from very short sections of data. Early exploration of the complex echoes reveals structure in the phase information similar to that noted for Envisat’s IEs.

Metocean Comparisons of Jason-2 and AltiKa—A Method to Develop a New Wind Speed Algorithm

As well as range, the AltiKa altimeter provides estimates of wave height, Hs and normalized backscatter, s0, that need to be assessed prior to statistics based on them being included in climate databases. An analysis of crossovers with the Jason-2 altimeter shows AltiKa Hs values to be biased high by only »0.05m, with a standard deviation (s.d.) of »0.1m for seven-point averages. AltiKa’s s 0 values are 2.5–3 dB less than those from Jason-2, with a s.d. of »0.3 dB, with these relatively large mismatches to be expected as AltiKa measures a different part of the spectrum of sea surface roughness. A new wind speed algorithm is developed through matchinghistogram of s0 values to that for Jason-2 wind speeds. The algorithm is robust to the use of short durations of data, with a consistency at roughly the 0.1 m/s level. Incorporation of Hs as a secondary input reduces the assessed error at crossovers from 0.82 m/s to 0.71 m/s. A comparison across all altimeter frequencies used to date demonstrates that the lowest wind speeds preferentially develop the shortest scales of roughness.

Applications of Satellite Earth Observations section - NEODAAS: Providing satellite data for efficient research

The NERC Earth Observation Data Acquisition and Analysis Service (NEODAAS) provides a central point of Earth Observation (EO) satellite data access and expertise for UK researchers. The service is tailored to individual users’ requirements to ensure that researchers can focus effort on their science, rather than struggling with correct use of unfamiliar satellite data.

Identifying predictable foraging habitats for a wide-ranging marine predator using ensemble ecological niche models

Aim: Ecological niche modelling can provide valuable insight into species' environmental preferences and aid the identification of key habitats for populations of conservation concern. Here, we integrate biologging, satellite remote-sensing and ensemble ecological niche models (EENMs) to identify predictable foraging habitats for a globally important population of the grey-headed albatross (GHA) Thalassarche chrysostoma. Location: Bird Island, South Georgia; Southern Atlantic Ocean. Methods: GPS and geolocation-immersion loggers were used to track at-sea movements and activity patterns of GHA over two breeding seasons (n = 55; brood-guard). Immersion frequency (landings per 10-min interval) was used to define foraging events. EENM combining Generalized Additive Models (GAM), MaxEnt, Random Forest (RF) and Boosted Regression Trees (BRT) identified the biophysical conditions characterizing the locations of foraging events, using time-matched oceanographic predictors (Sea Surface Temperature, SST; chlorophyll a, chl-a; thermal front frequency, TFreq; depth). Model performance was assessed through iterative cross-validation and extrapolative performance through cross-validation among years. Results: Predictable foraging habitats identified by EENM spanned neritic (<500 m), shelf break and oceanic waters, coinciding with a set of persistent biophysical conditions characterized by particular thermal ranges (3–8 °C, 12–13 °C), elevated primary productivity (chl-a > 0.5 mg m−3) and frequent manifestation of mesoscale thermal fronts. Our results confirm previous indications that GHA exploit enhanced foraging opportunities associated with frontal systems and objectively identify the APFZ as a region of high foraging habitat suitability. Moreover, at the spatial and temporal scales investigated here, the performance of multi-model ensembles was superior to that of single-algorithm models, and cross-validation among years indicated reasonable extrapolative performance. Main conclusions: EENM techniques are useful for integrating the predictions of several single-algorithm models, reducing potential bias and increasing confidence in predictions. Our analysis highlights the value of EENM for use with movement data in identifying at-sea habitats of wide-ranging marine predators, with clear implications for conservation and management.

Basking sharks and oceanographic fronts: quantifying associations in the north-east Atlantic

Understanding the mechanisms linking oceanographic processes and marine vertebrate habitat use is critical to effective management of populations of conservation concern. The basking shark Cetorhinus maximus has been shown to associate with oceanographic fronts – physical interfaces at the transitions between water masses – to exploit foraging opportunities resulting from aggregation of zooplankton. However, the scale, significance and variability of these observed associations have not yet been established. Here, we quantify the influence of mesoscale (10s – 100s km) frontal activity on habitat use over timescales of weeks to months. We use animal-mounted archival tracking with composite front mapping via Earth Observation (EO) remote sensing to provide an oceanographic context to individual shark movements. We investigate levels of association with fronts occurring over two spatio-temporal scales, (i) broad-scale seasonally persistent frontal zones and (ii) contemporaneous mesoscale thermal and chl-a fronts. Using random walk simulations and logistic regression within an iterative generalised linear mixed modelling (GLMM) framework, we find that seasonal front frequency is a significant predictor of shark presence. Temporally-matched oceanographic metrics also indicate that sharks demonstrate a preference for productive regions, and associate with contemporaneous thermal and chl-a fronts more frequently than could be expected at random. Moreover, we highlight the importance of cross-frontal temperature change and persistence, which appear to interact to affect the degree of prey aggregation along thermal fronts. These insights have clear implications for understanding the preferred habitats of basking sharks in the context of anthropogenic threat management and marine spatial planning in the northeast Atlantic.

Satellite monitoring of harmful algal blooms (HABs) to protect the aquaculture industry

Harmful algal blooms (HABs) can cause sudden and considerable losses to fish farms, for example 500,000 salmon during one bloom in Shetland, and also present a threat to human health. Early warning allows the industry to take protective measures. PML's satellite monitoring of HABs is now funded by the Scottish aquaculture industry. The service involves processing EO ocean colour data from NASA and ESA in near-real time, and applying novel techniques for discriminating certain harmful blooms from harmless algae. Within the AQUA-USERS project we are extending this capability to further HAB species within several European countries.

NEODAAS: Providing satellite data for efficient research

The NERC Earth Observation Data Acquisition and Analysis Service (NEODAAS) provides a central point of Earth Observation (EO) satellite data access and expertise for UK researchers. The service is tailored to individual users’ requirements to ensure that researchers can focus effort on their science, rather than struggling with correct use of unfamiliar satellite data.

The effect of the North Atlantic Subpolar Front as a boundary in pelagic biogeography decreases with increasing depth and organism size

Broad-scale patterns in the distribution of deep-sea pelagic species and communities are poorly known. An important question is whether biogeographic boundaries identified from surface features are important in the deep mesopelagic and bathypelagic. We present community analyses of discrete-depth samples of mesozooplankton and micronekton to full-ocean depth collected in the area where the Mid-Atlantic Ridge is crossed by the Subpolar Front. The results show that the distributional discontinuity associated with the front, which is strong near the surface, decreases with increasing depth. Both the frontal separation near the surface and the community convergence at increasing depths were clearer for mesozooplankton than for micronekton.

Identifying predictable foraging habitats for a wide-ranging marine predator using ensemble ecological niche models

Aim: Ecological niche modelling can provide valuable insight into species' environmental preferences and aid the identification of key habitats for populations of conservation concern. Here, we integrate biologging, satellite remote-sensing and ensemble ecological niche models (EENMs) to identify predictable foraging habitats for a globally important population of the grey-headed albatross (GHA) Thalassarche chrysostoma. Location: Bird Island, South Georgia; Southern Atlantic Ocean. Methods: GPS and geolocation-immersion loggers were used to track at-sea movements and activity patterns of GHA over two breeding seasons (n = 55; brood-guard). Immersion frequency (landings per 10-min interval) was used to define foraging events. EENM combining Generalized Additive Models (GAM), MaxEnt, Random Forest (RF) and Boosted Regression Trees (BRT) identified the biophysical conditions characterizing the locations of foraging events, using time-matched oceanographic predictors (Sea Surface Temperature, SST; chlorophyll a, chl-a; thermal front frequency, TFreq; depth). Model performance was assessed through iterative cross-validation and extrapolative performance through cross-validation among years. Results: Predictable foraging habitats identified by EENM spanned neritic (<500 m), shelf break and oceanic waters, coinciding with a set of persistent biophysical conditions characterized by particular thermal ranges (3–8 °C, 12–13 °C), elevated primary productivity (chl-a > 0.5 mg m−3) and frequent manifestation of mesoscale thermal fronts. Our results confirm previous indications that GHA exploit enhanced foraging opportunities associated with frontal systems and objectively identify the APFZ as a region of high foraging habitat suitability. Moreover, at the spatial and temporal scales investigated here, the performance of multi-model ensembles was superior to that of single-algorithm models, and cross-validation among years indicated reasonable extrapolative performance. Main conclusions: EENM techniques are useful for integrating the predictions of several single-algorithm models, reducing potential bias and increasing confidence in predictions. Our analysis highlights the value of EENM for use with movement data in identifying at-sea habitats of wide-ranging marine predators, with clear implications for conservation and management.

A method for semi-automated objective quantification of linear bedforms from multi-scale Digital Elevation Models

The increasing availability of large, detailed digital representations of the Earth’s surface demands the application of objective and quantitative analyses. Given recent advances in the understanding of the mechanisms of formation of linear bedform features from a range of environments, objective measurement of their wavelength, orientation, crest and trough positions, height and asymmetry is highly desirable. These parameters are also of use when determining observation-based parameters for use in many applications such as numerical modelling, surface classification and sediment transport pathway analysis. Here, we (i) adapt and extend extant techniques to provide a suite of semi-automatic tools which calculate crest orientation, wavelength, height, asymmetry direction and asymmetry ratios of bedforms, and then (ii) undertake sensitivity tests on synthetic data, increasingly complex seabeds and a very large-scale (39 000km2) aeolian dune system. The automated results are compared with traditional, manually derived,measurements at each stage. This new approach successfully analyses different types of topographic data (from aeolian and marine environments) from a range of sources, with tens of millions of data points being processed in a semi-automated and objective manner within minutes rather than hours or days. The results from these analyses show there is significant variability in all measurable parameters in what might otherwise be considered uniform bedform fields. For example, the dunes of the Rub’ al Khali on the Arabian peninsula are shown to exhibit deviations in dimensions from global trends. Morphological and dune asymmetry analysis of the Rub’ al Khali suggests parts of the sand sea may be adjusting to a changed wind regime from that during their formation 100 to 10 ka BP.

Decadal reanalysis of biogeochemical indicators and fluxes in the North West European shelf-sea ecosystem

In this paper we present the first decadal reanalysis simulation of the biogeochemistry of the North West European shelf, along with a full evaluation of its skill and value. An error-characterized satellite product for chlorophyll was assimilated into a physical-biogeochemical model of the North East Atlantic, applying a localized Ensemble Kalman filter. The results showed that the reanalysis improved the model predictions of assimilated chlorophyll in 60% of the study region. Model validation metrics showed that the reanalysis had skill in matching a large dataset of in situ observations for ten ecosystem variables. Spearman rank correlations were significant and higher than 0.7 for physical-chemical variables (temperature, salinity, oxygen), ∼0.6 for chlorophyll and nutrients (phosphate, nitrate, silicate), and significant, though lower in value, for partial pressure of dissolved carbon dioxide (∼0.4). The reanalysis captured the magnitude of pH and ammonia observations, but not their variability. The value of the reanalysis for assessing environmental status and variability has been exemplified in two case studies. The first shows that between 340,000-380,000 km2 of shelf bottom waters were oxygen deficient potentially threatening bottom fishes and benthos. The second application confirmed that the shelf is a net sink of atmospheric carbon dioxide, but the total amount of uptake varies between 36-46 Tg C yr−1 at a 90% confidence level. These results indicate that the reanalysis output dataset can inform the management of the North West European shelf ecosystem, in relation to eutrophication, fishery, and variability of the carbon cycle.