The Combined Effects of SST and the North Atlantic Subtropical High-Pressure System on the Atlantic Basin Tropical Cyclone Interannual Variability

: The combined effect of the sea surface temperature (SST) and the North Atlantic subtropical high-pressure system (NASH) in the interannual variability of the genesis of tropical cyclones (TCs) and landfalling in the period 1980–2019 is explored in this study. The SST was extracted from the Centennial Time Scale dataset from the National Oceanic and Atmospheric Administration (NOAA), and TC records were obtained from the Atlantic Hurricane Database of the NOAA/National Hurricane Center. The genesis and landfalling regions were objectively clustered for this analysis. Seven regions of TC genesis and ﬁve for landfalling were identiﬁed. Intercluster differences were observed in the monthly frequency distribution and annual variability, both for genesis and landfalling. From the generalized least square multiple regression model, SST and NASH (intensity and position) covariates can explain 22.7% of the variance of the frequency of TC genesis, but it is only statistically signiﬁcant ( p < 0.1) for the NASH center latitude. The SST mostly modulates the frequency of TCs formed near the West African coast, and the NASH latitudinal variation affects those originated in the Lesser Antilles arc. For landfalling, both covariates explain 38.7% of the variance; however, signiﬁcant differences are observed in the comparison between each region. With a statistical signiﬁcance higher than 90%, SST and NASH explain 33.4% of the landfalling variability in the archipelago of the Bahamas and central–eastern region of Cuba. Besides, landfalls in the Gulf of Mexico and Central America seem to be modulated by SST. It was also found there was no statistically signiﬁcant relationship between the frequency of genesis and landfalling with the NASH intensity. However, the NASH structure modulates the probability density of the TCs trajectory that make landfall once or several times in their lifetime. Thus, the NASH variability throughout a hurricane season affects the TCs trajectory in the North Atlantic basin. Moreover, we found that the landfalling frequency of TCs formed near the West Africa coast and the central North Atlantic is relatively low. Furthermore, the SST and NASH longitude center explains 31.6% ( p < 0.05) of the variance of the landfalling intensity in the archipelago of the Bahamas, while the SST explains 26.4% ( p < 0.05) in Central America. Furthermore, the 5-year moving average ﬁlter revealed decadal and multidecadal variability in both genesis and landfalling by region. Our ﬁndings conﬁrm the complexity of the atmospheric processes involved in the TC genesis and landfalling.


Introduction
Tropical cyclones (TCs) can be defined as non-frontal synoptic low-pressure systems, which are formed over tropical or subtropical waters at a synoptic scale, with a warm core, organized deep convection, and a cyclonic circulation of wind on the surface. These systems mainly form from June to November (hereafter "TC season" or "hurricane season"). Besides, they are one of the most destructive natural phenomena in the world [1]. Therefore, any trend in TC activity is relevant to society, and particularly for coastal populations where the effects of TCs are most damaging. Historically, but particularly in recent years, the devastating impacts of TCs formed in the North Atlantic (NATL) basin have been a major concern for the economies and societies of affected countries [2]. In 2001, Goldenberg et al. [3] showed that the frequency of TCs did not exhibit a clear trend from 1971 to 2000, while TC intensity shows an increasing trend. Later, Trenberth [4] proposed that TC activity in the NATL basin had increased significantly by the end of the 20th century. Observations of TCs are considered to be more accurate since the introduction of meteorological satellites in the late 1970s, and since then, there has been no clear trend in the frequency of TCs on a global scale [5][6][7]. Nevertheless, the low quality of historical TC records prevents an accurate analysis of trends in intensity [8,9] and frequency [10][11][12]. Consequently, there are doubts as to whether changes in tropical cyclonic activity have exceeded the expected natural variability [13].
Several authors [14][15][16] have shown that enthalpy fluxes play an important role in the TCs evolution and are strongly modulated by sea surface temperature (SST). Surface fluxes of latent and sensible heat from the oceans provide the potential energy to TCs [14]. Rapid changes in SST are directly linked to changes in TCs intensity and trajectory [17]. Indeed, it has been related to the poleward migration in the average latitude at which TCs reach maximum intensity [7]. Furthermore, the changes in the cloud pattern and the rainfall associated with TCs are influenced by the SST through the atmospheric water vapor [13,18]. These relationships contribute to the air-sea enthalpy disequilibrium and therefore influence the TCs intensification and weakening process. The TCs over warm waters tend to intensify [19,20], while dissipating or weakening over cold waters [21].
The close relationship between TCs and SST has opened a debate on the possibility of an increase in the total number of TCs that form in a warmer ocean [22,23]. By 2006, Santer et al. [24] showed that SST in most regions of TC formation increased by several tenths of a degree Celsius during the past decades. Moreover, many scientists have found that increased SST [24,25] affects not only the hurricane frequency [26,27] but also hurricane intensity [28][29][30][31] and TC size [32]. Likewise, some studies have proposed that the increase in SST due to global warming has caused an increase in destructive TCs in recent decades [30,31,33]. Vecchi and Knutson [12] also showed that multidecadal variations in TC counts in the NATL basin appear to be strongly correlated with SST; however, the upward trend in NATL TCs counts since 1851 contains a spurious trend component due to missing storms before the satellite era [10][11][12].
One of the most important aspects to take into account in the analysis of TC activity is the landfalling events. When a TC makes landfall, it can produce high winds, storm surges, and severe flooding that can result in loss of life and widespread damage to coastal populations. The intensity of the North Atlantic subtropical high-pressure system (NASH) can also influence the TC activity in the NATL [34,35]. A westward extension and a strengthening of the NASH favors the non-curving moving tracks of TCs, while a weakening of the NASH leads to a recurving trajectory over the Atlantic Ocean, decreasing landfalling frequency [36]. Fudeyasu et al. [37] determined that most of the TCs landfalling on the United States (U.S.) form in the Caribbean Sea and the Gulf of Mexico and move north-westward toward the southern and eastern coasts. In agreement, previous findings of Kossin et al. [35] revealed that most landfalling storms form in the Gulf of Mexico, and have not experienced discernible trends.
Dailey et al. [38] examined the relationship between NATL SST and U.S. landfalling hurricanes, suggesting that landfall risk is sensitive to SST conditions, but that sensitivity Atmosphere 2021, 12, 329 3 of 18 varies by region and intensity. Some authors have even used a Poisson regression model to examine the relationship between U.S. landfalling hurricane events and principal climate variability modes, e.g., [38,39].
Despite the great number and variety of studies, most have evaluated individually the influence of the SST, wind shear, or teleconnection modes on the cyclonic activity in the NATL, while less attention has been given to landfall frequency and geographic distribution. However, this, beside the changes in their intensity, and frequency, are crucial factors to assess the human, environmental, and economic impacts of TCs. Indeed, a key question that has motivated the development of many research studies on this topic, e.g., [2,6,9,23,24,32], is whether the activity of TCs has changed and will change due to global warming. Additionally, previous works have shown a strengthening of the NASH in a warmer climate, but the cause of its long-term variability, as well as the impact of this variability on the TC activity, has been poorly investigated [36,40]. Thus, the goal of this study is to examine the influence of SST and NASH intensity and position on the interannual variability observed in TCs genesis and landfalling frequency from 1980 to 2018 in the whole NATL basin. To investigate further, the aim is also to consider different regions of the genesis of TCs and landfall.
This paper is organized into four sections. In Section 2 the hurricane data, the NASH, and the SST dataset are described. Additionally, it contains the methodology used, including the description of the generalized least square model and the Kernel density estimation. Section (Section 3.2). Section 4 contains conclusions.

Dataset
Historical records of TCs in the NATL basin from 1851 to 2019 were taken from the Atlantic hurricane database (HURDAT2) of the NOAA/National Hurricane Center. This dataset comprised a text format file of TC parameters, including location and intensity (maximum wind speed and minimum central pressure), taken every six hours [41][42][43]. The HURDAT2 database also includes the TCs landfall events.
The possibility of detecting TCs that did not approach land or ships during their trajectories was minimal before the beginning of the meteorological satellites era in the mid-1970s. Therefore, historical records in the HURDAT2 database may contain physically unsupported trends due to the uncertain storm detection at this time (e.g., [8][9][10][11][12]). Kossin et al. [8] point out that the intensity estimation methods based on satellite were improved with the introduction of the Dvorak technique [44,45].
Considering these limitations of TC records in the NATL basin, we followed the criteria of Kossin et al. [8] and Bathia et al. [46] for the selection of the study period. Several authors (e.g., [8,[46][47][48]) pointed out that after 1980 the TCs records in the NATL present higher quality for climate analysis. Thus, we focus the following analysis on the period after 1980.
The June-November monthly mean SST data was obtained from the Centennial Time Scale (COBE SST2) dataset [49,50] of the National Oceanic and Atmospheric Administration (NOAA). Daily SST of this dataset was constructed as the sum of a trend, interannual variations, and daily changes, using in situ SST and sea ice concentration observations. The annual series of SST anomalies shown in Figure 1a was computed for the region confined between 5 • N-50 • N in latitude and 10 • W-100 • W in longitude. This region encompasses the genesis and peak of maximum intensity portions of all TCs formed from 1980 to 2019 in this area, as shown in Figure 1b.
To obtain a coherent physical insight of TCs activity in the NATL basin during the study period, we used the June-November monthly averages of the mean sea level pressure (MSLP) from the National Center for Environmental Prediction-National Center for Atmospheric Research (NCEP-NCAR) reanalysis fields on a 2.5 • latitude by 2.5 • longitude grid [51]. To obtain a coherent physical insight of TCs activity in the NATL basin dur study period, we used the June-November monthly averages of the mean sea leve sure (MSLP) from the National Center for Environmental Prediction-National Cen Atmospheric Research (NCEP-NCAR) reanalysis fields on a 2.5° latitude by 2.5° tude grid [51].

Methodology
To find inconsistencies in the HURDAT2 database records, we compared the tasets with the International Best Track Archive for Climate Stewardship da (IBTrACS) [52]. The IBTrCAS database collects the historical data of the best mon of tropical cyclones from all Regional Specialized Meteorological Centers and othe cies, and combines the different datasets in a single product. In our study period databases show similar records. Therefore, the analysis presented in this research is on the HURDAT2 database.
The historical records of TC genesis (defined as the first reported location of t and the landfalling events are separated into clusters applying statistical methods, permitted to investigate the whole basin and each cluster separately. The cluster tec applied is similar to that described by Corporal-Lodangco et al. [53].
The center of the NASH was determined based on the nearest neighbor tec [54][55][56][57][58][59] to find the highest pressure value using the MSLP in the area between 15° N and 60° W-0°. Each grid point was compared to the 24 neighboring nodes. A grid is a candidate to be the core of the NASH if its MSLP is higher than or equal to th

Methodology
To find inconsistencies in the HURDAT2 database records, we compared these datasets with the International Best Track Archive for Climate Stewardship database (IBTrACS) [52]. The IBTrCAS database collects the historical data of the best monitoring of tropical cyclones from all Regional Specialized Meteorological Centers and other agencies, and combines the different datasets in a single product. In our study period, both databases show similar records. Therefore, the analysis presented in this research is based on the HURDAT2 database.
The historical records of TC genesis (defined as the first reported location of the TC) and the landfalling events are separated into clusters applying statistical methods, which permitted to investigate the whole basin and each cluster separately. The cluster technique applied is similar to that described by Corporal-Lodangco et al. [53].
The center of the NASH was determined based on the nearest neighbor tech nique [54][55][56][57][58][59] to find the highest pressure value using the MSLP in the area between 15 • N-60 • N and 60 • W-0 • . Each grid point was compared to the 24 neighboring nodes. A grid point is a candidate to be the core of the NASH if its MSLP is higher than or equal to the eight near neighbors and higher than the other 16 nodes surrounding the eight neighbors around the single center grid point. In the case that several grid points have the same maximum pressure value, the grid point with the lowest latitude was considered as NASH center [54]. The Pearson's correlation coefficient (hereafter r) was determined by applying a twotailed test, using the effective number of degrees of freedom [60]. Furthermore, a t-test [61] was applied to determine the statistical significance of the correlations.
To understand the combined effect of the SST and NASH intensity and position on the interannual variability of cyclonic activity in the NATL basin, we applied the generalized least squares (GLS) estimation model. Similar to other multiple linear regression models, GLS allows us to estimate the unknown parameters of the linear relationship according to Equation (1): where Y is the dependent variable, a is a vector containing the regression coefficients, ] is a matrix (n_data x n_features) of the independent variables, and b is the intercept. Further information of the GLS model is provided by Safi and Saif [62] and Davidson and MacKinnon [63].

Kernel Density Estimation
To estimate the probability density of the TCs trajectory and intensity at landfall, the nonparametric kernel density estimation (KDE) was applied, following Equation (2): where x is the independent variable of the kernel function K(x) and x i are the data points; f(x) is the smooth estimate; h is the selected bandwidth; and n is the length of the dataset. The KDE smooths each data point x i into a small density bumps and then sums all these small bumps together to obtain the final density estimate. The kernel function K(x) must commonly satisfy the condition K(x)dx = 1, but may be considered arbitrary. The Gaussian kernel (Equation (3)) is a common choice [64]. Here, we used the classical KDE [65].

TC Genesis
From the cluster technique, seven regions, in which the genesis of TCs in the NATL frequently occurs, were identified. These are the tropical Central North Atlantic (G 1 , purple) that represents 14.9% of the total genesis counts in the study period, the Caribbean Sea (G 2 , red, 11.7%), the Lesser Antilles arc (G 3 , green, 13.1%), the coast of West Africa (G 4 , orange, 18.2%), the Central North Atlantic (G 5 , yellow, 6.2%), the northeast of the Bahamian archipelago (G 6 , brown, 21.3%), and the Gulf of Mexico (G 7 , blue, 14.6%), as shown in Figure 2 (top). The monthly seasonality of TCs formed on each cluster is represented in Figure 2 (bottom). Genesis at the beginning of the hurricane season generally occurs at high latitudes due to favorable thermodynamic conditions and in the Gulf of Mexico, while at the end of the season they are more frequent in the western Caribbean Sea and the central North Atlantic. The regions G 1 and G 4 exhibit the highest frequency of genesis in August and September, coinciding with the peak of the cyclonic season. These G 1,4 TCs are frequently originated from African easterly waves.
The interannual variability of the cyclonic activity is deduced from Figure 3; however, the application of the 5-year moving average filter reveals the existence of a decadal time scale variability, more frequent in the TCs formed in G 5 y G 7 regions.
The relationship between the mean annual SST during the TC season in the NATL basin and the number of genesis was first analyzed using Pearson's correlation. A strong statistically significant correlation was observed at 95% significance level between SST and TC genesis throughout the basin, as well as in the G 1 and G 4 regions, as shown in Table 1. In G 2 , G 3 , and G 5 , the correlation coefficient ranged between 0.13 and 0.30, however, it was not statistically significant. The inverse correlation (p > 0.05) observed in the G 6

and G 7 clusters
Atmosphere 2021, 12, 329 6 of 18 is curious, and constitutes a strong evidence of the complexity of the physical processes for TC genesis. In these two regions, the other preconditions described by Gray [66] seem to be more important. The NASH latitude center is significantly correlated with TC genesis frequency (r = 0.36) in the whole basin, while G 3 (r = 0.36) is significantly correlated, too. No significant correlations were found between TC genesis and NASH intensity and NASH longitude center (Table 1) The interannual variability of the cyclonic activity is deduced from F ever, the application of the 5-year moving average filter reveals the existen time scale variability, more frequent in the TCs formed in G5 y G7 regions.  Table 1. Pearson's correlation coefficients between annual time series of mean SST, the North Atlantic subtropical high-pressure system (NASH) intensity and NASH position (June to November) with the number of tropical cyclones (TCs) genesis in the whole NATL and each cluster. Statistical significance is marked with the bold text (p < 0.05). NASH latitude and NASH longitude represent the NASH center latitude and NASH center longitude, respectively.  The interannual variability of the cyclonic activity is deduced from Figure 3; however, the application of the 5-year moving average filter reveals the existence of a decadal time scale variability, more frequent in the TCs formed in G5 y G7 regions.  When the frequency of genesis is modeled with the SST as a single covariate from the GLS linear regression model, both for the whole basin and for the G N regions (N = 1, 2, 3, 4, 5, 7), the genesis frequency increases with increasing SST, while an inverse pattern is observed in G 6 , as shown in Figure 4. However, SST only explains 12% (p < 0.05) of the genesis frequency variance throughout the basin and 9%, 15.1%, and 9.1% in G 1 , G 4 , and G 5 regions, respectively. The northward latitudinal movement of the NASH center favors the genesis processes, explaining 12.9% (p < 0.05) of the variance in the entire basin. By regions, the latitude center of the NASH explains only the variability of the G 3 region, with 13.1% (p < 0.05). It is notable that the intensity of NASH individually does not play an important role in the genesis of TCs. These results fully support the linear relationship obtained from Pearson's correlation coefficient.
As a logical next step, we consider how SST and NASH (intensity and position) modulate the frequency of genesis both throughout the basin and within each cluster (the four covariates do not contain significant linear trends over the study period). The results of the regression are shown in Table 2. The GLS multiple linear regression suggests that the SST and the NASH (intensity and position) covariates explain 22.7% of the genesis frequency variance over the whole basin; however, the relationship is only statistically significant (p < 0.05) for NASH center latitude. The GLS regression reveals that the G 4 region is more dependent on the mean SST, while G 3 is dependent on the latitudinal variation of the NASH center. As a logical next step, we consider how SST and NASH (intensity and position) modulate the frequency of genesis both throughout the basin and within each cluster (the four covariates do not contain significant linear trends over the study period). The results of the regression are shown in Table 2. The GLS multiple linear regression suggests that the SST and the NASH (intensity and position) covariates explain 22.7% of the genesis frequency variance over the whole basin; however, the relationship is only statistically significant (p < 0.05) for NASH center latitude. The GLS regression reveals that the G4 region is more dependent on the mean SST, while G3 is dependent on the latitudinal variation of the NASH center. The landfalling monthly variation shows the maximum in August and September, although in regions L 2 and L 3 , secondary maximums are observed in June and July. The secondary maximums are a consequence of the TCs genesis in the Gulf of Mexico at the beginning of the TC season, while in the region L 4 , a secondary maximum is observed in October, mainly because of the TCs that originate in the Western Caribbean Sea at the end of the hurricane season. For the entire NATL basin, the highest frequency of landfalling is observed in August and September. The landfalling frequency shows a similar monthly distribution to genesis frequency. Table 2. Generalized least squares (GLS) multiple linear regression of annual NATL genesis frequency (for the whole basin and for each cluster) onto the SST and NASH (intensity and position). The SST and NASH values are based on June-November. Statistical significance is marked with bold text (p < 0.05). SST is in • C and NASH intensity in hPa. NASH latitude and NASH longitude represent the NASH center latitude and NASH center longitude, respectively.  The landfalling monthly variation shows the maximum in August and September, although in regions L2 and L3, secondary maximums are observed in June and July. The secondary maximums are a consequence of the TCs genesis in the Gulf of Mexico at the beginning of the TC season, while in the region L4, a secondary maximum is observed in October, mainly because of the TCs that originate in the Western Caribbean Sea at the end of the hurricane season. For the entire NATL basin, the highest frequency of landfalling is observed in August and September. The landfalling frequency shows a similar monthly As appreciated in Figure 6, the annual evolution of landfalling events experienced a major frequency during the 2000 decade for all the clusters. A visual analysis of the time series filtered with a centered 5-year moving window reveals that great similarities exist between the landfall occurrence in L 1 and L 3 , L 2 and L 4 , and are greatly different in L 5 . The increase of annual landfalling in the latter part of the record is also evident, particularly in L 1 and L 2 . TC landfalling seems to have a direct relationship with the increase of SST. The son correlation between the mean annual SST (June-November) and the NATL wide (L1 +…+ L5) landfalling count of 0.52 is statistically significant at 95% of the s cance level, as well as in regions L1, L3, and L4 with r = 0.44, r = 0.46, and r = 0.33, r tively. However, it is notable that the Pearson correlation between SST and landfal low in L2 and L5 (r = 0.30), but statistically significant at 90% of the significance lev the whole basin, SST and landfalling counts exhibit a strong correlation (r = 0.52), s cally significant at 95% of the significance level. No statistically significant correlatio found between NASH intensity and TC landfalling counts. However, the correlati tween the frequency of landfalling and the latitude NASH center in a whole basin tistically significant (r = 0.42), as shown in Table 3. A significant correlation is al served in L1 (r = 0.34) and L4 (r = 0.38) landfalling regions with the latitudinal variat the NASH center. Moreover, L1 exhibits a strong Pearson correlation (0.40, p < 0.05 the longitude of the NASH center. Table 3. Pearson's correlation coefficients between annual time series of mean SST, the NAS intensity, and NASH position (center latitude and longitude) (June to November) with the l falling counts in the whole NATL (L1 + ... + L5) and each cluster. Statistical significance is ma with bold text (p < 0.05). NASH latitude and NASH longitude represent the NASH center la and NASH center longitude, respectively.  TC landfalling seems to have a direct relationship with the increase of SST. The Pearson correlation between the mean annual SST (June-November) and the NATL basin-wide (L 1 + . . . + L 5 ) landfalling count of 0.52 is statistically significant at 95% of the significance level, as well as in regions L 1 , L 3 , and L 4 with r = 0.44, r = 0.46, and r = 0.33, respectively. However, it is notable that the Pearson correlation between SST and landfalling is low in L 2 and L 5 (r = 0.30), but statistically significant at 90% of the significance level. For the whole basin, SST and landfalling counts exhibit a strong correlation (r = 0.52), statistically significant at 95% of the significance level. No statistically significant correlation was found between NASH intensity and TC landfalling counts. However, the correlation between the frequency of landfalling and the latitude NASH center in a whole basin is statistically significant (r = 0.42), as shown in Table 3. A significant correlation is also observed in L 1 (r = 0.34) and L 4 (r = 0.38) landfalling regions with the latitudinal variations of the NASH center. Moreover, L 1 exhibits a strong Pearson correlation (0.40, p < 0.05) with the longitude of the NASH center. From the GLS model, the SST and the NASH latitude covariates explain 38.7% of the variance of the landfalling events over the NATL basin with a statistically significant dependence (p < 0.05) on the SST. Similar results are obtained for cluster L 1 , but explaining only 33.4% of the variance. The SST mostly explains the variance (p < 0.05) of the L 3 and L 4 regions with 26.4% and 24.1%, respectively. In the case of the L 2 and L 5 regions, the landfalling events do not appear to have a significant modulation of the SST and the position of the NASH. In the same way, as for genesis, no dependence was observed between the landfalling events and the intensity of the NASH, as shown in Table 4.

All Clusters
The landfalling of the TCs formed in G 1 are distributed among all the L N regions (N = 1, 2, 3, 4, 5), although the highest frequency is observed in L 5 . Similarly, G 3 and G 4 also exhibit high frequencies of landfalling in L 5 , as shown in Figure 7. We assume that these distributions are mainly due to the wide range of trajectories for these TCs. Additionally, G 3 has the highest frequency of landfalling in the central−eastern region of Cuba, the archipelago of the Bahamas, and the west of La Española (L 1 ), and another secondary maximum in L 2 . TCs formed in the western Caribbean Sea frequently make landfall on the coasts of Central America, the Yucatan Peninsula, and the western region of Cuba (L 4 ). In L 2 , TCs that form in G 4 and G 6 frequently make landfall, while those originating in the Gulf of Mexico (G 7 ) generally struck the coastline in L 3 . Table 3. Pearson's correlation coefficients between annual time series of mean SST, the NASH intensity, and NASH position (center latitude and longitude) (June to November) with the landfalling counts in the whole NATL (L 1 + ... + L 5 ) and each cluster. Statistical significance is marked with bold text (p < 0.05). NASH latitude and NASH longitude represent the NASH center latitude and NASH center longitude, respectively.  Table 4. Generalized least squares (GLS) multiple linear regression of annual NATL landfalling frequency (for the whole basin and each cluster) onto the SST and NASH (intensity and position). The SST and NASH values are based on June-November. Statistical significance is marked with the bold text (p < 0.05). SST is in • C and NASH intensity in hPa. NASH latitude and NASH longitude represent the NASH center latitude and NASH center longitude, respectively.  The landfalling events mostly depend on the TC trajectory, which is influenc the large-scale steering flow. Changes in the NASH position and intensity induce scale flow changes that affect the TCs track and, consequently, the probability of lan ing. Figure 8 shows the TCs track density of TCs that made landfall in each region a mean position of the NASH (composite using the specific month and specific year o falling event). As can be seen, there is no significant variation in the position of N center, however, there are appreciable differences in the NASH structure, which mo the preferential trajectory followed by the TCs. For L1 and L2, the anticylonic ridge pronounced, while for L3 and L4, it penetrates to the Florida Peninsula. For the L5 r the structure of the NASH shows a pattern similar to that observed for L1 and L2, alt more contracted, suggesting that TCs that make landfall in L5 have a high probabi make landfall again in L1 or L2. The NASH structure directly influences the distribu the probability density tracks of TC that make landfall. The effect is different depe on the genesis location. TCs that develop in the eastern part of the main develop region (G4) tend to have a low probability of landfalling than those that develop western part (G1 and G3), due to the longer time over the tropical ocean. These result support the findings of Colbert and Soden [36], who pointed out that the evolut NASH throughout a hurricane season affects the TCs trajectory in the NATL basin. the NASH is strong and centered to the west, TCs are directed south toward the G Mexico, increasing the landfalling probability. The landfalling events mostly depend on the TC trajectory, which is influenced by the large-scale steering flow. Changes in the NASH position and intensity induce large-scale flow changes that affect the TCs track and, consequently, the probability of landfalling. Figure 8 shows the TCs track density of TCs that made landfall in each region and the mean position of the NASH (composite using the specific month and specific year of landfalling event). As can be seen, there is no significant variation in the position of NASH center, however, there are appreciable differences in the NASH structure, which modulate the preferential trajectory followed by the TCs. For L 1 and L 2 , the anticylonic ridge is less pronounced, while for L 3 and L 4 , it penetrates to the Florida Peninsula. For the L 5 region, the structure of the NASH shows a pattern similar to that observed for L 1 and L 2 , although more contracted, suggesting that TCs that make landfall in L 5 have a high probability to make landfall again in L 1 or L 2 . The NASH structure directly influences the distribution of the probability density tracks of TC that make landfall. The effect is different depending on the genesis location. TCs that develop in the eastern part of the main development region (G 4 ) tend to have a low probability of landfalling than those that develop in the western part (G 1 and G 3 ), due to the longer time over the tropical ocean. These results fully support the findings of Colbert and Soden [36], who pointed out that the evolution of NASH throughout a hurricane season affects the TCs trajectory in the NATL basin. When the NASH is strong and centered to the west, TCs are directed south toward the Gulf of Mexico, increasing the landfalling probability.

All
These results are confirmed in Figure 9, which shows the position and mean structure of the NASH (considering only the months where landfalling occurred) for the year with the highest frequency of landfalling, both for the entire basin and each region. The steering flow imposed by the NASH circulation modulates the frequency of landfalling in the different regions. No significant variations are observed in the intensity of the NASH; however, the frequency of landfalling depends on the latitudinal variations of the NASH center. When the NASH center is shifted to the north, it favors that TCs make landfall in the L 1 , L 2 , and L 5 regions, while an NASH center low in latitude modulates the movement of the TCs through the Caribbean Sea towards the Gulf of Mexico, favoring landfalling in Central America, the western region of Cuba, and the Mexican coast. These results are confirmed in Figure 9, which shows the position and mean structure of the NASH (considering only the months where landfalling occurred) for the year with the highest frequency of landfalling, both for the entire basin and each region. The steering flow imposed by the NASH circulation modulates the frequency of landfalling in the different regions. No significant variations are observed in the intensity of the NASH; however, the frequency of landfalling depends on the latitudinal variations of the NASH center. When the NASH center is shifted to the north, it favors that TCs make landfall in the L1, L2, and L5 regions, while an NASH center low in latitude modulates the movement of the TCs through the Caribbean Sea towards the Gulf of Mexico, favoring landfalling in Central America, the western region of Cuba, and the Mexican coast. From Figure 9, we can infer that landfalling occurs more frequently in regions where positive SST anomalies predominate in the surrounding waters. However, in L 2 in 2004 (Table 5), there were 10 landfallings with negative SST anomalies in the northeast of the Gulf of Mexico, the Straits of Florida, and the seas surrounding the archipelago of the Bahamas, which is evidence of the strong modulating role of steering flow imposed by the NASH. As shown in Figure 10, the probability density of TC landfalling intensity shows the most probable intensity ranging between 18 and 35 m/s. Despite the similarity of the patterns, the regions slightly differ from another in the probability density of the TC intensity spectrum at landfall. In all cases, the mean landfall intensity is approximately 35 m/s, although region L1 shows a mean intensity of 42.5 m/s, with two maximums of the probability density. Figure 9. Composite of NASH position and structure (contour) for months with landfalling events in the year of maximum landfalling frequency for the whole basin and each region. The shaded plots represent the mean SST anomalies for the same months and years of NASH position and structure composite. The black star point represents the NASH center, and cyan points represent landfalling events. "All clusters" represents the whole basin and LN (N = 1, 2, 3, 4, 5) represents each landfalling cluster.
From Figure 9, we can infer that landfalling occurs more frequently in regions where positive SST anomalies predominate in the surrounding waters. However, in L2 in 2004 (Table 5), there were 10 landfallings with negative SST anomalies in the northeast of the Gulf of Mexico, the Straits of Florida, and the seas surrounding the archipelago of the Bahamas, which is evidence of the strong modulating role of steering flow imposed by the NASH. Table 5. TC season of maximum landfalling frequency for whole basin and each region. L1  L2  L3  L4  L5  Year  2008  2008  2004  2005  2005  2017  Landfall counts  29  10  10  10  8  9 As shown in Figure 10, the probability density of TC landfalling intensity shows the most probable intensity ranging between 18 and 35 m/s. Despite the similarity of the patterns, the regions slightly differ from another in the probability density of the TC intensity Figure 9. Composite of NASH position and structure (contour) for months with landfalling events in the year of maximum landfalling frequency for the whole basin and each region. The shaded plots represent the mean SST anomalies for the same months and years of NASH position and structure composite. The black star point represents the NASH center, and cyan points represent landfalling events. "All clusters" represents the whole basin and L N (N = 1, 2, 3, 4, 5) represents each landfalling cluster.  Moreover, the SST and NASH (intensity and position) covariates explain 31.6% of the variance of the intensity at which the TCs make landfall in the L1 region; however, the dependence is only significant (p < 0.05) for the SST and the NASH center longitude. In the case of region L4, the covariates explain 25.1%, but it is statistically significant for the Moreover, the SST and NASH (intensity and position) covariates explain 31.6% of the variance of the intensity at which the TCs make landfall in the L 1 region; however, the dependence is only significant (p < 0.05) for the SST and the NASH center longitude. In the case of region L 4 , the covariates explain 25.1%, but it is statistically significant for the NASH center latitude. No statistically significant relationships were found for the other regions.

Conclusions
The genesis of tropical cyclones (TCs) in the North Atlantic basin from 1980 to 2019 were objectively separated into groups and analyzed to identify the combined effects of sea surface temperature (SST) over the domain previously defined in Figure 1, and the North Atlantic subtropical high-pressure system (NASH) in its interannual variability. A strong Pearson correlation was identified between the frequency of genesis throughout the whole basin and the SST, as well as between the genesis frequency and NASH position. Nevertheless, the generalized least squares (GLS) multiple linear regression model showed that the SST and NASH covariates can only explain 22.7% of the variance in the genesis frequency, but the relationship is statistically significant with the NASH latitude center (p < 0.1). Furthermore, the multiple regression revealed that TC genesis near the coast of West Africa is more dependent on SST, while TCs formed near the Lesser Antilles arc are more influenced by the NASH position than by the SST. In the other genesis clusters, no statistically significant relationship was identified between SST and NASH. These findings show the complexity of the atmospheric processes involved in the TC genesis. Moreover, our results also reveal that the NATL TC genesis exhibits clear interbasin differences in frequency, but similar interannual variability.
Differences between each landfalling region, as well as the interannual variability of landfalling events, were also identified by Pearson's correlation and multiple regression analysis. The frequency of landfalling in the Florida Peninsula and the east coast of the United States show a marked decadal variability, while the other regions exhibit a multidecadal variability. Similar to the frequency of genesis, the frequency of landfalling events throughout the basin shows a strong Pearson correlation with the SST and the NASH latitude center. However, the SST and NASH covariates only explain 38.7% (p < 0.05 for SST) of the variance in the study period. Regionally, both covariates explain 33.4% of the variance of the frequency of landfalling in the archipelago of the Bahamas and the central−eastern region of Cuba, while landfalling events in the coasts of the Gulf of Mexico and Central America seem to have a statistically significant dependence on SST. Similar to the frequency of genesis, no direct modulation of the NASH intensity and the landfalling frequency was found. However, the structure of NASH modulates the probability density of TC trajectories that made landfall one or more times in their lifetime. The TCs formed near the West Africa coasts and central North Atlantic show the lowest frequency of landfalls in the study period. The latitudinal shift of the NASH modulates the frequency of landfalling in each region. Furthermore, the SST explains only 24.4% of the variance of the intensity at which TCs make landfall in the Bahamas archipelago, while the NASH position explains 24.6% of those who make landfall in Central America and the western region of Cuba.
The results discussed above support the observation that it is not optimal to consider the genesis and landfalling regions as a whole when attempting to study the interannual variability of TC activity. The differences observed in both genesis and landfalling regions show that the TC activity is also dependent on local-scale factors and processes such as the SST, but also of the synoptic-scale influence of the NASH, which is an important feature to consider when making future projections of TC activity.