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.

Sea-level control on the connection between shelf-edge deltas and the Bourcart canyon head (western Mediterranean) during the last glacial/interglacial cycle

A dense grid of high- and very high resolution seismic data, together with piston cores and borehole data providing time constraints, enables us to reconstruct the history of the Bourcart canyon head in the western Mediterranean Sea during the last glacial/interglacial cycle. The canyon fill is composed of confined channel–levee systems fed by a series of successively active shelf fluvial systems, originating from the west and north. Most of the preserved infill corresponds to the interval between Marine Isotope Stage (MIS) 3 and the early deglacial (19 cal ka BP). Its deposition was strongly controlled by a relative sea level that impacted the direct fluvial/canyon connection. During a period of around 100 kyr between MIS 6 and MIS 2, the canyon “prograded” by about 3 km. More precisely, several parasequences can be identified within the canyon fill. They correspond to forced-regressed parasequences (linked to punctuated sea-level falls) topped by a progradational-aggradational parasequence (linked to a hypothetical 19-ka meltwater pulse (MWP)). The bounding surfaces between forced-regressed parasequences are condensed intervals formed during intervals of relative sediment starvation due to flooding episodes. The meandering pattern of the axial incision visible within the canyon head, which can be traced landward up to the Agly paleo-river, is interpreted as the result of hyperpycnal flows initiated in the river mouth in a context of increased rainfall and mountain glacier flushing during the early deglacial.