{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Regional tropical cyclone impact functions\n",
"\n",
"### CITATION\n",
"Eberenz S, Lüthi S, Bresch DN. Regional tropical cyclone impact functions for globally consistent risk assessments. Natural Hazards and Earth System Sciences. 2021 Jan 29;21(1):393-415.\n",
"\n",
"### LINK\n",
"https://nhess.copernicus.org/articles/21/393/2021/\n",
"\n",
"### NOTES\n",
"The data is used to provide a 'generic asset' vulnerability distribution. The fitted results are copied into Table_A2_Impact_Function_Slope.csv. \n",
"(a) the global default impact function (uncalibrated), (b) calibrated by optimizing RMSF, and (c) calibrated by optimizing TDR. The regions NA1 to WP4 are defined in Table A1. The row “combined” summarizes results for all regions combined based on the regionalized calibration; the row “global calibration” is based on one unified global calibration based on all matched TC 473 events. RMSF: root-mean-squared fraction; TDR: total damage ratio.\n",
"\n",
"The nine calibration regions are the Caribbean with Central America and Mexico (NA1), the USA and Canada (NA2), North Indian Ocean (NI), Oceania with Australia (OC), South Indian Ocean without Australia (SI), South East Asia (WP1), the Philippines (WP2), mainland China (WP3), and the north West Pacific (WP4) (see Fig. 3 and Table A1).\n",
"\n",
"The functional form of the damage function is:\n",
"\n",
"$$\n",
"d(v) = v_n^3 / (1 + v_n^3) \\\\\n",
"v_n = \\frac{\\max(v - v_t, 0)}{v_h - v_t}\n",
"$$\n",
"\n",
"$v_t$ is the threshold wind speed of 25.7 ms$^-1$. $v_h$ is the 'half' speed at which a damage ratio of 50% occurs.\n",
"\n",
"A simple use of the data is provided below. For China Mainland, North America and South America, the $v_h$ parameter is taken to be the mean of the RMSF and TDR fits. Europe and Generic curves are mapped to North America. For Asia the global default of 74.7 m/s is retained."
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import plotly.io\n",
"from plotly.subplots import make_subplots\n",
"from physrisk.vulnerability_models.configuration.config_builders import (\n",
" ConfigBuilderEberenzWind,\n",
")\n",
"\n",
"plotly.io.renderers.default = \"notebook\"\n",
"\n",
"builder = ConfigBuilderEberenzWind()\n",
"config_items = builder.build_config()"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
" \n",
" "
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/html": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"fig1 = make_subplots(rows=1, cols=1)\n",
"for item in config_items:\n",
" fig1.add_scatter(\n",
" x=item.points_x, y=item.points_y, row=1, col=1, name=item.asset_identifier\n",
" )\n",
"fig1.update_xaxes(title=\"Wind speed (m/s)\", title_font={\"size\": 14}, row=1, col=1)\n",
"fig1.update_yaxes(\n",
" title=\"Damage as fraction of TIV\", title_font={\"size\": 14}, row=1, col=1\n",
")\n",
"fig1.update_layout(legend=dict(orientation=\"h\", y=-0.25))\n",
"fig1.update_layout(margin=dict(l=20, r=20, t=20, b=20))"
]
}
],
"metadata": {
"kernelspec": {
"display_name": ".venv",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.10.15"
},
"orig_nbformat": 4
},
"nbformat": 4,
"nbformat_minor": 2
}