Merge pull request #2 from MAMoises/main

flooding

wikiIA/index.md

@@ -21,7 +21,7 @@ In each section you will find a list of developments, projects, explorations mad
 TITLE | AUTHOR | GROUP |  MORE INFO
 ------|--------|-------|-----------
 Tsusy | Albert Gallego | Costas (A3) |  [description sheet](sheets/tsusy.md)
-
+Machine learning for coastal flooding | Moisés Álvarez | ADAPTA (A2) | [description sheet](sheets/flooding.md)
 
 ### Deep Learning
 TITLE | AUTHOR | GROUP |  MORE INFO

wikiIA/sheets/flooding.md

@@ -0,0 +1,135 @@
+# Machine learning methods for coastal flooding
+## Coastal flooding modelling strategies
+The choice of the modelling strategy depends on the spatio-temporal scale to be solved. For spatial scales of few hundreds of meters, the direct option is to solve the surfzone hydrodynamics considering waves and water levels and subsequent overland flow using direct numerical simulation. However, for larger spatial scales, the wave and water level transformation problem is uncoupled from the overland flooding simulation to speed up the computations.
+
+<div align="center">
+  <figure>
+    <img src="../_static/images/flood_ms.png" width="50%">
+    <figcaption style="text-align: center;">Figure 1: Coastal flooding modelling strategies.</figcaption>
+  </figure>
+</div>
+<br/>
+
+For regional-scale applications, it is usual to solve the surfzone hydrodynamics using specific numerical tools to derive the coastal forcings for a reduced-complexity overland flood model as shown in Figure.
+<div align="center">
+  <figure>
+    <img src="../_static/images/uncoupled_methodology.png" width="50%">
+    <figcaption style="text-align: center;">Figure 2: Flow diagram for a regional-scale flood modeling applications.</figcaption>
+  </figure>
+</div>
+<br/>
+
+## Limitations of numerical flooding prediction
+However, in spite of the computational advantages of reduced-complexity process modeling, numerical modeling of coastal flooding has some limitations:
+* Computational cost
+* Compatibility with operational forecasting systems
+* Probabilistic analysis
+* Uncertainty analyses
+
+## Problem statement
+In order to overcome these limitations, machine learning is a great tool. However, the predictability problem implies that a spatial output (millions of points) needs to be predicted based on a reduced set of parameters describing the coastal forcings of the storm (wave height, period, direction and still water level).
+
+$$
+Y = F(X) \text{, where } Y \text{ is the spatial map and } X \text{ represents Hs, Tp, Dir, SWL.}
+$$
+
+<div align="center">
+  <figure>
+    <img src="../_static/images/multipoints.png" width="50%">
+    <figcaption style="text-align: center;">Figure 3: Flood prediction points.</figcaption>
+  </figure>
+</div>
+<br/>
+
+## Methodology
+The methodology consists of creating a surrogate statistical model based on a set of numerical flood simulations. It is divided in several steps:
+1. Numerical model setup
+2. Extreme event selection
+3. Numerical modeling
+4. Statistical model training
+5. Statistical model evaluation
+## 1 Numerical model setup
+A 2DH XBeach simulation in a 10x10m grid is setup considering both topobathymetry of the San Lorenzo beach in Gijón consisting of 45k numerical cells.
+<div align="center">
+  <figure>
+    <img src="../_static/images/mallatbati.png" width="50%">
+    <figcaption style="text-align: center;">Figure 4: Numerical model grid.</figcaption>
+  </figure>
+</div>
+<br/>
+
+## 2 Event selection
+20 historical extreme events are selected considering a peaks over treshold (POT) method. 100 events are evaluated by combining the selected historical cases with 5 SLR scenarios.
+<div align="center">
+  <figure>
+    <img src="../_static/images/eventselection.png" width="50%">
+    <figcaption style="text-align: center;">Figure 5: Event selection.</figcaption>
+  </figure>
+</div>
+<br/>
+
+## 3 Numerical simulation
+The numerical simulation of the 100 cases yields a training dataset consisting of 100 realizations of the 45k grid (representing the Y predictand variable) against the 100 wave and water level parameters representing the storms (X predictor variable). 
+<div align="center">
+  <figure>
+    <img src="../_static/images/floodsimulations.png" width="50%">
+    <figcaption style="text-align: center;">Figure 6: Some simulated flood events.</figcaption>
+  </figure>
+</div>
+<br/>
+
+
+## 4 Training the statistical model
+The statistical model consists of projecting the training predictand dataset (N=45000xM=100) into a reduced subset using principal component analysis (PCA). 
+
+$$
+Y_{N,M}=U_{N,N} \Delta_{N,M} V_{M,M}
+$$
+
+<div align="center">
+  <figure>
+    <img src="../_static/images/EOFS.png" width="50%">
+    <figcaption style="text-align: center;">Figure 7: Some EOFs.</figcaption>
+  </figure>
+</div>
+<br/>
+
+Then, a given flood map (Y) can be solved as a linear summation of fixed EOFs (U) multiplied by storm-dependent latent variables (PCs) 
+
+$$
+Y_{j}=\sum_{i=1}^N \alpha_{ji}U_j
+$$
+
+Then, the spatial characteristics of the problem is represented by invariang EOFs and the inference problem is restricted to finding the storm-dependent latent variables that weights every EOF as a function of the storm characteristics which much more doable than the initial problem. The latent variables are inferred using Gaussian processes:
+
+$$
+\alpha_{ji}=GP(\underbrace{X}_{Hs, Tp, Dir, SWL})
+$$
+
+A balanced solution in terms of accuracy is obtained when choosing the EOFs that represent 90% of the total variance:
+
+<div align="center">
+  <figure>
+    <img src="../_static/images/modos90.png" width="50%">
+    <figcaption style="text-align: center;">Figure 8:Metamodel results when considering the EOFs that capture 90% of the total variance.</figcaption>
+  </figure>
+</div>
+<br/>
+
+## 5 Model evaluation
+The surrogate GP model is evaluated against the true numerical solution and results highlight a slight underprediction of the total flooded area. 
+<div align="center">
+  <figure>
+    <img src="../_static/images/testmodelo.png" width="50%">
+    <figcaption style="text-align: center;">Figure 9: Model evaluation. </figcaption>
+  </figure>
+</div>
+<br/>
+
+## 6 Way forward
+In order to improve the results we are exploring the following:
+* Crop the training dataset to the inland area (main interest is overland flooding)
+* Improve the GP surrogate model
+* Test deep learning nonlinear projection techniques
+* Test deep learning multioutput techniques
+* Feel free to chat if you are curious or want to help!