Changeset 6954
 Timestamp:
 May 5, 2009, 10:49:44 AM (15 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

anuga_work/publications/boxing_day_validation_2008/patong_validation.tex
r6953 r6954 25 25 %Abstract 26 26 \begin{abstract} 27 In this paper a new benchmark for tsunami model validation is proposed. The benchmark is based upon the 2004 Indian Ocean tsunami, which provides a uniquely large amount of observational data for model comparison. Unlike the small number of existing benchmarks, the proposed test validates all three stages of tsunami evolution  generation, propagation and inundation. Specifically we use geodetic measurements of the SumatraAndaman earthquake to validate the tsunami source, altimetry data from the JASONsatellite to test open ocean propagation, eyewitness accounts to assess near shore propagation and a detailed inundation survey of Patong Bay, Thailand to compare model and observed inundation. Furthermore we utilise this benchmark to further validate the hydrodynamic modelling tool \textsc{anuga} which is used to simulate the tsunami inundation and run raininduced floods.27 In this paper a new benchmark for tsunami model validation is proposed. The benchmark is based upon the 2004 Indian Ocean tsunami, which provides a uniquely large amount of observational data for model comparison. Unlike the small number of existing benchmarks, the proposed test validates all three stages of tsunami evolution  generation, propagation and inundation. Specifically we use geodetic measurements of the SumatraAndaman earthquake to validate the tsunami source, altimetry data from the \textsc{jason} satellite to test open ocean propagation, eyewitness accounts to assess near shore propagation and a detailed inundation survey of Patong Bay, Thailand to compare model and observed inundation. Furthermore we utilise this benchmark to further validate the hydrodynamic modelling tool \textsc{anuga} which is used to simulate the tsunami inundation and run raininduced floods. 28 28 \end{abstract} 29 29 … … 134 134 135 135 \subsubsection{Buildings and Other Structures} 136 FIXME(John): Complete 137 Where did elevation data come from??? 136 Human made build and structures can significantly effect tsunamni inundation. The location and size and number of floors of the buildings in Patong Bay were extracted from a GIS data set~\cite{}. The heights of these buildings was estimated assuming that each floor was 3m high. 138 137 139 138 \subsubsection{Eyewitness Accounts} … … 168 167 \subsection{Generation}\label{sec:modelGeneration} 169 168 170 There are various approaches to modelling the expected crustal deformation from an earthquake at depth. Most approaches model the earthquake as a dislocation in linear, elastic medium. Here we use the method of~\cite{wang03}. One of the main advantages of their method is that it allows the dislocation to be located in a stratified linear elastic halfspace with an arbitary number of layers. Other methods (such as those based on Okada's equations) can only model the dislocation in a homogeneous elastic half space, or can only include a limited number of layers, and thus cannot model the effect of the depth dependence of the elasticity of the Earth~\cite{wang03}. The original versions of the codes described here are available from http://www.iamg.org/CGEditor/index.htm. The first program, \textsc{EDGRN}, calculates elastic Green's function for a set of point sources at a regular set of depths out to a specified distance. The equations controlling the deformation are solved by using a combination of Hankel's transform and Wang et al's implementation of the ThomsonHaskell propagator algorithm~\cite{wang03}. Once the Green's functions are calculated we use a slightly modified version of \textsc{EDCMP} to calculate the sea floor deformation for a specific subfault. This second code discretises the subfault into a set of unit sources and sums the elastic Green's functions calculated from \textsc{EDGRN} for all the unit sources on the fault plane in order to calculate the final static deformation caused by a two dimensional dislocation along the subfault. This step is possible because of the linearity of the governing equations. For this study, we have made minor modifications to \textsc{EDCMP} in order for it to output in a file format compatible with the propagation code in the following section but it is otherwise the similar to the original code.169 There are various approaches to modelling the expected crustal deformation from an earthquake at depth. Most approaches model the earthquake as a dislocation in linear, elastic medium. Here we use the method of~\cite{wang03}. One of the main advantages of their method is that it allows the dislocation to be located in a stratified linear elastic halfspace with an arbitary number of layers. Other methods (such as those based on Okada's equations) can only model the dislocation in a homogeneous elastic half space, or can only include a limited number of layers, and thus cannot model the effect of the depth dependence of the elasticity of the Earth~\cite{wang03}. The original versions of the codes described here are available from \url{http://www.iamg.org/CGEditor/index.htm}. The first program, \textsc{EDGRN}, calculates elastic Green's function for a set of point sources at a regular set of depths out to a specified distance. The equations controlling the deformation are solved by using a combination of Hankel's transform and Wang et al's implementation of the ThomsonHaskell propagator algorithm~\cite{wang03}. Once the Green's functions are calculated we use a slightly modified version of \textsc{EDCMP} to calculate the sea floor deformation for a specific subfault. This second code discretises the subfault into a set of unit sources and sums the elastic Green's functions calculated from \textsc{EDGRN} for all the unit sources on the fault plane in order to calculate the final static deformation caused by a two dimensional dislocation along the subfault. This step is possible because of the linearity of the governing equations. For this study, we have made minor modifications to \textsc{EDCMP} in order for it to output in a file format compatible with the propagation code in the following section but it is otherwise the similar to the original code. 171 170 172 171 In order to calculate the crustal deformation using these codes we thus need to have a model describing the variation in elastic properties with depth and a slip model of the earthquake to describe the dislocation. The elastic parameters used for this study are the same as those in Table 2 of \cite{burbidge08}. For the slip model, there are many possible models for the 2004 AndamanSumatran earthquake to choose from ~\cite{chlieh07,asavanant08,arcas06,grilli07,ioualalen07}. Some are determined from various geological surveys of the site, others solve an inverse problem which calibrates the source based upon the tsunami wave signal, the seismic signal and/or the runup. The source parameters used here to simulate the 2004 Indian Ocean Tsunami were taken from model GM9.15 from~\cite{chlieh07}. This model was created by inversion of wide range of geodetic and seismic data. The slip model consists of 686 20km x 20km subsegments each with a different slip, strike and dip angle. The dip subfaults go from $17.5^0$ in the north and $12^0$ in the south. Refer to~\cite{chlieh07} for a detailed discussion of this model and its derivation. Note that the geodetic data used in the validation was also included by~\cite{chlieh07} in the inversion used to find GM9.15, thus the validation is not completely independent. However, a successful validation would still show that the crustal deformation and elastic properties model used here is at least as valid as the one used in~\cite{chlieh07} and can reproduce the observations just as accurately. … … 221 220 \includegraphics[width=5.0cm,keepaspectratio=true]{ursgaDomain.jpg} 222 221 \includegraphics[width=5.0cm,keepaspectratio=true]{extent_of_ANUGA_model.jpg} 223 \caption{Computational domain of the ursga simulation (left) and the \textsc{anuga} simulation (right s). FIXME: Add lat longs to anuga and make fig for ursga.Show where ANUGA domain fits in ursgaDomain}222 \caption{Computational domain of the ursga simulation (left) and the \textsc{anuga} simulation (right). FIXME: Show where ANUGA domain fits in ursgaDomain} 224 223 \label{fig:computational_domain} 225 224 \end{center} 226 225 \end{figure} 227 226 228 The domain was discretised into approximately ...,000 triangles. The resolution of the grid was increased in certain regions to efficiently increase the accuracy of the simulation. The grid resolution ranged between a maximum triangle area of $...\times 10^5$ m$^2$ near the Western ocean boundary to $...$ m$^2$ in the small regions surrounding the inundation region in Patong Bay. Due to a lack of available data, friction was set to a constant throughout the computational domain. For the reference simulation a Manning's coefficient of 0.01 was chosen based upon previous numerical experiments conducted by the authors (FIXME: Citation Tom Baldock?? Or Duncan??).227 The domain was discretised into approximately 350,000 triangles. The resolution of the grid was increased in certain regions to efficiently increase the accuracy of the simulation. The grid resolution ranged between a maximum triangle area of $1\times 10^5$ m$^2$ near the Western ocean boundary to $3$ m$^2$ in the small regions surrounding the inundation region in Patong Bay. Due to a lack of available data, friction was set to a constant throughout the computational domain. For the reference simulation a Manning's coefficient of 0.01 was chosen based upon previous numerical experiments conducted by the authors (FIXME: Citation Tom Baldock?? Or Duncan??). 229 228 230 229 The boundary condition at each side of the domain towards the south and the north where no data was available was chosen as a transmissive boundary condition effectively replicating the time dependent wave height present just inside the computational domain. Momentum was set to zero. Other choices include applying the mean tide value as a Dirichlet type boundary condition but experiments as well as the result of the verification reported here showed that this approach tends to under estimate the tsunami impact due to the tempering of the wave near the side boundaries. 231 230 231 During the \textsc{anuga} simulation the tide was kept constant at $0.80$m. This value was chosen to correspond to the tidal height specified by the Thai Navy tide charts (\url{http://www.navy.mi.th/hydro/}) at the time the tsunami arrived at Patong Bay. Although the tsunami propagated for approximately 3 hours before it reach Patong Bay, the period of time during which the wave propagated through the \textsc{anuga} domain is much smaller. Consequently the assumption of constant tide height is reasonable 232 232 233 FIXME(John): Need a commentary on the dynamics of what is being observed and whether it aligns with eye witness observations. 233 Both the URS model and the \textsc{anuga} inundation model shows that the event comprises a train of waves some with preceding drawdown effects (ADD details of waveform with a graph from URL and a gauge from \textsc{anuga} and discuss).234 Both the URS model and the \textsc{anuga} inundation model shows that the event comprises a train of waves some with preceding drawdown effects Add details of waveform with a graph from URL and a gauge from \textsc{anuga} and discuss. This will come from the work of Kristy (plots of water depth and speed at two locations on beach and one of shore) and Richard (approximating water depth and wave speeds from videos). 234 235 235 236 Maximum onshore inundation elevation was simulated throughout the entire Patong Bay region. Figure~\ref{fig:inundationcomparison1cm} shows very good agreement between the measured and simulated inundation. The \textsc{anuga} simulation determines a region to be inundated if at some point in time it was covered by at least 1cm of water. This precision in field measurements is impossible to obtain. The inundation boundary is determined by observing water marks and other signs left by the receding waters. The precision of the observed inundation map is, most likely, at least an order of magnitude worse than the \textsc{anuga} simulation. The simulated inundation based upon a 10cm threshold is shown in Figure~\ref{fig:inundationcomparison1cm}. An inundation threshold of 10cm was selected for all future simulations to reflect the likely accuracy of the survey and subsequently facilitate a more appropriate comparison between the modelled and observed inundation area. … … 248 249 Here we introduce the measure 249 250 \begin{equation} 250 A _{in}=\frac{A_m\cap A_o}{A_o}251 A(I_{in})=\frac{A(I_m\cap I_o)}{A(I_o)} 251 252 \end{equation} 252 to quantify the fraction of the observed inundation area $A_o$ captured by the model $A_m$. Another useful measure is the fraction of the modelled inundation area that falls outside the observed inundation area given by the formula253 to quantify the fraction of the area $A(I_{in})$ of observed inundation region $I_o$ captured by the model $I_m$. Another useful measure is the fraction of the modelled inundation area that falls outside the observed inundation area given by the formula 253 254 \begin{equation} 254 A _{out}=\frac{A_m\setminus (A_m\cap A_o)}{A_o}255 A(I_{out})=\frac{A(I_m\setminus (I_m\cap I_o))}{A(I_o)} 255 256 \end{equation} 256 257 These values for the two aforementioned simulations are given in Table~\ref{table:inundationAreas} … … 281 282 \begin{center} 282 283 \label{table:inundationAreas} 283 \caption{$A _{in}$ and $A_{out}$ of the reference simulation and all sensitivity studies}284 \caption{$A(I_{in})$ and $A(I_{out})$ of the reference simulation and all sensitivity studies} 284 285 \begin{tabular}{ccc} 285 286 \hline 286 & $A _{in}$ & $A_{out}$ \\287 & $A(I_{in})$ & $A(I_{out})$ \\ 287 288 \hline\hline 288 Reference & Ã &\\289 Reference & 0.76 & 0.22\\ 289 290 Min. Friction & Ã & \\ 290 291 Max. Friction & Ã & \\ … … 296 297 \end{center} 297 298 \end{table} 298 299 FIXME(Ole): It would be nice if we could be a little more quantitative  e.g. along the lines of the MISG study that John and Jane participated in. Thoughts anyone?300 299 301 300 %================Section===========================
Note: See TracChangeset
for help on using the changeset viewer.