Coastal Upwelling Front Detection off Central Chile ( 36 . 5 – 37 ◦ S ) and Spatio-Temporal Variability of Frontal Characteristics

In Eastern Boundary Upwelling Systems, cold coastal waters are separated from offshore by a strong cross-shore Sea Surface Temperature (SST) gradient zone. This upwelling front plays a major role for the coastal ecosystem. This paper proposes a method to automatically identify the front and define its main characteristics (position, width, and intensity) from high resolution data. The spatio-temporal variability of the front characteristics is then analyzed in a region off Central Chile (37◦S), from 2003 to 2016. The front is defined on daily 1 km-resolution SST maps by isotherm T0 with T0 computed from mean SST with respect to the distance from the coast. The probability of detecting a front, as well as the front width and intensity are driven by coastal wind conditions and increased over the 2007–2016 period compared to the 2003–2006 period. The front position, highly variable, is related to the coastal jet configuration and does not depend on the atmospheric forcing. This study shows an increase by 14% in the probability of detecting a front and also an intensification by 17% of the cross-front SST difference over the last 14 years. No trend was found in the front position.


Introduction
In the main Eastern Boundary Current Systems (EBUS: California, Humboldt, Canarias, and Benguela), upwelling forced by alongshore equatorward winds usually creates a surface cold tongue in the coastal band.This cold tongue is limited offshore by a strong cross-shore density gradient zone, referred to as the upwelling front (e.g., [1][2][3][4]).Such front separates the nutrient-rich shelf waters nearshore from the lower nutrient water masses in the oceanic area or Coastal Transition Zone (CTZ) (e.g., [5][6][7]).In addition, this front acts as a retention and accumulation zone, and usually represents a barrier for the offshore distribution of coastal planktonic components [8][9][10][11].Altogether, upwelling dynamics and retention of coastal waters are important processes in the understanding of the dispersal of coastal waters, nutrients and plankton and, thus, the dynamics of EBUS.
The upwelling cold tongue frequently presents strong seasonality in association with the annual cycle of upwelling-favorable winds and it can be easily observed from high resolution satellite Sea Surface Temperature (SST) distribution (e.g., [12,13]).However, the upwelling front characteristics (i.e., position, magnitude and width) and their spatio-temporal variability have been infrequently documented for the main EBUS, especially at the high frequency scales of variability.Studies of upwelling fronts using time-averaged SST maps underestimate the SST gradients intensity.For example, off Central Chile, the SST gradient intensity over the front computed from seasonal means was found to be ∼0.01 • C km −1 , whereas weekly SST gradients can be one order of magnitude higher, ∼0.1-0.25 • C km −1 [7].
Mesoscale and submesoscale structures (i.e., eddies, filaments, and meanders), common in EBUS, also contribute to create additional SST gradients in the CTZ.They are here referred as secondary fronts as opposed to the main upwelling front.In the early stages of their development, some of these structures have been observed interacting with the main upwelling front [6,7] and allowing cross-shore exchanges of coastal and oceanic diatom species [14].In such cases, the upwelling front is locally deformed by the mesoscale or submesoscale feature.
Several methods have been developed that can automatically identify fronts from SST maps (e.g., [15][16][17][18]).Most of them are based on the detection of horizontal gradients in satellite SST (e.g., [19,20]) or on the statistical analysis of the SST distribution (e.g., [16,[21][22][23]).Statistical methods assess whether a temperature threshold can be defined to separate the SST distribution into two distinct classes.When this condition is fulfilled, the isotherm corresponding to the threshold temperature defines the front in the SST map.Such methods are designed to identify all kinds of SST fronts in a given region and not specifically the main upwelling front.Previous studies in the upwelling region off Central Chile have used a fixed isotherm to define the upwelling front [7,14].However, this has limitations because this value is likely to vary in time, especially under strong seasonality in upwelling at mid-latitudes.Thus, a new method is required to automatically determine for each SST map the most adequate isotherm defining the main upwelling front in EBUS.
Here, we focus on the upwelling front in a mid-latitude region of the Humboldt Current System (HCS), in the area off Point Lavapie (36.5-37 • S, Central Chile).In the HCS, SST is affected by both the wind-induced upwelling and the heat fluxes at the air-sea interface (e.g., [24][25][26][27]).The propagation of coastal-trapped waves (CTW) also plays a role in the SST variability [24,28], in particular during El Niño Southern Oscillation (ENSO) events [29,30].The imprint of CTWs is strong in the Northern HCS, but off Central Chile (30-40 • S) their amplitude is largely attenuated so that the SST can be inferred from a one-dimensional model forced by the local wind stress and solar heating [24].The region off Central-South Chile (35-40 • S) is characterized by a strong seasonality in the coastal wind stress and the wind stress curl, with upwelling-favorable wind conditions from September to April (e.g., [13,31]).The upwelling-induced coastal cold tongue is present from November to May, suggesting a lag in the seasonal response of the SST to the wind conditions [13].This delay has been explained by the concurrent effects of increasing solar heating and upwelling-favorable wind, which result in a warming of the surface layer and a cooling in the 20-80 m layer during austral spring Indeed, the temperature difference between the surface and 80 m depth is only ∼1-2 • C during this season, indicating that the effect of the upwelling on the SST is low.During summer and fall, the water column stratification increases (up to ∼4 • C difference between the same layers) and even weak upwelling events can strongly affect the SST [13].
In central-southern Chile, Point Lavapie (37 • S) has been recognized as a main upwelling center [32].This location has also been associated with the formation of an offshore meander of the coastal equatorward jet and of recurrent mesoscale features [7,[31][32][33].Studying a few cases in this area, Letelier et al. [7] suggested that the upwelling front position corresponds to that of the coastal jet meander.Just north of Point Lavapie, relative high mean frontal probabilities and an enhancement of frontal activity over the shelf and up to 200 km offshore have been observed [34].In this area, the characteristics of the main upwelling front and their frequency of variability remain poorly described.This issue is nevertheless of high relevance in terms of the front influence on coastal planktonic communities and their contribution to high primary production levels, which support particularly large fish landings in this area [35].
The objective of this paper is first to propose a new method for the automatic detection of the main upwelling front using high resolution SST maps in EBUS.This method is then applied to describe and understand the variability of the front characteristics (probability of detecting a front, magnitude, and position) in a small area off Point Lavapie.The extent to which this variability can be related to local atmospheric forcing and oceanic circulation is also explored.The front characteristics are analyzed in Section 3, before the results are discussed in Section 4. Finally, we conclude in Section 5.

Observational Datasets
To detect the main upwelling front, the Group for High Resolution Sea Surface Temperature (GHRSST) Multi-scale Ultra-high Resolution (MUR) SST dataset [36] is used.MUR is a merged product of satellite and in situ observations from the NASA Advanced Microwave Scanning Radiometer-EOS (AMSRE), the Moderate Resolution Imaging Spectroradiometer (MODIS), the US Navy microwave WindSat radiometer, the Advanced Very High Resolution Radiometer (AVHRR), and in situ SST observations from the NOAA iQuam project.This dataset provides daily SST maps with 1 km resolution from 2003 to 2016.
During the PHYTO-FRONT short-term ship survey [14], two cross-shore transects (36.75 • S and 36.5 • S) were carried out during early February 2014.The shallowest CTD temperature measures (3.5 m depth) are used here to compare the upwelling front identified using in situ data and MUR SST.
The oceanic circulation is studied using surface geostrophic velocities computed from SSALTO/DUACS Sea Surface Height (SSH) satellite observations [37] distributed by the Copernicus Marine Environment Monitoring Service (CMEMS).Previously distributed by Aviso+, this product merges data from all altimeter missions: Jason-3, Sentinel-3A, HY-2A, Saral/AltiKa, Cryosat-2, Jason-2, Jason-1, T/P, ENVISAT, GFO, and ERS1.The L4 version provides daily SSH anomaly maps with 25 km resolution from 2003 to 2016 that allows us to study the coastal jet meanders with respect to its mean position.Coastal jets are defined here by alongshore velocities V along above 0.1 m s −1 within 300 km nearshore.Consecutive pixels presenting such values of V along are considered part of the same coastal jet.Two pixels with V along > 0.1 m s −1 and separated by at least one pixel with V along 0.1 m s −1 are considered part of two different jets.All daily maps are analyzed and classified according to the number of coastal jets detected in the area.
Coastal wind stress conditions can be inferred from scatterometer observations.The 2003-2016 period was continuously monitored by several instruments deployed over different periods.However, when analyzing raw data, differences among the missions calibrations introduce bias in the wind stress tendency [38].To solve this issue, wind data from sparse scatterometer fields (ERS, QuikSCAT, and ASCAT) have been processed by the IFREMER and compared to other data sources, such as radiometer data (SSM/I) and atmospheric wind reanalyses (ERA-Interim), to produce a new long time wind product described by Desbiolles et al. [39].This dataset provides 6 h wind stress maps at 1/4 • spatial resolution over the 2003-2016 period.To avoid data contamination by the land-ocean transition, the upwelling-favorable wind conditions are estimated here by averaging the daily alongshore wind stress intensity between 75 and 150 km from the coast and from 37 to 36.5 • S.
Short-wave heat forcing at the air-sea interface is studied using the National Oceanography Centre Southampton (NOCS) Version 2.0 Surface Flux Dataset [40] over the 2003-2014 period.NOCS surface fluxes are computed from daily fields of the surface meteorological parameters (ICOADS Release 2.4 ship data for the years 2003-2006, and ICOADS Release 2.5 for the years 2007-2014) using bulk parameterizations [41].The released dataset provides monthly mean values with a 1 • resolution.Short-wave heat flux conditions are estimated here by averaging the solar heat flux over the 300 km nearshore (area used to detect and analyze the upwelling front; see Section 2.2).

General Methodology
Here, the presence of an upwelling front is detected from the SST daily maps on a reduced latitude range (0.5 • ).The methodology, described in the following, is divided into two steps: In Section 2.2.2, a "frontal zone" (i.e, the zone with high mean cross-shore SST gradient) is first detected.The upwelling front is defined afterwards as the isotherm corresponding to the mean SST over the frontal zone.Then, in Section 2.2.3, the front characteristics are computed from the alongshore mean SST in the frontal zone.Figure 1 shows some examples of front detections.The presence of intense and large mesoscale structures can sometimes create a mean cross-shore SST gradient outside the frontal zone and of comparable intensity (Figure 1b).In such cases, which represent ∼30% of the maps, both the main frontal zone and this secondary frontal zone are identified.The method is summarized on Figure 2.

Front Detection
To identify the frontal zone, SST are latitude-averaged in the 300 km nearshore to compute the mean SST with respect to the distance from the coast (Figure 1).The mean SST at 1 km (T nearshore ) and at 300 km (T o f f shore ) from the coast are computed.When the difference T o f f shore − T nearshore is below 1 • C, we consider that there is no upwelling, and, thus no upwelling front is detected.When T o f f shore − T nearshore is stronger than 1 • C, the mean SST is smoothed using a 30 km running mean.The goal of the latitude average and the cross-shore smoothing is to filter small scales variations associated with mesoscale structures.The smoothed SST is derived to compute the cross-shore SST gradient (grad) and its minimum (grad min ).grad min values are always negative because T nearshore < T o f f shore .We consider that there is no upwelling front that can be detected when grad min is weak (grad min grad LI M with grad LI M = −1.5 × 10 −2 • C km −1 ).Thus, only when grad min < grad LI M , the frontal zone is identified (see Appendix A).
The SST gradient minimum outside the frontal zone (grad min2 ) is also computed.When grad min2 is lower than grad LI M , we look for a second frontal zone applying the methodology described in Appendix A to the SST gradient outside the first frontal zone.After this step, grad outside the two frontal zones is always lower than grad LI M .Then, the two frontal zones are compared to identify the main upwelling front and the secondary front.To this goal, the temperature differences across the first (DT1) and the second (DT2) frontal zone are calculated as the SST differences between the extremities of the frontal zones.The mean SST gradient in the first (grad1) and the second (grad2) frontal zones are also computed.If , then the first (second) front is considered more intense than the second (first) one.In this case, the first (second) front is the main upwelling front and the second (first, respectively) front is the secondary front.For example, in Figure 1c,d Instead of presenting a frontal zone with intense cross-shore SST gradient and low grad outside of it, ∼3% of the analyzed maps present an almost homogeneous grad over the 300 km nearshore.This configuration does not correspond to a well-defined upwelling front separating coastal cold waters from warm offshore waters.Those cases correspond to the detection of a very wide "frontal zone" ( 150 km) and they are discarded from the analysis.For the other maps where a frontal zone is detected, the front is defined on the SST map as the T 0 isotherm, with T 0 being the mean SST over the frontal zone.

Determination of the Front Characteristics
The front parameters are inferred from the smoothed SST profile over the (main) frontal zone.The front width (dl) is computed as the distance between the frontal zone extremities, while the front intensity can be estimated from the SST difference between the frontal zone extremities (DT) and from the cross-front SST gradient ( DT dl ).The position of the front (DX) is computed as the distance from the coast to the center of the frontal zone.
To study the temporal variability of the front characteristics (Section 3), the upwelling front detection method is applied to the 37 • S-36.5 • S area over the 2003-2016 upwelling period.First, for each month, we compute the percentage of daily maps where an upwelling front is actually detected ("probability of detecting a front", Section 3.2.1).Second, to study the front intensity, T o f f shore − T nearshore , DT, dl and DT dl are computed for each daily map over the 2003-2016 upwelling months.Monthly and yearly means, as well as a climatology are then obtained by time-averaging these daily values (Section 3.2.2).Finally, to analyze the impact of the jet meanders on the front position, an histogram of DX is computed for each jet configuration (Section 3.2.3).This histogram is obtained using the front positions detected for all days presenting such jet configuration.

Statistical Validation of the Detected Fronts
To verify whether isotherm T 0 on the daily map actually separates properly the nearshore cold waters from the offshore region (Section 4.2.1), a statistical test is performed following the methodology described by Cayula and Cornillon [16].Pixels included in a box, centered on the front mean position and covering the 0.5 • latitude range with a zonal width of dx + 50 km, are extracted from the daily map.A SST histogram is computed (Figure 3) and used to define two SST classes (SST colder than T 0 and SST warmer than T 0 ).To assess whether this segmentation is statistically relevant, θ and the signal-to-noise (σ) ratios defined by Cayula and Cornillon [16] are computed (Table 1).These ratios compare the difference between the mean SST in each class with the standard deviations around these means: where N 1 (N 2 ) is the number of pixels, µ 1 (µ 2 ) is the mean SST and S 1 (S 2 ) is the standard deviation in the first (second, respectively) class.The segmentation is considered relevant when θ 0.7 and σ 4.0 [16].To evaluate the robustness of our method for the Chile EBUS, we performed this statistical  0.9 5.8 13 December 2010 0.9 5.0 14 November 2016 0.9 6.1

Atmospheric and Oceanic Forcing
Here, we first describe the oceanic and atmospheric conditions over the studied area that could potentially affect the upwelling front.

Local Atmospheric Forcing and Its Impact on the SST
The front characteristics depend on the SST field within 300 km from the coast that is primarily forced by the solar heat flux and the wind stress [24].Figure 4 shows the temporal variability of the SST conditions during the first and last two years of the study period.Variations in the area are dominated by the seasonal cycle.The short-wave flux climatology presents a maximum in January and a minimum in June while the alongshore coastal wind is upwelling-favorable from September to April, with stronger intensity (>0.08 N m −2 ) from November to March and a maximum in January (Figure 5a).Offshore, the SST is little affected by coastal upwelling and the seasonal cycle of T o f f shore is strongly linearly correlated (R = 0.98) to solar heating with a 1.5 month lag (Figure 5b).The seasonal cycle of T nearshore is also correlated (R = 0.98) to the solar heating with the same lag, but its amplitude is reduced by 40% compared to T o f f shore because upwelling is slowing heating down when the solar flux is the strongest.The effect of upwelling can be evidenced from the seasonal cycle of T o f f shore − T nearshore that is maximum in March and minimum in September.This seasonal cycle is correlated to the alongshore wind stress seasonal cycle with a lag of two months (R = 0.97).Thus, although wind conditions are upwelling-favorable during September-April, the imprint of this forcing on the SST is from November to June.This is why we now refer to the November-June period as the "upwelling months".Note that this delay has already been evidenced and explained by [13] (see Section 1).When removing the climatology, the daily nearshore SST anomalies are negatively correlated (R = −0.42)during the upwelling month to the alongshore wind stress intensity averaged over the previous 10 days.This underlines the importance of wind stress forcing on T nearshore for shorter timescales than the seasonal one.
Monthly values of T nearshore anomalies are only weakly correlated to the Oceanic Niño Index (ONI), [42].Indeed, the highest correlations are obtained with a three-month delay and do not exceed 0.24.For example, differences in nearshore SST anomalies between January-June 2003 and November 2003-June 2004 (Figure 4) cannot be explained by the ENSO cycle.Indeed, both periods correspond to neutral ENSO conditions (not shown) while monthly anomalies of T nearshore are weak ( 0. The yearly values of T nearshore are correlated to those of T o f f shore (R = 0.83) but the correlations with the local forcing (both the alongshore wind and the solar heating) are not significant (p-value > 0.05).This suggests that the interannual variations of T nearshore and of T o f f shore are dominated by large-scale variations instead of local atmospheric conditions.However, neither the yearly values of T nearshore nor those of T o f f shore are significantly correlated to the ONI (not shown).
Although correlations between the yearly values of T o f f shore − T nearshore and of the alongshore wind intensity are not significant, on a longer time scale, the trend of the wind conditions and the cross-shore SST difference are consistent.Indeed, in the studied area, the alongshore wind stress intensity presents an intensification after 2007 (Figure 5c), as described by [43], and explained by a poleward displacement of the South East Pacific anticyclone.The alongshore wind stress intensity averaged over the 2007-2016 period is 25% stronger than over the years 2003-2006 while T o f f shore − T nearshore is enhanced by 24% (Section 3.2.2).This 0.5 • C increase of T o f f shore − T nearshore is mainly due to a 0.6 • C cooling of T nearshore (T o f f shore remaining almost unchanged, not shown) due to the coastal upwelling intensification.

Coastal Jet Circulation
It has been suggested that the position of the SST upwelling front is related to the position of the coastal jet [7].Here, the coastal circulation is studied using surface geostrophic currents (see Section 2.1) to identify the most frequent coastal jet configurations.First, the number of coastal jets present within 300 km from the coast are computed.No jets are detected for 21% of the daily maps, only one jets is detected for 64% of cases, and 15% of maps present two jets.From all maps where only one jet is detected, a distance from the coast Probability Density Function (PDF) is computed using all pixels included in the jet (Figure 6a).The PDF presents two maxima at ∼25-75 and 200-225 km from the coast, separated by a minimum at ∼100 km.This shows that two positions for the jet are most probable (within 100 km from the coast, centered ∼50 km, or offshore from 100 km from the coast, centered ∼225 km).For all maps where two jets are detected, pixels with alongshore velocities above 0.1 m s −1 are classified whether they belong to the most nearshore or the most offshore jet.A first distance from the coast PDF is computed using all pixels located in the most nearshore jets, and a second one is computed using all pixels located in the most offshore jets (Figure 6b).These PDFs show that, when two jets are detected, the most nearshore is preferentially located at ∼50 km from the coast and the most offshore is at ∼250 km.Almost all most nearshore (offshore) jets are located within (offshore from, respectively) ∼150 km from the coast.We now call "nearshore jet" ("offshore jet") a jet with a center within (offshore from, respectively) 125 km from the coast.The distance from the coast to the jet center is defined as the mean distance from the coast of all pixels included in the jet.Climatologies of the numbers of maps with no jet, with one nearshore jet, with one offshore jet, and with two jets ("jet configurations") are computed for the upwelling months (Figure 6c).Maps with no jet are more frequent during spring (November-December) and late fall (June), when they represent ∼30% of the daily maps.From January to May, maps without any jet only represent ∼10-20% of cases.The presence of only one nearshore jet represents ∼55% of the maps during November.This proportion decreases during the upwelling months down to 4% in June.On the contrary, the presence of only one offshore jet is rare (10% of the cases) during November and then increases during the next months up to ∼70% of maps in April.Finally, two jets are detected for ∼20% of cases from January to April (with a maximum of ∼30% in May) and for ∼5-10% of cases during spring and late Fall.and (b) two jets in the 300 km nearshore.Black (red) line represents the PDF for pixels associated with the nearshore (offshore, respectively) jet.(c) Numbers of daily maps (in %) with no jet in the 300 km nearshore (black line), with only one jet in the 300 km nearshore, located within (red line) and offshore from (magenta line) the 125 km nearshore, and with two jets in the 300 km nearshore (cyan line).Numbers of maps where normalized by the number of analyzed maps during the selected month.

Probability of Detecting a Front
The percentage of daily maps where an upwelling front is actually detected is above 50% from November to June (Figure 7a), period that corresponds to the "upwelling months" (Section 3.1.1).This is why the front characteristics presented in this paper are computed only during those months.The number of daily maps with an upwelling front is maximum in March, where almost 100% of the maps exhibit an upwelling front.This is consistent with the alongshore wind stress maximum in January and the two-month lag between the wind forcing and its impact on the SST field (Section 3.1.1).
The yearly values of the number of days with an upwelling front evidence an increase after 2007 (Figure 7b) and the mean value over the 2007-2016 period is 14% stronger than over the 2003-2006 period.This is consistent with the increase in upwelling-favorable wind conditions over the same period (Section 3.1.1).However, the standard deviations around these means are large and a longer time series is necessary to confirm this tendency.

Front Intensity
Monthly values of the front intensity show that DT is highly correlated with T o f f shore − T nearshore (Table 2), suggesting a driving of DT by T o f f shore − T nearshore .Not surprisingly, dl is strongly correlated to DT and so is DT dl .These results evidence that a stronger (weaker) T o f f shore − T nearshore is associated with enhanced (reduced) DT, a wider (narrower) front and also a stronger (weaker, respectively) DT dl .The seasonal cycles of DT, DT dl , and dl show a maximum during February and March (Figure 8).Indeed, DT is increased by 81% in March compared to November.This is consistent with the T o f f shore − T nearshore intensification during summer (Figure 8a) that results from the SST response to the increased upwelling-favorable wind conditions (Section 3.1.1).Note, however, that the increase in T o f f shore − T nearshore between March and November (128%) is stronger than the increase in DT.This is likely due to the seasonal cycle of SST gradients associated with mesoscale activity.Indeed, mesoscale structures can create cross-shore SST gradients out of the main upwelling front e.g., (Figure 1b).In those cases, T o f f shore − T nearshore can be approximated by the sum of DT and the SST differences across all these secondary fronts.The occurrence of such secondary SST gradients is enhanced during summer [34], increasing the difference between DT and T o f f shore − T nearshore .Between March and November, the front is also broadened (Figure 8b).However, the increase in dl (40%) is reduced compared to DT, so that the mean cross-shore gradient across the front DT dl is also slightly enhanced (31%, Figure 8c).Yearly values show the increase of T o f f shore − T nearshore discussed in Section 3.1.1(Figure 9a).Resulting from the stronger T o f f shore − T nearshore , DT is also enhanced by 17% during 2007-2016 compared to the 2003-2006 period (Figure 9b).A slight increase (12%) in the front width can also be observed (Figure 9c), although the high interannual variability makes it hardly distinguishable.As a result, DT  dl is almost unchanged and a longer time series would be necessary to identify any increase.

Front Position
The position of the front (DX) presents high temporal variability, and can vary from 100 km to 200 km from the coast within a few days (Figure 4).The standard deviations around the yearly values are large (∼100-150 km) and data do not show any trend over the 2003-2016 period (Figure 10).There was no direct relation between the front position and the coastal wind or heat fluxes (the correlation between the two time series are not significant).We now study the relation between the coastal jet configuration and the front position.Figure 11a shows that the presence of only one nearshore jet is associated with a front located very nearshore, between 50 and 100 km from the coast for ∼80% of the cases.In this jet configuration, DX > 125 km represents ∼10% of the cases.When two jets are detected, the front is also preferentially located between 50 and 100 km from the coast (∼65% of the cases), but more offshore fronts (>125 km) are more frequent (∼30%) than when only one nearshore jet is detected.When no jet are detected, the front position is located from 50 to 125 km offshore for 80% of days.Finally, when only one offshore jet is identified, the front position is much less restrained nearshore, and the proportion of fronts is ∼10% for all positions between 50 and 250 km (with a maximum of 15% at 50 km from the coast).During the "upwelling months", the probability of detecting only one nearshore jet (that is associated with more nearshore fronts) decreases and the probability of detecting one offshore jet (associated with more offshore fronts) increases (Figure 6c).This explains the fact that the mean distance from the coast to the front (and the standard deviation) increases during the "upwelling months", from 80 km in November to 150 km in May (Figure 11b).In EBUS, a multiplicity of intense SST gradients are present near the frontal zone.This make difficult the upwelling front identification among all the secondary fronts.Our method is able to detect this main front, despite the large diversity of situations that can be observed.For example, in Figure 12a,c, the geostrophic circulations are very similar, with a large anticyclone centered on 36.3 • S, 77.3 • W in both cases.Nevertheless, the SST fields are different.Indeed, in Figure 12a,b, the coastal upwelled waters advected by the eddy result in a relatively homogeneous (cold) SST from the coast to the western eddy side, and the upwelling front is located offshore from the anticyclone.In Figure 12c, although cold waters are also advected by the eddy, SST on the eastern eddy side is still warmer than closer to the coast.As a result, two fronts are detected by the algorithm on each eddy side.The nearshore front, which is more intense, is the main upwelling front (Figure 12d).On 11 January 2005, a cyclone advects coastal waters offshore, creating a cold filament (Figure 12e).Although the SST inside and outside the filament differs by ∼1 • C, the reduced size of this feature leads to a very weak imprint on the mean SST (Figure 12f) and no secondary front is detected by our method.On 10 March 2009, the mean SST (Figure 12g) shows the presence of two strong SST gradients ∼50 km and 150 km from the coast.The offshore frontal zone presents a stronger SST difference and a stronger mean SST gradient (Figure 12h).Thus, our method considers that the coastal SST water extends from the coast to the offshore front, which is identified as the main upwelling front.In some cases, a mesoscale structure can deform the upwelling front, so that the front in this area is no longer oriented in the alongshore direction and can even be locally cross-shore oriented (Figure 12i).Our method is able to capture such features and to define a statistically consistent front (Figure 12i,j and Table 1).This is because the front is defined as an isotherm (thus, there is no condition over the front direction) and because, even in such cases, there are two very low cross-shore gradient zones, nearshore from the filament base (in the cold tongue), and offshore from the filament extension.Thus, the algorithm is able to compute the mean SST between these two zones and to define the front.The algorithm has proved to be robust off Central Chile, where it was able to distinguish the main upwelling front from this strong mesoscale field.Mesoscale SST gradients in this region are particularly strong [34].Thus, the method can all the more be applied to SST fields presenting less strong SST mesoscale field such as other areas, time-averaged SST maps, or lower resolution products.In addition, our method analyzes the SST daily conditions to compute case-by-case the most adequate isotherm T 0 that defines the front.This is why this methodology can be applied to different time-periods and other EBUS regions without major code modifications.

Limitations
Nevertheless, our method includes some arbitrary parameters (e.g., grad LI M ) that were empirically defined by applying the algorithm to HCS areas included between 38 • S and 35 • S.These parameters may need to be adjusted for other areas and make the method somewhat subjective.For this reason, we included in the algorithm a step that analyzes the frontal zone interval identified to verify afterward whether the interval extremities actually correspond to edges in the cross-shore SST gradient (see Appendix A).
The method described here only detects the main upwelling front in a 0.5 • latitude range.This latitude range was also empirically determined.Indeed, the meridional extension of the studied region must be large enough to filter noise when computing the mean SST with respect to the distance from the coast.However, this extension must also be reduced enough so that one isotherm can actually define the front on the daily SST map (this isotherm must be located in the strong cross-shore gradient area for all latitudes in the studied region), and 0.5 • was the largest meridional extension that fulfill this requirement.To obtain a front defined over a larger region, the following method could be applied: for each latitude lat 1 , the algorithm is applied to the area between lat 1 −0.25 • and lat 1 + 0.25 • and the value of the front temperature T 0 (lat 1 ) is computed.The upwelling front is then defined on the SST map for each latitude lat 1 as the point where SST = T 0 (lat 1 ).
When clouds are present over a pixel, the high resolution infrared measurement are unavailable and the available micro-wave satellite measurements present lower resolution (25 km).This could affect the front detection and the measures of dl and DT dl .However, considering its typical size (∼50-100 km, Figure 8b), the upwelling frontal zone is large enough to be well defined even on 25 km-resolution observations.In addition, Morales et al. [44] found that missing data in the case of satellite surface Chl-a off central-southern Chile, which can be mostly attributed to cloud coverage, represent 55-60% of the data, with seasonality not influencing these results.Thus, we assume that the impact of clouds on the front characteristics and their seasonality is weak.Cayula and Cornillon (1992) The front detection method from Cayula and Cornillon [16] was designed to detect all fronts present in a given area and not just the main upwelling front.However, it provides a method to assess whether a given isotherm on a SST map separates two water regions with statistically different SST characteristics.To evaluate whether the isotherm T 0 (used in our method to define the upwelling front) do represent a threshold separating the coastal cold water from offshore, this statistical test is applied to the fronts detected over 180 maps (Section 2.3).In 55 maps, no front is detected (31% of the total).A frontal zone larger than 150 km is detected in five maps (3% of the total), and, thus, discarded from the test.Over the 120 remaining maps where a relevant front has been detected, 19 maps (15% of 120) do not pass the test (because θ < 0.7 or because the signal-to-noise ratio is <4.0).However, among these 19 cases, 12 of them correspond to a front with a nearshore extremity closer than 25 km from the coast.In those cases (e.g., Figure 13 and Table 1), there are too few pixels between the front and the coast to form a population statistically comparable to the class of pixels warmer than T 0 .The seven other maps (6% of the 120 maps) correspond to cases where the algorithm fails to properly identify an upwelling front.

Consistency with In Situ Measurement
The front position computed with the new method is now compared to the front identified using in situ observations (Figure 14).During the PHYTO-FRONT cruise, in both cross-shore transects, ship measurement evidenced a narrow zone of strong surface temperature and salinity gradients between two consecutive stations [14].These zones correspond to SST values ∼17 • C for both transects and were identified as the upwelling front.Although in situ measurements at 3.5 m depth and satellite SST observations differ by up to 1 • C in the offshore region, the front automatically detected by our method is located between the two stations delimiting the front extension in Morales et al. [14].

Errors and Uncertainties
The front characteristics show high variability at all temporal scales and the standard deviations around the climatological and yearly values presented here are large.As a result, the relatively small variations in the seasonal cycle (Figure 8c) and in the trend (Figure 9d) during 2003-2016 for DT dl are located within the standard deviation, and, thus, are not statistically significant.In the next years, the inclusion of new yearly values may reduce the uncertainty in the evolution of DT dl so that a significant trend may be identified out of a longer timeseries.The front position DX also shows large standard deviations and no trend can be identified (Figure 10).On the contrary, the variations in the seasonal cycles and the trends of the probability of detecting a front (Figure 7), DT (Figures 8a and 9b) and dl (Figures 8b and 9c) are significant.

Jet Configuration and Front Position
A striking result of this article is the statistical analysis of the coastal jet position and its relation with the position of the upwelling front.Four most common configurations have been identified (no jet, one nearshore jet, one offshore jet, and two jets), which is consistent with the configurations observed in the studied cases from previous studies [7,32,33].Our study shows that the most common configuration is the presence of one offshore jet.We also evidence a seasonal variability of the proportion of each configuration, consistent with the seasonal current maps from Aguirre et al. [31], that shows a nearshore jet in spring, two jets during summer, and an offshore jet in fall.This westward propagation of the coastal jet north of Point Lavapie during summer has been related to the negative wind stress curl conditions [31].Coastal Ekman pumping facilitates the jet separation from the coast near a cape (the separation occurs earlier and the jet propagates farther) [45].The offshore jet propagation during the upwelling season is also a known feature in the California EBUS, where the importance of wind stress curl [46,47], interaction with alongshore topography [48] and beta-effect (Rossby wave dynamics) [5,46] have been underlined.The main upwelling front is commonly described as co-located with the coastal jet (e.g., [7]).We show here that this is true when the jet is nearshore, nevertheless the presence of only one jet located offshore is associated with very different front positions.Nearshore jets (and, thus, colocated fronts) are more frequent during spring/early summer.At the beginning of the upwelling season, the coastal jet formation has been directly linked to the upwelling-induced cross-shore density gradient [45,[47][48][49].The coastal jet is created by geostrophic adjustment to the SST front, and thus is co-located with it.Offshore jets are more frequent at the end of the upwelling period, and the jet meanders are no longer related to the local SST gradient.

Reinforcement of the Upwelling Front in EBUS
Our study shows a reinforcement of the upwelling front (more frequent front, with a more intense cross-front SST difference) linked to the intensification of alongshore coastal winds at 37 • S after 2007.However, this results may not be extrapolated to other regions of the HCS or other EBUS.Indeed, the wind intensification has been linked to a poleward displacement of the South East Pacific anticyclone.This displacement favors upwelling conditions off Central Chile but may decrease the coastal wind in the northern part of the HCS [50], and the future of coastal winds in the HCS under climate change conditions remains an open question [50][51][52][53].

Conclusions
An algorithm to identify the main upwelling front from high resolution SST daily maps is used here to automatically compute the front magnitude, width and position.The main asset of this method is its ability to automatically identify the main front among all strong SST gradients created by mesoscale and submesoscale features.In addition, although some arbitrary parameters may need to be adjusted, the method is general enough to be applied to different Eastern Boundary Upwelling Systems.However, clouds making unavailable infrared high resolution measurements, the cloud coverage may affect the results, in particular when studying very narrow (<25 km) upwelling fronts.
This algorithm has been applied to an area in the highly productive region off Central Chile (37 • S-36.5 • S) during 2003-2016 to analyze the front characteristics and their temporal variability.We show that the intensity of upwelling-favorable wind conditions is the main driver of the probability of detecting a front, and also of the front magnitude and its width, for both the seasonal cycle and the long-term trend.The coastal wind intensification after 2007 (+25%) is responsible for an increase in the percentage of days presenting a main upwelling front (probability of detecting a front increased by 14%), in the SST difference across the front (by 17%) and in the front width (by 12%), while changes in the cross-front SST gradient are very low.On the contrary, the position of the front is highly variable within a few days and does not show any trend over the studied period.In addition, this position is not related to the wind conditions but to the coastal jet meanders.The coastal jet(s) location(s) during the upwelling months has (have) been detected to obtain for the first time a statistics over the years 2003-2016 and to identify the most frequent configurations.We show that the presence of a very nearshore jet (∼25-75 km from the coast) is associated with a front located very close while when the coastal jet meander is located offshore from 125 km from the coast, the front is not co-located with the jet.
Similar front analysis could be performed for other Eastern Boundary Upwelling System areas to monitor the long-term variability associated with climate change and to identify similarities and differences in their evolution.The analysis of other oceanic fields (from satellite or in situ observations) in the Coastal Transition Zone can now be related to the front characteristics defined here.This would be useful to study the impact of the upwelling front fine-scale variations on the biological components and activity.
we restrain the interval around x min : the value of r is increased by 0.1 and we determine the new interval where grad > r • grad LI M .At this step, the "frontal zone" identified corresponds to the interval around grad min where the mean grad is below a threshold (grad LI M if dx 50 km, and 0.7 dx if dx < 50 km).Note that this threshold is in any case lower than grad LI M .The value of grad LI M corresponds to a strong cross-shore grad for the HCS, it is negative enough to make sure that the interval is not too wide and includes too weak values of grad near its extremities.However, x 1 and x 2 may need to be adjusted so that they actually delimit the strong grad extension.This is necessary for two typical cases.First, when grad is weak (grad min > −0.025 • C km −1 , e.g., Figure A2a-c), the threshold is sometimes too low and the interval does not include all the area of strong grad.In this case, the frontal zone should be widened; this represents ∼2% of the analyzed maps.Second, when the peak around grad min is asymmetric (e.g., Figure A2d-f), the frontal zone should be reduced on one side, so that the interval only includes the strongest values of grad.This case corresponds to ∼10% of the analyzed maps.To adjust x 1 and x 2 , the algorithm first analyzes whether the frontal zone should be widened.grad is averaged over x 1 − dx 2 , x 1 and over x 2 , x 2 + dx 2 .If one of these values is lower than 1 3 grad min , the front is widened.To widen the frontal zone, grad is smoothed using a running mean of size dx.The frontal zone is now determined as the interval where dx-smoothed grad is lower than 30% of its minimum value.The algorithm also evaluates the peak symmetry by averaging grad over x 1 , x 1 + dx 2 and over x 1 + dx 2 , x 2 .If the difference between the two averages is larger than 40% of the peak intensity (grad min − grad(x 1 )+grad(x 2 ) 2 ), then the interval is reduced on the interval side with the weakest averaged grad.To reduce the frontal zone, grad is smoothed using a running mean of size dx 2 .The interval extremity that needs to be reduced is now determined by the point where the dx 2 -smoothed grad is

Figure 1 .Figure 2 .
Figure 1.Daily SST (colors, in • C) observations and geostrophic currents (arrows) derived from SSH observations for: (a) 13 December 2010; (c) 24 March 2007; and (e) 14 November 2016.Black box indicates the studied area and the main (secondary) upwelling front is marked with plain (dotted, respectively) magenta line.Magenta dot indicates Point Lavapie.Mean daily SST (in • C, averaged over the studied latitude range) with respect to the distance from the coast (black line) for: (b) 13 December 2010; (d) 24 March 2007; and (f) 14 November 2016.Dotted blue line represents the smoothed mean SST and magenta (red) circles mark the extremity (center, respectively) of the main frontal zone.Green circles mark the secondary frontal zone extremities and center.
, two fronts are identified around 75 km and 200 km from the coast.The offshore front corresponds to a filament created by a deflection of the coastal jet.Over the nearshore frontal zone, DT1 = 2.1 • C and grad1 = 0.03 • C km −1 , while over the offshore one, DT2 = 2.0 • C and grad2 = 0.02 • C km −1 .The nearshore front being more intense ( 86), it is identified by the algorithm as the main upwelling front.

Figure 3 .
Figure 3. SST histograms of the region within 50 km from the front detected by the algorithm for: (a) 13 December 2010; (b) 23 March 2007; and (c) 14 November 2016.Red line represents T 0 .

Figure 4 .Figure 5 .
Figure 4. Mean daily SST (in • C) between 37 • S and 36.5 • S with respect to the distance from the coast and time during the upwelling months (November-June) for years: (a) 2003-2004; and (b) 2015-2016.Black line marks the center of the stronger gradient zone ("upwelling front") identified by the algorithm.
5 • C) during the first period and warm during the second one ( 0.5 • C during Spring and even 1 • C during Fall 2004).The ENSO influence is more notable during 2014 and 2015, when the ONI shows warm (ONI 0.5 • C) and increasing ENSO conditions from November 2014 (not shown).Consistently with the three-month delay between the ONI and T nearshore anomalies, the latter increases from February 2015.The ONI values 1.5 • C from August 2015 to March 2016 (not shown) indicates a very strong ENSO phase that would affect T nearshore from November to June.T nearshore anomalies are indeed warm and 0.5 • C, except in December, February and March (|T nearshore | 0.5 • C, neutral conditions).

Figure 6 .
Figure 6.(a,b) Distance from the coast Probability Density Functions (PDF) of pixels with alongshore surface currents above 0.1 m•s −1 in the 300 km nearshore.PDFs are computed using the 2003-2016 geostrophic currents daily maps during upwelling months with: (a) only one jet in the 300 km nearshore;and (b) two jets in the 300 km nearshore.Black (red) line represents the PDF for pixels associated with the nearshore (offshore, respectively) jet.(c) Numbers of daily maps (in %) with no jet in the 300 km nearshore (black line), with only one jet in the 300 km nearshore, located within (red line) and offshore from (magenta line) the 125 km nearshore, and with two jets in the 300 km nearshore (cyan line).Numbers of maps where normalized by the number of analyzed maps during the selected month.

Figure 7 .
Figure 7. Number of daily SST maps presenting an upwelling front compared to the total number of maps (in %): (a) seasonal cycle (climatology over the 2003-2016 period, error bars represent the standard deviations); and (b) yearly values (red line).Blue lines represent the averages (horizontal lines) and standard deviations (vertical lines) over the 2003-2006 and the 2007-2016 periods.

Figure 8 .
Figure 8. Climatology of mean front characteristics using daily SST maps during the 2003-2016 upwelling months: (a) SST difference (in o C) across the front (DT, plain blue line) and SST difference between the coast and 300 km offshore (T o f f shore − T nearshore , dotted green line); (b) front width (dl, in km); and (c) cross-front SST gradient ( DT dl , • C km −1 ).Error bars represent the standard deviations.

Figure 9 .
Figure 9. Mean front characteristics (yearly averages for the fronts detected on daily SST maps, blue line).Black lines represent the averages (horizontal lines) and standard deviations (vertical lines) over the 2003-2006 and the 2007-2016 periods: (a) SST difference (in • C) between the coast and 300 km offshore; (b) SST difference (in • C) across the front; (c) Front width (in km); and (d) cross-front SST gradient (in • C km −1 ).

Figure 10 .
Figure 10.Seasonal cycle of the distance between the coast and the upwelling front (in km).Error bars represent the standard deviations.

Figure 11 .
Figure 11.(a) Histograms of DX computed using all daily maps with the same jet configuration.DX is normalized by the total number of maps with the considered jet configuration.Black line corresponds to the maps with no jet in the 300 km nearshore.Red (magenta) line corresponds to the maps with only one jet in the 300 km nearshore, located within (offshore from, respectively) 125 km from the coast.Cyan line corresponds to the maps with two jets in the 300 km nearshore.(b) Seasonal cycle of the mean front position.Error bars represent the standard deviations.

Figure 13 .
Figure 13.(a,b) Same as Figure 1a,b; and (c) same as Figure 3, for 15 January 2007 and the 38 • S-37.5 • S zone as studied area.

Figure 14 .
Figure 14.(a) Mean satellite SST (colors, in • C) during 4-6 February 2014.Black and magenta crosses mark the PHYTO-FRONT sampling stations positions.Magenta line represents the front detected on the mean SST map for the 36.875 • S-36.375 • S zone.(b) Mean satellite SST at 36.5 • S during 5-6 February 2014 period (blue line).Black and magenta crosses indicate the 3.5 m depth in situ temperature along the 36.5 • S transect.Red dot marks the value of T 0 computed using the 05-06 February 2014 SST included in the 0.5 • latitude-large zone around 36.5 • S. (c) Same as (b) but for latitude 36.75 • S and 5-6 February2014 period.Magenta crosses in (a-c) indicate the stations delimiting the front according to Morales et al. [14].

Figure A1 .
Figure A1.Scheme of the frontal zone identification from the smoothed SST profile.

Table 1 .
θ and σ ratios computed for several daily SST maps (area within 50 km from the detected upwelling front) and using T 0 as the threshold temperature.

Table 2 .
Correlations between the time-averaged (monthly mean, yearly mean and seasonal cycle) front characteristics, computed from the daily maps.