Future changes in intense tropical cyclone hazards in the Pearl River Delta region: An air-wave-ocean coupled model study

Future changes in intense tropical cyclone hazards in the Pearl River Delta region: An air-wave-ocean coupled model study | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Future changes in intense tropical cyclone hazards in the Pearl River Delta region: An air-wave-ocean coupled model study Zhenning Li, Jimmy C.H. Fung, Mau Fung Wong, Shangfei Lin, Fenying Cai, and 2 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-3625663/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 09 Mar, 2024 Read the published version in Natural Hazards → Version 1 posted 5 You are reading this latest preprint version Abstract The Pearl River Delta (PRD) region is highly vulnerable to tropical cyclone (TC)-caused coastal hazards due to its long and meandering shoreline and well-developed economy. With global warming expected to continue or worsen in the rest of the 21st century, this study examines the TC impact on the PRD coastal regions by reproducing three intense landfalling TCs using a sophisticated air-wave-ocean coupled model of high spatial resolution (1-km atmosphere and 500-m wave and ocean). The simulations are conducted using present-day reanalysis data and the same TCs occurring in a pseudo-global warming scenario projected for the 2090s. Results indicate that the coupled model accurately reproduces the air-wave-ocean status during the TC episodes. The 2090s thermodynamic status effectively increases the intensity of intense TCs, leading to more severe coastal hazards including gale, rainstorm, and storm surges and waves. On average, the maximum surface wind speed within 50–200 km to the right of the TC center can increase by 4.3 m/s (+ 22%). The 99th and the 99.9th percentile of accumulated rainfall will increase from 405 mm to 475 mm (+ 17.3%), and from 619 mm to 735 mm (+ 18.6%), respectively. The maximum significant wave height at the ocean is lifted by an average of 57 cm (+ 13.8%), and the coastline typically faces a 40–80 cm increase. The maximum storm surges are lifted by 30–80 cm over the open sea but aggravate much higher along the coastline, especially for narrowing estuaries. For Typhoon Vicente (2012), there is more than a 200 cm wave height increase observed both at open sea and along the coastline. In the 2090s context, a combination of mean sea level rise, storm surge, and wave height can reach more than 300 cm increase in total water level at certain hot-spot coastlines, without considering the superposition of spring tides. Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Introduction Landfalling tropical cyclones (TCs) are intense weather systems that can cause a range of coastal hazards, including storm surges and waves, coastal erosion, wind induced damage, rainstorms, and secondary hazards like landslides and inland flooding (Tran et al. 2022 , Wu et al. 2022 ). The impact of landfalling TCs can vary depending on the severity and location of the storm, as well as the vulnerability of the affected communities (Wang et al. 2008 , Knutson et al. 2010 , Knutson et al. 2019 , Chu et al. 2020 ). The Pearl River Delta (PRD) region, with its long and meandering coastline and well-developed economy, is highly vulnerable to TC-caused coastal disasters (Chan 2008 , Chen et al. 2020 ). Recent TCs, such as Typhoon Hato (2017) and Mangkhut (2018), have demonstrated the devastating impact these storms can have, causing widespread damage and losses (Liu and Chan 2020 ). The estimated direct economic loss due to Mangkhut (2018) and Hato (2017) in Hong Kong is about HK $ 4.60 billion (~ US $ 590 million) and HK $ 1.21 billion (~ US $ 155 million), respectively (Choy et al. 2020 ). Meanwhile, the impacts of climate change on TCs have the potential to significantly alter their influence on future coastal residents (Knutson et al. 2010 , Tsuboki et al. 2015 , Tran et al. 2022 ). The global warming trend is expected to continue or even worsen in the rest of the 21st century (Knutson et al. 2010 ). A recent comprehensive assessment by Knutson et al ( 2020 ) shows that a 2°C anthropogenic global warming is projected to impact TC activity as follows: 1) Medium-to-high confidence in the increased precipitation rates by around + 14%, or close to the rate of tropical water vapor increase with warming; 2) Medium-to-high confidence in the + 5% stronger lifetime maximum surface wind speeds; 3) Medium-to-high confidence in a + 13% increase of the proportion of TCs that reach very intense (category 4–5) levels; and 4) The most confident TC-related projection is that sea level rise accompanying the warming will lead to higher storm inundation levels, assuming all other factors are unchanged (Knutson et al. 2020 ). However, research conclusions are more mixed, and confidence levels are lower for the basin-scale or regional projections (Knutson et al. 2020 , Bloemendaal et al. 2022 ). Due to the limitation in spatial resolution, Global Climate Model (GCM) studies also introduce large biases when reproducing TC activities on a regional scale (Bloemendaal et al. 2022 ). These factors fall far short of satisfying the demand to guide regional climate change adaptation and resilience. On a regional scale, air-sea interaction plays a crucial role in the development of TCs. As noted by Emanuel ( 1986 ), the ocean surface enthalpy fluxes transferred to the atmosphere are the main energy sources of TCs. The Coupled Model Intercomparison Project Phase 6 (CMIP6) model results indicate that in the future, there will be a warmer and fresher surface ocean. This implies an increased thermal stratification in a warming climate, which in turn will feedback on the TC development (Huang et al. 2015 , Balaguru et al. 2016 ). Additionally, studies have shown that at extremely high surface wind speeds (30 m/s or higher), the roughness length, determined by the sea state, tends to plateau. This can maintain or even enhance the inner-core intensity of the TC (Moon et al. 2004 , Chen et al. 2013 ). Numerical studies incorporating such effects can construct more realistic TC structures than atmosphere-only or even atmosphere–ocean coupled models (Davis et al. 2008, Chen et al. 2013 , Wu et al. 2019 ). Consequently, operational centers such as the European Centre for Medium Range Weather Forecasts (ECMWF) are now using atmosphere-wave-ocean (AWO) coupled models to generate routine global products for better performance in TC forecasts (Mogensen et al. 2017 , Magnusson et al. 2019 ). Incorporating a dynamically evolving 3-dimensional ocean model and a wave model at convection-permitting scales (i.e., 4 km or finer; see Zhang et al. 2011 ) is crucial to better capture the TC-related coastal disasters on a regional scale. By generating crucial products related to coastal hazards across the atmosphere (wind and rainfall), sea wave (significant wave height), and ocean (storm surge level) from the synergistic framework, valuable supporting evidence can be provided to guide regional adaptation and resilience against TCs in a changing climate. The objective of this study is to investigate the impact of future TC-related coastal disasters on the PRD region in a changing climate, using a sophisticated AWO coupled model. The PRD region is home to numerous large populations and is a pivot of manufacturing, finance, and marine industries such as shipping and aquaculture in China. To achieve this objective, we examine three intense landfalling TCs since 2010 in “parallel digital worlds”: a present-day world and a warming future world. Specifically, we conduct simulations using (a) present-day reanalysis data and (b) the same TCs occurring in a pseudo-global warming (PGW) scenario projected for the end of the 21st century. Conventional PGW experiments have been extensively used in atmospheric-only regional models such as the WRF (Chen et al. 2020 , Jyoteeshkumar et al. 2021, Tran et al. 2022 ). To our knowledge, this study is the first to implement PGW in an AWO fully coupled system over the PRD region, achieving such high horizontal resolutions (1-km atmosphere and 500-m ocean and sea wave). The perturbated warming initial and boundary conditions are derived from both atmosphere and ocean GCMs projected by CMIP6 and enforced in both the atmosphere and the ocean components of the fully coupled system. The results of this study will provide valuable insights into the potential impacts of future TC-related coastal disasters on the PRD region and guide regional adaptation and resilience against such disasters in the changing climate. Modeling system and experimental design 2.1 The Coupled Ocean–Atmosphere–Wave–Sediment Transport (COAWST) modeling system The COAWST system serves as the fundamental framework for this study. It is an agglomeration of open-source modeling components tailored to investigate coupled processes of the atmosphere, ocean, and waves in coastal regions. Since its release by the United States Geological Survey in 2010, it has been widely used in various studies (Warner et al., 2010 ). COAWST has been applied to examine several Atlantic TCs that caused considerable damage in North America, such as Ivan (2004), Ida (2009), and Sandy (2012) (Olabarrieta et al. 2012 ,Zambon et al. 2014a , Zambon et al. 2014b ). In this study, we use COAWST v3.6, which comprises the following modules: (1) an atmospheric component named Weather Research and Forecasting Model (WRF v4.1.2); (2) a Regional Ocean Modeling System (ROMS) v3.7 as an oceanic component; (3) a SWAN v41.31 as a sea wave component; and (4) a Model Coupling Toolkit (MCT) v2.6.0 as a coupler. In previous work by Li et al ( 2022 ), the COAWST has been customized and tuned extensively to suit local weather and climate studies in the PRD region. It was used to investigate the impact of sea wave-wind interaction on TC development, including typhoon tracks and intensities. The current system builds upon this framework and is organized into a hierarchical architecture with two tiers, each tailored to simulate dominant physical processes at corresponding spatial scales (see Fig. 1 ). Tier 1 is designed for open sea scale modeling. The WRF model is configured with a horizontal resolution of 9 km × 9 km in its outermost domain (WRF D01). The ROMS and SWAN models (ROMS + SWAN D01) are nested within WRF D01, with an ocean eddy-resolving resolution of 2.2 km × 2.2 km. The intermediate domain for WRF (WRF D02) is embedded in ROMS + SWAN D01 at sea, with a convection-permitting resolution of 3 km × 3 km. Tier 2 is designed to provide high spatial resolution for modeling the PRD coastal region. The innermost WRF domain (WRF D03) has a resolution of 1 km × 1 km and covers the entire PRD region. The innermost ROMS + SWAN domain (ROMS + SWAN D02) has a resolution of 500 m × 500 m and covers the Pearl River estuary and the nearshore sea. The configurations for each of the three component models follow Li et al. ( 2022 ), with minor modifications made for computational efficiency. These modifications include a smaller mesh size in WRF D01/02 and ROMS + SWAN D01. In addition, the number of vertical layers in the WRF model is reduced to 39 layers, but the number of vertical layers within 1,500 m above the surface is increased to 19 for improved accuracy in representing physical processes and structures in the planetary boundary layer. For air-wave interaction, we implemented Lin and Sheng ( 2020 ) scheme, which is embedded in the Mellor-Yamada-Nakanishi-Niino (MYNN; Nakanishi and Niino 2004 ) planetary boundary layer scheme and associated surface layer schemes. Detailed configurations are summarized in Table S1 . Tide information was obtained from the Oregon State University Tidal Prediction Software (OTPS; Egbertand and Erofeeva 2002), which provides data on eight primary harmonic constituents (M2, S2, N2, K2, K1, O1, P1, Q1). 2.2 Experimental design In this study, we selected three intense landfalling TCs - Mangkhut (2018), Hato (2017), and Vicente (2012) - which caused significant impacts on the PRD coastal areas. These selected TCs were then simulated using the COAWST framework. Two groups of experiments were conducted: present-day control (CTRL) experiments and pseudo-global warming (PGW) experiments. The CTRL experiments aimed to reproduce the selected TCs and related coastal disasters with high fidelity, including rainfall, wind, storm surges, and storm waves. The PGW experiments aimed to simulate the same storms occurring in a global warming scenario projected for the end of the 21st century. We utilized the high-resolution version of the Max Planck Institute Earth System Model (MPI-ESM-1-2-HR) from CMIP6 datasets to extract PGW signals (Gutjahr et al. 2019 ). Thermo-variables, including air temperature, humidity, sea water temperature, and sea water salinity, were extracted from 2090–2099 in the Shared Socioeconomic Pathways 5-8.5 (SSP5-8.5) scenario and 2005–2014 in the historical simulation. The differences between the two epochs were then added to the initial and boundary conditions of the WRF and ROMS models. Figure S2 and figure S3 show the changes in July-August-September 2090s for the SSP5-8.5 scenario, including a projected 2.4°C warming of Sea Surface Temperature (SST) in the outermost domain and a 2.6°C warming of the 25–50 m subsurface ocean. Specific humidity near the surface level is expected to increase by 3.0–4.0 g/kg, and the surface air temperature will increase by 3.2°C, with more than 6.0°C warming in the upper level (near 250 hPa, see figure S3). We compared the results with the multi-model ensemble mean and a previous study by Chen et al. ( 2020 ) based on CMIP5 output and found similar changes in variable magnitudes and large-scale patterns. To our knowledge, this is the first time the PGW method has been applied in a regional AWO coupled model to study future TC activities. Results 3.1 Simulation basics and model evaluation We begin by comparing the tracks and intensities of the selected TCs in both the CTRL and PGW experiments. The successful reproduction of TC evolutions in the CTRL experiments provides a higher level of confidence for future projections. As depicted in Figs. 2 a–c, the COAWST model performed quite well in capturing the tracks of all selected TCs, including the stable steering flow caused northwest-southeast tracks of Mangkhut (2018) and Hato (2017) from the western Philippine waters to the southern China coast. Notably, the model effectively simulated the subtle but essential episode of Vicente (2012), where it slowed down and remained stagnant over the northern South China Sea on July 22, 2012, as the subtropical ridge northeast of the storm suddenly weakened (Chen et al. 2017 ). While the track for TC Vicente (2012) exhibited an evident southward track shift compared with the observation and future run near 111°E, 22°N. This southward track deflection is a typical behavior regarding the terrain effect over the western Guangdong Province, similar effect in Taiwan has been well archived by Wu et al. ( 2015 ). Figures 2 d–f and 2 g–i display the central minimum sea level pressures (SLPs) and maximum 10-m wind speeds, respectively. It is noteworthy that the COAWST underestimates the central SLPs in all cases, which is a known bias of the WRF model, even in cloud-resolving resolution without vortex initialization (Wu et al. 2006 , Tong et al. 2018 , Chang et al. 2020 , Chen et al. 2020 ). To address this issue, we tested the insertion of the Rankine vortex to the initial condition based on the real TC intensity and location through the WRF TC “bogus” scheme in Hato (2017). However, this technique resulted in a remarkable overestimation of surface wind and unreasonably high storm surge levels when compared with in-situ observations (discussed later). As a trade-off, experiments were conducted without “bogus” TC to better simulate more direct coastal disaster indicators, despite the underestimation of central SLPs. We compared the changes in tracks and intensities between the PGW and CTRL experiments. To ensure the constancy of ensuing analysis, we used large-scale spectral nudging of steering wind with appropriate relaxation coefficients, and all TCs in the PGW experiments kept similar tracks to those in CTRL. Additionally, since spectral nudging tends to suppress the TC intensity (Moon et al. 2018 ), we conducted experiments without nudging, and the results indicated that the tracks could deflect by hundreds of kilometers, although the minimum SLP deepened by around 5 hPa (Figure S4). For intensity comparison between the PGW and CTRL, both Mangkhut (2018) and Hato (2017) exhibited a slight drop in the minimum SLP by 3–7 hPa, which is similar to recent studies on TCs impacting Southeast Asia (Chen et al. 2020 , Jyoteeshkumar et al. 2021, Tran et al. 2022 ). However, for Vicente (2012), we found a sharp drop in the central SLP in the PGW experiment, which was 15 hPa lower than that in CTRL. The unusual stationary track over the northern South China Sea allowed Vicente (2012) to absorb excess energy from the warmer and stabler ocean in a warming future. In the CTRL simulation, the mean SST beneath the stationary location of Vicente (2012) dropped from around 29.3°C to 27.5°C, whereas in PGW, it dropped from around 32.1°C to 30.0°C. Even after strong upper-ocean mixing by Vicente (2012), the SST is still warmer than the initial SST over that location in CTRL (Figure not shown). Regarding the maximum surface wind speed, the time series show similar evolution, but intense peak wind strengthening can be found in the PGW. Mangkhut (2018) and Hato (2017) exhibited a peak surface wind speed that was approximately 4 m/s (~ 10%) higher, which is consistent with the results from Chen et al. ( 2020 ) study on three different TCs affecting the PRD in the far future (2075–2099). Particularly, Vicente (2012) even attained nearly 10 m/s (~ 30%) increase in maximum 10-meter wind speed due to further intensification. In addition to evaluating performance in tracks and intensities, evaluating the performance of the COAWST using in-situ station-based observations are also important, as major coastal disasters caused by TCs are fierce wind, torrential rainstorms, and formidable storm surges and waves. Consequently, 10-m wind speed and accumulated rainfall are compared with 2739 surface automatic weather stations in Guangdong province (Figure S5a). Quality control is implemented, and stations with accumulated rainfall less than 50 mm in each TC exposure are eliminated in rainfall evaluation. Table 1 shows that the model overestimates the wind speed by 1.2 m/s (+ 7.4%) with 5.7 m/s root mean square error (RMSE) during the TC exposure. For accumulated rainfall, the COAWST underestimated accumulated rainfall by 4.8 mm (-2.8%) across all qualified stations. While the RMSE is 70.6 mm, relatively small mean bias and large RMSE indicate the model performs okay in reproducing the mean rainfall amount but poorly captures the fine-scale spatial distribution. We then compared storm surge levels using the Hong Kong Observatory tide gauges in Hong Kong (Fig. S4b). Tide-induced water level fluctuations relative to the Chart Datum were linearly combined with the pure storm surge level produced by COAWST to represent the total surge. Results indicate that the simulated surge level is overestimated by 0.13m (+ 5.5%) with RMSE 0.35m. The two metrics show that the model can reproduce the surge variability for the three cases with reasonable accuracy (order of 10 cm) in this region, similar to the magnitude compared with previous studies using various models (Wang et al. 2008 , Yin et al. 2020 , Chen et al. 2020 ). Following Li et al. ( 2022 ), significant wave height is compared with in-situ wave observations at four buoys (QF303, QF304, QF306, and QF307) located in the continental shelf sea in the northern South China Sea (see Figure S4c), within 100-km distances to the nearest coast. Due to limited data availability, the significant wave height is only evaluated for Mangkhut (2018), with − 0.56 m (-9.3%) underestimation in mean bias and RMSE 1.12 m. Still, it shows the ability of the wave model to produce reliable wave height, given that the wind can be reasonably simulated by the atmospheric model. Table 1 Station observations in comparison with model simulations Mean bias Relative bias Root mean square error 10-m wind speed (m/s) + 1.2 + 7.4% 5.7 Accumulated rainfall (mm) -4.8 -2.8% 70.6 Storm surge level (m) + 0.13 + 5.5% 0.35 Significant wave height (m) -0.56 -9.3% 1.12 3.2 Changes in wind and precipitation Severe wind can cause dramatic damage in coastal regions, especially for metropolises such as Shenzhen and Hong Kong in the PRD region. During Mangkhut (2018), the collapse of trees caused extensive road blockage, paralyzing land transportation. Broken windows and falling debris due to building damages seriously threaten the safety of residents (Choy et al. 2020 ). Figure 3 estimates the grid-wise maximum 10-m wind speed change using the PGW experiment minus the CTRL. Negative values are masked out for clarity. All three cases result in an increase in the 10-m maximum wind speed in the PGW experiments for the areas exposed by the TC circulation, especially for the coastal area where it is exposed to the TC right-of-track. Please note that the difference in 10-m wind are probably coming from the variation (drift) of TC’s center in the two experiments rather than the change in intensity. To eliminate this effect, quantitative analysis based on the distance to TC’s center has been implemented. For example, the maximum 10-m wind speed within 50–200-km right to the TC center increases by 4.3 m/s (22%). Tran et al. ( 2022 ) attributed the increase in the TC right-of-track surface maximum wind speed to the faster transition speed of the TCs. While our study's three cases cannot definitively conclude that TCs are moving faster, the stronger intensity of these storms will certainly contribute to higher wind speed exposure for vulnerable coastal regions. In the case of Vicente (2012), which experienced the largest intensity increment in the PGW experiment, large areas of increased wind speed were found within the coverage of the TC's circulation. On average, the maximum 10-m wind speed within 50–200 km to the right of the TC center increased by 6.9 m/s (39%). In a warming future, special attention should be paid to TCs with stationary tracks, like Vicente (2012), as they are likely to experience abrupt and acute rapid intensification near the shore, resulting in much greater wind impacts on coastal areas (Fig. 3 c). TCs often bring torrential rainfall, which can cause severe flooding and landslides, adding to the list of coastal disasters associated with these storms. The resulting damage to infrastructure, property, and loss of life can be catastrophic, particularly in low-lying or densely populated areas. Figure 4 shows the grid-wise histogram of accumulated rainfall (mm) and hourly rain rate within a radius of 500 km of the TC center over the lifetime. The blue solid curve represents the CTRL experiment, and the red dashed curve represents the PGW experiment. A robust accumulated rainfall increase was observed in all three cases, with an average of 96.9 mm in CTRL and 107.9 mm in PGW, representing an 11.7% increase. Previous studies have shown that the global TC precipitation increases by slightly more than 7% per 1°C of tropical SST warming (the Clausius–Clapeyron relation) in multiple climate model projections (Knutson et al. 2020 ). Given that the domain-averaged additional SST perturbation is about 2.4°C, an 11.7% increase in accumulated rainfall is slightly lower than the Clausius–Clapeyron relation. However, the extreme tails exhibit a slightly super-Clausius–Clapeyron relation, suggesting that special attention should be paid to extreme rainfall events. As shown in Fig. 4 a–c, the PGW experiment significantly increased the percentage of accumulated rainfall by more than 200 mm in the three cases, from 11.0–13.9% (+ 26.4%). The 99th and the 99.9th percentile of accumulated rainfall threshold increased, from 405 mm to 475 mm (+ 17.3%), and from 619 mm to 735 mm (+ 18.6%), respectively. For the hourly rain rate (≥0.1mm/h), the 99th and the 99.9th percentile of hourly rain rate increased from 40.9 mm/h to 45.3 mm/h (+ 10.6%), and from 82.2 mm/h to 85.9 mm/h (+ 4.5%), respectively (Fig. 4 d–f). 3.3 Sea waves and storm surges TCs can cause devastating storm surges, which can be further exacerbated by spring tides, leading to catastrophic consequences such as flooding in low-lying areas, infrastructure damage, and loss of life, particularly in areas with poor drainage or inadequate flood protection. Local coastline geometry can also enhance flooding (Lowe et al. 2005, Wang et al. 2008 , Wu et al. 2019 , Lin et al. 2019 ). Wave-surge and wave-current interactions are important processes during storm surges and can be resolved by coupling a 3-dimensional ocean model and a spectral sea wave model, such as COAWST (Olabarrieta et al. 2011 , Benetazzo et al. 2013 ). Grid-level maximum significant wave height and storm surge height difference using PGW minus CTRL in ROMS + SWAN D02 are shown in Fig. 5 . The changes in maximum significant wave height exhibit a complex pattern due to sensitivity to TC tracks (Fig. 5 a/c/e). We found that, by masking out the ocean grid approximating the shoreline (see Figure S5), the maximum significant wave height at the ocean is lifted by 57 cm (+ 13.8%) on average. Notably, there is a large difference between Hato (2017) and Vicente (2012).In the PGW experiment, the surface wind speed on the left-hand side of the TC does not increase for Hato (2017), limiting wave generation. In contrast, more than a 200 cm wave height increase is observed at the open sea near the TC eyewall and along the coastline for Vicente (2012). Storm surge exhibits a relatively uniform increment pattern (Fig. 5 b/d/f). On average, maximum pure surge levels are lifted by 30–80 cm over the open sea but increase much higher along the coastline, especially for narrowing estuaries where the maximum surge level can be lifted by up to 200 cm. To better illustrate the effect on coastal regions, averaged changes for three cases on the oceanic coastal grids (see Figure S5) are shown in Fig. 5 g–h. More than a 200 cm increase in significant wave height is observed over the offshore islands facing the open ocean, while the mainland coastline typically faced a 40–80 cm increase. The lee side of the islands generally did not face significant wave height increases more than 60 cm, and for most inner riverbanks, only a lower than 20 cm wave height increase can be found. Conversely, storm surge level is sensitive to coastline geometry (Fig. 5 h), with the inner riverbanks facing more than a 60 cm increase in surge level. Most coastlines would face a 40–60 cm increase in storm surge, while only the sheltered Hong Kong central waters evaded the large surge level changes, with only a 20–40 cm additional water elevation. According to multi-model mean results from CMIP6, compared to the present-day climatology (1995–2014), global mean sea level rise will be 0.65 m at the end of the 21st century (2081–2100) in SSP5-8.5 (Hermans et al. 2021 ). In the context of the 2090s, a combination of mean sea level rise, storm surge, and wave height can reach more than a 300 cm increase in total water level at certain hot-spot coastlines, without considering the superposition of spring tides. Conclusions and discussion The PGW experiments, carried out by a novel AWO coupled model, mimic the thermodynamic conditions expected in the 2090s. The results suggest that changes in intense tropical cyclones (TCs) could have a significant impact on the PRD region, home to large populations, and the pivot of manufacturing, finance, and marine industries in China. Several coastal disaster indicators, including surface maximum wind speed, accumulated rainfall, hourly rain rate, storm surge height, and storm wave height, become more severe in the PGW experiments as the intensity of the selected TC cases increases. Our study highlights the particular concern of typhoons like Vicente (2012), which have stationary tracks and remain in warm waters for a long time. In a warming future, the warmer and stabler surface ocean, even after a long TC-induced "cold wake", which has an important negative feedback on TC intensity (Bender and Ginis 2000 , Knutson et al. 2001 ), cannot weaken these types of TCs. They tend to experience abrupt and acute rapid intensification nearshore, resulting in a much shorter response time and greater impacts on coastal areas. Our analysis focused on the extreme values in the wind speed distribution, and we found that the maximum 10-m wind speed within 50–200 km to the right of the TC center increased by 4.3 m/s (or 22%). For Vicente (2012), the increase was even greater at 6.9 m/s (or 39%). The increase in the TC right-of-track surface maximum wind speed may be attributed to the faster transition speed of the TCs, while the stronger intensity of these TCs tends to bring wind impacts further inland (Tran et al. 2022 ). Regarding the rainfall impact, we observed an average 11.7% increase in accumulated rainfall, corresponding to a slightly less than 7% per 1°C increase relative to tropical SST warming (the Clausius-Clapeyron relation). It is noteworthy that the percentage of accumulated rainfall by more than 200 mm will increase by 26.4%, and the 99th and 99.9th percentile of accumulated rainfall threshold will increase by + 17.3% and + 18.6%, respectively. The greater than 7% per 1°C increase of the right tail of the rainfall distribution indicates much more severe extreme precipitation caused by TCs, significantly increasing the risks of coastal landslides and flooding in a warming future. Previous studies have shown that, assuming all other factors are unchanged, the most confident TC-related projection is that sea level rise accompanying the warming will lead to higher storm inundation levels (Knutson et al. 2020 ). In this study, we went beyond previous research by considering the surface wave component and its coupling with the 3-dimensional ocean and atmosphere. Simulation results show that even without considering the mean sea level rise and spring tides, stronger TC intensity in a warming climate can already push more than 100 cm of additional water level to most PRD coastlines by combining surge level and wave height changes. Given the significant uncertainties in mean sea level rise, especially for the mass loss in Greenland and Antarctic ice sheets (Hermans et al. 2021 ), from a risk management perspective, we suggest a stress test considering a total water level increase of more than 300 cm in the PRD region due to storm surge in the context of the 2090s, even without considering the overlap with spring tides. Such an approach would help to ensure that the PRD coastal areas and their populations are better prepared for future climate risks and can adapt to the potential impacts of TCs on coastal flooding. It is important to note that our current study has several limitations. We have only examined three intense landfalling TCs since 2010s. The small sampling size would cause potential bias in future projections. Further studies, especially with a much larger sample size, are needed to confirm the current findings and explore other potential impacts of future changes in intense TCs on the PRD coastal regions. Meanwhile, it is important to acknowledge the limitations of our approach in studying the storm surge impact, as we did not use an inundation model to estimate the extent and depth of coastal flooding during storm surge events. This limitation may have restricted our ability to fully quantify the potential impacts of climate change on storm surges in the PRD region. Many other important parameters, e.g. the potential shift of the TC tracks in the warming climate and potential rapid intensification or rapid weakening of TCs over the study region are needed to investigated by ensemble simulation technique. Ensemble technique also support robust analysis based on impact indices such as the power dissipation index (PDI, Emanuel 2005 ) and cyclone damage potential (CDP, Holland et al. 2019 ), which can be further incorporated to study the TC changes. Declarations Funding We appreciate the assistance of the Hong Kong Observatory (HKO) and Environmental Central Facility (HKUST), which provided the meteorological and model input data. This work is supported by grants from the Research Grants Council of the Hong Kong Special Administrative Region, China (Project No. AoE/E-603/18 and T31-603/21-N), the Environment and Conservation Fund (2022-02), the CORE Project 2022/23, and the CORE is a joint research center for ocean research between QNLM and HKUST. Financial and non-financial interests The authors have no relevant financial or non-financial interests to disclose. References Balaguru, K., Foltz, G.R., Leung, L.R. and Emanuel, K.A., 2016. Global warming-induced upper-ocean freshening and the intensification of super typhoons. Nature Communications, 7(1), p.13670. Bender, M.A. and Ginis, I., 2000. Real-case simulations of hurricane–ocean interaction using a high-resolution coupled model: Effects on hurricane intensity. Monthly Weather Review, 128(4), pp.917-946. Benetazzo, A., Carniel, S., Sclavo, M. and Bergamasco, A., 2013. Wave–current interaction: Effect on the wave field in a semi-enclosed basin. Ocean Modelling, 70, pp.152-165. Bloemendaal, N., de Moel, H., Martinez, A.B., Muis, S., Haigh, I.D., van der Wiel, K., Haarsma, R.J., Ward, P.J., Roberts, M.J., Dullaart, J.C. and Aerts, J.C., 2022. A globally consistent local-scale assessment of future tropical cyclone risk. Science Advances, 8(17), p.eabm8438. Chan JCL (2008) Decadal variations of intense typhoon occurrence in the western North Pacific. Proc R Soc Lond A 464:249–272. doihttps://doi.org/10.1098/rspa.2007.0183 Chang, Y.P., Yang, S.C., Lin, K.J., Lien, G.Y. and Wu, C.M., 2020. Impact of tropical cyclone initialization on its convection development and intensity: A case study of Typhoon Megi (2010). Journal of the Atmospheric Sciences, 77(2), pp.443-464. Chen, J., Wang, Z., Tam, C.Y., Lau, N.C., Lau, D.S.D. and Mok, H.Y., 2020. Impacts of climate change on tropical cyclones and induced storm surges in the Pearl River Delta region using pseudo-global-warming method. Scientific Reports, 10(1), pp.1-10. Chen, S. S., Zhao, W., Donelan, M. A., & Tolman, H. L. 2013. Directional wind-wave coupling in fully coupled atmosphere-wave-ocean models: Results from CBLAST-Hurricane. Journal of the Atmospheric Sciences, 70(10), 3198–3215. Chen, X., Wang, Y., Zhao, K. and Wu, D., 2017. A numerical study on rapid intensification of Typhoon Vicente (2012) in the South China Sea. Part I: Verification of simulation, storm-scale evolution, and environmental contribution. Monthly Weather Review, 145(3), pp.877-898. Choy, C.W., of Insurers, T.H.K.F., Wu, M.C. and Lee, T.C., 2020. Assessment of the damages and direct economic loss in Hong Kong due to Super Typhoon Mangkhut in 2018. Tropical Cyclone Research and Review, 9(4), pp.193-205. Chu, J.E., Lee, S.S., Timmermann, A., Wengel, C., Stuecker, M.F. and Yamaguchi, R., 2020. Reduced tropical cyclone densities and ocean effects due to anthropogenic greenhouse warming. Science Advances, 6(51), p.eabd5109. Cummings, J.A. and Smedstad, O.M., 2013. Variational data assimilation for the global ocean. In Data Assimilation for Atmospheric, Oceanic and Hydrologic Applications (Vol. II) (pp. 303-343). Springer, Berlin, Heidelberg. Egbert, G.D. and Erofeeva, S.Y., 2002. Efficient inverse modeling of barotropic ocean tides. Journal of Atmospheric and Oceanic Technology, 19(2), pp.183-204. Emanuel, K.A., 1986. An air-sea interaction theory for tropical cyclones. Part I: Steady-state maintenance. Journal of Atmospheric Sciences, 43(6), pp.585-605. Emanuel, K. A., 2005: Increasing destructiveness of tropical cyclones over the past 30 years. Nature, 436, 686–688. https://doi.org/10.1038/nature03906 Gutjahr, O., Putrasahan, D., Lohmann, K., Jungclaus, J.H., von Storch, J.S., Brüggemann, N., Haak, H. and Stössel, A., 2019. Max Planck Institute Earth System Model (MPI-ESM1. 2) for the high-resolution model intercomparison project (HighResMIP). Geoscientific Model Development, 12(7), pp.3241-3281. Hermans, T.H., Gregory, J.M., Palmer, M.D., Ringer, M.A., Katsman, C.A. and Slangen, A.B., 2021. Projecting global mean sea‐level change using CMIP6 models. Geophysical Research Letters, 48(5), p.e2020GL092064. Holland, G.J., Done, J.M., Douglas, R., Saville, G.R. and Ge, M., 2019. Global tropical cyclone damage potential. Hurricane risk, pp.23-42. Huang, P., Lin, I.I., Chou, C. and Huang, R.H., 2015. Change in ocean subsurface environment to suppress tropical cyclone intensification under global warming. Nature Communications, 6(1), p.7188. Jyoteeshkumar Reddy, P., Sriram, D., Gunthe, S.S. and Balaji, C., 2021. Impact of climate change on intense Bay of Bengal tropical cyclones of the post-monsoon season: a pseudo global warming approach. Climate Dynamics, 56, pp.2855-2879. Knutson, T., Camargo, S.J., Chan, J.C., Emanuel, K., Ho, C.H., Kossin, J., Mohapatra, M., Satoh, M., Sugi, M., Walsh, K. and Wu, L., 2019. Tropical cyclones and climate change assessment: Part I: Detection and attribution. Bulletin of the American Meteorological Society, 100(10), pp.1987-2007. Knutson, T., Camargo, S.J., Chan, J.C., Emanuel, K., Ho, C.H., Kossin, J., Mohapatra, M., Satoh, M., Sugi, M., Walsh, K. and Wu, L., 2020. Tropical cyclones and climate change assessment: Part II: Projected response to anthropogenic warming. Bulletin of the American Meteorological Society, 101(3), pp.E303-E322. Knutson, T.R., McBride, J.L., Chan, J., Emanuel, K., Holland, G., Landsea, C., Held, I., Kossin, J.P., Srivastava, A.K. and Sugi, M., 2010. Tropical cyclones and climate change. Nature Geoscience, 3(3), pp.157-163. Knutson, T.R., Tuleya, R.E., Shen, W. and Ginis, I., 2001. Impact of CO2-induced warming on hurricane intensities as simulated in a hurricane model with ocean coupling. Journal of Climate, 14(11), pp.2458-2468. Li, Z., Tam, C.Y., Li, Y., Lau, N.C., Chen, J., Chan, S.T., Dickson Lau, D.S. and Huang, Y., 2022. How Does Air‐Sea Wave Interaction Affect Tropical Cyclone Intensity? An Atmosphere‐Wave‐Ocean Coupled Model Study Based on Super Typhoon Mangkhut (2018). Earth and Space Science, 9(3), p.e2021EA002136. Lin, N., Marsooli, R. and Colle, B.A., 2019. Storm surge return levels induced by mid-to-late-twenty-first-century extratropical cyclones in the Northeastern United States. Climatic Change, 154, pp.143-158. Lin, S. and Sheng, J., 2020. Revisiting dependences of the drag coefficient at the sea surface on wind speed and sea state. Continental Shelf Research, 207, p.104188. Liu, K.S. and Chan, J.C., 2020. Recent increase in extreme intensity of tropical cyclones making landfall in South China. Climate Dynamics, 55(5-6), pp.1059-1074. Lowe, J.A. and Gregory, J.M., 2005. The effects of climate change on storm surges around the United Kingdom. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 363(1831), pp.1313-1328. Magnusson, L., Bidlot, J. R., Bonavita, M., Brown, A. R., Browne, P. A., De Chiara, G., et al. (2019). ECMWF activities for improved hurricane forecasts. Bulletin of the American Meteorological Society, 100(3), 445–458. Mogensen, K. S., Magnusson, L., & Bidlot, J. R. (2017). Tropical cyclone sensitivity to ocean coupling in the ECMWF coupled model. Journal of Geophysical Research: Oceans, 122, 4392–4412. https://doi.org/10.1002/2017JC012753 Moon, I. J., Ginis, I., & Hara, T. (2004). Effect of surface waves on Charnock coefficient under tropical cyclones. Geophysical Research Letters, 31, L20302. https://doi.org/10.1029/2004GL020988 Moon, J., Cha, D.H., Lee, M. and Kim, J., 2018. Impact of spectral nudging on real‐time tropical cyclone forecast. Journal of Geophysical Research: Atmospheres, 123(22), pp.12-647. Nakanishi, M. and Niino, H., 2004. An improved Mellor–Yamada level-3 model with condensation physics: Its design and verification. Boundary-layer Meteorology, 112, pp.1-31. Olabarrieta, M., Warner, J.C. and Kumar, N., 2011. Wave‐current interaction in Willapa Bay. Journal of Geophysical Research: Oceans, 116(C12). Olabarrieta, M., Warner, J.C., Armstrong, B., Zambon, J.B. and He, R. (2012) Ocean– atmosphere dynamics during Hurricane Ida and Nor’Ida: An application of the coupled ocean–atmosphere–wave–sediment transport (COAWST) modeling system. Ocean Modelling, 43, pp.112-137. Taylor, P. K., & Yelland, M. J. (2001). The dependence of sea surface roughness on the height and steepness of the waves. Journal of Physical Oceanography, 31(2), 572–590. Tong, M., Sippel, J.A., Tallapragada, V., Liu, E., Kieu, C., Kwon, I.H., Wang, W., Liu, Q., Ling, Y. and Zhang, B., 2018. Impact of assimilating aircraft reconnaissance observations on tropical cyclone initialization and prediction using operational HWRF and GSI ensemble–variational hybrid data assimilation. Monthly Weather Review, 146(12), pp.4155-4177. Tran, T.L., Ritchie, E.A., Perkins‐Kirkpatrick, S.E., Bui, H. and Luong, T.M., 2022. Future Changes in Tropical Cyclone Exposure and Impacts in Southeast Asia From CMIP6 Pseudo‐Global Warming Simulations. Earth's Future, 10(12), p.e2022EF003118. Tsuboki, K., Yoshioka, M.K., Shinoda, T., Kato, M., Kanada, S. and Kitoh, A., 2015. Future increase of supertyphoon intensity associated with climate change. Geophysical Research Letters, 42(2), pp.646-652. Wang, S., McGrath, R., Hanafin, J., Lynch, P., Semmler, T. and Nolan, P., 2008. The impact of climate change on storm surges over Irish waters. Ocean Modelling, 25(1-2), pp.83-94. Warner, J.C., Armstrong, B., He, R. and Zambon, J.B. (2010) Development of a coupled ocean–atmosphere–wave–sediment transport (COAWST) modeling system. Ocean Modelling, 35(3), pp.230-244. Wu, C.C., Chou, K.H., Wang, Y. and Kuo, Y.H., 2006. Tropical cyclone initialization and prediction based on four-dimensional variational data assimilation. Journal of the Atmospheric Sciences, 63(9), pp.2383-2395. Wu, C.C., Li, T.H. and Huang, Y.H., 2015. Influence of mesoscale topography on tropical cyclone tracks: Further examination of the channeling effect. Journal of the Atmospheric Sciences, 72(8), pp.3032-3050. Wu, J., Chen, Y., Liao, Z., Gao, X., Zhai, P. and Hu, Y., 2022. Increasing risk from landfalling tropical cyclone-heatwave compound events to coastal and inland China. Environmental Research Letters, 17(10), p.105007. Wu, L., Breivik, Ø. and Rutgersson, A., 2019. Ocean‐wave‐atmosphere interaction processes in a fully coupled modeling system. Journal of Advances in Modeling Earth Systems, 11(11), pp.3852-3874. Wu, W. and Leonard, M., 2019. Impact of ENSO on dependence between extreme rainfall and storm surge. Environmental Research Letters, 14(12), p.124043. Yin, K., Xu, S., Zhao, Q., Huang, W., Yang, K. and Guo, M., 2020. Effects of land cover change on atmospheric and storm surge modeling during typhoon event. Ocean Engineering, 199, p.106971. Zambon, J.B., He, R. and Warner, J.C. (2014a) Investigation of hurricane Ivan using the coupled ocean–atmosphere–wave–sediment transport (COAWST) model. Ocean Dynamics, 64(11), pp.1535-1554. Zambon, J.B., He, R. and Warner, J.C. (2014b) Tropical to extratropical: Marine environmental changes associated with Superstorm Sandy prior to its landfall. Geophysical Research Letters, 41(24), pp.8935-8943. Zhang, F., Weng, Y., Gamache, J. F., & Marks, F. D. (2011). Performance of convection-permitting hurricane initialization and prediction during 2008–2010 with ensemble data assimilation of inner-core airborne Doppler radar observations. Geophysical Research Letters, 38, L15810. Supplementary Files SupplementaryInfo.docx Cite Share Download PDF Status: Published Journal Publication published 09 Mar, 2024 Read the published version in Natural Hazards → Version 1 posted Reviewers agreed at journal 26 Nov, 2023 Reviewers invited by journal 24 Nov, 2023 Editor invited by journal 18 Nov, 2023 Editor assigned by journal 16 Nov, 2023 First submitted to journal 15 Nov, 2023 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-3625663","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":253064806,"identity":"9e3920cf-906e-4e70-9929-c58bea96ddfb","order_by":0,"name":"Zhenning Li","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA+UlEQVRIiWNgGAWjYDACCRBhwMDAx8x84MAHIJuNnVgtbMxsiQdngLQwE6UFpJKBx/gwD4hFSIv87OZjjwsK7ti1sfMYHLb5tU2ej5mB8cPHHNxaGOccSzeeYfAsuY2ZreBwbt9twzZmBmbJmdtwa2GWyDGTBpqfzMbMvOFwbs9tRqAWNmZePFrYJPK/QbUwGBy27LltT1ALj0QOG0iLHRszi8Fhhh+3EwlqkZBIAzssARjICQd7G24DPcXYjNcv8jOSn0nz/Dlsz89/+PCHH39u285vbz744SMeLTCQ2AAiGdvAZANh9UBgD6H+EKV4FIyCUTAKRhgAAHfISW6ecmikAAAAAElFTkSuQmCC","orcid":"https://orcid.org/0000-0002-7770-6867","institution":"The Hong Kong University of Science and Technology","correspondingAuthor":true,"submittingAuthor":false,"prefix":"","firstName":"Zhenning","middleName":"","lastName":"Li","suffix":""},{"id":253064807,"identity":"2e706c0e-c11a-4107-9ef4-929d2769c55a","order_by":1,"name":"Jimmy C.H. Fung","email":"","orcid":"","institution":"The Hong Kong University of Science and Technology","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Jimmy","middleName":"C.H.","lastName":"Fung","suffix":""},{"id":253064808,"identity":"fdc921da-8033-4b71-b22a-9e3ad2be32ab","order_by":2,"name":"Mau Fung Wong","email":"","orcid":"","institution":"The Hong Kong University of Science and Technology","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Mau","middleName":"Fung","lastName":"Wong","suffix":""},{"id":253064809,"identity":"afe04a60-eb46-428e-af00-b9d28841a76d","order_by":3,"name":"Shangfei Lin","email":"","orcid":"","institution":"The Hong Kong University of Science and Technology","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Shangfei","middleName":"","lastName":"Lin","suffix":""},{"id":253064810,"identity":"dd14d981-3980-4455-9bcb-de2bff23b87a","order_by":4,"name":"Fenying Cai","email":"","orcid":"","institution":"Potsdam Institute for Climate Impact Research (PIK) e V: Potsdam-Institut fur Klimafolgenforschung (PIK) e V","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Fenying","middleName":"","lastName":"Cai","suffix":""},{"id":253064811,"identity":"9868f2a6-1a3b-464e-b23f-49d2276df2dd","order_by":5,"name":"Wenfeng Lai","email":"","orcid":"","institution":"The Hong Kong University of Science and Technology","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Wenfeng","middleName":"","lastName":"Lai","suffix":""},{"id":253064812,"identity":"eaa0f236-0609-43dd-9065-943f6e74bc9d","order_by":6,"name":"Alexis K.H. Lau","email":"","orcid":"","institution":"The Hong Kong University of Science and Technology","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Alexis","middleName":"K.H.","lastName":"Lau","suffix":""}],"badges":[],"createdAt":"2023-11-17 12:40:56","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-3625663/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-3625663/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1007/s11069-024-06510-7","type":"published","date":"2024-03-09T15:01:49+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":47307624,"identity":"ef653c05-6336-4b80-b4bc-f769f149f648","added_by":"auto","created_at":"2023-11-29 16:10:42","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":1506522,"visible":true,"origin":"","legend":"\u003cp\u003eDomain configurations for the COAWST. (a) Tier 1 consists of the outermost 9-km WRF D01, with the intermediate 3-km WRF D02 enclosed in the red box. The blue dashed box represents the 2.2-km ROMS/SWAN D01. Tier 2 includes (b) the innermost WRF 1-km D02 that covers the entire PRD region, and (c) the innermost ROMS/SWAN D02, with a horizontal resolution of 500 m, covering the Pearl River estuary. The terrain height and bathymetry are shown using colored shading.\u003c/p\u003e","description":"","filename":"1.png","url":"https://assets-eu.researchsquare.com/files/rs-3625663/v1/c312df2023796a8c82e7e4ce.png"},{"id":47307625,"identity":"8d8d0701-0ce7-4d31-9951-c9a212d1d573","added_by":"auto","created_at":"2023-11-29 16:10:42","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":211013,"visible":true,"origin":"","legend":"\u003cp\u003e(a) Tracks of Mangkhut (2018) represented by the black solid curve for CMA best Track data, the blue dashed curve for the CTRL experiment, and the red dashed curve for the PGW experiment. (d) same as (a) but for the TC’s minimum sea level pressure (hPa) evolution. (g) is same as (d) but for the maximum surface wind speed (m/s). (b) and (c), (e) and (f), and (h) and (i) show the same for Hato (2017) and Vicente (2012), respectively.\u003c/p\u003e","description":"","filename":"2.png","url":"https://assets-eu.researchsquare.com/files/rs-3625663/v1/4a3499cb90573ebfe83a71fe.png"},{"id":47307626,"identity":"1363d871-57c7-4357-a0f3-1ce957cf700e","added_by":"auto","created_at":"2023-11-29 16:10:42","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":409356,"visible":true,"origin":"","legend":"\u003cp\u003e(a) Grid-wise maximum 10-m wind speed difference (in m/s) during the simulation period between PGW and CTRL simulations for Mangkhut (2018) in WRF d02. Negative values are masked out for clarity. (b–c) Same as (a), but for Hato (2017) and Vicente (2012), respectively.\u003c/p\u003e","description":"","filename":"3.png","url":"https://assets-eu.researchsquare.com/files/rs-3625663/v1/4c04dc14c8a11e3f03863887.png"},{"id":47307623,"identity":"f09dfada-b072-482b-a55e-b2447ad63ba0","added_by":"auto","created_at":"2023-11-29 16:10:42","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":137792,"visible":true,"origin":"","legend":"\u003cp\u003eHistograms of accumulated rainfall (in mm) larger than 50 mm within a radius of 500 km from the TC center in WRF D02. The blue curve represents the CTRL experiment, while the red dashed curve represents the PGW experiment. (a) Mangkhut (2018), (b) Hato (2017), and (c) Vicente (2012) for accumulated rainfall. (d) Mangkhut (2018), (e) Hato (2017), and (f) Vicente (2012) for hourly precipitation larger than 5 mm/hr.\u003c/p\u003e","description":"","filename":"4.png","url":"https://assets-eu.researchsquare.com/files/rs-3625663/v1/25575aebb22d01ea37d74734.png"},{"id":47307628,"identity":"9015515d-2e3f-46ee-b9f8-d8dc2b3af219","added_by":"auto","created_at":"2023-11-29 16:10:43","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":271925,"visible":true,"origin":"","legend":"\u003cp\u003e(a) Grid-level maximum significant wave height difference (cm) and (b) storm surge level difference (cm) using PGW minus CTRL for Mangkhut (2018) in ROMS+SWAN d02, (c-d) and (e-f), same as (a-b) but for Hato (2017) and Vicente (2012), respectively. (g) and (h), three-case averaged coastal grid maximum significant wave height difference (cm) and storm surge level difference (cm).\u003c/p\u003e","description":"","filename":"5.png","url":"https://assets-eu.researchsquare.com/files/rs-3625663/v1/dfd636329ecebb4ae71c7608.png"},{"id":52431990,"identity":"e57e5dd9-04cd-4adc-9a4a-050321378217","added_by":"auto","created_at":"2024-03-11 15:10:14","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":2205730,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-3625663/v1/3d068c6d-4906-40bb-8296-ee50b5b3a399.pdf"},{"id":47308882,"identity":"d25b9848-87e7-4526-a516-8bafd3ded153","added_by":"auto","created_at":"2023-11-29 16:18:43","extension":"docx","order_by":4,"title":"","display":"","copyAsset":false,"role":"supplement","size":2773545,"visible":true,"origin":"","legend":"","description":"","filename":"SupplementaryInfo.docx","url":"https://assets-eu.researchsquare.com/files/rs-3625663/v1/d7ccf82937d2d6b34416e18a.docx"}],"financialInterests":"","formattedTitle":"Future changes in intense tropical cyclone hazards in the Pearl River Delta region: An air-wave-ocean coupled model study","fulltext":[{"header":"Introduction","content":"\u003cp\u003eLandfalling tropical cyclones (TCs) are intense weather systems that can cause a range of coastal hazards, including storm surges and waves, coastal erosion, wind induced damage, rainstorms, and secondary hazards like landslides and inland flooding (Tran et al. \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2022\u003c/span\u003e, Wu et al. \u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). The impact of landfalling TCs can vary depending on the severity and location of the storm, as well as the vulnerability of the affected communities (Wang et al. \u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e2008\u003c/span\u003e, Knutson et al. \u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e2010\u003c/span\u003e, Knutson et al. \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e2019\u003c/span\u003e, Chu et al. \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). The Pearl River Delta (PRD) region, with its long and meandering coastline and well-developed economy, is highly vulnerable to TC-caused coastal disasters (Chan \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e2008\u003c/span\u003e, Chen et al. \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). Recent TCs, such as Typhoon Hato (2017) and Mangkhut (2018), have demonstrated the devastating impact these storms can have, causing widespread damage and losses (Liu and Chan \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). The estimated direct economic loss due to Mangkhut (2018) and Hato (2017) in Hong Kong is about HK\u003cspan\u003e$\u003c/span\u003e4.60\u0026nbsp;billion (~\u0026thinsp;US\u003cspan\u003e$\u003c/span\u003e 590\u0026nbsp;million) and HK\u003cspan\u003e$\u003c/span\u003e1.21\u0026nbsp;billion (~\u0026thinsp;US\u003cspan\u003e$\u003c/span\u003e 155\u0026nbsp;million), respectively (Choy et al. \u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). Meanwhile, the impacts of climate change on TCs have the potential to significantly alter their influence on future coastal residents (Knutson et al. \u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e2010\u003c/span\u003e, Tsuboki et al. \u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e2015\u003c/span\u003e, Tran et al. \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). The global warming trend is expected to continue or even worsen in the rest of the 21st century (Knutson et al. \u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e2010\u003c/span\u003e). A recent comprehensive assessment by Knutson et al (\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e2020\u003c/span\u003e) shows that a 2\u0026deg;C anthropogenic global warming is projected to impact TC activity as follows: 1) Medium-to-high confidence in the increased precipitation rates by around +\u0026thinsp;14%, or close to the rate of tropical water vapor increase with warming; 2) Medium-to-high confidence in the +\u0026thinsp;5% stronger lifetime maximum surface wind speeds; 3) Medium-to-high confidence in a\u0026thinsp;+\u0026thinsp;13% increase of the proportion of TCs that reach very intense (category 4\u0026ndash;5) levels; and 4) The most confident TC-related projection is that sea level rise accompanying the warming will lead to higher storm inundation levels, assuming all other factors are unchanged (Knutson et al. \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). However, research conclusions are more mixed, and confidence levels are lower for the basin-scale or regional projections (Knutson et al. \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e2020\u003c/span\u003e, Bloemendaal et al. \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Due to the limitation in spatial resolution, Global Climate Model (GCM) studies also introduce large biases when reproducing TC activities on a regional scale (Bloemendaal et al. \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). These factors fall far short of satisfying the demand to guide regional climate change adaptation and resilience.\u003c/p\u003e \u003cp\u003eOn a regional scale, air-sea interaction plays a crucial role in the development of TCs. As noted by Emanuel (\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e1986\u003c/span\u003e), the ocean surface enthalpy fluxes transferred to the atmosphere are the main energy sources of TCs. The Coupled Model Intercomparison Project Phase 6 (CMIP6) model results indicate that in the future, there will be a warmer and fresher surface ocean. This implies an increased thermal stratification in a warming climate, which in turn will feedback on the TC development (Huang et al. \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e2015\u003c/span\u003e, Balaguru et al. \u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). Additionally, studies have shown that at extremely high surface wind speeds (30 m/s or higher), the roughness length, determined by the sea state, tends to plateau. This can maintain or even enhance the inner-core intensity of the TC (Moon et al. \u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e2004\u003c/span\u003e, Chen et al. \u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e2013\u003c/span\u003e). Numerical studies incorporating such effects can construct more realistic TC structures than atmosphere-only or even atmosphere\u0026ndash;ocean coupled models (Davis et al. 2008, Chen et al. \u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e2013\u003c/span\u003e, Wu et al. \u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). Consequently, operational centers such as the European Centre for Medium Range Weather Forecasts (ECMWF) are now using atmosphere-wave-ocean (AWO) coupled models to generate routine global products for better performance in TC forecasts (Mogensen et al. \u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e2017\u003c/span\u003e, Magnusson et al. \u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). Incorporating a dynamically evolving 3-dimensional ocean model and a wave model at convection-permitting scales (i.e., 4 km or finer; see Zhang et al. \u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e2011\u003c/span\u003e) is crucial to better capture the TC-related coastal disasters on a regional scale. By generating crucial products related to coastal hazards across the atmosphere (wind and rainfall), sea wave (significant wave height), and ocean (storm surge level) from the synergistic framework, valuable supporting evidence can be provided to guide regional adaptation and resilience against TCs in a changing climate.\u003c/p\u003e \u003cp\u003eThe objective of this study is to investigate the impact of future TC-related coastal disasters on the PRD region in a changing climate, using a sophisticated AWO coupled model. The PRD region is home to numerous large populations and is a pivot of manufacturing, finance, and marine industries such as shipping and aquaculture in China. To achieve this objective, we examine three intense landfalling TCs since 2010 in \u0026ldquo;parallel digital worlds\u0026rdquo;: a present-day world and a warming future world. Specifically, we conduct simulations using (a) present-day reanalysis data and (b) the same TCs occurring in a pseudo-global warming (PGW) scenario projected for the end of the 21st century. Conventional PGW experiments have been extensively used in atmospheric-only regional models such as the WRF (Chen et al. \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2020\u003c/span\u003e, Jyoteeshkumar et al. 2021, Tran et al. \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). To our knowledge, this study is the first to implement PGW in an AWO fully coupled system over the PRD region, achieving such high horizontal resolutions (1-km atmosphere and 500-m ocean and sea wave). The perturbated warming initial and boundary conditions are derived from both atmosphere and ocean GCMs projected by CMIP6 and enforced in both the atmosphere and the ocean components of the fully coupled system. The results of this study will provide valuable insights into the potential impacts of future TC-related coastal disasters on the PRD region and guide regional adaptation and resilience against such disasters in the changing climate.\u003c/p\u003e"},{"header":"Modeling system and experimental design","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003e2.1 The Coupled Ocean\u0026ndash;Atmosphere\u0026ndash;Wave\u0026ndash;Sediment Transport (COAWST) modeling system\u003c/h2\u003e \u003cp\u003eThe COAWST system serves as the fundamental framework for this study. It is an agglomeration of open-source modeling components tailored to investigate coupled processes of the atmosphere, ocean, and waves in coastal regions. Since its release by the United States Geological Survey in 2010, it has been widely used in various studies (Warner et al., \u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e2010\u003c/span\u003e). COAWST has been applied to examine several Atlantic TCs that caused considerable damage in North America, such as Ivan (2004), Ida (2009), and Sandy (2012) (Olabarrieta et al. \u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e2012\u003c/span\u003e,Zambon et al. \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2014a\u003c/span\u003e, Zambon et al. \u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e2014b\u003c/span\u003e). In this study, we use COAWST v3.6, which comprises the following modules: (1) an atmospheric component named Weather Research and Forecasting Model (WRF v4.1.2); (2) a Regional Ocean Modeling System (ROMS) v3.7 as an oceanic component; (3) a SWAN v41.31 as a sea wave component; and (4) a Model Coupling Toolkit (MCT) v2.6.0 as a coupler.\u003c/p\u003e \u003cp\u003eIn previous work by Li et al (\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2022\u003c/span\u003e), the COAWST has been customized and tuned extensively to suit local weather and climate studies in the PRD region. It was used to investigate the impact of sea wave-wind interaction on TC development, including typhoon tracks and intensities. The current system builds upon this framework and is organized into a hierarchical architecture with two tiers, each tailored to simulate dominant physical processes at corresponding spatial scales (see Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). Tier 1 is designed for open sea scale modeling. The WRF model is configured with a horizontal resolution of 9 km \u0026times; 9 km in its outermost domain (WRF D01). The ROMS and SWAN models (ROMS\u0026thinsp;+\u0026thinsp;SWAN D01) are nested within WRF D01, with an ocean eddy-resolving resolution of 2.2 km \u0026times; 2.2 km. The intermediate domain for WRF (WRF D02) is embedded in ROMS\u0026thinsp;+\u0026thinsp;SWAN D01 at sea, with a convection-permitting resolution of 3 km \u0026times; 3 km. Tier 2 is designed to provide high spatial resolution for modeling the PRD coastal region. The innermost WRF domain (WRF D03) has a resolution of 1 km \u0026times; 1 km and covers the entire PRD region. The innermost ROMS\u0026thinsp;+\u0026thinsp;SWAN domain (ROMS\u0026thinsp;+\u0026thinsp;SWAN D02) has a resolution of 500 m \u0026times; 500 m and covers the Pearl River estuary and the nearshore sea.\u003c/p\u003e \u003cp\u003eThe configurations for each of the three component models follow Li et al. (\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2022\u003c/span\u003e), with minor modifications made for computational efficiency. These modifications include a smaller mesh size in WRF D01/02 and ROMS\u0026thinsp;+\u0026thinsp;SWAN D01. In addition, the number of vertical layers in the WRF model is reduced to 39 layers, but the number of vertical layers within 1,500 m above the surface is increased to 19 for improved accuracy in representing physical processes and structures in the planetary boundary layer. For air-wave interaction, we implemented Lin and Sheng (\u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e2020\u003c/span\u003e) scheme, which is embedded in the Mellor-Yamada-Nakanishi-Niino (MYNN; Nakanishi and Niino \u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e2004\u003c/span\u003e) planetary boundary layer scheme and associated surface layer schemes. Detailed configurations are summarized in Table \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e. Tide information was obtained from the Oregon State University Tidal Prediction Software (OTPS; Egbertand and Erofeeva 2002), which provides data on eight primary harmonic constituents (M2, S2, N2, K2, K1, O1, P1, Q1).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec4\" class=\"Section2\"\u003e \u003ch2\u003e2.2 Experimental design\u003c/h2\u003e \u003cp\u003eIn this study, we selected three intense landfalling TCs - Mangkhut (2018), Hato (2017), and Vicente (2012) - which caused significant impacts on the PRD coastal areas. These selected TCs were then simulated using the COAWST framework. Two groups of experiments were conducted: present-day control (CTRL) experiments and pseudo-global warming (PGW) experiments. The CTRL experiments aimed to reproduce the selected TCs and related coastal disasters with high fidelity, including rainfall, wind, storm surges, and storm waves. The PGW experiments aimed to simulate the same storms occurring in a global warming scenario projected for the end of the 21st century.\u003c/p\u003e \u003cp\u003eWe utilized the high-resolution version of the Max Planck Institute Earth System Model (MPI-ESM-1-2-HR) from CMIP6 datasets to extract PGW signals (Gutjahr et al. \u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). Thermo-variables, including air temperature, humidity, sea water temperature, and sea water salinity, were extracted from 2090\u0026ndash;2099 in the Shared Socioeconomic Pathways 5-8.5 (SSP5-8.5) scenario and 2005\u0026ndash;2014 in the historical simulation. The differences between the two epochs were then added to the initial and boundary conditions of the WRF and ROMS models. Figure S2 and figure S3 show the changes in July-August-September 2090s for the SSP5-8.5 scenario, including a projected 2.4\u0026deg;C warming of Sea Surface Temperature (SST) in the outermost domain and a 2.6\u0026deg;C warming of the 25\u0026ndash;50 m subsurface ocean. Specific humidity near the surface level is expected to increase by 3.0\u0026ndash;4.0 g/kg, and the surface air temperature will increase by 3.2\u0026deg;C, with more than 6.0\u0026deg;C warming in the upper level (near 250 hPa, see figure S3). We compared the results with the multi-model ensemble mean and a previous study by Chen et al. (\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2020\u003c/span\u003e) based on CMIP5 output and found similar changes in variable magnitudes and large-scale patterns. To our knowledge, this is the first time the PGW method has been applied in a regional AWO coupled model to study future TC activities.\u003c/p\u003e \u003c/div\u003e"},{"header":"Results","content":"\u003cdiv id=\"Sec6\" class=\"Section2\"\u003e \u003ch2\u003e3.1 Simulation basics and model evaluation\u003c/h2\u003e \u003cp\u003eWe begin by comparing the tracks and intensities of the selected TCs in both the CTRL and PGW experiments. The successful reproduction of TC evolutions in the CTRL experiments provides a higher level of confidence for future projections. As depicted in Figs.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ea\u0026ndash;c, the COAWST model performed quite well in capturing the tracks of all selected TCs, including the stable steering flow caused northwest-southeast tracks of Mangkhut (2018) and Hato (2017) from the western Philippine waters to the southern China coast. Notably, the model effectively simulated the subtle but essential episode of Vicente (2012), where it slowed down and remained stagnant over the northern South China Sea on July 22, 2012, as the subtropical ridge northeast of the storm suddenly weakened (Chen et al. \u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). While the track for TC Vicente (2012) exhibited an evident southward track shift compared with the observation and future run near 111\u0026deg;E, 22\u0026deg;N. This southward track deflection is a typical behavior regarding the terrain effect over the western Guangdong Province, similar effect in Taiwan has been well archived by Wu et al. (\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e2015\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eFigures \u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ed\u0026ndash;f and \u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eg\u0026ndash;i display the central minimum sea level pressures (SLPs) and maximum 10-m wind speeds, respectively. It is noteworthy that the COAWST underestimates the central SLPs in all cases, which is a known bias of the WRF model, even in cloud-resolving resolution without vortex initialization (Wu et al. \u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e2006\u003c/span\u003e, Tong et al. \u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e2018\u003c/span\u003e, Chang et al. \u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e2020\u003c/span\u003e, Chen et al. \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). To address this issue, we tested the insertion of the Rankine vortex to the initial condition based on the real TC intensity and location through the WRF TC \u0026ldquo;bogus\u0026rdquo; scheme in Hato (2017). However, this technique resulted in a remarkable overestimation of surface wind and unreasonably high storm surge levels when compared with in-situ observations (discussed later). As a trade-off, experiments were conducted without \u0026ldquo;bogus\u0026rdquo; TC to better simulate more direct coastal disaster indicators, despite the underestimation of central SLPs.\u003c/p\u003e \u003cp\u003eWe compared the changes in tracks and intensities between the PGW and CTRL experiments. To ensure the constancy of ensuing analysis, we used large-scale spectral nudging of steering wind with appropriate relaxation coefficients, and all TCs in the PGW experiments kept similar tracks to those in CTRL. Additionally, since spectral nudging tends to suppress the TC intensity (Moon et al. \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e2018\u003c/span\u003e), we conducted experiments without nudging, and the results indicated that the tracks could deflect by hundreds of kilometers, although the minimum SLP deepened by around 5 hPa (Figure S4).\u003c/p\u003e \u003cp\u003eFor intensity comparison between the PGW and CTRL, both Mangkhut (2018) and Hato (2017) exhibited a slight drop in the minimum SLP by 3\u0026ndash;7 hPa, which is similar to recent studies on TCs impacting Southeast Asia (Chen et al. \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2020\u003c/span\u003e, Jyoteeshkumar et al. 2021, Tran et al. \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). However, for Vicente (2012), we found a sharp drop in the central SLP in the PGW experiment, which was 15 hPa lower than that in CTRL. The unusual stationary track over the northern South China Sea allowed Vicente (2012) to absorb excess energy from the warmer and stabler ocean in a warming future. In the CTRL simulation, the mean SST beneath the stationary location of Vicente (2012) dropped from around 29.3\u0026deg;C to 27.5\u0026deg;C, whereas in PGW, it dropped from around 32.1\u0026deg;C to 30.0\u0026deg;C. Even after strong upper-ocean mixing by Vicente (2012), the SST is still warmer than the initial SST over that location in CTRL (Figure not shown). Regarding the maximum surface wind speed, the time series show similar evolution, but intense peak wind strengthening can be found in the PGW. Mangkhut (2018) and Hato (2017) exhibited a peak surface wind speed that was approximately 4 m/s (~\u0026thinsp;10%) higher, which is consistent with the results from Chen et al. (\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2020\u003c/span\u003e) study on three different TCs affecting the PRD in the far future (2075\u0026ndash;2099). Particularly, Vicente (2012) even attained nearly 10 m/s (~\u0026thinsp;30%) increase in maximum 10-meter wind speed due to further intensification.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eIn addition to evaluating performance in tracks and intensities, evaluating the performance of the COAWST using in-situ station-based observations are also important, as major coastal disasters caused by TCs are fierce wind, torrential rainstorms, and formidable storm surges and waves. Consequently, 10-m wind speed and accumulated rainfall are compared with 2739 surface automatic weather stations in Guangdong province (Figure S5a). Quality control is implemented, and stations with accumulated rainfall less than 50 mm in each TC exposure are eliminated in rainfall evaluation. Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e shows that the model overestimates the wind speed by 1.2 m/s (+\u0026thinsp;7.4%) with 5.7 m/s root mean square error (RMSE) during the TC exposure. For accumulated rainfall, the COAWST underestimated accumulated rainfall by 4.8 mm (-2.8%) across all qualified stations. While the RMSE is 70.6 mm, relatively small mean bias and large RMSE indicate the model performs okay in reproducing the mean rainfall amount but poorly captures the fine-scale spatial distribution. We then compared storm surge levels using the Hong Kong Observatory tide gauges in Hong Kong (Fig. S4b). Tide-induced water level fluctuations relative to the Chart Datum were linearly combined with the pure storm surge level produced by COAWST to represent the total surge. Results indicate that the simulated surge level is overestimated by 0.13m (+\u0026thinsp;5.5%) with RMSE 0.35m. The two metrics show that the model can reproduce the surge variability for the three cases with reasonable accuracy (order of 10 cm) in this region, similar to the magnitude compared with previous studies using various models (Wang et al. \u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e2008\u003c/span\u003e, Yin et al. \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e2020\u003c/span\u003e, Chen et al. \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). Following Li et al. (\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2022\u003c/span\u003e), significant wave height is compared with in-situ wave observations at four buoys (QF303, QF304, QF306, and QF307) located in the continental shelf sea in the northern South China Sea (see Figure S4c), within 100-km distances to the nearest coast. Due to limited data availability, the significant wave height is only evaluated for Mangkhut (2018), with \u0026minus;\u0026thinsp;0.56 m (-9.3%) underestimation in mean bias and RMSE 1.12 m. Still, it shows the ability of the wave model to produce reliable wave height, given that the wind can be reasonably simulated by the atmospheric model.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eStation observations in comparison with model simulations\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"4\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMean bias\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eRelative bias\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eRoot mean square error\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e10-m wind speed (m/s)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e+\u0026thinsp;1.2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e+\u0026thinsp;7.4%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e5.7\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eAccumulated rainfall (mm)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e-4.8\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-2.8%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e70.6\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eStorm surge level (m)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e+\u0026thinsp;0.13\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e+\u0026thinsp;5.5%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e0.35\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eSignificant wave height (m)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e-0.56\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e-9.3%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e1.12\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec7\" class=\"Section2\"\u003e \u003ch2\u003e3.2 Changes in wind and precipitation\u003c/h2\u003e \u003cp\u003eSevere wind can cause dramatic damage in coastal regions, especially for metropolises such as Shenzhen and Hong Kong in the PRD region. During Mangkhut (2018), the collapse of trees caused extensive road blockage, paralyzing land transportation. Broken windows and falling debris due to building damages seriously threaten the safety of residents (Choy et al. \u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). Figure\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e estimates the grid-wise maximum 10-m wind speed change using the PGW experiment minus the CTRL. Negative values are masked out for clarity.\u003c/p\u003e \u003cp\u003eAll three cases result in an increase in the 10-m maximum wind speed in the PGW experiments for the areas exposed by the TC circulation, especially for the coastal area where it is exposed to the TC right-of-track. Please note that the difference in 10-m wind are probably coming from the variation (drift) of TC\u0026rsquo;s center in the two experiments rather than the change in intensity. To eliminate this effect, quantitative analysis based on the distance to TC\u0026rsquo;s center has been implemented. For example, the maximum 10-m wind speed within 50\u0026ndash;200-km right to the TC center increases by 4.3 m/s (22%). Tran et al. (\u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2022\u003c/span\u003e) attributed the increase in the TC right-of-track surface maximum wind speed to the faster transition speed of the TCs. While our study's three cases cannot definitively conclude that TCs are moving faster, the stronger intensity of these storms will certainly contribute to higher wind speed exposure for vulnerable coastal regions. In the case of Vicente (2012), which experienced the largest intensity increment in the PGW experiment, large areas of increased wind speed were found within the coverage of the TC's circulation. On average, the maximum 10-m wind speed within 50\u0026ndash;200 km to the right of the TC center increased by 6.9 m/s (39%). In a warming future, special attention should be paid to TCs with stationary tracks, like Vicente (2012), as they are likely to experience abrupt and acute rapid intensification near the shore, resulting in much greater wind impacts on coastal areas (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ec).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eTCs often bring torrential rainfall, which can cause severe flooding and landslides, adding to the list of coastal disasters associated with these storms. The resulting damage to infrastructure, property, and loss of life can be catastrophic, particularly in low-lying or densely populated areas. Figure\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e shows the grid-wise histogram of accumulated rainfall (mm) and hourly rain rate within a radius of 500 km of the TC center over the lifetime. The blue solid curve represents the CTRL experiment, and the red dashed curve represents the PGW experiment. A robust accumulated rainfall increase was observed in all three cases, with an average of 96.9 mm in CTRL and 107.9 mm in PGW, representing an 11.7% increase. Previous studies have shown that the global TC precipitation increases by slightly more than 7% per 1\u0026deg;C of tropical SST warming (the Clausius\u0026ndash;Clapeyron relation) in multiple climate model projections (Knutson et al. \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). Given that the domain-averaged additional SST perturbation is about 2.4\u0026deg;C, an 11.7% increase in accumulated rainfall is slightly lower than the Clausius\u0026ndash;Clapeyron relation. However, the extreme tails exhibit a slightly super-Clausius\u0026ndash;Clapeyron relation, suggesting that special attention should be paid to extreme rainfall events. As shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ea\u0026ndash;c, the PGW experiment significantly increased the percentage of accumulated rainfall by more than 200 mm in the three cases, from 11.0\u0026ndash;13.9% (+\u0026thinsp;26.4%). The 99th and the 99.9th percentile of accumulated rainfall threshold increased, from 405 mm to 475 mm (+\u0026thinsp;17.3%), and from 619 mm to 735 mm (+\u0026thinsp;18.6%), respectively. For the hourly rain rate (\u0026ge;0.1mm/h), the 99th and the 99.9th percentile of hourly rain rate increased from 40.9 mm/h to 45.3 mm/h (+\u0026thinsp;10.6%), and from 82.2 mm/h to 85.9 mm/h (+\u0026thinsp;4.5%), respectively (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ed\u0026ndash;f).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec8\" class=\"Section2\"\u003e \u003ch2\u003e3.3 Sea waves and storm surges\u003c/h2\u003e \u003cp\u003eTCs can cause devastating storm surges, which can be further exacerbated by spring tides, leading to catastrophic consequences such as flooding in low-lying areas, infrastructure damage, and loss of life, particularly in areas with poor drainage or inadequate flood protection. Local coastline geometry can also enhance flooding (Lowe et al. 2005, Wang et al. \u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e2008\u003c/span\u003e, Wu et al. \u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e2019\u003c/span\u003e, Lin et al. \u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). Wave-surge and wave-current interactions are important processes during storm surges and can be resolved by coupling a 3-dimensional ocean model and a spectral sea wave model, such as COAWST (Olabarrieta et al. \u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e2011\u003c/span\u003e, Benetazzo et al. \u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e2013\u003c/span\u003e). Grid-level maximum significant wave height and storm surge height difference using PGW minus CTRL in ROMS\u0026thinsp;+\u0026thinsp;SWAN D02 are shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e. The changes in maximum significant wave height exhibit a complex pattern due to sensitivity to TC tracks (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003ea/c/e). We found that, by masking out the ocean grid approximating the shoreline (see Figure S5), the maximum significant wave height at the ocean is lifted by 57 cm (+\u0026thinsp;13.8%) on average. Notably, there is a large difference between Hato (2017) and Vicente (2012).In the PGW experiment, the surface wind speed on the left-hand side of the TC does not increase for Hato (2017), limiting wave generation. In contrast, more than a 200 cm wave height increase is observed at the open sea near the TC eyewall and along the coastline for Vicente (2012). Storm surge exhibits a relatively uniform increment pattern (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eb/d/f). On average, maximum pure surge levels are lifted by 30\u0026ndash;80 cm over the open sea but increase much higher along the coastline, especially for narrowing estuaries where the maximum surge level can be lifted by up to 200 cm.\u003c/p\u003e \u003cp\u003eTo better illustrate the effect on coastal regions, averaged changes for three cases on the oceanic coastal grids (see Figure S5) are shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eg\u0026ndash;h. More than a 200 cm increase in significant wave height is observed over the offshore islands facing the open ocean, while the mainland coastline typically faced a 40\u0026ndash;80 cm increase. The lee side of the islands generally did not face significant wave height increases more than 60 cm, and for most inner riverbanks, only a lower than 20 cm wave height increase can be found. Conversely, storm surge level is sensitive to coastline geometry (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eh), with the inner riverbanks facing more than a 60 cm increase in surge level. Most coastlines would face a 40\u0026ndash;60 cm increase in storm surge, while only the sheltered Hong Kong central waters evaded the large surge level changes, with only a 20\u0026ndash;40 cm additional water elevation. According to multi-model mean results from CMIP6, compared to the present-day climatology (1995\u0026ndash;2014), global mean sea level rise will be 0.65 m at the end of the 21st century (2081\u0026ndash;2100) in SSP5-8.5 (Hermans et al. \u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). In the context of the 2090s, a combination of mean sea level rise, storm surge, and wave height can reach more than a 300 cm increase in total water level at certain hot-spot coastlines, without considering the superposition of spring tides.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e"},{"header":"Conclusions and discussion","content":"\u003cp\u003eThe PGW experiments, carried out by a novel AWO coupled model, mimic the thermodynamic conditions expected in the 2090s. The results suggest that changes in intense tropical cyclones (TCs) could have a significant impact on the PRD region, home to large populations, and the pivot of manufacturing, finance, and marine industries in China. Several coastal disaster indicators, including surface maximum wind speed, accumulated rainfall, hourly rain rate, storm surge height, and storm wave height, become more severe in the PGW experiments as the intensity of the selected TC cases increases. Our study highlights the particular concern of typhoons like Vicente (2012), which have stationary tracks and remain in warm waters for a long time. In a warming future, the warmer and stabler surface ocean, even after a long TC-induced \"cold wake\", which has an important negative feedback on TC intensity (Bender and Ginis \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2000\u003c/span\u003e, Knutson et al. \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2001\u003c/span\u003e), cannot weaken these types of TCs. They tend to experience abrupt and acute rapid intensification nearshore, resulting in a much shorter response time and greater impacts on coastal areas.\u003c/p\u003e \u003cp\u003eOur analysis focused on the extreme values in the wind speed distribution, and we found that the maximum 10-m wind speed within 50\u0026ndash;200 km to the right of the TC center increased by 4.3 m/s (or 22%). For Vicente (2012), the increase was even greater at 6.9 m/s (or 39%). The increase in the TC right-of-track surface maximum wind speed may be attributed to the faster transition speed of the TCs, while the stronger intensity of these TCs tends to bring wind impacts further inland (Tran et al. \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Regarding the rainfall impact, we observed an average 11.7% increase in accumulated rainfall, corresponding to a slightly less than 7% per 1\u0026deg;C increase relative to tropical SST warming (the Clausius-Clapeyron relation). It is noteworthy that the percentage of accumulated rainfall by more than 200 mm will increase by 26.4%, and the 99th and 99.9th percentile of accumulated rainfall threshold will increase by +\u0026thinsp;17.3% and +\u0026thinsp;18.6%, respectively. The greater than 7% per 1\u0026deg;C increase of the right tail of the rainfall distribution indicates much more severe extreme precipitation caused by TCs, significantly increasing the risks of coastal landslides and flooding in a warming future.\u003c/p\u003e \u003cp\u003ePrevious studies have shown that, assuming all other factors are unchanged, the most confident TC-related projection is that sea level rise accompanying the warming will lead to higher storm inundation levels (Knutson et al. \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). In this study, we went beyond previous research by considering the surface wave component and its coupling with the 3-dimensional ocean and atmosphere. Simulation results show that even without considering the mean sea level rise and spring tides, stronger TC intensity in a warming climate can already push more than 100 cm of additional water level to most PRD coastlines by combining surge level and wave height changes. Given the significant uncertainties in mean sea level rise, especially for the mass loss in Greenland and Antarctic ice sheets (Hermans et al. \u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e2021\u003c/span\u003e), from a risk management perspective, we suggest a stress test considering a total water level increase of more than 300 cm in the PRD region due to storm surge in the context of the 2090s, even without considering the overlap with spring tides. Such an approach would help to ensure that the PRD coastal areas and their populations are better prepared for future climate risks and can adapt to the potential impacts of TCs on coastal flooding.\u003c/p\u003e \u003cp\u003eIt is important to note that our current study has several limitations. We have only examined three intense landfalling TCs since 2010s. The small sampling size would cause potential bias in future projections. Further studies, especially with a much larger sample size, are needed to confirm the current findings and explore other potential impacts of future changes in intense TCs on the PRD coastal regions. Meanwhile, it is important to acknowledge the limitations of our approach in studying the storm surge impact, as we did not use an inundation model to estimate the extent and depth of coastal flooding during storm surge events. This limitation may have restricted our ability to fully quantify the potential impacts of climate change on storm surges in the PRD region. Many other important parameters, e.g. the potential shift of the TC tracks in the warming climate and potential rapid intensification or rapid weakening of TCs over the study region are needed to investigated by ensemble simulation technique. Ensemble technique also support robust analysis based on impact indices such as the power dissipation index (PDI, Emanuel \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e2005\u003c/span\u003e) and cyclone damage potential (CDP, Holland et al. \u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e2019\u003c/span\u003e), which can be further incorporated to study the TC changes.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003e\u003cem\u003eFunding\u003c/em\u003e\u003c/strong\u003e We appreciate the assistance of the Hong Kong Observatory (HKO) and Environmental Central Facility (HKUST), which provided the meteorological and model input data. This work is supported by grants from the Research Grants Council of the Hong Kong Special Administrative Region, China (Project No. AoE/E-603/18 and T31-603/21-N), the Environment and Conservation Fund (2022-02), the CORE Project 2022/23, and the CORE is a joint research center for ocean research between QNLM and HKUST.\u003cstrong\u003e\u003cem\u003e\u0026nbsp;\u003c/em\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e\u003cem\u003eFinancial and non-financial interests\u003c/em\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors have no relevant financial or non-financial interests to disclose.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eBalaguru, K., Foltz, G.R., Leung, L.R. and Emanuel, K.A., 2016. Global warming-induced upper-ocean freshening and the intensification of super typhoons. Nature Communications, 7(1), p.13670.\u003c/li\u003e\n\u003cli\u003eBender, M.A. and Ginis, I., 2000. Real-case simulations of hurricane\u0026ndash;ocean interaction using a high-resolution coupled model: Effects on hurricane intensity. Monthly Weather Review, 128(4), pp.917-946.\u003c/li\u003e\n\u003cli\u003eBenetazzo, A., Carniel, S., Sclavo, M. and Bergamasco, A., 2013. Wave\u0026ndash;current interaction: Effect on the wave field in a semi-enclosed basin. Ocean Modelling, 70, pp.152-165.\u003c/li\u003e\n\u003cli\u003eBloemendaal, N., de Moel, H., Martinez, A.B., Muis, S., Haigh, I.D., van der Wiel, K., Haarsma, R.J., Ward, P.J., Roberts, M.J., Dullaart, J.C. and Aerts, J.C., 2022. A globally consistent local-scale assessment of future tropical cyclone risk. Science Advances, 8(17), p.eabm8438.\u003c/li\u003e\n\u003cli\u003eChan JCL (2008) Decadal variations of intense typhoon occurrence in the western North Pacific. Proc R Soc Lond A 464:249\u0026ndash;272. doihttps://doi.org/10.1098/rspa.2007.0183\u003c/li\u003e\n\u003cli\u003eChang, Y.P., Yang, S.C., Lin, K.J., Lien, G.Y. and Wu, C.M., 2020. Impact of tropical cyclone initialization on its convection development and intensity: A case study of Typhoon Megi (2010). Journal of the Atmospheric Sciences, 77(2), pp.443-464.\u003c/li\u003e\n\u003cli\u003eChen, J., Wang, Z., Tam, C.Y., Lau, N.C., Lau, D.S.D. and Mok, H.Y., 2020. Impacts of climate change on tropical cyclones and induced storm surges in the Pearl River Delta region using pseudo-global-warming method. Scientific Reports, 10(1), pp.1-10.\u003c/li\u003e\n\u003cli\u003eChen, S. S., Zhao, W., Donelan, M. A., \u0026amp; Tolman, H. L. 2013. Directional wind-wave coupling in fully coupled atmosphere-wave-ocean models: Results from CBLAST-Hurricane. Journal of the Atmospheric Sciences, 70(10), 3198\u0026ndash;3215.\u003c/li\u003e\n\u003cli\u003eChen, X., Wang, Y., Zhao, K. and Wu, D., 2017. A numerical study on rapid intensification of Typhoon Vicente (2012) in the South China Sea. Part I: Verification of simulation, storm-scale evolution, and environmental contribution. Monthly Weather Review, 145(3), pp.877-898.\u003c/li\u003e\n\u003cli\u003eChoy, C.W., of Insurers, T.H.K.F., Wu, M.C. and Lee, T.C., 2020. Assessment of the damages and direct economic loss in Hong Kong due to Super Typhoon Mangkhut in 2018. Tropical Cyclone Research and Review, 9(4), pp.193-205.\u003c/li\u003e\n\u003cli\u003eChu, J.E., Lee, S.S., Timmermann, A., Wengel, C., Stuecker, M.F. and Yamaguchi, R., 2020. Reduced tropical cyclone densities and ocean effects due to anthropogenic greenhouse warming. Science Advances, 6(51), p.eabd5109.\u003c/li\u003e\n\u003cli\u003eCummings, J.A. and Smedstad, O.M., 2013. Variational data assimilation for the global ocean. In Data Assimilation for Atmospheric, Oceanic and Hydrologic Applications (Vol. II) (pp. 303-343). Springer, Berlin, Heidelberg.\u003c/li\u003e\n\u003cli\u003eEgbert, G.D. and Erofeeva, S.Y., 2002. Efficient inverse modeling of barotropic ocean tides. Journal of Atmospheric and Oceanic Technology, 19(2), pp.183-204.\u003c/li\u003e\n\u003cli\u003eEmanuel, K.A., 1986. An air-sea interaction theory for tropical cyclones. Part I: Steady-state maintenance. Journal of Atmospheric Sciences, 43(6), pp.585-605.\u003c/li\u003e\n\u003cli\u003eEmanuel, K. A., 2005: Increasing destructiveness of tropical cyclones over the past 30 years. Nature, 436, 686\u0026ndash;688. https://doi.org/10.1038/nature03906\u003c/li\u003e\n\u003cli\u003eGutjahr, O., Putrasahan, D., Lohmann, K., Jungclaus, J.H., von Storch, J.S., Br\u0026uuml;ggemann, N., Haak, H. and St\u0026ouml;ssel, A., 2019. Max Planck Institute Earth System Model (MPI-ESM1. 2) for the high-resolution model intercomparison project (HighResMIP). Geoscientific Model Development, 12(7), pp.3241-3281.\u003c/li\u003e\n\u003cli\u003eHermans, T.H., Gregory, J.M., Palmer, M.D., Ringer, M.A., Katsman, C.A. and Slangen, A.B., 2021. Projecting global mean sea‐level change using CMIP6 models. Geophysical Research Letters, 48(5), p.e2020GL092064.\u003c/li\u003e\n\u003cli\u003eHolland, G.J., Done, J.M., Douglas, R., Saville, G.R. and Ge, M., 2019. Global tropical cyclone damage potential. Hurricane risk, pp.23-42.\u003c/li\u003e\n\u003cli\u003eHuang, P., Lin, I.I., Chou, C. and Huang, R.H., 2015. Change in ocean subsurface environment to suppress tropical cyclone intensification under global warming. Nature Communications, 6(1), p.7188.\u003c/li\u003e\n\u003cli\u003eJyoteeshkumar Reddy, P., Sriram, D., Gunthe, S.S. and Balaji, C., 2021. Impact of climate change on intense Bay of Bengal tropical cyclones of the post-monsoon season: a pseudo global warming approach. Climate Dynamics, 56, pp.2855-2879.\u003c/li\u003e\n\u003cli\u003eKnutson, T., Camargo, S.J., Chan, J.C., Emanuel, K., Ho, C.H., Kossin, J., Mohapatra, M., Satoh, M., Sugi, M., Walsh, K. and Wu, L., 2019. Tropical cyclones and climate change assessment: Part I: Detection and attribution. Bulletin of the American Meteorological Society, 100(10), pp.1987-2007.\u003c/li\u003e\n\u003cli\u003eKnutson, T., Camargo, S.J., Chan, J.C., Emanuel, K., Ho, C.H., Kossin, J., Mohapatra, M., Satoh, M., Sugi, M., Walsh, K. and Wu, L., 2020. Tropical cyclones and climate change assessment: Part II: Projected response to anthropogenic warming. Bulletin of the American Meteorological Society, 101(3), pp.E303-E322.\u003c/li\u003e\n\u003cli\u003eKnutson, T.R., McBride, J.L., Chan, J., Emanuel, K., Holland, G., Landsea, C., Held, I., Kossin, J.P., Srivastava, A.K. and Sugi, M., 2010. Tropical cyclones and climate change. Nature Geoscience, 3(3), pp.157-163.\u003c/li\u003e\n\u003cli\u003eKnutson, T.R., Tuleya, R.E., Shen, W. and Ginis, I., 2001. Impact of CO2-induced warming on hurricane intensities as simulated in a hurricane model with ocean coupling. Journal of Climate, 14(11), pp.2458-2468.\u003c/li\u003e\n\u003cli\u003eLi, Z., Tam, C.Y., Li, Y., Lau, N.C., Chen, J., Chan, S.T., Dickson Lau, D.S. and Huang, Y., 2022. How Does Air‐Sea Wave Interaction Affect Tropical Cyclone Intensity? An Atmosphere‐Wave‐Ocean Coupled Model Study Based on Super Typhoon Mangkhut (2018). Earth and Space Science, 9(3), p.e2021EA002136.\u003c/li\u003e\n\u003cli\u003eLin, N., Marsooli, R. and Colle, B.A., 2019. Storm surge return levels induced by mid-to-late-twenty-first-century extratropical cyclones in the Northeastern United States. Climatic Change, 154, pp.143-158.\u003c/li\u003e\n\u003cli\u003eLin, S. and Sheng, J., 2020. Revisiting dependences of the drag coefficient at the sea surface on wind speed and sea state. Continental Shelf Research, 207, p.104188.\u003c/li\u003e\n\u003cli\u003eLiu, K.S. and Chan, J.C., 2020. Recent increase in extreme intensity of tropical cyclones making landfall in South China. Climate Dynamics, 55(5-6), pp.1059-1074.\u003c/li\u003e\n\u003cli\u003eLowe, J.A. and Gregory, J.M., 2005. The effects of climate change on storm surges around the United Kingdom. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 363(1831), pp.1313-1328.\u003c/li\u003e\n\u003cli\u003eMagnusson, L., Bidlot, J. R., Bonavita, M., Brown, A. R., Browne, P. A., De Chiara, G., et al. (2019). ECMWF activities for improved hurricane forecasts. Bulletin of the American Meteorological Society, 100(3), 445\u0026ndash;458.\u003c/li\u003e\n\u003cli\u003eMogensen, K. S., Magnusson, L., \u0026amp; Bidlot, J. R. (2017). Tropical cyclone sensitivity to ocean coupling in the ECMWF coupled model. Journal of Geophysical Research: Oceans, 122, 4392\u0026ndash;4412. https://doi.org/10.1002/2017JC012753\u003c/li\u003e\n\u003cli\u003eMoon, I. J., Ginis, I., \u0026amp; Hara, T. (2004). Effect of surface waves on Charnock coefficient under tropical cyclones. Geophysical Research Letters, 31, L20302. https://doi.org/10.1029/2004GL020988\u003c/li\u003e\n\u003cli\u003eMoon, J., Cha, D.H., Lee, M. and Kim, J., 2018. Impact of spectral nudging on real‐time tropical cyclone forecast. Journal of Geophysical Research: Atmospheres, 123(22), pp.12-647.\u003c/li\u003e\n\u003cli\u003eNakanishi, M. and Niino, H., 2004. An improved Mellor\u0026ndash;Yamada level-3 model with condensation physics: Its design and verification. Boundary-layer Meteorology, 112, pp.1-31.\u003c/li\u003e\n\u003cli\u003eOlabarrieta, M., Warner, J.C. and Kumar, N., 2011. Wave‐current interaction in Willapa Bay. Journal of Geophysical Research: Oceans, 116(C12).\u003c/li\u003e\n\u003cli\u003eOlabarrieta, M., Warner, J.C., Armstrong, B., Zambon, J.B. and He, R. (2012) Ocean\u0026ndash; atmosphere dynamics during Hurricane Ida and Nor\u0026rsquo;Ida: An application of the coupled ocean\u0026ndash;atmosphere\u0026ndash;wave\u0026ndash;sediment transport (COAWST) modeling system. Ocean Modelling, 43, pp.112-137.\u003c/li\u003e\n\u003cli\u003eTaylor, P. K., \u0026amp; Yelland, M. J. (2001). The dependence of sea surface roughness on the height and steepness of the waves. Journal of Physical Oceanography, 31(2), 572\u0026ndash;590.\u003c/li\u003e\n\u003cli\u003eTong, M., Sippel, J.A., Tallapragada, V., Liu, E., Kieu, C., Kwon, I.H., Wang, W., Liu, Q., Ling, Y. and Zhang, B., 2018. Impact of assimilating aircraft reconnaissance observations on tropical cyclone initialization and prediction using operational HWRF and GSI ensemble\u0026ndash;variational hybrid data assimilation. Monthly Weather Review, 146(12), pp.4155-4177.\u003c/li\u003e\n\u003cli\u003eTran, T.L., Ritchie, E.A., Perkins‐Kirkpatrick, S.E., Bui, H. and Luong, T.M., 2022. Future Changes in Tropical Cyclone Exposure and Impacts in Southeast Asia From CMIP6 Pseudo‐Global Warming Simulations. Earth\u0026apos;s Future, 10(12), p.e2022EF003118.\u003c/li\u003e\n\u003cli\u003eTsuboki, K., Yoshioka, M.K., Shinoda, T., Kato, M., Kanada, S. and Kitoh, A., 2015. Future increase of supertyphoon intensity associated with climate change. Geophysical Research Letters, 42(2), pp.646-652.\u003c/li\u003e\n\u003cli\u003eWang, S., McGrath, R., Hanafin, J., Lynch, P., Semmler, T. and Nolan, P., 2008. The impact of climate change on storm surges over Irish waters. Ocean Modelling, 25(1-2), pp.83-94.\u003c/li\u003e\n\u003cli\u003eWarner, J.C., Armstrong, B., He, R. and Zambon, J.B. (2010) Development of a coupled ocean\u0026ndash;atmosphere\u0026ndash;wave\u0026ndash;sediment transport (COAWST) modeling system. Ocean Modelling, 35(3), pp.230-244.\u003c/li\u003e\n\u003cli\u003eWu, C.C., Chou, K.H., Wang, Y. and Kuo, Y.H., 2006. Tropical cyclone initialization and prediction based on four-dimensional variational data assimilation. Journal of the Atmospheric Sciences, 63(9), pp.2383-2395.\u003c/li\u003e\n\u003cli\u003eWu, C.C., Li, T.H. and Huang, Y.H., 2015. Influence of mesoscale topography on tropical cyclone tracks: Further examination of the channeling effect. Journal of the Atmospheric Sciences, 72(8), pp.3032-3050.\u003c/li\u003e\n\u003cli\u003eWu, J., Chen, Y., Liao, Z., Gao, X., Zhai, P. and Hu, Y., 2022. Increasing risk from landfalling tropical cyclone-heatwave compound events to coastal and inland China. Environmental Research Letters, 17(10), p.105007.\u003c/li\u003e\n\u003cli\u003eWu, L., Breivik, \u0026Oslash;. and Rutgersson, A., 2019. Ocean‐wave‐atmosphere interaction processes in a fully coupled modeling system. Journal of Advances in Modeling Earth Systems, 11(11), pp.3852-3874.\u003c/li\u003e\n\u003cli\u003eWu, W. and Leonard, M., 2019. Impact of ENSO on dependence between extreme rainfall and storm surge. Environmental Research Letters, 14(12), p.124043.\u003c/li\u003e\n\u003cli\u003eYin, K., Xu, S., Zhao, Q., Huang, W., Yang, K. and Guo, M., 2020. Effects of land cover change on atmospheric and storm surge modeling during typhoon event. Ocean Engineering, 199, p.106971.\u003c/li\u003e\n\u003cli\u003eZambon, J.B., He, R. and Warner, J.C. (2014a) Investigation of hurricane Ivan using the coupled ocean\u0026ndash;atmosphere\u0026ndash;wave\u0026ndash;sediment transport (COAWST) model. Ocean Dynamics, 64(11), pp.1535-1554.\u003c/li\u003e\n\u003cli\u003eZambon, J.B., He, R. and Warner, J.C. (2014b) Tropical to extratropical: Marine environmental changes associated with Superstorm Sandy prior to its landfall. Geophysical Research Letters, 41(24), pp.8935-8943.\u003c/li\u003e\n\u003cli\u003eZhang, F., Weng, Y., Gamache, J. F., \u0026amp; Marks, F. D. (2011). Performance of convection-permitting hurricane initialization and prediction during 2008\u0026ndash;2010 with ensemble data assimilation of inner-core airborne Doppler radar observations. Geophysical Research Letters, 38, L15810.\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":true,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"natural-hazards","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"nhaz","sideBox":"Learn more about [Natural Hazards](https://www.springer.com/journal/11069)","snPcode":"11069","submissionUrl":"https://submission.nature.com/new-submission/11069/3","title":"Natural Hazards","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false},"keywords":"","lastPublishedDoi":"10.21203/rs.3.rs-3625663/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-3625663/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eThe Pearl River Delta (PRD) region is highly vulnerable to tropical cyclone (TC)-caused coastal hazards due to its long and meandering shoreline and well-developed economy. With global warming expected to continue or worsen in the rest of the 21st century, this study examines the TC impact on the PRD coastal regions by reproducing three intense landfalling TCs using a sophisticated air-wave-ocean coupled model of high spatial resolution (1-km atmosphere and 500-m wave and ocean). The simulations are conducted using present-day reanalysis data and the same TCs occurring in a pseudo-global warming scenario projected for the 2090s. Results indicate that the coupled model accurately reproduces the air-wave-ocean status during the TC episodes. The 2090s thermodynamic status effectively increases the intensity of intense TCs, leading to more severe coastal hazards including gale, rainstorm, and storm surges and waves. On average, the maximum surface wind speed within 50\u0026ndash;200 km to the right of the TC center can increase by 4.3 m/s (+\u0026thinsp;22%). The 99th and the 99.9th percentile of accumulated rainfall will increase from 405 mm to 475 mm (+\u0026thinsp;17.3%), and from 619 mm to 735 mm (+\u0026thinsp;18.6%), respectively. The maximum significant wave height at the ocean is lifted by an average of 57 cm (+\u0026thinsp;13.8%), and the coastline typically faces a 40\u0026ndash;80 cm increase. The maximum storm surges are lifted by 30\u0026ndash;80 cm over the open sea but aggravate much higher along the coastline, especially for narrowing estuaries. For Typhoon Vicente (2012), there is more than a 200 cm wave height increase observed both at open sea and along the coastline. In the 2090s context, a combination of mean sea level rise, storm surge, and wave height can reach more than 300 cm increase in total water level at certain hot-spot coastlines, without considering the superposition of spring tides.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e","manuscriptTitle":"Future changes in intense tropical cyclone hazards in the Pearl River Delta region: An air-wave-ocean coupled model study","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2023-11-29 16:10:38","doi":"10.21203/rs.3.rs-3625663/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"reviewerAgreed","content":"","date":"2023-11-26T11:45:48+00:00","index":0,"fulltext":""},{"type":"reviewersInvited","content":"","date":"2023-11-24T10:33:14+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"Natural Hazards","date":"2023-11-18T19:22:48+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2023-11-16T08:36:49+00:00","index":"","fulltext":""},{"type":"submitted","content":"Natural Hazards","date":"2023-11-16T03:21:52+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"natural-hazards","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"nhaz","sideBox":"Learn more about [Natural Hazards](https://www.springer.com/journal/11069)","snPcode":"11069","submissionUrl":"https://submission.nature.com/new-submission/11069/3","title":"Natural Hazards","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false}}],"origin":"","ownerIdentity":"3f1587ad-d063-4893-ae9e-5eff76b0a49c","owner":[],"postedDate":"November 29th, 2023","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[],"tags":[],"updatedAt":"2024-03-11T15:06:13+00:00","versionOfRecord":{"articleIdentity":"rs-3625663","link":"https://doi.org/10.1007/s11069-024-06510-7","journal":{"identity":"natural-hazards","isVorOnly":false,"title":"Natural Hazards"},"publishedOn":"2024-03-09 15:01:49","publishedOnDateReadable":"March 9th, 2024"},"versionCreatedAt":"2023-11-29 16:10:38","video":"","vorDoi":"10.1007/s11069-024-06510-7","vorDoiUrl":"https://doi.org/10.1007/s11069-024-06510-7","workflowStages":[]},"version":"v1","identity":"rs-3625663","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-3625663","identity":"rs-3625663","version":["v1"]},"buildId":"7rjqhiLT3MXkJMwkYKINL","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}