Detection of mesoscale thermal fronts from 4km data using smoothing techniques: Gradient-based fronts classification and basin scale application

In order to optimize frontal detection in sea surface temperature fields at 4 km resolution, a combined statistical and expert-based approach is applied to test different spatial smoothing of the data prior to the detection process. Fronts are usually detected at 1 km resolution using the histogram-based, single image edge detection (SIED) algorithm developed by Cayula and Cornillon in 1992, with a standard preliminary smoothing using a median filter and a 3 × 3 pixel kernel. Here, detections are performed in three study regions (off Morocco, the Mozambique Channel, and north-western Australia) and across the Indian Ocean basin using the combination of multiple windows (CMW) method developed by Nieto, Demarcq and McClatchie in 2012 which improves on the original Cayula and Cornillon algorithm. Detections at 4 km and 1 km of resolution are compared. Fronts are divided in two intensity classes (“weak” and “strong”) according to their thermal gradient. A preliminary smoothing is applied prior to the detection using different convolutions: three type of filters (median, average and Gaussian) combined with four kernel sizes (3 × 3, 5 × 5, 7 × 7, and 9 × 9 pixels) and three detection window sizes (16 × 16, 24 × 24 and 32 × 32 pixels) to test the effect of these smoothing combinations on reducing the background noise of the data and therefore on improving the frontal detection. The performance of the combinations on 4 km data are evaluated using two criteria: detection efficiency and front length. We find that the optimal combination of preliminary smoothing parameters in enhancing detection efficiency and preserving front length includes a median filter, a 16 × 16 pixel window size, and a 5 × 5 pixel kernel for strong fronts and a 7 × 7 pixel kernel for weak fronts. Results show an improvement in detection performance (from largest to smallest window size) of 71% for strong fronts and 120% for weak fronts. Despite the small window used (16 × 16 pixels), the length of the fronts has been preserved relative to that found with 1 km data. This optimal preliminary smoothing and the CMW detection algorithm on 4 km sea surface temperature data are then used to describe the spatial distribution of the monthly frequencies of occurrence for both strong and weak fronts across the Indian Ocean basin. In general strong fronts are observed in coastal areas whereas weak fronts, with some seasonal exceptions, are mainly located in the open ocean. This study shows that adequate noise reduction done by a preliminary smoothing of the data considerably improves the frontal detection efficiency as well as the global quality of the results. Consequently, the use of 4 km data enables frontal detections similar to 1 km data (using a standard median 3 × 3 convolution) in terms of detectability, length and location. This method, using 4 km data is easily applicable to large regions or at the global scale with far less constraints of data manipulation and processing time relative to 1 km data.

Stable isotope ratios in bentho-demersal biota along a depth gradient in the bay of biscay: A multitrophic study

Although stable isotope ratios are increasingly used to investigate the trophic ecology of marine organisms, their spatial variations are still poorly understood in the coastal environment. In this study, we measured the stable isotope composition (δ13C, δ15N) of suspended particulate organic matter (SPOM) (primary producer), a suspension feeder, the great scallop Pecten maximus (primary consumer), megabenthic decapods and benthic fishes (secondary consumers) along a depth gradient (from 5m to 155m depth) across the continental shelf of the Bay of Biscay. Although the three trophic levels exhibited similar δ13C patterns along the gradient, the δ15N patterns varied between SPOM, scallops and carnivores. The δ15N difference between SPOM and scallops decreased with increasing depth, suggesting that non trophic factors may affect the stable isotope composition of scallops at deepest sampling stations. An opposed trend was found between scallops and carnivores, suggesting that the trophic level of these carnivores increased at higher depth, possibly as an adaptation to lower prey abundances. Although our results suggest that primary consumers are suitable to establish isotopic baselines in coastal environments, we stress the need for further studies aiming at characterizing the variability of stable isotopes in coastal biota, and the respective effects of baseline, trophic and metabolic factors in their isotopic composition.

The Southeast Indian Ridge between 88°E and 118°E: Gravity anomalies and crustal accretion at intermediate spreading rates

Although slow spreading ridges characterized by a deep axial valley and fast spreading ridges characterized by an axial bathymetric high have been extensively studied, the transition between these two modes of axial morphology is not well understood. We conducted a geophysical-survey of the intermediate spreading rate Southeast Indian Ridge between 88 degrees E and 118 degrees E, a 2300-km-long section of the ridge located between the Amsterdam hot spot and the Australian-Antarctic Discordance where satellite gravity data suggest that the Southeast Indian Ridge (SEIR) undergoes a change from an axial high in the west to an axial valley in the east. A basic change in axial morphology is found near 103 degrees 30'E in the shipboard data; the axis to the west is marked by an axial high, while a valley is found to the east. Although a well-developed axial high, characteristic of the East Pacific Rise (EPR), is occasionally present, the more common observation is a rifted high that is lower and pervasively faulted, sometimes with significant (> 50 m throw) faults within a kilometer of the axis. A shallow axial valley (< 700 m deep) is observed from 104 degrees E to 114 degrees E with a sudden change to a deep (>1200 m deep) valley across a transform at 114 degrees E. The changes in axial morphology along the SEIR are accompanied by a 500 m increase in near-axis ridge flank depth from 2800 m near 88 degrees E to 3300 m near 114 degrees E and by a 50 mGal increase in the regional level of mantle Bouguer gravity anomalies over the same distance, The regional changes in depth and mantle Bouguer anomaly (MBA) gravity can be both explained by a 1.7-2.4 km change in crustal thickness or by a mantle temperature change of 50 degrees C-90 degrees C. In reality, melt supply (crustal thickness) and mantle temperature are linked, so that changes in both may occur simultaneously and these estimates serve as upper bounds. The along-axis MBA gradient is not uniform. Pronounced steps in the regional level of the MBA gravity occur at 103 degrees 30'E-104 degrees E and at 114 degrees E-116 degrees E and correspond to the changes in the nature of the axial morphology and in the amplitude of abyssal hill morphology suggesting that the different forms of morphology do not grade into each other but rather represent distinctly different forms of axial (s)tructure and tectonics with a sharp transition between them. The change from an axial high to an axial valley requires a threshold effect in which the strength of the lithosphere changes quickly. The presence or absence of a quasi-steady state magma chamber may provide such a mechanism. The different forms of axial morphology are also associated with different intrasegment MBA gravity patterns. Segments with an axial high have an MBA low located at a depth minimum near the center of the segment, At EPR-like segments, the MBA low is about 10 mGal with along-axis gradients of 0.15-0.25 mGal/km, similar to those observed at the EPR, Rifted highs have a shallower low and lower gradients suggesting an attenuated composite magma chamber and a reduced and perhaps episodic melt supply. Segments with a shallow axial valley have very flat along-axis MBA profiles with little correspondence between axial depth and axial MBA gravity.