
Impact-based forecasting of tropical cyclone-related human displacement to support anticipatory action
- Select a language for the TTS:
- UK English Female
- UK English Male
- US English Female
- US English Male
- Australian Female
- Australian Male
- Language selected: (auto detect) - EN
Play all audios:

ABSTRACT Tropical cyclones (TCs) displace millions every year. While TCs pose hardships and threaten lives, their negative impacts can be reduced by anticipatory actions like evacuation and
humanitarian aid coordination. In addition to weather forecasts, impact forecast enables more effective response by providing richer information on the numbers and locations of people at
risk of displacement. We introduce a fully open-source implementation of a globally consistent and regionally calibrated TC-related displacement forecast at low computational costs,
combining meteorological forecast with population exposure and respective vulnerability. We present a case study of TC Yasa which hit Fiji in December 2020. We emphasise the importance of
considering the uncertainties associated with hazard, exposure, and vulnerability in a global uncertainty analysis, which reveals a considerable spread of possible outcomes. Additionally, we
perform a sensitivity analysis on all recorded TC displacement events from 2017 to 2020 to understand how the forecast outcomes depend on these uncertain inputs. Our findings suggest that
for longer forecast lead times, decision-making should focus more on meteorological uncertainty, while greater emphasis should be placed on the vulnerability of the local community shortly
before TC landfall. Our open-source codes and implementations are readily transferable to other users, hazards, and impact types. SIMILAR CONTENT BEING VIEWED BY OTHERS RECENT NORTHWARD
SHIFT OF TROPICAL CYCLONE ECONOMIC RISK IN CHINA Article Open access 31 May 2024 UNCERTAINTIES AND SENSITIVITIES IN THE QUANTIFICATION OF FUTURE TROPICAL CYCLONE RISK Article Open access 13
October 2023 INCREASING TYPHOON IMPACT AND ECONOMIC LOSSES DUE TO ANTHROPOGENIC WARMING IN SOUTHEAST CHINA Article Open access 08 September 2022 INTRODUCTION Human displacement occurs when
people are forced to leave their homes or places of habitual residence due to external events such as natural extreme weather hazards1. Every year, weather extremes cause millions of people
around the world to be displaced. Tropical cyclones (TCs) account for the second largest contribution to these human displacements after floods, with an average of 9.3 million people being
displaced every year between 2017 and 20202. The duration of displacement and its humanitarian impact vary widely: from short-term pre-emptive evacuations to long-term displacement if houses
or community infrastructures are significantly damaged, or people lose access to their economic activities. To cope with the post-disaster recovery, international and national assistance
for humanitarian response and relief funds are often required. The range of assistance includes providing emergency shelter, clean water and food, health care, psychological support, and
long-term community and livelihood recovery and restoration3,4. Anticipatory actions can help to reduce the negative impacts of extreme weather events. Examples of anticipatory actions
include evacuation planning, emergency protection, and humanitarian aid coordination. The World Meteorological Organization (WMO) proposed the WMO Coordination Mechanism (WCM)5 and a pilot
project Weather4UN5 that enable access to weather and climate information and the provision of expert advice from WMO Members to the member states and other humanitarian agencies. Another
similar scheme is the Forecast-based financing (FbF) proposed by the International Federation of Red Cross and Red Crescent Societies (IFRC)6. It is a specific financial scheme that allows
access to financial resources based on scientific forecasts and risk analysis for humanitarian actions agreed in advance. The allocation of the funds is automatically released when a certain
forecast threshold is reached, which allows mitigation measures to take place prior to the hazard and thus reducing impacts7. Currently, the planning for anticipatory action is mostly based
on weather forecasts, which may potentially miss crucial nuances required to minimise negative impacts. For instance, knowing whether the same intensity of storm hits a populated city or a
scattered rural area, or being aware of the resilience of the directly affected community, will in general lead to different preemptive measures. At the moment, this auxiliary information is
indirectly taken into account by decision-makers who use expert knowledge and past experiences to guide their assessments, but there is rarely a systematic quantification available. An
impact forecast moves a step forward from conventional weather forecasts to give information on how the weather will affect people. It systematically translates the weather information into
risk by combining it with social variables8. There have been multiple emergency evacuation decisions support tools for tropical cyclones in the US that combine weather forecast with traffic
information to support evacuation plannings (e.g. Harris et al.9, Davidson et al.10, and Blanton et al.11), and establishing platforms that provides weather risk communication with
integrates societal information flow (e.g. CHIME; Morss et al.12). Currently, there are however no tools that provide globally consistent impact information for human displacement. Here we
introduce a proof-of-concept implementation of a global TC-related displacement impact forecast, that could provide more comparable, standardised, and less singular information to support
decision-making for anticipatory action. We use the open-source probabilistic natural catastrophe risk assessment platform CLIMADA (CLIMate ADAptation)13 for the displacement impact
forecast. The platform is widely used in long-term weather and climate impact assessments for building insurance and other socioeconomic impacts14,15,16,17,18,19,20, and has also been used
for establishing operational impact forecasts for building damages from winter storms21. CLIMADA generates risk information and quantifies socioeconomic impacts by integrating hazard,
exposure, and vulnerability data. In our study, we forecast the number of people who are at risk of being displaced due to upcoming TC events by combining datasets from i) TC track forecasts
and their associated windfields as hazard, ii) the global population distribution as exposures, and iii) the vulnerability of people being displaced represented by a univariate impact
function relating the TC wind speed to the probability of displacement at a given location, see Fig. 1. In typical disaster risk assessments, risks are quantified by considering long-term
aggregated average impacts from past events13. However, when it comes to forecasting a single event, we argue that it is important to not only consider the most probable outcome but also
look at the overall distribution of plausible forecast outcomes when deciding on the anticipatory action. The uncertainties in the forecasted impact arise from the complex interplay between
the meteorological forecast variability based on the ensemble forecast system22,23, and uncertainties associated with all other input data. We will show that it is crucial to account for
uncertainties in all inputs to provide a more comprehensive perspective on the forecasted outcomes. Furthermore, to make this uncertainty practical in the context of decision-making, we
conduct sensitivity analyses to gain a deeper understanding of how the relevant outcomes depend on these uncertain inputs. Here we demonstrate the TC impact forecast for displacement by
performing an analysis for the TC Yasa that hit Fiji in December 2020 and caused the displacement of 23,414 people2,24. We chose Fiji as a demonstrating case as the Pacific Islands are an
under-studied region, and while the island characteristics enable a showcase of hit-or-miss scenarios, the probability of whether there will be impacts are important when providing the
forecast information. We first forecast risk of displacement by considering only the meteorological forecast variability. Then we quantify the full uncertainty range of forecast outcomes
that includes uncertainty from the exposures and vulnerability. We repeat the calculation using a quasi-Monte Carlo sampling method to cover a wide range of possible input variations25,26,
and a sensitivity analysis to attribute which input variations contribute most to the overall uncertainty. We also generalise the main findings by performing the same displacement forecast
analysis at different lead times prior to the landfall for all past TC events which have displacement records from 2017 to 20202. Our results provide decision-makers with a wide range of
possible outcomes and guidance on which input variations should be considered for better-informed decisions on anticipatory actions. RESULTS IMPACT FORECAST FOR TC YASA Here we first
demonstrate the displacement forecast by performing a case study of the TC Yasa, which passed through Fiji on 17 December 2020. TC Yasa reached category 5 of the Saffir-Simpson wind scale,
with maximum sustained winds reaching 146 kts (75.1 m/s) before making its landfall in Fiji27. The cyclone was the strongest since TC Winston in 2016 and was the most destructive cyclone in
Fiji during the cyclone season 2020–2021, damaging infrastructure, buildings and agricultural areas24. IDMC has recorded 23,414 people on the Fijian islands as being displaced2. Figure 2a
shows a spatially explicit map of the forecast-ensemble-averaged displacement in Fiji, overlayed with the 51 ensemble forecast TC tracks from the Integrated Forecasting System maintained by
the European Centre for Medium-Range Weather Forecasts (ECMWF-IFS)22,23. The exposure is the static estimation of population distribution as described in subsection “Population exposures” of
section “Methods”, and the impact function (or vulnerability curve) is a calibrated function that minimises the relative errors over all the recorded displacement events in the Pacific
Island region, detailed in subsection “Calibration of impact functions for human displacement” of section “Methods”. The average forecasted number of people at risk of displacement in Fiji
is 172,000 (orange dashed line in Fig. 2b), whilst the total number of displaced people ranges from 3500 to 450,000 based on the maximum 1-minute sustained wind fields calculated from the 51
ensemble forecast TC tracks. The distribution of the number of displaced people per ensemble member, as shown in Fig. 2b, is right skewed with a long tail of high impacts. The large
displacement uncertainty here only comes from the spread of the TC tracks. The magnitude of the impact mainly depends on whether the TC hits the main island and passes through the populated
cities, or remains mostly over the ocean. The case that gives the highest number of people who could be potentially displaced features a forecasted TC track that crosses Fiji’s main island
of Viti Levu near the populated cities such as Suva, Lautoka and Nadi, whereas the lowest numbers come from tracks that do not make landfall and pass the Northeastern side of Fiji. GLOBAL
UNCERTAINTY AND SENSITIVITY ANALYSIS An effective impact forecast should not only include the variability of the meteorological forecast, but also the uncertainties within the exposure
information and impact functions21,25,28. Here we include the uncertainties from the total population count and replace the single impact function with an ensemble of functions, with each
member calibrated against a past event to represent the range of plausible functions (detailed in subsections “Calibration of impact functions for human displacement” and “Uncertainty and
sensitivity” of section “Methods”). Note that this calibration method differs from the one used for Fig. 2, as discussed in section “Discussion”. We ran the impact forecast model more than
8000 times with each run randomly sampling parameters for each model component (meteorological forecast, exposures, and impact functions). The computational time of this uncertainty
distribution starting from the TC forecast tracks extraction is around 10 minutes on a MacBook Pro equipped with 2.3 GHz Quad-Core Intel Core i7 and 16 GB random-access memory (RAM). Figure
3a shows the probability distribution of the global uncertainty analysis. The resulting distribution is even more right-skewed compared to 2b) (meteorological variability plot), yielding an
average of 123,391 people in Fiji being displaced by TC Yasa, but with the average and peak probability closer to the reported number of displacements by the IDMC. Furthermore, the
corresponding first-order sensitivity analysis shown in Fig. 3 indicates that the overall impact forecast uncertainty is most sensitive to the uncertainty from impact functions (sensitivity
index 0.411), with the meteorological forecast uncertainty also making a substantial contribution (0.320) at two days lead time of the forecast. The sensitivity index for the total
population, on the other hand, is very small (0.005). DISPLACEMENT IMPACT FORECAST AT DIFFERENT LEAD TIME We repeat the preceding analysis for lead times ranging from 3.5 days until just
before landfall at intervals of 0.5 days. Figure 4a shows the kernel density fit to the probability distribution of the global uncertainty of the displacement impact forecast at different
lead times. The distribution is more spread out for longer lead times, but narrows for the forecast close before landfall. Figure 4b shows the first-order Sobel sensitivity index29 of all
the uncertain input parameters. We observe a general decreasing trend of the sensitivity index for the meteorological forecast, while an increasing trend for the impact function. In the
meantime, the index for the change in total population remains close to zero at all forecast lead time. We remark that the increase in sensitivity to the vulnerability uncertainty at shorter
lead times cannot be equated with an increase in uncertainty from vulnerability, but only that the relative contribution to the total uncertainty from vulnerability increases, and from
forecast decreases. UNCERTAINTY AND SENSITIVITY ANALYSIS FOR PAST TC DISPLACEMENT EVENTS FROM 2017 TO 2020 We ran the same uncertainty and sensitivity analysis on the total number of
displaced people for worldwide recorded TC events from 2017 to 2020 at lead times from 3 to 0.5 days before landfall. We find the impact forecast for displacement show a bias of
overestimation when comparing to the reported number of displacement from IDMC (see Fig. S6.) Details of the impact forecast validation and biases are included in the supplementary
information S4. Figure 5 summarises the distributions of first-order sensitivity indices for all TC events. Similar to the trends identified in Fig. 4b, the median of the first-order
sensitivity index for meteorological forecast is the largest when the forecast lead time is longer, but decreases closer to the TC landfall. The impact function’s sensitivity index, in
contrast, is smaller with longer lead times but increases over time. Hence, on average over all considered events, we find that the relative contribution to the overall uncertainty from the
forecast is largest at long lead times, and from vulnerability at short lead times. This overall trend does however not hold for all individual cases as exemplified with TC Harold in the
Supplementary Figs. S4 and S5 for which also at long lead times the sensitivity to the vulnerability uncertainty is largest. Note that the forecast uncertainty for TC Harold is comparable to
TC Yasa even though the sensitivities are different. DISCUSSION Our work shows that it is possible to provide spatially explicit impact forecasts for tropical cyclones-related displacement
from publicly available data in near-real time. We combine the tropical cyclone (TC) track forecasts from ECMWF with the population exposure and vulnerability of displacement in the
open-source probabilistic risk assessment platform CLIMADA to forecast the displacement risk by TCs. We demonstrate the impact forecast for displacement with TC Yasa that hits Fiji in
December 2020. We argue that our probabilistic approach is more useful than a deterministic model. The spread of the displacement predictions reflects the statistical uncertainty in the
model, and it is important to consider the full uncertainty distribution when making decisions for anticipatory action, rather than only employing the often-used ensemble-averaged metrics.
We show that the modelled impact follows a strongly skewed distribution, and thus the averaged impact cannot reflect the full range of possible outcomes. Importantly, the tail risk of
high-impact scenarios, which can be the triggers for major disasters, is inadequately represented by average measures and must be derived from a quantified uncertainty distribution.
Uncertainty around the impact forecast stems from the interplay between meteorological forecast variability, uncertainties from the population number estimation, and the spread of the impact
functions. In general, it is difficult to exhaustively characterise all uncertainties of input data25,28, and one should always carefully consider which elements are needed for the purpose
of the modelling task30. For impact-based forecast for anticipatory action planning, it is important to understand what individual events are plausible (e.g., will the TC hit the Eastern or
the Western part of the area?) instead of summary values (e.g., the median outcome of the TC forecast impacts is in both the Eastern and Western areas), where the latter is more common in
risk assessment. The chosen uncertainties in this study reflect this purpose where data is available. For the hazard, we consider the individual outcomes from each forecast ensemble member
as a representation of the meteorological uncertainty. For the vulnerability, qualitative impact data are often not available. To the best of our knowledge, IDMC provides the most
comprehensive and consistent source of global displacement data, and hence we use the IDMC displacement data for impact function calibration for the global consistency. At the same time,
despite many efforts from IDMC dedicated to systematically collecting and harmonising displacement data, data gaps and inconsistencies remain31. Some of these gaps stem from the variability
in reporting, for instance, displacement can be described as “homeless” or “moved”, or included in the category “directly affected” 31. In order to capture all the plausible displacement
outcomes from the TC events, we consider the ensemble of individual impact functions calibrated to best represent each recorded displacement event in the used IDMC database (per region).
This is in contrast to calibrating one optimal impact function that minimises errors for all past events at once (which is the impact function used in estimating displacement with
meteorological uncertainty only in subsection “Impact forecast for TC Yasa” of section “Results” and Fig. 2), and subsequently sample the uncertainty from the confidence interval of the
optimisation process itself. The latter is typically done in climate risk assessment studies and is best suited to estimate the average impact over many events and its uncertainty. For the
exposures, due to a lack of uncertainty information from the population data layer provider, we only vary homogeneously the total population, which can be understood as a proxy for the
movement of people in the area. Importantly, our analysis demonstrates the importance of performing a global uncertainty analysis (i.e., varying all input data simultaneously)28 since it is
the interplay of hazard, exposure and vulnerability that can lead to extreme scenarios (e.g., a strong TC hitting directly highly populated cities with large displacement vulnerability). To
understand how the relevant outputs depend on uncertain inputs we also perform a sensitivity analysis. Conducted at various lead times, we observe that in general meteorological
uncertainties play a dominant role at larger forecast lead times, whereas local vulnerability becomes more significant as TCs approach landfall and forecast uncertainty reduces. This finding
offers valuable insights for decision-makers, guiding their considerations when devising anticipatory actions at different points in time. In situations with longer lead times,
decision-makers may find it beneficial to seek expert input from meteorologists regarding the TC’s expected path and further cyclogenesis. On the other hand, as the TC approaches landfall,
involving individuals with local community knowledge becomes crucial for effective planning of anticipatory action. However, we also point out that there are large variations of the
sensitivity indices for both meteorological forecast and impact function (error bars in Fig. 5). There can be large lead times where the uncertainty contribution from the meteorological
forecast is little but the contribution from the uncertainty of impact functions is large, and vice versa for shorter lead times. For example, this was the case with TC Harold, which
impacted Vanuatu and displaced 80,000 people in April 2020 (further details in Supplementary Information S3). Such instances could be related to the fact that some meteorological situations
are easier to predict than others. Several sources of uncertainty are not represented in our model setup which if included could add more nuances to the previous conclusions from the
uncertainty and sensitivity analysis at different lead times. Here we model population displacement as a direct function of TC wind speeds, but we acknowledge that the relationship between
displacements and TC hazards is much more complex. On the physical hazard side, displacement occurs not only due to damages from wind, but also due to storm surge, flood and torrential rain.
However, since both surges and rain correlate to the TC wind speed32, the latter is often taken as a proxy of the overall TC hazard intensity33,34. In addition, our model does not resolve
the compounding risks if multiple hazards occur at the same time or in close succession. We choose the simpler single-hazard approach because TC wind speed can be provided by computationally
inexpensive models, and considering multiple and compound hazards strongly increases model complexity35. In estimating population distribution, it is crucial to recognise that locational
uncertainty stems from the accuracy of estimating densely populated clusters, which could be a significant parameter in uncertainty analysis. However, due to the lack of uncertainty
information, we only apply a uniform scaling uncertainty of the total population. This likely explains the low sensitivity to the exposures. It is imperative to understand that this does not
mean that the exposures’ uncertainties are not important. It only means that the impact forecast is not sensitive to homogeneous scaling of the population. The uncertainty and sensitivity
analysis can only represent what is used as input for the model28. Furthermore, our model only considers displacement as a direct impact of hazards. But people can also be displaced by
indirect impacts such as the loss of access to basic services (e.g., water provision) due to cascading failures of critical infrastructures18. Our impact functions are calibrated using the
IDMC recorded displacement attributed to TC events which include people who have been displaced due to indirect impacts. Maximum sustained wind speed is generally considered a good proxy for
compound damages afflicted by TCs32,34. Thus, the modelled total impact in the area affected by TC wind should reflect displacement from direct and indirect causes, but only in areas of
high wind speeds. The model could further be improved by incorporating TC sub-hazard footprints like rainfall and storm surge, or explicitly model displacement triggered by housing damages,
loss of livelihoods or loss of access to critical infrastructures36. The modelling of the impacts of compound hazards is an ongoing scientific effort (e.g. Rossi et al.36 and Stalhandske et
al.37). We further remark that for the planning of anticipatory action, not only the displaced people, but also the people that may decide to remain in the impacted area38, or even become
trapped there if escape routes are blocked39,40 may be of relevance. Those trapped are often the most vulnerable, but this effect is not captured in our analysis. Additionally, displacement
is only one of the numerous socio-economic impacts related to TCs. While forecasts of asset damages and other impacts are possible with approaches similar to ours, we focus on displacement
because it is arguably the most relevant impact for humanitarian agencies. Future research is encouraged to disentangle the complex nature of disaster displacement and incorporate them in
the next iterations of impact forecast systems, and to integrate additional types of impact for a more complete picture on impending disasters. Impact information provided on top of the
weather forecasts can support humanitarian communities to prepare properly before hazards strike and to act more efficiently. Our presented disaster impact forecast system is open-source and
based on solely open-access data which allows any actor even with limited resources to run the model at a low operational cost. Our code and implementation are also transferable to other
users (e.g. meteorological services, businesses, governments, individuals), hazards (e.g., floods, heatwaves, storm surges, drought), and impact types (e.g. displacement, mortality). We hope
the versatility of our work encourages the ongoing development of impact forecasting systems and fosters collaborative efforts between the risk modelling community and relevant
stakeholders. METHODS RISK ASSESSMENT PLATFORM CLIMADA We use the open-source probabilistic risk assessment platform CLIMADA v3.3 for our displacement impact forecast written in
Python13,41,42. CLIMADA is designed to simulate the interactions of weather- and climate-related hazards, the exposure of people or assets to the hazard, and the vulnerability of those
exposed people or assets in a globally consistent fashion. We use the tropical cyclone weather forecast as hazard with the population as exposure and their displacement vulnerability as an
impact function to predict the number of people at risk of displacement. We forecast the spatially explicit displacement number at 150 arc second resolution on land (~4 km) globally. TC
TRACK ENSEMBLE FORECASTS AND DERIVATION OF TC WINDSPEED AS HAZARD We take the maximum 1-minute sustained wind speed at 10 m above the ground as a proxy of the TC hazard intensity. The TC
wind speed is derived from the TC track forecast from the European Centre for Medium-Range Weather Forecasts Integrated Forecasting System (ECMWF-IFS)22,23, updated in real-time at 00:00 and
12:00 UTC every day. The ECMWF-IFS forecast products consist of a high-resolution deterministic run and 51 ensemble members. We consider only the 51 ensemble members for the displacement
risk impact forecast. Each TC track forecast consists of the forecasted positions, central pressure, ambient pressure, and maximum wind speed, available at 6-hour intervals for 240 h, which
we further interpolate to 1-hour intervals. The real-time TC forecast tracks are openly accessible from the ECMWF servers via file transfer protocol (FTP)43, and are piped to CLIMADA through
its tropical cyclone forecast modules. Past TC track forecasts from 2017 to 2020 are available from the THORPEX Interactive Grand Global Ensemble (TIGGE) project44. In total, there are 161
displacement events around the globe from the IDMC database with matched ECMWF TC archived forecast tracks2. The TC wind speed is then computed at a horizontal resolution of 150 arc seconds
on land and 1800 arc seconds on sea from the TC ensemble forecast tracks, based on the revised hurricane pressure–wind model by refs. 13,45. POPULATION EXPOSURE For the exposure dataset, we
take the spatially explicit representation of the world’s population from the Gridded Population of the World (GPW v4) dataset46, openly accessible from the Socioeconomic Data and
Application Center (SEDAC). The dataset is created by uniformly distributing the numbers of people from census data or population figures provided by national statistics offices at the
smallest administrative unit, without considering any ancillary sources. CALIBRATION OF IMPACT FUNCTIONS FOR HUMAN DISPLACEMENT The vulnerability of people to displacement by TCs is
represented by mathematical functions (named impact functions in the CLIMADA terminology, and often called vulnerability curves) relating the TC wind speed to the percentage of people
displaced at a location. Prior studies have focused on defining impact functions for TC-related building asset damages in terms of monetary values33,47,48. To the best of our knowledge,
there is no research yet on impact functions for displacement due to TCs. In general, the impact functions must always be calibrated for the data (exposure and hazard) at hand. This can be
done via numerical optimisation and may be complemented by empirical knowledge49 or expert knowledge14. Of central importance is choosing a calibration method suitable for the purpose of the
model. Here we assume a sigmoid curve functional form and calibrate a set of impact functions using the reported displacement data from the Internal Displacement Monitoring Centre (IDMC)
database (accessed in September 2022). This approach differs from previous approaches that used housing damage as a proxy to estimate displacement50 and which have been suspected to provide
a conservative estimate likely to underestimate displacement3. Our approach in particular strives to reproduce the total reported displacement per TC event. It thus implicitly includes the
effects of all TC sub-hazards (wind, flood, surge), preemptive measures such as evacuation, as well as regional and cultural differences. We calibrate the displacement impact functions to
the reported displacement data from the IDMC database for events recorded from 2008 to 2020. These reported numbers of displacements can vary from a few persons to millions in some
high-impact events2. IDMC systematically collects displacement data from governmental and non-governmental institutions post-disasters, and dedicated efforts to verify and harmonise data to
ensure their interoperability. All the data are curated in the global database aggregated at country level2, and the analyses are presented in the IDMC yearly Global Report on Internal
Displacement (GRID) (e.g., GRID 202351). In this study, we take the 394 events for which IDMC has reported new displacement due to TCs as the impact data for the impact function calibration.
The hazard events are represented by the 1-min maximum sustained wind speed at 10 m from surface derived from the corresponding historical TC tracks from the International Best Track
Archive for Climate Stewardship (IBTrACS)27, and the associated wind speed computed from the45 model as implemented in CLIMADA42. We note that there might exist biases in terms of the TC
hazard intensity between these historical tracks and the ECMWF forecast tracks52,53 (c.f., subsection “TC track ensemble forecasts and derivation of TC windspeed as hazard” of section
“Methods”) used for the impact forecast. The displacement obtained with ECMWF tracks (0.5-days lead-time) shows an overestimation as compared to both the reported displacement from IDMC and
the modelled displacement using IBTrACS, even though ECMWF tracks have a low-intensity bias. This is likely because the tracks from the ECMWF numerical weather predictions are simulated for
longer times than is recorded by IBTrACS (see detailed discussion in the supplementary information S4). Currently no bias correction is made, but this could be addressed in future iterations
of the model. As exposures, we use the global static population layer described in subsection “Population exposure” of section “Methods”. To account for the local adaptation to TC around
the globe, one would ideally calibrate the vulnerability at high resolution. However, due to the rare nature of TCs, reported displacement data is only available for a few locations
aggregated at country levels. To ensure sufficiently robust statistics, we group countries into ten different regions with at least 19 reported TCs and derive separate impact functions for
each one of them (see Fig. S1 in the supplementary material). We use the 9 regions defined by Eberenz et al.48, and we additionally separate Oceania into two regions: Australia and New
Zealand, and the Pacific Islands. The Pacific Islands region impact functions are used for the displacement estimation in Fiji due to TC Yasa in Figs. 2–4. We consider two optimisation
options for the calibration. In both cases, we use the same third-order sigmoid-type candidate function. For the results shown in Fig. 2 where we use a single deterministic impact function,
we calibrate one impact function per region which minimises the root mean square fraction (RMSF) error over all events in the region. This is designed to minimise the spread of relative
error of single events. For all the uncertainty results shown in Figs. 3–5 where we use a family of functions to capture uncertainty, we instead calibrate one impact function per historical
event by minimising the root mean square error. Thus, instead of one impact function per region, one obtains a bundle of impact functions per region. The details of the calibration and the
resulting set of functions are shown in the supplementary information. UNCERTAINTY AND SENSITIVITY To characterise the weather forecast uncertainty alone (c.f. Fig. 2), the impact is
computed separately for each of the forecast ensemble members with exposures and impact function remaining constant. In general, however, an impact forecast depends non-linearly on the three
input variables exposures, impact function and hazard. Therefore, uncertainties in each component may interact non-linearly. To capture these effects, we perform a global (as opposed to
one-at-the-time) uncertainty and sensitivity analysis28. The uncertainty analysis gives information on the spread of the output variables (here the number of displaced people). The
sensitivity analysis provides indices that subsume the sensitivity of a model output variable to the uncertainty of each input parameter25,28. To compute the global uncertainty distribution
and sensitivity indices (c.f. Figs. 2, and 4) we use a standard quasi-Monte-Carlo numerical simulation approach25,54. We consider in addition to the forecast ensemble members one uncertainty
parameter for the exposures and impact function each. For the exposures, we varied the total value of the population estimate uniformly between [80%, 120%]. For the uncertainty distribution
of the impact functions per region, we use all the impact functions in a region that lie within the 80% confidence interval around the calibrated best estimate (the light blue lines in Fig.
S2). We remark that this discrete uncertainty distribution covers a larger interval than the 80% interval of the error from the RMSF optimisation over all events in a region (c.f.
supplementary material). We argue that this larger uncertainty range better represents the variability in vulnerability and is thus well-suited for impact forecasting. Using the Sobol
sampling algorithm26, we generate a total of more than 8000 samples of the three input parameters. For each sample, we compute the impact at each exposure location, and from this the total
number of displaced people. For the sensitivity analysis, the same samples are used to compute the first-order Sobol sensitivity indices as defined in Saltelli & Annoni29 for the total
number of displaced people, which characterise how much each individual input parameter contributes to the output metrics’ variance. REPORTING SUMMARY Further information on research design
is available in the Nature Portfolio Reporting Summary linked to this article. DATA AVAILABILITY The real-time updated TC tracks forecast are distributed under the Creative Commons CC-4.0-BY
licence, accessible through the ECMWF website (https://www.ecmwf.int/en/forecasts/datasets/wmo-essential) and can be retrieved through the CLIMADA platform. The archived TC forecast tracks
are available through the TIGGE platform (https://www.ecmwf.int/en/forecasts/datasets/wmo-essential). For the observed historical TC tracks used for impact function calibration, the tracks
are distributed under the IBTrACS website (https://www.ncei.noaa.gov/products/international-best-track-archive). All tracks can be imported to the CLIMADA platform for calculating impacts.
The recorded displacement data are publicly accessible via IDMC displacement database (https://www.internal-displacement.org/database/displacement-data). The exposure population data is
obtained from the Gridded Population of the World (GPW v4) dataset, openly accessible from the Socioeconomic Data and Application Center (SEDAC) (https://doi.org/10.7927/H4JW8BX5). CODE
AVAILABILITY All code necessary to reproduce the analysis is made available on https://github.com/manniepmkam/TC_displacement_forecast and permanently stored at
https://zenodo.org/doi/10.5281/zenodo.13342825. REFERENCES * UN. _Guiding Principles on Internal Displacement_. 14.
https://www.ohchr.org/en/special-procedures/sr-internally-displaced-persons/international-standards (1998). * IDMC. Global Internal Displacement Database. IDMC.
https://www.internal-displacement.org/database. * Desai, B. et al. Addressing the human cost in a changing climate. _Science_ 372, 1284–1287 (2021). Article ADS PubMed CAS Google Scholar
* IFRC. _Forecast-Based Financing and Disaster Displacement: Acting Early to Reduce the Humanitarian Impacts of Displacement_.
https://reliefweb.int/report/world/forecast-based-financing-and-disaster-displacement-acting-early-reduce-humanitarian (2020). * WMO. _WMO Coordination Mechanism Implementation Plan_.
https://meetings.wmo.int/EC-76/_layouts/15/WopiFrame.aspx?sourcedoc=%7b66c26e13-2fa7-4dfe-a191-647a9e54a70f%7d&action=default (2023). * Coughlan de Perez, E. et al. Forecast-based
financing: an approach for catalyzing humanitarian action based on extreme weather and climate forecasts. _Nat. Hazards Earth Syst. Sci._ 15, 895–904 (2015). Article ADS Google Scholar *
Thalheimer, L., Simperingham, E. & Jjemba, E. W. The role of anticipatory humanitarian action to reduce disaster displacement. _Environ. Res. Lett._ 17, 014043 (2022). Article ADS
Google Scholar * WMO. _WMO Guidelines on Multi-Hazard Impact-Based Forecast and Warning Services_ (WMO, 2015). * Harris, A., Roebber, P. & Morss, R. A new verification approach? Using
coupled natural–human models to evaluate the impact of forecast errors on evacuations. _Bull. Am. Meteorol. Soc._ 104, E1166–E1178 (2023). Article ADS Google Scholar * Davidson, R. A. et
al. An integrated scenario ensemble-based framework for hurricane evacuation modeling: part 1—decision support system. _Risk Anal._ 40, 97–116 (2020). Article PubMed Google Scholar *
Blanton, B. et al. An integrated scenario ensemble-based framework for hurricane evacuation modeling:part 2—hazard modeling. _Risk Anal._ 40, 117–133 (2020). Article PubMed Google Scholar
* Morss, R. E. et al. Hazardous weather prediction and communication in the modern information environment. _Bull. Am. Meteorol. Soc._ 98, 2653–2674 (2017). Article ADS Google Scholar *
Aznar-Siguan, G. & Bresch, D. N. CLIMADA v1: a global weather and climate risk assessment platform. _Geosci. Model Dev._ 12, 3085–3097 (2019). Article ADS Google Scholar * Kam, P. M.
et al. Global warming and population change both heighten future risk of human displacement due to river floods. _Environ. Res. Lett._ 16, 044026 (2021). Article ADS Google Scholar *
Sauer, I. J. et al. Climate signals in river flood damages emerge under sound regional disaggregation. _Nat. Commun._ 12, 2128 (2021). Article ADS PubMed PubMed Central CAS Google
Scholar * Lüthi, S. et al. Rapid increase in the risk of heat-related mortality. _Nat. Commun._ 14, 4894 (2023). Article ADS PubMed PubMed Central Google Scholar * Stalhandske, Z. et
al. Projected impact of heat on mortality and labour productivity under climate change in Switzerland. _Nat. Hazards Earth Syst. Sci. Discuss_. 1–20. https://doi.org/10.5194/nhess-2021-361
(2021). * Mühlhofer, E., Koks, E. E., Kropf, C. M., Sansavini, G. & Bresch, D. N. A generalized natural hazard risk modelling framework for infrastructure failure cascades. _Reliab. Eng.
Syst. Saf._ 234, 109194 (2023). Article Google Scholar * Meiler, S. et al. Intercomparison of regional loss estimates from global synthetic tropical cyclone models. _Nat. Commun._ 13,
6156 (2022). Article ADS PubMed PubMed Central CAS Google Scholar * Hülsen, S. et al. Global protection from tropical cyclones by coastal ecosystems—past, present, and under climate
change. _Environ. Res. Lett._ 18, 124023 (2023). Article ADS Google Scholar * Röösli, T., Appenzeller, C. & Bresch, D. N. Towards operational impact forecasting of building damage
from winter windstorms in Switzerland. _Meteorol. Appl._ 28, e2035 (2021). Article ADS Google Scholar * Buizza, R. The ECMWF ensemble prediction system. in _Predictability of Weather and
Climate_ (eds Hagedorn, R. & Palmer, T.) 459–488 (Cambridge University Press, Cambridge, 2006). * Isaksen, L. et al. Ensemble of data assimilations at ECMWF (ECMWF, 2010). * IFRC.
_Fiji/Pacific: Tropical Cyclone Yasa/Ana - DREF Final Report (MDRFJ005)_. https://reliefweb.int/report/fiji/fijipacific-tropical-cyclone-yasaana-dref-final-report-mdrfj005 (2022). * Kropf,
C. M. et al. Uncertainty and sensitivity analysis for probabilistic weather and climate-risk modelling: an implementation in CLIMADA v.3.1.0. _Geosci. Model Dev._ 15, 7177–7201 (2022).
Article ADS Google Scholar * Sobol′, I. M. Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates. _Math. Comput. Simul._ 55, 271–280 (2001). Article
MathSciNet Google Scholar * Knapp, K. R., Kruk, M. C., Levinson, D. H., Diamond, H. J. & Neumann, C. J. The international best track archive for climate stewardship (IBTrACS):
unifying tropical cyclone data. _Bull. Am. Meteorol. Soc._ 91, 363–376 (2010). Article ADS Google Scholar * Pianosi, F. et al. Sensitivity analysis of environmental models: a systematic
review with practical workflow. _Environ. Model. Softw._ 79, 214–232 (2016). Article Google Scholar * Saltelli, A. & Annoni, P. How to avoid a perfunctory sensitivity analysis.
_Environ. Model. Softw._ 25, 1508–1517 (2010). Article Google Scholar * Knüsel, B., Baumberger, C., Zumwald, M., Bresch, D. N. & Knutti, R. Argument-based assessment of predictive
uncertainty of data-driven environmental models. _Environ. Model. Softw._ 134, 104754 (2020). Article Google Scholar * IDMC. _Disaster Displacement - A Global Review, 2008-2018_. 54
https://www.internal-displacement.org/sites/default/files/publications/documents/201905-disaster-displacement-global-review-2008-2018.pdf (2019). * Czajkowski, J. & Done, J. As the wind
blows? Understanding hurricane damages at the local level through a case study analysis. _Weather Clim. Soc._ 6, 202–217 (2014). Article Google Scholar * Emanuel, K. Global warming effects
on U.S. hurricane damage. _Weather Clim. Soc._ 3, 261–268 (2011). Article Google Scholar * Gettelman, A., Bresch, D. N., Chen, C. C., Truesdale, J. E. & Bacmeister, J. T. Projections
of future tropical cyclone damage with a high-resolution global climate model. _Clim. Change_ 146, 575–585 (2017). Article ADS Google Scholar * Zscheischler, J. et al. Future climate risk
from compound events. _Nat. Clim. Change_ 8, 469–477 (2018). Article ADS Google Scholar * Rossi, L. et al. A new methodology for probabilistic flood displacement risk assessment: the
case of Fiji and Vanuatu. _Front. Clim_. 6. https://doi.org/10.3389/fclim.2024.1345258 (2024). * Stalhandske, Z. et al. Global multi-hazard risk assessment in a changing climate. _Sci. Rep._
14, 5875 (2024). Article ADS PubMed CAS Google Scholar * Farbotko, C. & McMichael, C. Voluntary immobility and existential security in a changing climate in the Pacific. _Asia Pac.
Viewp._ 60, 148–162 (2019). Article Google Scholar * Black, R., Arnell, N. W., Adger, W. N., Thomas, D. & Geddes, A. Migration, immobility and displacement outcomes following extreme
events. _Environ. Sci. Policy_ 27, S32–S43 (2013). Article Google Scholar * Zickgraf, C. Theorizing (im)mobility in the face of environmental change. _Reg. Environ. Change_ 21, 126 (2021).
Article Google Scholar * Bresch, D. N. & Aznar-Siguan, G. CLIMADA v1.4.1: towards a globally consistent adaptation options appraisal tool. _Geosci. Model Dev._ 14, 351–363 (2021).
Article ADS Google Scholar * Kropf, C. M. et al. CLIMADA-project/climada_python: v3.1.0. Zenodo https://doi.org/10.5281/zenodo.5947271 (2022). * ECMWF. _WMO Essential_ (ECMWF, 2018). *
Swinbank, R. et al. The TIGGE project and its achievements in: bulletin of the American Meteorological Society Volume 97 Issue 1 (2016). _Bull. Am. Meteorol. Soc._ 97, 49–67 (2016). Article
ADS Google Scholar * Holland, G. A Revised Hurricane pressure–wind model. _Mon. Weather Rev._ 136, 3432–3445 (2008). Article ADS Google Scholar * Doxsey-Whitfield, E. et al. Taking
advantage of the improved availability of census data: a first look at the gridded population of the world, version 4. _Pap. Appl. Geogr._ 1, 226–234 (2015). Article Google Scholar *
UNDRR. _Global Assessment Report on Disaster Risk Reduction 2015_. http://www.undrr.org/publication/global-assessment-report-disaster-risk-reduction-2015 (2015). * Eberenz, S., Lüthi, S.
& Bresch, D. N. Regional tropical cyclone impact functions for globally consistent risk assessments. _Nat. Hazards Earth Syst. Sci._ 21, 393–415 (2021). Article ADS Google Scholar *
Severino, L. G. et al. Projections and uncertainties of future winter windstorm damage in Europe. _EGUsphere_ 1–31. https://doi.org/10.5194/egusphere-2023-205 (2023). * IDMC. _Global
Disaster Displacement Risk - A Baseline for Future Work_. https://www.internal-displacement.org/publications/global-disaster-displacement-risk-a-baseline-for-future-work (2017). * IDMC.
_GRID 2023 | 2023 Global Report on Internal Displacement_. https://www.internal-displacement.org/global-report/grid2023/ (2023). * Chan, M. H. K., Wong, W. K. & Au-Yeung, K. C. Machine
learning in calibrating tropical cyclone intensity forecast of ECMWF EPS. _Meteorol. Appl._ 28, e2041 (2021). Article ADS Google Scholar * Aijaz, S., Kepert, J. D., Ye, H., Huang, Z.
& Hawksford, A. Bias correction of tropical cyclone parameters in the ecmwf ensemble prediction system in Australia. _Mon. Weather Rev._ 147, 4261–4285 (2019). Article ADS Google
Scholar * Herman, J. & Usher, W. SALib: an open-source python library for sensitivity analysis. _J. Open Source Softw._ 2, 97 (2017). Article ADS Google Scholar * Met Office.
_Cartopy: a Cartographic Python Library with a Matplotlib Interface_ (Met Office, 2010). * Hunter, J. D. Matplotlib: A 2D graphics environment. _Comput. Sci. Eng._ 9, 90–95 (2007). Article
Google Scholar * Saltelli, A. Making best use of model evaluations to compute sensitivity indices. _Comput. Phys. Commun._ 145, 280–297 (2002). Article ADS CAS Google Scholar Download
references ACKNOWLEDGEMENTS C.M.K. acknowledges funding from the European Union’s Horizon 2020 research and innovation program (grant agreement No. 820712). All authors would like to
acknowledge Sylvain Ponserre, Maria Teresa Miranda Espinosa and all colleagues from the Internal Displacement Monitoring Centre (IDMC) for providing guidance and insights on the displacement
data. FUNDING Open access funding provided by Swiss Federal Institute of Technology Zurich. AUTHOR INFORMATION AUTHORS AND AFFILIATIONS * Institute for Environmental Decisions, ETH Zürich,
Zurich, Switzerland Pui Man Kam, Fabio Ciccone, Chahan M. Kropf, Lukas Riedel, Christopher Fairless & David N. Bresch * Internal Displacement Monitoring Centre, Geneva, Switzerland Pui
Man Kam * Federal Office of Meteorology and Climatology MeteoSwiss, Zurich, Switzerland Chahan M. Kropf, Lukas Riedel & David N. Bresch Authors * Pui Man Kam View author publications You
can also search for this author inPubMed Google Scholar * Fabio Ciccone View author publications You can also search for this author inPubMed Google Scholar * Chahan M. Kropf View author
publications You can also search for this author inPubMed Google Scholar * Lukas Riedel View author publications You can also search for this author inPubMed Google Scholar * Christopher
Fairless View author publications You can also search for this author inPubMed Google Scholar * David N. Bresch View author publications You can also search for this author inPubMed Google
Scholar CONTRIBUTIONS P.M.K., C.F., C.M.K. and D.N.B. conceptualised the study. P.M.K. curated the TC forecast tracks from ECMWF to CLIMADA, conducted the formal analysis, undertook
visualisation, and wrote the manuscript. P.M.K. and F.C. contributed to the historical TC displacement impact forecast analysis from 2017 to 2020. C.M.K., L.R., C.F. and D.N.B. supervised
the analysis and contributed to the interpretation of the results. D.N.B. acquired funding. All authors reviewed and edited the manuscript. CORRESPONDING AUTHOR Correspondence to Pui Man
Kam. ETHICS DECLARATIONS COMPETING INTERESTS The authors declare no competing interests. PEER REVIEW PEER REVIEW INFORMATION _Nature Communications_ thanks the anonymous reviewers for their
contribution to the peer review of this work. A peer review file is available. ADDITIONAL INFORMATION PUBLISHER’S NOTE Springer Nature remains neutral with regard to jurisdictional claims in
published maps and institutional affiliations. SUPPLEMENTARY INFORMATION SUPPLEMENTARY INFORMATION PEER REVIEW FILE REPORTING SUMMARY RIGHTS AND PERMISSIONS OPEN ACCESS This article is
licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give
appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in
this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative
Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a
copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. Reprints and permissions ABOUT THIS ARTICLE CITE THIS ARTICLE Kam, P.M., Ciccone, F., Kropf, C.M. _et al._
Impact-based forecasting of tropical cyclone-related human displacement to support anticipatory action. _Nat Commun_ 15, 8795 (2024). https://doi.org/10.1038/s41467-024-53200-w Download
citation * Received: 29 November 2023 * Accepted: 04 October 2024 * Published: 10 October 2024 * DOI: https://doi.org/10.1038/s41467-024-53200-w SHARE THIS ARTICLE Anyone you share the
following link with will be able to read this content: Get shareable link Sorry, a shareable link is not currently available for this article. Copy to clipboard Provided by the Springer
Nature SharedIt content-sharing initiative