Recent Changes in the Low-Level Jet along the Subtropical West Coast of South America

: Surface winds along the subtropical west coast of South America are characterized by the quasi-weekly occurrences of low-level jet events. These short lived but intense wind events impact the coastal ocean environment. Hence, identifying long-term trends in the coastal low-level jet (CLLJ) is essential for understanding changes in marine ecosystems. Here we use ERA5 reanalysis (1979–2019) and an objective algorithm to track anticyclones to investigate recent changes in CLLJ events off central Chile (25–43 ◦ S). Results present evidence that the number of days with intense wind ( ≥ 10 ms − 1 ), and the number and duration of CLLJ events have signiﬁcantly changed off central Chile in recent decades. There is an increase in the number of CLLJ events in the whole study area during winter (June-July-August; JJA), while during summer (December–January–February; DJF) a decrease is observed at lower latitudes (29–34 ◦ S), and an increase is found at the southern boundary of the Humboldt system. We suggest that changes in the central pressures and frequency of extratropical, migratory anticyclones that reach the coast of South America, which force CLLJs, have played an important role in the recent CLLJ changes observed in this region. data and reanalysis output Higher correlations were observed austral they were lower in austral These highlight the ability of representing the strong wind our area, and, therefore, southerly


Introduction
Low-level winds along the west coast of subtropical South America are strongly influenced by the Southeast Pacific Anticyclone (SPA), a quasi-stationary atmospheric high-pressure cell that drives equatorward winds off Perú and Chile ( Figure 1). In this vast region the seasonal cycle of these low-level winds differs between its northern and southern portions, mainly driven by the latitudinal migration of the SPA through the year. During austral winter (JJA) the SPA moves northward, producing mean westerly winds in south-central Chile (35-42 • S). During austral summer (DJF) the SPA moves southward, producing mean equatorward alongshore winds as far south as 42 • S [1,2]. Off northern Chile (18-28 • S) southerly winds prevail through the year, but they are rather weak and stable, with low synoptic and seasonal variability [3]. During austral spring and summer there is a quasi-weekly occurrence of a southerly atmospheric coastal low-level jet (CLLJ) alternating with periods of relaxation or even weak northerly winds [4,5]. This CLLJ is characterized by a band of southerly winds in excess of 10 ms −1 extending up to 1000 km along the coast but only 200-300 km in the cross-shore direction [6]. Vertically, the strong winds encompass the whole marine boundary layer (MBL), peaking just below the base of the capping subsidence temperature inversion at about 1 km above the surface. Studies in other ocean eastern boundaries indicate that CLLJs are associated with the quasi-stationary subtropical high-pressure cells concomitant with an inland thermal low, which enhance the cross-shore sea level pressure gradient (e.g., [7]). The CLLJ off central Chile, however, is driven by the alongshore surface pressure gradient between a coastal low over north-central Chile and a migratory anticyclone farther south, as illustrated for a case of study in January 2007 ( Figure  2). This is because the very steep coastal topography breaks down the geostrophic balance (that otherwise would result in an easterly flow), thus accelerating the equatorward alongshore low-level flow, only constrained by surface friction [1,5,6]. During austral spring and summer there is a quasi-weekly occurrence of a southerly atmospheric coastal low-level jet (CLLJ) alternating with periods of relaxation or even weak northerly winds [4,5]. This CLLJ is characterized by a band of southerly winds in excess of 10 ms −1 extending up to 1000 km along the coast but only 200-300 km in the cross-shore direction [6]. Vertically, the strong winds encompass the whole marine boundary layer (MBL), peaking just below the base of the capping subsidence temperature inversion at about 1 km above the surface. Studies in other ocean eastern boundaries indicate that CLLJs are associated with the quasi-stationary subtropical high-pressure cells concomitant with an inland thermal low, which enhance the cross-shore sea level pressure gradient (e.g., [7]). The CLLJ off central Chile, however, is driven by the alongshore surface pressure gradient between a coastal low over north-central Chile and a migratory anticyclone farther south, as illustrated for a case of study in January 2007 ( Figure 2). This is because the very steep coastal topography breaks down the geostrophic balance (that otherwise would result in an easterly flow), thus accelerating the equatorward alongshore low-level flow, only constrained by surface friction [1,5,6]. The CLLJ plays an important role in the regional climate of the midlatitude and su tropical western continental coastal areas (e.g., [8]), impacting directly on the coastal oce [9][10][11][12][13][14]. The CLLJ prompts intense coastal upwelling, introducing high synoptic-scale v iability in sea surface temperature (SST), ocean diffusivities and surface heat fluxes alo the coast [10][11][12]15] in the upper-ocean circulation and biogeochemical cycles [16,17]. D ing austral spring and summer, a well-defined near-surface ocean equatorward flow w a jet-like structure develops, following the climatological position of the CLLJ axis [9,1 showing that during summer, in ocean eastern boundaries, wind sea wave heights comparable to the swell wave heights. Particularly at the Humboldt system, authors [14] have found that the CLLJ strongly modulates the seasonal cycle of significant wa heights.
Physical mechanisms associated with the generation of the CLLJ include comp regional-to large-scale atmosphere-ocean-land interactions. Therefore, how clim The CLLJ plays an important role in the regional climate of the midlatitude and subtropical western continental coastal areas (e.g., [8]), impacting directly on the coastal ocean [9][10][11][12][13][14]. The CLLJ prompts intense coastal upwelling, introducing high synoptic-scale variability in sea surface temperature (SST), ocean diffusivities and surface heat fluxes along the coast [10][11][12]15] in the upper-ocean circulation and biogeochemical cycles [16,17]. During austral spring and summer, a well-defined near-surface ocean equatorward flow with a jet-like structure develops, following the climatological position of the CLLJ axis [9,13], showing that during summer, in ocean eastern boundaries, wind sea wave heights are comparable to the swell wave heights. Particularly at the Humboldt system, authors in [14] have found that the CLLJ strongly modulates the seasonal cycle of significant wave heights.
Physical mechanisms associated with the generation of the CLLJ include complex regional-to large-scale atmosphere-ocean-land interactions. Therefore, how climate change will impact the future behavior of the CLLJ is a difficult task and a crucial environmental issue. Climate projections have been used for assessing the responses of CLLJs to climate change [18][19][20]. Although the CLLJ is a regional-scale process not resolved in coarse resolution global circulation models (GCMs), relative changes in the synoptic-scale patterns that force CLLJ events, well represented in GCMs, show a poleward shift of midlatitude migratory anticyclones in the Canary and Humboldt systems [21]. Several efforts have been conducted in order to characterize trends in mean southerly winds during the last decades [22][23][24][25]. However, knowledge about recent trends in the CLLJ in the southern Humboldt system is still limited (e.g., [8]).
Here we use ERA5 reanalysis data of sea level pressure and 10 m above-the-surface winds to characterize the evolution of intense equatorward coastal winds (as a proxy of the CLLJ) during the last four decades (1979-2019) off central Chile. Changes in frequency and duration of the CLLJ are reported as a function of latitude and related to changes in central pressures and density of extratropical migratory anticyclones. Documenting the contemporaneous change in the CLLJ is important in its own, as it can assist interpretation of other coastal variables and serve as a stringent benchmark for model validation. For comparison purposes, we consider available wind and sea level pressure records at coastal stations. The rest of the paper is organized as follows: In Section 2 we describe the data used and analyses performed. The main results are presented in Section 3, which is followed by a discussion in Section 4. Finally, the main findings of this work are summarized in Section 5.

Data Sets
We focused on the region covering the coast of Chile between 25 and 44 • S, where the CLLJ is found more frequently [4]. We used hourly mean fields of meridional (V10) and zonal (U10) winds 10 m above the surface, and sea level pressure (SLP) data from the ERA5 reanalysis for the period 1979-2019. This reanalysis is produced by the European Centre for Medium-Range Weather Forecast (ECMWF) and consists of a global-scale regular grid with a 0.25 × 0.25 lat-lon resolution [26]. We used the ERA5 because it offers the best performance for surface wind speed and variability [27]. Data were downloaded from the Copernicus Climate Data Store (https://cds.climate.copernicus.eu/).
In addition, SLP and winds from 7 meteorological stations along the Chilean coast were also analyzed ( Figure 3). Meteorological data were obtained from the Meteorological Service of the Chilean Navy (Dirección General del Territorio Marítimo y de Marina Mercante, DIRECTEMAR) except data from Punta Lengua de Vaca station that were obtained from the Center for Advanced Studies in Arid Zones (Centro de Estudios Avanzados en Zonas Áridas, CEAZA). All the stations have a temporal resolution of 3 h, except Punta Lengua de Vaca, where 30-min records are available (Table 1). In the latter case, the 30-min time series was transformed into an hourly mean time series. Data-void periods and outliers were discarded. latitude migratory anticyclones in the Canary and Humboldt systems [21]. Several efforts have been conducted in order to characterize trends in mean southerly winds during the last decades [22][23][24][25]. However, knowledge about recent trends in the CLLJ in the southern Humboldt system is still limited (e.g., [8]).
Here we use ERA5 reanalysis data of sea level pressure and 10 m above-the-surface winds to characterize the evolution of intense equatorward coastal winds (as a proxy of the CLLJ) during the last four decades (1979-2019) off central Chile. Changes in frequency and duration of the CLLJ are reported as a function of latitude and related to changes in central pressures and density of extratropical migratory anticyclones. Documenting the contemporaneous change in the CLLJ is important in its own, as it can assist interpretation of other coastal variables and serve as a stringent benchmark for model validation. For comparison purposes, we consider available wind and sea level pressure records at coastal stations. The rest of the paper is organized as follows: In Section 2 we describe the data used and analyses performed. The main results are presented in Section 3, which is followed by a discussion in Section 4. Finally, the main findings of this work are summarized in Section 5.

Data Sets
We focused on the region covering the coast of Chile between 25 and 44 °S, where the CLLJ is found more frequently [4]. We used hourly mean fields of meridional (V10) and zonal (U10) winds 10 m above the surface, and sea level pressure (SLP) data from the ERA5 reanalysis for the period 1979-2019. This reanalysis is produced by the European Centre for Medium-Range Weather Forecast (ECMWF) and consists of a global-scale regular grid with a 0.25 × 0.25 lat-lon resolution [26]. We used the ERA5 because it offers the best performance for surface wind speed and variability [27]. Data were downloaded from the Copernicus Climate Data Store (https://cds.climate.copernicus.eu/).
In addition, SLP and winds from 7 meteorological stations along the Chilean coast were also analyzed ( Figure 3). Meteorological data were obtained from the Meteorological Service of the Chilean Navy (Dirección General del Territorio Marítimo y de Marina Mercante, DIRECTEMAR) except data from Punta Lengua de Vaca station that were obtained from the Center for Advanced Studies in Arid Zones (Centro de Estudios Avanzados en Zonas Áridas, CEAZA). All the stations have a temporal resolution of 3 h, except Punta Lengua de Vaca, where 30-min records are available (Table 1). In the latter case, the 30min time series was transformed into an hourly mean time series. Data-void periods and outliers were discarded.   Furthermore, we used wind data from the SeaWinds scatterometer onboard QuikSCAT between July 1999 and November 2009. Daily winds data with a spatial resolution of 0.5 • Were obtained from Centre d'Exploitation et de Recherche Satellitaire d'Archivage et de Traitement (CERSAT), at Institut Francais de Recherche pour l'Exploitation de la Mer (IFREMER). This dataset is available from ftp://ftp.ifremer.fr/ifremer/cersat/products/ gridded/mwf-quikscat/, and for more information the reader is referred to the user manual at ftp://ftp.ifremer.fr/ifremer/cersat/products/gridded/mwf-quikscat/documentation/ mutwqscat.pdf.

Data Processing
The long-time frame of the ERA5 data allowed us to describe changes of the coastal winds over the last four decades. First, hourly time series of U10 and V10 at 75 points ( Figure 1) along the coast of central Chile were used to calculate the alongshore winds (V a ) with respect to the adjacent coastal orientation. These alongshore winds were defined with positive values when directed equatorward. Then hourly values of V a and SLP were daily averaged. To describe changes in CLLJs during the last decades, our approach included changes in windy days (WDs), here defined as those with a mean V a above a threshold of 10 ms −1 . Off central Chile, this value approximately indicates the limit of the upper 20% of the daily frequency distribution of the meridional wind speed [18]. First, we constructed seasonal time series of number of WDs and number of CLLJ events, defined as sets of consecutive WDs; and duration of CLLJ events time-series trends for the period 1979-2019 were calculated through least-square linear fits.
Six-hour SLP data from ERA5 were interpolated onto a 2.5 • × 2.5 • horizontal grid and used as input to a cyclone-tracking scheme, originally developed and implemented at the University of Melbourne [28]. A resolution of 2.5 • allows for the representation of anticyclones since they have an average diameter in the order of 20 • in the southern hemisphere [21]. This tracking algorithm has been and still is widely used for tracking both cyclones and anticyclones (e.g., [21,[29][30][31][32]). Here we located and tracked anticyclones to assess relationships with the occurrence of WDs. From the total number of anticyclones identified, we separated stationary from migratory anticyclones. Following [21], migratory anticyclones were defined as anticyclones that last at least 3 days within the selected domain and move eastward faster than about 3 deg. lon day −1 . Using these data, we calculated the probability that WDs occur along the coast of central Chile when a migratory anticyclone (identified by the tracking algorithm) resides inside the target area shown in Figure 1. Within our target area, stationary anticyclones, particularly the SPA, frequently merge with extratropical migratory anticyclones, producing an SLP spatial pattern that seems to be a poleward extension of the SPA. This pattern is not detected as a migratory anticyclone by the procedure after running the tracking algorithm. However, the relationship between migratory anticyclones and intense southerly winds found in this region using the Climate Forecast System Reanalysis [21] encourage us to use this methodology.
To test the statistical significance of trends we used a Monte Carlo method in which one-thousand synthetic randomized time series were constructed. Finally, if the trend of the original time series was larger than the 90th percentile in the distribution of randomly generated trends, we considered the trend to be significant. To assess any possible seasonality, we split the analyses into austral summer (December-January-February; DJF), fall (March-April-May; MAM), winter (June-July-August; JJA) and spring (September-October-November; SON).

Wind Quality Assessments
To evaluate the quality of the ERA5 coastal V a and SLP, we compared them with observations at coastal meteorological stations and from satellite observations described in Section 2.1. To construct error metrics, common periods of data availability between ERA5 and those for each meteorological station were considered ( Figure 3). The reanalysis data were linearly interpolated to match the location of each meteorological station. The following statistics were used as error metrics to assess the difference between the observed and ERA5 data: bias, root mean squared error (RMSE) and temporal linear correlation (R). This statistical analysis was performed using daily means time series.
Statistics error metrics are presented in Tables 2 and 3 for V a and SLP, respectively. Figure 4 shows scatter plots for selected meteorological stations located at the northern, center and southern regions of our study area.  The comparison between Va data from ERA5 and measurements indicated a high level of agreement for both time series, with R ranging between 0.76 and 0.84 (Table 2).  The comparison between V a data from ERA5 and measurements indicated a high level of agreement for both time series, with R ranging between 0.76 and 0.84 (Table 2). Higher correlations were generally observed during austral winter (mean R-JJA = 0.84) along the entire coast. Lower correlation values appear in austral summer with a mean R-DJF = 0.76 (Table 2).

PEER
Mean values of bias indicated that V a in ERA5 data had, in general, higher values than in situ measurements. Particularly during austral winter and fall, when northerly winds increase in frequency, positive bias represented stronger (negative) winds in observations than reanalysis. In contrast, during spring and summer positive bias was associated with more intense (positive) winds reanalysis. On the other hand, RMSE showed similar values at different seasons, but in general, the values increased southward.
The good performance of ERA5 in representing surface wind variability was also ratified by the comparison of V a with satellite data from QuikSCAT (R = 0.94). A heat scatter plot ( Figure 5) shows that reanalysis and satellite V a data ranged mostly between 3 and 12 ms −1 , the mean ERA5 winds values being only slightly stronger than satellite winds (Table 2). Furthermore, lower values of bias and RMSE were observed in comparison with those obtained from the comparison of reanalysis with meteorological stations. On the other hand, as it can be seen in Table 3, the comparison between measurements and reanalysis of SLP time series indicated a higher degree of agreement than the Va ones, with mean R values between 0.91 and 0.95. Seasonally, the strongest correlations were generally seen during austral winter (mean R-JJA = 0.95), while the opposite was observed during austral summer and fall, both seasons with mean R values of 0.91.
These metrics indicate that reanalysis data and in situ measurements are in close agreement for Va (mean R-annual >0.82) as well as for SLP (mean R annual >0.91), despite the spurious trends observed in reanalysis output derived of changing the source of information over time, among others. Higher correlations were observed during austral winter but they were lower in austral summer. These results highlight the ability of ERA5 for representing the strong wind events in our study area, and, therefore, the southerly low-level coastal jet along the southwestern coast of South America.

Seasonal Trends of the Southerly Coastal Jets
Now we describe the mean seasonal characteristics of the daily-averaged Va off central Chile, with emphasis on the CLLJ events. Figure 6 shows maps of the seasonal mean On the other hand, as it can be seen in Table 3, the comparison between measurements and reanalysis of SLP time series indicated a higher degree of agreement than the V a ones, with mean R values between 0.91 and 0.95. Seasonally, the strongest correlations were generally seen during austral winter (mean R-JJA = 0.95), while the opposite was observed during austral summer and fall, both seasons with mean R values of 0.91.
These metrics indicate that reanalysis data and in situ measurements are in close agreement for V a (mean R-annual >0.82) as well as for SLP (mean R annual >0.91), despite the spurious trends observed in reanalysis output derived of changing the source of information over time, among others. Higher correlations were observed during austral winter but they were lower in austral summer. These results highlight the ability of ERA5 for representing the strong wind events in our study area, and, therefore, the southerly low-level coastal jet along the southwestern coast of South America.

Seasonal Trends of the Southerly Coastal Jets
Now we describe the mean seasonal characteristics of the daily-averaged V a off central Chile, with emphasis on the CLLJ events. Figure 6 shows maps of the seasonal mean wind speeds between 25 • S and 44 • S. Weak mean wind speeds with a mean S-SW direction prevailed during fall and winter north of~35 • S, while the mid-latitude westerlies dominated to the south of it. During spring and summer, wind speeds reached maximum intensities within an area of maximum meridional winds. The wind maxima at 30 • S in spring and 37 • S in summer were anchored to the major capes, Punta Lengua de Vaca and Punta Lavapie, respectively, marking a break in the otherwise straight northsouth coastline [1]. This seasonal wind speed variability was consistent with the SPA poleward displacement during summer, delimiting the latitudinal span of the area where CLLJs develop. Southerly CLLJ events off central Chile occur year-round, but they are more frequent during spring-summer when they are present over 60% of the time [4]. Furthermore, CLLJ events are intermixed with periods of relaxed southerlies or weak northerlies that result in a pronounced day-to-day variability that reaches higher values to the south of the wind maxima in each season (white contours in Figure 6). During fall and spring, the core of this synoptic-scale variability is rooted at Punta Lavapie (~37 • S). During winter maximum V a variability is centered around 35 • S, but during summer, when the core of the wind maxima is located near 37 • S, the maximum of V a variability shifts southward.
OR PEER REVIEW 9 of 19 the core of this synoptic-scale variability is rooted at Punta Lavapie (~37 °S). During winter maximum Va variability is centered around 35 °S, but during summer, when the core of the wind maxima is located near 37 °S, the maximum of Va variability shifts southward. Figure 6e shows the seasonal percent frequency of WDs (Va ≥ 10 ms −1 ) as a function of latitude. These windy days were present more than 20% of the time during winter between 27 °S and 31 °S but decreased rapidly with latitude. During spring this value expanded southward to 37 °S. In contrast, during summer, WDs occurred less than 10% of the time in the northern zone between 25 °S and 30 °S, being more frequent south of 33 °S with a peak frequency of 30% at 37 °S. Overall, these results are in agreement with previous climatologies based on scatterometer data [4], the Climate Forecast System Reanalysis [1] and regional climate simulations [18], lending further support for the use of ERA5 reanalysis.   Figure 6e shows the seasonal percent frequency of WDs (V a ≥ 10 ms −1 ) as a function of latitude. These windy days were present more than 20% of the time during winter between 27 • S and 31 • S but decreased rapidly with latitude. During spring this value expanded southward to 37 • S. In contrast, during summer, WDs occurred less than 10% of the time in the northern zone between 25 • S and 30 • S, being more frequent south of 33 • S with a peak frequency of 30% at 37 • S. Overall, these results are in agreement with previous climatologies based on scatterometer data [4], the Climate Forecast System Reanalysis [1] and regional climate simulations [18], lending further support for the use of ERA5 reanalysis.
During summer, positive trends in WDs were found off south-central Chile, particularly between 37 • S and 44 • S, contrasting with negative values at lower latitudes (29-34 • S) (red line Figure 7a). While a corresponding increase in the number of CLLJ events (consecutive windy days) was also observed between 37 • S and 44 • S, the negative trend at lower latitudes was not clear (red line, Figure 7b) suggesting a shortening of CLLJ event durations between 28 • S and 33 • S (red line, Figure 7c). In contrast, the duration of CLLJ events showed a significant increase only south of 42 • S. During austral winter, significant positive trends in WDs and CLLJ events were obtained in the whole domain, although these trends were higher at the northern part of the domain (gray line, Figure 7a,b), Furthermore, during this season the duration of the CLLJ events exhibited a significant increase south of 37 • S. In fall, a significant increase in the number of WDs was observed between 32 • S and 36 • S (blue line, Figure 7a), but the increase in the number of CLLJ events was observed between 34 • S and 39 • S. This pattern was consistent with the trends found in the duration of CLLJ events, which showed an increase around 33 • S and a decrease around 37 • S. These spatial patterns in changes in CLLJ were, in general, reproduced by other reanalysis products, such as CFSR ( Figure S1).

OR PEER REVIEW
10 of 19 thermore, during this season the duration of the CLLJ events exhibited a significant increase south of 37 °S. In fall, a significant increase in the number of WDs was observed between 32 °S and 36 °S (blue line, Figure 7a), but the increase in the number of CLLJ events was observed between 34 °S and 39 °S. This pattern was consistent with the trends found in the duration of CLLJ events, which showed an increase around 33 °S and a decrease around 37 °S. These spatial patterns in changes in CLLJ were, in general, reproduced by other reanalysis products, such as CFSR ( Figure S1). The time evolution of the number of WDs and CLLJ events are shown at selected locations as time series of seasonal means (austral summer and winter) offshore northcentral Chile at 30 °S and south-central Chile at 38 °S (Figure 8). Off 30 °S positive trends were observed during winter, while negative trends were found in summer. In contrast, positive trends were found at 38 °S, more markedly during summer than winter. Trends that were significant at p = 0.1 are highlighted in Figure 8. The time series show that these long-term trends in the number of windy days and CLLJ events were linear in a first ap-  (Figure 8). Off 30 • S positive trends were observed during winter, while negative trends were found in summer. In contrast, positive trends were found at 38 • S, more markedly during summer than winter. Trends that were significant at p = 0.1 are highlighted in Figure 8. The time series show that these long-term trends in the number of windy days and CLLJ events were linear in a first approximation, and they were embedded within high variability at interannual time scales.

Relationship with Migratory Anticyclones
The average density of all anticyclones and the split between stationary and mig tory anticyclones are shown in Figure 9 for the southern hemisphere. There was a we defined belt of maximum density between 30 °S and 50 °S. Most of the tracked anti clones were stationary, and the density increased with the position of the quasi-stationa subtropical highs at the southeast of each ocean basin. As expected, migratory anti clones, embedded in the westerly wind belt, tracked preferentially to the south of the s tionary anticyclones and exhibited a more homogeneous pattern in longitude. These sults are consistent with a similar analysis using the Climate Forecast System Reanaly [21].
Since strong southerly wind events along the coast of central Chile have been link

Relationship with Migratory Anticyclones
The average density of all anticyclones and the split between stationary and migratory anticyclones are shown in Figure 9 for the southern hemisphere. There was a well-defined belt of maximum density between 30 • S and 50 • S. Most of the tracked anticyclones were stationary, and the density increased with the position of the quasi-stationary subtropical highs at the southeast of each ocean basin. As expected, migratory anticyclones, embedded in the westerly wind belt, tracked preferentially to the south of the stationary anticyclones and exhibited a more homogeneous pattern in longitude. These results are consistent with a similar analysis using the Climate Forecast System Reanalysis [21].
ing fall, increasing to 70% during spring and summer. In this region where a high probability of strong wind in the presence of a migratory anticyclone was found, a positive trend in CLLJ events has been observed during the last four decades (Section 3.2).   Since strong southerly wind events along the coast of central Chile have been linked to the enhanced alongshore pressure gradient associated with the passage of a migratory anticyclone through southern Chile [1,4,6], we evaluated such a link in the ERA5 reanalysis. We calculated the probability that a V a that exceeds 10 ms −1 occurs along the coast of central Chile when the center of a migratory anticyclone, identified by the tracking software, is located inside the box shown in Figure 1. Results showed that when a migratory anticyclone approaches the coast of Chile, there is a high probability (>90%) that a WDs occurs. To refine this relationship, the same analysis was repeated for winds at different latitudes ( Figure 10). As expected, the overall relationship between migratory anticyclones and intense wind events at a particular latitude exhibited lower values. During winter, results showed that when a migratory anticyclone is located within the box there is a probability between 40% and 50% that a strong V a occurs along the coast off central Chile (27-36 • S), this probability being decreasing strongly south of 37 • S. In other seasons, the highest probability was found at 37 • S (Punta Lavapié Cape) with values of 50% during fall, increasing to 70% during spring and summer. In this region where a high probability of strong wind in the presence of a migratory anticyclone was found, a positive trend in CLLJ events has been observed during the last four decades (Section 3.2).
To provide a regional scale picture of changes in migratory anticyclones, Figure 11 shows the difference between the density of the last decade (2000-2019) and the first decade of the data (1979-1988) separately for each season and averaged in longitude in the box shown in Figure 1. An increase in migratory anticyclone density was noticeable during summer around 40 • S, while a slight decrease was observed between 27 • S and 38 • S. Minor changes in density of migratory anticyclones were found during winter without any spatial pattern. Furthermore, a comparison in the distribution of the central SLP in migratory anticyclones between the last decade (2000-2019) and the first decade (1979)(1980)(1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988) of the data is shown in Figure 12. Changes in the central pressures of migratory anticyclones that reach the coast of the Southeast Pacific (box in Figure 1) showed an increase during winter and spring and no major changes during summer and fall. Hence, the observed trends in strong southerly winds and CLLJ events seem to be linked to changes in pressures and frequency of migratory anticyclones at the southwestern coast of South America. The increased coastal winds in the southern part of the study area during summer is then consistent with an increase in migratory anticyclone density, while during winter the observed increase in WDs of intense southerly winds and CLLJ events at the northern region would be more related to an increase in central pressures of migratory anticyclones. These results are consistent with SLP and wind anomalies during the CLLJ events, which show that during the last decade (2000-2019) they were increased compared to the first decade of the data (1979)(1980)(1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988). Figures S2 and S3 show examples of mean anomalies associated with CLLJ centered at Punta Lengua de Vaca (30 • S) and Punta Lavapie, respectively.   To provide a regional scale picture of changes in migratory anticyclones, Figure 11 shows the difference between the density of the last decade (2000-2019) and the first decade of the data (1979)(1980)(1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988) separately for each season and averaged in longitude in the box shown in Figure 1. An increase in migratory anticyclone density was noticeable during summer around 40 °S, while a slight decrease was observed between 27 °S and 38 °S. Minor changes in density of migratory anticyclones were found during winter without any spatial pattern. Furthermore, a comparison in the distribution of the central SLP in migratory anticyclones between the last decade (2000-2019) and the first decade (1979)(1980)(1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988) of the data is shown in Figure 12. Changes in the central pressures of migratory anticyclones that reach the coast of the Southeast Pacific (box in Figure 1) showed an increase during winter and spring and no major changes during summer and fall. Hence, the observed trends in strong southerly winds and CLLJ events seem to be linked to changes in pressures and frequency of migratory anticyclones at the southwestern coast of South America. The increased coastal winds in the southern part of the study area during summer is then consistent with an increase in migratory anticyclone density, while during winter the observed increase in WDs of intense southerly winds and CLLJ events at the northern region would be more related to an increase in central pressures of migratory anticyclones. These results are consistent with SLP and wind anomalies during the CLLJ events, which show that during the last decade (2000-2019) they were increased compared to the first decade of the data (1979)(1980)(1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988). Figures S2 and S3 show examples of mean anomalies associated with CLLJ centered at Punta Lengua de Vaca (30 °S) and Punta Lavapie, respectively.

Discussion
Studies in eastern boundary current systems provide examples of the association between CLLJ and subtropical high-pressure systems, frequently enhanced by the presence of an inland thermal low (e.g., [7]) and that CLLJ events are more frequent during summer [8]. Nevertheless, the most intense southerly wind events observed off central Chile are driven by the alongshore pressure gradient between a coastal low developing in central Chile and a migratory anticyclone drifting eastward farther south [1,5,6]. Therefore, the annual cycle of the strong equatorward winds along the coast is clearly different in the northern and southern part of our study area. Our results are consistent with previous research that have shown the seasonality in the frequency of CLLJ using satellite data [4] and characterized the intense southerly winds at Punta Lengua de Vaca (30 °S) and Punta Lavapié (37 °S) using the Climate Forecast System Reanalysis data [1]. Here, we have shown that strong southerly winds in these two areas have evolved differently in the last four decades. In the northern area the intense alongshore winds have strengthened, particularly in the region (~30 °S) where the maximum winds intensity occurs, and during austral winter. Conversely, the southern area shows no change where the maximum winds occur (~36.5 °S) in summer but shows a decrease in southerly wind events to the north and an increase to the south of 37 °S. These seasonal and spatial patterns of change have also been observed in trends of mean alongshore winds [23,25] and have been mainly attributed to an intensification and poleward expansion of the SPA [22] and to the Bakun's effect in the northern region [25]. Furthermore, this spatial configuration of the mean and intense alongshore wind trends is like that projected for the Humboldt system under global warming scenarios. Regional projections show an increase in mean upwelling-favorable winds off the Chilean coast south of 35 °S [18,33,34]. This is consistent with results from GCMs that show increased mean alongshore winds at the poleward boundaries of the major coastal upwelling systems. There, the region off south-central Chile stands out as one where major and more robust changes in timing, intensity and spatial distribution of mean alongshore winds occur, when compared to other eastern ocean boundaries [35,36]. Similarly, projections for the Humboldt system for austral summer show a clear increase in the intensity of southerly wind events in the southern region and decrease in the northern area [21,37]. A poleward shift of the SPA explains the strengthening (weakening) of mean southerly winds at high (low) latitudes of the Humboldt upwelling system

Discussion
Studies in eastern boundary current systems provide examples of the association between CLLJ and subtropical high-pressure systems, frequently enhanced by the presence of an inland thermal low (e.g., [7]) and that CLLJ events are more frequent during summer [8]. Nevertheless, the most intense southerly wind events observed off central Chile are driven by the alongshore pressure gradient between a coastal low developing in central Chile and a migratory anticyclone drifting eastward farther south [1,5,6]. Therefore, the annual cycle of the strong equatorward winds along the coast is clearly different in the northern and southern part of our study area. Our results are consistent with previous research that have shown the seasonality in the frequency of CLLJ using satellite data [4] and characterized the intense southerly winds at Punta Lengua de Vaca (30 • S) and Punta Lavapié (37 • S) using the Climate Forecast System Reanalysis data [1]. Here, we have shown that strong southerly winds in these two areas have evolved differently in the last four decades. In the northern area the intense alongshore winds have strengthened, particularly in the region (~30 • S) where the maximum winds intensity occurs, and during austral winter. Conversely, the southern area shows no change where the maximum winds occur (~36.5 • S) in summer but shows a decrease in southerly wind events to the north and an increase to the south of 37 • S. These seasonal and spatial patterns of change have also been observed in trends of mean alongshore winds [23,25] and have been mainly attributed to an intensification and poleward expansion of the SPA [22] and to the Bakun's effect in the northern region [25]. Furthermore, this spatial configuration of the mean and intense alongshore wind trends is like that projected for the Humboldt system under global warming scenarios. Regional projections show an increase in mean upwelling-favorable winds off the Chilean coast south of 35 • S [18,33,34]. This is consistent with results from GCMs that show increased mean alongshore winds at the poleward boundaries of the major coastal upwelling systems. There, the region off south-central Chile stands out as one where major and more robust changes in timing, intensity and spatial distribution of mean alongshore winds occur, when compared to other eastern ocean boundaries [35,36]. Similarly, projections for the Humboldt system for austral summer show a clear increase in the intensity of southerly wind events in the southern region and decrease in the northern area [21,37]. A poleward shift of the SPA explains the strengthening (weakening) of mean southerly winds at high (low) latitudes of the Humboldt upwelling system in climate change projections [34,35], whereas a link between intense southerly events and the passage of migratory anticyclones is observed in the present-day climate and projections, as they become more/less frequent at higher/lower latitudes [1,21]. The trends seen during the last four decades for austral summer strong equatorward winds are expected to be continued throughout the XXI century.
Understanding the response of the coastal ocean to CLLJ events has extensive implications for coastal physics, biogeochemistry, ecology and fisheries. For example, it has been broadly recognized that CLLJ produces intense upwelling-favorable wind events that result in coastal cooling [10][11][12]15]. Thus, the increase in frequency in CLLJ reported here could be, at least partially, related to the observed coastal cooling during the last decades in the Humboldt system [22,24,38]. Furthermore, there are intense high-frequency biogeochemical fluctuations in the surface coastal ocean strongly coupled with CLLJ and associated relaxation wind events. During intense upwelling conditions, the dissolved oxygen and pH decrease, while nutrients and air-sea fluxes of heat and CO 2 increase [15,39,40]. The decreasing dissolved oxygen and pH conditions during these intense upwelling events respond mainly to the vertical advection of equatorial subsurface water with dissolved organic matter remineralization products [15,41]. This high-frequency biogeochemical variability can cause changes in the microbiological community and primary production rates, impacting the lower trophic level dynamics [42][43][44]. Hence, the analysis of changes in CLLJ duration is essential to understand the environmental modification along the Humboldt system, which exhibits much more marked and conspicuous spatial and temporal heterogeneity in upwelling intensity and coastal water productivity than previously thought [25].
A long-term study of anticyclone variability through reanalysis products reveals their dependency on the southern annular mode (SAM) in the southern hemisphere. According to [30], strong trends over recent decades in anticyclone frequency have been observed in several reanalysis datasets at the southern hemisphere, and ERA5 in particular shows an increase in migratory anticyclone density ( Figure S4). These changes can be partially attributed to changes in SAM, particularly during summer [30]. This spatial pattern is similar to the differences we found in the density of migratory anticyclones during summer, suggesting that latitudinal changes in migratory anticyclone density (and then in CLLJ) would be related to positive trends of SAM observed during the last decades, particularly during austral summer. This SAM behavior has been attributed to stratospheric ozone depletion and increased greenhouse gas concentrations [45,46]. It is relevant to note that the increase in wind speed in the southern portion of the Humboldt system shows a robust trend in several reanalysis datasets [27]. Particularly, the summer changes in upwellingfavorable wind occurrences provide a key mechanism for changes in the coastal waters of western Patagonia. For example, during the austral summer of 2016, a harmful algae bloom (HAB) developed close to Chiloé Island, generating economic losses and sanitary problems [47]. This HAB coincided with a strong El Niño event and an extreme positive phase of the SAM that altered the atmospheric circulation in southern South America, leading to positive anomalies in upwelling-favorable winds south of 40 • S and higher than normal solar radiation reaching the surface [48,49]. Hence, the positive trends in upwelling-favorable wind events found here in the southern Humboldt system could be contributing to the increase of HAB events in this region.
Understanding the origin of the present-day changes in southerly wind intensity is also especially relevant in the interpretation of paleoclimate records of eolian particles-a direct proxy for wind intensity (e.g., [50]). Along this line, several high-resolution paleoclimate studies have been carried out with the aim of reconstructing changes in the equatorward, alongshore winds in the coastal area of the Humboldt system (e.g., [51][52][53][54][55]). Two of these reconstructions [54,55], based on eolian lithic particles found on marine laminated sediments, emphasize the intensification of these alongshore winds in recent decades [37]. Our results provide new insights on the origin and variability of the alongshore wind intensity contributing to revisiting paleowind reconstructions in southwestern South America.

Conclusions
Based on 41 years of SLP and wind at a height of 10 m above sea level in the ERA5 reanalysis data along the southwestern coast of South America (25-43 • S), we have documented that, embedded within the strong variability at interannual time scales, significant long-term trends exist in the number of days with intense equatorward winds (i.e., > 10 ms −1 ) and in the number and duration of CLLJ events in recent decades. We have placed particular emphasis on a seasonal and latitudinal characterization of trends. An increase in the number of days with intense southerly winds and number of CLLJs events over the whole study area during winter contrasts during summer with a decrease in the number of days with intense southerly winds at lower latitudes (29-34 • S), and an increase farther south, particularly between 37 • S and 44 • S. During fall, an increase in the number of days with intense southerly winds and CLLJs events is observed between 32 • S and 39 • S.
Using an objective algorithm to track anticyclones, we have associated the passage of a migratory anticyclone through the southern portion of the study area with a high probability that an intense wind event occurs along the coast of central Chile on account of a temporary enhancement of the alongshore SLP gradient there, as stated originally in [1,5,6]. Thus, we suggest that changes in extratropical, synoptic-scale migratory anticyclones that reach the coast of South America and force CLLJs events could play an important role in recent trends observed in CLLJ in this region. In particular, the increase in the central surface pressure of migratory anticyclones observed during winter could explain the increase in the strength of CLLJs events. On the other hand, the increase in density of migratory anticyclones observed during summer at about 40 • S could explain the increase in the frequency of CLLJs events in the southern part of the study area. These results highlight the need to include extratropical processes to explain changes in the CLLJs in this coastal upwelling system, and that these changes complement large-scale tropical processes such as the poleward shift of the SPA during the last decades found in previous studies. Data Availability Statement: The ERA5 reanalysis analyzed in this study is available online at https://climate.copernicus.eu/climate-reanalysis.