The average latitude where tropical cyclones (TCs) reach their peak intensity has been observed to be shifting poleward in some regions over the past 30 years, apparently in concert with the independently observed expansion of the tropical belt. This poleward migration is particularly well observed and robust in the western North Pacific Ocean (WNP). Such a migration is expected to cause systematic changes, both increases and decreases, in regional hazard exposure and risk, particularly if it persists through the present century. Here, it is shown that the past poleward migration in the WNP has coincided with decreased TC exposure in the region of the Philippine and South China Seas, including the Marianas, the Philippines, Vietnam, and southern China, and increased exposure in the region of the East China Sea, including Japan and its Ryukyu Islands, the Korea Peninsula, and parts of eastern China. Additionally, it is shown that projections of WNP TCs simulated by, and downscaled from, an ensemble of numerical models from phase 5 of the Coupled Model Intercomparison Project (CMIP5) demonstrate a continuing poleward migration into the present century following the emissions projections of the representative concentration pathway 8.5 (RCP8.5). The projected migration causes a shift in regional TC exposure that is very similar in pattern and relative amplitude to the past observed shift. In terms of regional differences in vulnerability and resilience based on past TC exposure, the potential ramifications of these future changes are significant. Questions of attribution for the changes are discussed in terms of tropical belt expansion and Pacific decadal sea surface temperature variability.
1. Introduction and motivation
Over the past few decades, the average location where tropical cyclones (TCs) reach their peak lifetime intensity has systematically migrated poleward in some regions (Kossin et al. 2014). This has occurred concurrently with systematic changes in the mean tropical and subtropical environments (Kossin et al. 2014; Kang and Elsner 2016). Statistically significant poleward migration rates in the mean latitude where TCs reach their lifetime-maximum intensity (LMI; denoted herein as ϕLMI) are observed in the western North Pacific Ocean (WNP), southern Indian Ocean, and South Pacific Ocean, with statistically insignificant rates found in the eastern North Pacific and northern Indian Oceans and a statistically significant equatorward migration rate in the North Atlantic (Fig. 1). In a hemispheric sense, then, TCs in both the southern Indian and South Pacific Oceans apparently contribute to the observed poleward migration in the Southern Hemisphere, while WNP TCs contribute largely to the Northern Hemisphere trends identified by Kossin et al. (2014). The WNP is the most active ocean basin in terms of TC activity and also encompasses most regions with the highest human exposure and mortality risk from TCs (Peduzzi et al. 2012). In particular, there is a large amount of coastline poleward of the tropics in the WNP that is exposed to TC-related hazard. Contrarily, while the southern Indian and South Pacific Oceans experience a relatively high level of TC activity, there is generally low hazard and mortality risk and limited coastline poleward of the tropics (Peduzzi et al. 2012).
Any systematic regional changes in WNP TC exposure associated with decadal TC track migration would be highly problematic for a number of reasons. In addition to TC-related human risk factors brought by wind, storm surge, and extreme rainfall, many regions depend on TC passage for their freshwater supply (Jiang and Zipser 2010; Lam et al. 2012). In this case, the benefits of reduced TC exposure need to be tempered by the potential for potable water shortages. Regional changes in TC exposure also need to be considered in light of differing regional sensitivities to TC hazards. For example, the sensitivity to typhoons on Guam, in the Marianas, is relatively small because of local building codes (many homes in Guam are built similarly to concrete bunkers), warning systems, and public education based on the experience that climatologically high TC exposure has brought to the region. Contrarily, regions that experience comparatively low climatological TC exposure, such as Japan for example, can be more sensitive to TC hazards and impacts because of a comparative lack of experience-based adaptation strategies (Cardona et al. 2012). Increasing TC exposure in these regions has the potential to significantly increase vulnerability and, consequently, mortality risk. All of these aspects indicate a need for a more in-depth exploration of WNP TC track and exposure changes, both past and projected.
In this paper, we explore intrabasin changes in WNP TC track characteristics using a combination of observations and numerical model simulations from phase 5 of the Coupled Model Intercomparison Project (CMIP5; Taylor et al. 2012). Section 2 discusses the data and methods applied in the subsequent sections. In section 3, historical TC data from the four sources that provide best-track data in the WNP since 1980 are analyzed with a focus on regional changes in TC exposure patterns. Taking advantage of the fact that ϕLMI is comparatively insensitive to heterogeneities in the best-track data (Kossin et al. 2014), TC track variability is then considered over longer time periods in light of known modes of interannual and decadal variability in an attempt to disentangle natural and anthropogenic causes for the observed migration. In section 4, we explore twenty-first-century projections of WNP TC track and exposure variability simulated by, and downscaled from, an ensemble of numerical CMIP5 models. The paper concludes with a summary and discussion in section 5.
2. Data and methods
Observational TC data were taken from the International Best Track Archive for Climate Stewardship (IBTrACS), version 03r06, which provides estimates of TC center position and intensity every 6 h during the lifetime of each TC (Knapp et al. 2010). There are four best-track data sources from four distinct agencies in the WNP: the Joint Typhoon Warning Center (JTWC), Japan Meteorological Agency (JMA), China Meteorological Administration (CMA), and Hong Kong Observatory (HKO). We note here that these four sources often differ, sometimes substantially, in their reported values for the same TCs, but they are not likely to be completely independent and their ensemble does not necessarily represent four unique datasets. Lifetime-maximum intensity is the peak intensity reached by each TC during its lifetime (Elsner et al. 2008). Best-track data provide intensities in units of knots (kt; 1 kt ≈ 0.51 m s−1). The latitude where LMI is reached is denoted herein as ϕLMI. Only TCs with LMI ≥ 35 kt are used in the following analyses.
As introduced in Kossin et al. (2014), the metric ϕLMI is comparatively robust to the heterogeneities introduced into best-track data by both temporal inconsistency in data quality and analysis as well as interagency procedural differences. For example, ϕLMI is insensitive to the differences in wind-averaging techniques because it does not consider intensity in an absolute sense but only requires identifying when and where the maximum intensity occurs relative to a TC’s lifetime; that is, ϕLMI is unaffected by procedural heterogeneities across the different best-track data sources. The metric ϕLMI is also comparatively insensitive to the usual temporal heterogeneities of intensity estimates, again because the absolute value of intensity is not needed; rather, only knowledge that a maximum intensity was achieved for that particular TC is necessary (ϕLMI represents only the latitude at that time). In the case of ϕLMI, the well-known and substantial temporal heterogeneities in absolute intensity estimates that have been introduced into the best track by numerous changes in technology play a much smaller role. Whatever technology is applied in a given era, it will not systematically change during a TC’s lifetime, and the ability to identify the LMI, whatever its absolute value is, should remain similar as technology changes. The metric ϕLMI simply marks the latitude at that moment of a TC’s lifetime.
TC exposure is expressed in terms of track density on a 2° × 2° latitude–longitude grid in units of ensemble-average number of days per year that a TC center was in that grid box. In the model projection ensembles, the WNP is defined as the region north of the equator and within 100°E–180° longitude.
Monthly indices of El Niño–Southern Oscillation (ENSO) and the Pacific decadal oscillation (PDO) are provided by NOAA’s Earth System Research Laboratory Physical Sciences Division1 (Rayner et al. 2003; Mantua et al. 1997). Annual averages of the indices were taken over July–November (JASON), which comprise the months of greatest WNP TC activity (about 80% of annual TC genesis events occurs in these months). ENSO and PDO variability was removed from the time series of annual-mean ϕLMI by regressing ϕLMI onto their annual-mean index values and forming time series of the residuals of the regressions.
The TCs explicitly simulated by the CMIP5 models are the same as described in Camargo (2013), and the “synthetic” TCs that were downscaled from CMIP5 model output are the same as those described in Emanuel (2013, 2015). Here we have 100 synthetic TCs per year to sample from. Relative annual frequencies are determined by the percentage of random vortex seeds that reach a specified intensity threshold in the downscaling technique (Emanuel et al. 2008).
None of the time series analyzed here exhibits autocorrelation as determined by a Durbin–Watson test statistic, with the exception of the Wellington best-track data time series shown in Fig. 1. The Wellington time series confidence bounds were corrected to account for this. No other adjustments were stipulated.
3. Observed changes in WNP ϕLMI and associated changes in TC exposure
The poleward migration of ϕLMI in the WNP shown in Fig. 1 is consistent across the four best-track data sources available for that region (Fig. 2), which, when considered in concert with the temporal and interagency robustness of the metric ϕLMI, provides high confidence in the fidelity of the mean migration rate. On average, ϕLMI has migrated poleward in the WNP at a rate of 0.56° ± 0.40° decade−1 in the period 1980–20132. Associated with this change in mean ϕLMI is a pronounced and nearly uniform poleward shift in the distribution of ϕLMI (Fig. 3). While the mean ϕLMI increases by only about 2° latitude in the past 34 years, the changes in the tails of the distribution can be large and the probability of a TC reaching peak intensity at the lowest (highest) latitudes have proportionally become substantially smaller (greater) in this time period.
In the postgeostationary satellite period of highest-quality TC data, since around 1980, the overall frequency of WNP TCs has decreased by around 15% (from an average of about 27–28 yr−1 over 1980–96 to about 23–24 yr−1 over 1997–2013) leading to a general decrease in basinwide TC exposure (Liu and Chan 2013; Yokoi and Takayabu 2013; Lin and Chan 2015). Over this same period, however, the changing distribution of ϕLMI has been associated with a change in WNP TC track characteristics in general (e.g., Chu et al. 2012), and the combined effects of basinwide frequency changes and track changes have caused substantial changes in regional TC exposure through changes in track density (Fig. 4). In the more southern regions of the Philippine and South China Seas, the combined effects of the basinwide reduction in TC frequency and the poleward migration has led to a large reduction in TC exposure. In more northern regions, in and around the East China Sea, the basinwide reduction has been offset by the poleward migration. The TC exposure around the Marianas and the Philippines in particular has decreased substantially, by around 50% or more, in the latter half of the period 1980–2013. Contrarily, exposure in the regions of Japan and its Ryukyu Islands, the Korea Peninsula, and Taiwan has increased, by as much as 100% or more in some regions, despite the overall reduction in basinwide TC counts.
On interannual time scales, TC tracks are substantially modulated by the natural variability of ENSO. This modulation occurs in the North Atlantic (Kossin et al. 2010), the eastern North Pacific (Camargo et al. 2008), the Southern Hemisphere (Ramsay et al. 2012), and, of relevance to this study, the WNP (Camargo et al. 2007; Zhao et al. 2010; Goh and Chan 2012; Zhang et al. 2012; Wang and Wang 2013). This variability, however, does not project onto the decadal time scale of global TC migration (Kossin et al. 2014). For the particular case of WNP TCs, this is further demonstrated in Table 1, which shows the statistics of regressions of annual-mean ϕLMI from each data source onto the JASON average of a Niño-3.4 index and compares migration rates with and without ENSO variability included. The JASON-mean Niño-3.4 index explains 24%–43% of the interannual variance of ϕLMI and the regression coefficient is statistically significant for each data source. However, when ENSO variability is removed from the time series of annual-mean ϕLMI, the migration rates are only minimally affected while the uncertainty in the rates actually decreases substantially (cf. the migration rates of ϕLMI and the regression residuals εNiño-3.4 in Table 1). That is, ENSO variability simply adds noise to the ϕLMI time series but does not project onto the migration rate.
The tracks of WNP TCs are also systematically modulated on decadal time scales by various modes of Pacific decadal variability such as the Pacific decadal oscillation (Ho et al. 2004; Liu and Chan 2008; Zhao and Wu 2014; Mei et al. 2015), which raises questions about whether the observed migration is part of a natural decadal cycle or a response to anthropogenic forcing (or, more plausibly, some combination of the two). To explore this with much confidence requires a longer time period than 1980–2013, which is fortunately possible with the available WNP best-track data. The data from the JTWC, CMA, and HKO begin in years 1945, 1949, and 1961, respectively (JMA best-track data begin in 1979). Analysis of the earlier presatellite data should, in general, proceed with some caution because of the potential for data inaccuracies (e.g., Kossin et al. 2013), but we can proceed here with a reasonable level of confidence for three reasons: 1) prior to 1987, WNP TCs were routinely sampled by aircraft reconnaissance (Martin and Gray 1993), which significantly increases the accuracy of the earlier presatellite best-track data; 2) there are multiple best-track data sources, which provides some measure of robustness/consistency testing; and 3) as noted earlier and discussed in Kossin et al. (2014), ϕLMI is comparatively less sensitive to best-track data heterogeneity.
The relationship between the decadal behavior of ϕLMI and the PDO is broadly demonstrated in Fig. 5, which shows an inverse relationship between the two [i.e., cool (warm) PDO phase is related to poleward (equatorward) location of peak TC intensity]. This is consistent with the previous studies noted above. Of particular relevance here is what part of the observed poleward TC migration is explained by known modes of Pacific decadal variability such as the PDO and what part, if any, may be related to other (possibly anthropogenic) forms of slow decadal forcing. Before we proceed with this exploration, it should also be noted that the decadal variability described by the smoothed PDO index is not tacitly internal or natural and very likely describes a conflation of natural and anthropogenic factors (Dong et al. 2014; Meehl et al. 2013; Boo et al. 2015). Furthermore, the expansion of the tropical belt, which Kossin et al. (2014) showed is consistent with the observed poleward TC migration, has been linked to both PDO variability and anthropogenic aerosol forcing (Allen et al. 2014), as well as various other factors, both natural and anthropogenic (Lucas et al. 2014).
In an attempt to disentangle the contributions of these factors to WNP TC migration and exposure changes, Fig. 6 shows time series of the residuals of annual-mean ϕLMI regressed onto indices of the PDO and ENSO. The multiple regression coefficients are statistically significant, as expected, and the JASON-mean PDO and Niño-3.4 indices together explain 24%, 28%, and 34% of the variability of annual-mean ϕLMI derived from the JTWC, CMA, and HKO best-track data, respectively (Table 1). When ENSO and PDO variability are removed (regressed) from the annual-mean ϕLMI time series, the residuals of the regression (εNiño-3.4, PDO in Table 1) maintain a statistically significant poleward migration of 0.17°, 0.21°, and 0.33° decade−1, respectively. This result is highly suggestive that the observed WNP TC migration can be attributed to a combination of natural and anthropogenic factors with a substantial portion of the migration occurring independently of the known dominant modes of variability in the region.
To further explore attribution linkages between anthropogenic forcing and WNP TC migration, the following section considers WNP TC behavior projected into the twenty-first century. Before we enter into that topic, it should be explicitly acknowledged that we have limited our discussion and analyses to questions of TC track variability and have not considered TC intensity variability in an absolute sense, although absolute intensity is clearly an important factor in assessing TC hazard and mortality risk. Introducing absolute intensity data into these analyses, whether from historical best-track records or from numerically simulated TCs, introduces another layer of uncertainty (Kossin et al. 2013; Camargo and Wing 2016), and here we have chosen to focus only on the most robust metrics. Nonetheless, it should be noted that analysis of absolute intensity-based measures (e.g., TC power dissipation or simply mean intensity) shows similar patterns of change in both the best-track data and the model-based data discussed in the following sections. That is, there are no clear indications that absolute intensity changes have had or are projected to have any offsetting effects on the changes in track density and associated changes in TC exposure and risk explored here.
4. Projected changes in WNP ϕLMI and associated changes in TC exposure
If the observed poleward TC migration in the WNP is forced in part by anthropogenic factors, and consequently persists into the present century, then the associated changes in TC exposure would be expected to progressively impact various regions in the WNP. A number of previous studies have explored relationships between WNP TC track characteristics and anthropogenic climate change and have identified systematic shifts in TC genesis location and steering flow that combine to cause TC track shifts (Wu and Wang 2004; Li et al. 2010; Wang et al. 2011; Murakami et al. 2011, 2012; Yokoi et al. 2013; Bell et al. 2013; Mori et al. 2013; Manganello et al. 2014; Wu et al. 2014; Colbert et al. 2015; Mei et al. 2014, 2015; Roberts et al. 2015; Wang and Wu 2015). There is, however, some divergence among the model projections in some of these studies. While some high-resolution climate models (usually forced with fixed sea surface temperature) show a poleward shift that is generally attributed to the expansion of the tropics (e.g., Murakami et al. 2012; Wu et al. 2014), other modeling studies emphasize an eastward shift of the TC tracks toward the central Pacific (Li et al. 2010; Murakami et al. 2011; Yokoi et al. 2013; Mori et al. 2013). In some cases, northward and eastward shifts occur simultaneously as a northeastward migration (Murakami et al. 2012; Colbert et al. 2015; Roberts et al. 2015), and one model projects a southeastward shift (Manganello et al. 2014). While tempering our results with knowledge of this lack of firm consensus in model projections of WNP TC track variability, here we expand on previous work by analyzing the tracks of TCs formed in a suite of CMIP5 models (Camargo 2013) and downscaled from CMIP5 environmental data (Emanuel 2013, 2015), focusing on changes in ϕLMI and associated patterns of change of TC exposure, and comparing these patterns to observed patterns. Further tempering the interpretation of our results, a number of performance issues related to TCs generated by climate models or downscaled from climate model output have been documented (e.g., Camargo and Wing 2016; Emanuel 2015, and references therein). Here we will endeavor to present the model output as transparently as possible with the understanding that there are inherent uncertainties that exist.
a. TCs explicitly simulated by the CMIP5 models
Despite their relatively low spatial resolution, most of the CMIP5 models generate TC-like storms that can be tracked and analyzed. Here we use an ensemble of storms generated by 10 CMIP5 models, as described in Camargo (2013). The 10 models are described in Table 1 in Camargo (2013) with four models (INM-CM4.0, MIROC-ESM, MPI-ESM-LR, and NorESM1-M) omitted owing to a lack of TC generation in the WNP basin. The analyses here are applied first to CMIP5 historical simulations (1980–2005) for comparison with the observations and then to twenty-first-century projections under the representative concentration pathway 8.5 (RCP8.5) emissions scenario (Riahi et al. 2011).
In the historical simulations, the model ensemble tends to underestimate TC frequency, producing about 5.9 WNP TCs per year on average in the first half of the simulation period 1980–2005 and about 5.6 WNP TCs per year in the second half (roughly from one-fourth to one-fifth of the observed frequency). The underestimation of TC frequency is due, in part, to deficiencies in the models themselves, but there is also measurable sensitivity to the detection algorithm applied after the fact to the model output (Horn et al. 2014; Walsh et al. 2015). The rate of poleward migration and the poleward shift of the distribution of ϕLMI are also underestimated (cf. Figs. 7a and 8a to Figs. 2 and 3), although the general patterns of track density and track density change between the first and second halves of the period (Fig. 9) share some similarities with the observed changes (Fig. 4). If poleward TC migration is physically linked to the expansion of the tropics, the underestimation of the poleward migration in the historical simulations might be expected as the CMIP5 models have been independently shown to substantially underestimate the expansion of the tropical Hadley circulation (Hu et al. 2013; Nguyen et al. 2015).
In the twenty-first-century RCP8.5 projection simulations, the models again tend to underestimate TC frequency, producing an ensemble average of about 6.1 WNP TCs per year in the first 30 years of the simulation period 2006–99 and about 5.6 WNP TCs per year in the last 30 years. The projected reduction of WNP TC frequency is consistent with the general model consensus in that basin (Christensen et al. 2013, their Fig. 14.17). The projected rate of migration of annual-mean ϕLMI is about 0.2° ± 0.1° decade−1 (Fig. 7b), which is slower than the observed rate of about 0.56° ± 0.40° decade−1 in the period 1980–2013 (Fig. 2) but, interestingly, about the same as the observed poleward migration rate when ENSO and PDO variability are removed from the time series (Fig. 6). Similar to observations, the projected poleward migration of annual-mean ϕLMI manifests within a fairly uniform poleward shift in the distribution of ϕLMI (Fig. 8b) and is associated with substantial changes in projected regional TC exposure (Fig. 10). The projected pattern of change of TC exposure is similar to the past observed changes (cf. Figs. 4 and 10), and the projected decrease in basinwide TC frequency either compounds or offsets the TC exposure changes as a result of track changes. Similar to the observed changes, in the regions of the Philippines and Marianas and regions surrounding the South China Sea, there are substantial projected decreases in exposure, by as much as 50% or more, in the last 30 years of this century compared to the first 30 years. Contrarily, the exposure in regions around the East China Sea, such as Japan and its Ryukyu Islands and the Korea Peninsula, is projected to increase despite the overall projected reduction in basinwide TC counts.
b. TCs downscaled from CMIP5 model output
In addition to TCs that are generated within the CMIP5 models, the CMIP5 model environmental fields can be downscaled to create large samples of synthetic TCs (Emanuel et al. 2006, 2008; Emanuel 2010). Here we use the global TC set described in Emanuel (2013, 2015) and repeat the analyses performed above using the WNP TCs. The downscaling method allows the generation of large numbers of synthetic TCs (for this study, we have 100 TCs per year to sample from), and the annual frequency is determined by the percentage of random vortex seeds that reach a specified intensity threshold (Emanuel et al. 2008). Similar to section 4a above, we use an ensemble of TCs downscaled from output from multiple CMIP5 models. The eight models used here are listed in Table 1 of Emanuel (2015) with the addition of NCAR’s CCSM4.
Before analyzing the WNP TCs downscaled from the twenty-first-century CMIP5 model projections, we consider WNP TCs downscaled from past reanalysis data for comparison with the observations. The TCs are downscaled from the European Centre for Medium-Range Weather Forecasts (ECMWF) interim reanalysis (ERA-Interim; Dee et al. 2011) over the period 1980–2013. Here the observed mean annual frequency from the four best-track sources is used to sample the TCs from the larger sample of synthetic TCs. The time series of annual-mean ϕLMI from the downscaled WNP TCs (Fig. 11a) shows a significant poleward migration rate of about 0.44° ± 0.30° decade−1, which is comparable but somewhat less than the mean observed rate of about 0.56° ± 0.40° decade−1 within the same period. The poleward shift of the distribution of ϕLMI is captured well by the TCs downscaled from the reanalysis data (cf. Fig. 12a to Fig. 3), and the pattern of change in TC exposure (Fig. 13) is broadly similar to the observed pattern of change (Fig. 4), with increased exposure in regions of the East China Sea and generally decreased exposure in the region around the northern Philippines (cf. Fig. 13 and Fig. 4). Similar results are found using data from the Modern-Era Retrospective Analysis for Research and Applications (MERRA; Rienecker et al. 2011) and the National Centers for Environmental Prediction–National Center for Atmospheric Research (NCEP–NCAR) reanalyses (Kalnay et al. 1996).
When the same downscaling method is applied to the projected output from the eight CMIP5 models, a significant poleward migration of WNP ϕLMI emerges over the twenty-first century, with a rate of 0.13° ± 0.04° decade−1 (Fig. 11b), which is also associated with a poleward shift of the distribution of ϕLMI (Fig. 12b). The projected migration rate and interannual variability of ϕLMI in the downscaled TCs are smaller than found in the explicitly simulated TCs (for these analyses, the larger sample of available synthetic TCs are used, which is a likely factor in the reduced variance). As found in the explicitly simulated TCs, the poleward shift of ϕLMI is associated with regional changes in TC exposure. However, contrary to the explicitly simulated TCs, the projected frequency of the downscaled WNP TCs increases over the twenty-first century (Emanuel 2013), which causes a widespread and significant increase across the entire basin (Fig. 14). The increase in basinwide frequency is substantial (about 22% when comparing the first and last 30 years of the simulation) and dominates the TC exposure changes due to track changes. To better separate the effect of the track and frequency changes, Fig. 15 shows the projected changes in TC exposure with annual frequency held fixed. When this is done, the patterns of exposure change due to track migration are found to be qualitatively similar to those of the explicitly simulated TCs, although with some quantitative differences.
An observed characteristic of past WNP TC track density (Figs. 4a,b) that is also captured fairly well by the explicitly simulated TCs, both from the historical (Figs. 9a,b) and projection simulations (Figs. 10a,b), is the distinctly separated maxima on the eastern and western flanks of the northern Philippines. This characteristic is not as clearly captured in the track density downscaled from reanalysis (Figs. 13a,b) and is only marginally present in the projected track density (Figs. 15a,b). The maximum to the east of the northern Philippines is well defined, which demonstrates that the downscaled TCs tend to spend relatively less of their lifetimes in the South China Sea, both in the past and projected simulations, than the observed and explicitly simulated TCs. For the downscaled TCs, this is largely caused by the intensity model used by the method, which has difficulty simulating the recovery of storms that have passed over land, as often happens when storms pass over the Philippines en route to the South China Sea. As a side note, this stands in contrast to the simulated TCs of Mei et al. (2015), which display a similar minimum of track density in the South China Sea, but due to a substantial dearth of TC genesis in the region.
5. Summary and discussion
There is some general consistency in the past and projected changes in WNP ϕLMI and the associated changes in TC exposure, although differences exist in the various data considered here. Past observed changes in the period 1980–2013 show a fairly uniform poleward shift in the distribution of ϕLMI, which manifests in large changes in the tails of the distribution (i.e., in the deep tropics and subtropics) and large changes in regional patterns of TC exposure. Decreased exposure is observed in the regions of the Philippine and South China Seas, including the Marianas, the Philippines, Vietnam, and southern China, and increased exposure is observed in the region of the East China Sea, including Japan and its Ryukyu Islands, the Korea Peninsula, and parts of eastern China. Decreases of 50% or more and increases of 100% or more are observed in some regions. These changes result from a combination of basinwide TC frequency decreases (about a 15% decrease in this period) and a robust and significant poleward migration of ϕLMI of about 0.56° ± 0.40° decade−1. In particular, this combination of factors has a compounding effect at lower latitudes (frequency decreases in concert with track shifts away from the region) and an offsetting effect at higher latitudes (frequency decreases in concert with track shifts into the region). In the latter case, some regions have experienced an increase in exposure as the effect of the track shifts has dominated the effects of reduced basinwide frequency.
For TCs explicitly simulated in the CMIP5 historical simulations, the models substantially underestimate basinwide TC frequency (producing fewer than 6 yr−1 in an ensemble average, compared to about 25 yr−1 observed) and the poleward migration of ϕLMI is underestimated and statistically insignificant (0.3° ± 0.8° decade−1). Still, the general pattern and patterns of change of TC exposure in the historical simulations share a number of similarities to the observed patterns described above. For the TCs downscaled from reanalysis data, the poleward migration rate of ϕLMI is somewhat underestimated but significant (about 0.4° ± 0.3° decade−1). The patterns of TC exposure, however, do not adequately capture the substantial maximum observed in the South China Sea, in the region west of the Philippines, which makes the patterns of exposure change for the downscaled TCs in that region somewhat suspect.
Twenty-first-century projections of TCs, both explicitly generated within the CMIP5 models and downscaled from CMIP5 model output, identify a continued significant poleward migration of ϕLMI (0.2° ± 0.1° and 0.13° ± 0.04° decade−1, respectively) and associated changes in TC exposure patterns that share some similarities with past changes. However, while the explicitly generated TCs decrease in annual basinwide frequency over the twenty-first century, the frequency of downscaled TCs increases, which leads to substantially different projections of TC exposure. For projections based on the downscaled TCs, the exposure increases across the entire basin as the increased frequency dominates the effect of the track changes. The question, then, of how WNP TC frequency will change in the present century is critical to address from the standpoint of regional preparations for changes in TC exposure, particularly in the more southern regions of the WNP basin where the sign of the exposure changes can switch between negative and positive.
The twenty-first-century projections discussed here are based on the RCP8.5 emissions scenario and provide evidence that continued anthropogenic climate forcing can manifest in a continued poleward migration of WNP TC tracks. This is consistent with our analysis of extended time series of historical best-track data over the past ~70 years in which indices of ENSO and the PDO were removed via multiple regression. In this case, the poleward migration rates found in the best-track data from multiple sources remain statistically significant and on the order of about 0.2° decade−1. The residual migration rate after removing ENSO and PDO variability suggests that the migration has been occurring outside of the dominant known modes of variability in the region, which provides further evidence for an anthropogenic contribution.
While the poleward migration of WNP TCs is projected to occur in concert with basinwide frequency changes (of either sign), the associated changes in regional TC exposure should also be considered in concert with projected increases in TC intensity (Christensen et al. 2013; Peduzzi et al. 2012), which has the potential to further compound changes in regional human mortality risk. Another factor that should be considered is the potential for track changes to affect potable water supplies in some regions that depend on TC passage for their freshwater supplies. But perhaps the greatest potential threat from the changes in TC exposure discussed here lies in the potential for shifting patterns of exposure relative to patterns of sensitivity and resilience. The benefits of reduced TC exposure to regions that are comparatively insensitive and resilient to TC hazards because of past prolonged exposure (e.g., the Marianas) would likely not be as great as the consequences of increased exposure in regions that are more sensitive because of historically less frequent exposure to TC hazards (e.g., Japan and the Korea Peninsula). The known weak bias in the CMIP5 models’ ability to simulate the observed tropical expansion should also be considered, as the future migration rate and associated changes in TC exposure could be larger than shown here.
The representative concentration pathway used here (RCP8.5) is based on a presumption of little to no change from the present trajectories of greenhouse gas concentrations (sometimes colloquially referred to as “business as usual”). The potential for continued changes in regional TC exposure and human mortality risk demonstrated here could tenably be reduced, and even reversed, by changes in emissions standards.
JPK is supported by the National Oceanic and Atmospheric Administration National Centers for Environmental Information (NOAA/NCEI). KAE acknowledges support from National Science Foundation Grant 1461517 and Office of Naval Research Grant N000140910526. SJC acknowledges support from NOAA Grant NA110AR4310093 and National Science Foundation Grant AGS-1143959.
Available from http://www.esrl.noaa.gov/psd/gcos_wgsp/Timeseries/Nino34/ and http://www.esrl.noaa.gov/psd/gcos_wgsp/Timeseries/PDO/, respectively.