Changeset 6950


Ignore:
Timestamp:
May 5, 2009, 9:55:48 AM (16 years ago)
Author:
jakeman
Message:

John Jakeman: updated patong_validation paper with additional references

Location:
anuga_work/publications/boxing_day_validation_2008
Files:
1 added
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/publications/boxing_day_validation_2008/patong_validation.tex

    r6943 r6950  
    200200
    201201\subsection{Propagation}
    202 The deformation results described in Section~\ref{modelGeneration} was then used to provide an profile of the initial ocean surface displacement. This wave was used as an initial condition for \textsc{ursga} and was propagated the tsunami throughtout the Bay of Bengal. The rectangular computational domain of the largest grid extended from 90$^0$ to 100$^0$East and 0 to 15$^0$North and contained 1335$\times$1996 finite difference points. Inside this grid, a nested sequence of grids was used. The grid resolution of the nested grids went from 27 arc seconds in the coarsest grid, down to 9 arc seconds in the second grid, 3 arc seconds in the thrid grid and finally 1 arc second in the finest grid near Patong. The computational domain is shown in Figure\ref{gif:ursgaDomain}.
     202The deformation results described in Section~\ref{sec:modelGeneration} was then used to provide an profile of the initial ocean surface displacement. This wave was used as an initial condition for \textsc{ursga} and was propagated the tsunami throughtout the Bay of Bengal. The rectangular computational domain of the largest grid extended from 90$^0$ to 100$^0$East and 0 to 15$^0$North and contained 1335$\times$1996 finite difference points. Inside this grid, a nested sequence of grids was used. The grid resolution of the nested grids went from 27 arc seconds in the coarsest grid, down to 9 arc seconds in the second grid, 3 arc seconds in the thrid grid and finally 1 arc second in the finest grid near Patong. The computational domain is shown in Figure\ref{fig:computational_domain}.
    203203
    204204Figure \ref{fig:jasonComparison} provides a comparison of the \textsc{ursga} predicted surface elevation with the JASON satellite altimetry data. The \textsc{ursga} model replicates the amplitude and timing of the first peak and trough well. However the model does not resolve the double peak of the first wave. Also note that the \textsc{ursga} model prediction of the ocean surface elevation becomes out of phase with the JASON data at 3 to 7 degrees latitude. Chlieh et al~\cite{chlieh07} also observe these misfits and suggest it is caused by a reflected wave from the Aceh Peninsula that is not resolved in the model due to insufficient resolution of the computational mesh and bathymetry data. This is also a limitation of the model presented here, but probably could be improved by nesting grids near Aceh.
     
    218218\begin{center}
    219219%\includegraphics[width=5.0cm,keepaspectratio=true]{extent_of_ursga_model.jpg}
     220\includegraphics[width=5.0cm,keepaspectratio=true]{ursgaDomain.jpg}
    220221\includegraphics[width=5.0cm,keepaspectratio=true]{extent_of_ANUGA_model.jpg}
    221 \includegraphics[width=5.0cm,keepaspectratio=true]{extent_of_ANUGA_model.jpg}
    222 \caption{Computational domain of the ursga simulation (left) and the \textsc{anuga} simulation (rights). FIXME: Add lat longs to anuga and make fig for ursga}
     222\caption{Computational domain of the ursga simulation (left) and the \textsc{anuga} simulation (rights). FIXME: Add lat longs to anuga and make fig for ursga. Show where ANUGA domain fits in ursgaDomain}
    223223\label{fig:computational_domain}
    224224\end{center}
     
    232232Both 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).
    233233
    234 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:inundationcomparison10cm}. 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.
     234Maximum 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.
    235235
    236236\begin{figure}[ht]
  • anuga_work/publications/boxing_day_validation_2008/tsunami07.bib

    r6917 r6950  
    838838}
    839839
    840 @Article{amnon05,
     840@Article{ammon05,
    841841author = {C.J., Ammon and C. Ji and H. Thio and D. Robinson and Sidao Ni and V. Hjorleifsdottir and H. and T. Lay and S. Das and D. Helmberger and G. Ichinose and J. Polet and D. Wald},
    842842title = {Rupture Process of the 2004 Sumatra-Andaman Earthquake},
     
    10411041  pages =       "37-38",
    10421042  year =        2005,}
     1043
     1044@Article{burbidge08,
     1045author = {Burbidge, D. and Cummins, P.R. and Mleczko, R. and Thio, H.K.},
     1046title = {A Probabilistic Tsunami Hazard Assessment for Western Australia},
     1047journal = {Pure appl. geophys.},
     1048year = {2008},
     1049volume = {165},
     1050pages = {2059--2088},
     1051doi = {10.1007/s00024-008-0421-x},
     1052}
     1053
     1054@Article{thio08,
     1055author = {Thio, H.K. and Somerville, P. and Inchinose, G.},
     1056title = {Probabilistic analysis of tsunami hazards in southeast Asia},
     1057journal = {J. Earthquakes and Tsunami},
     1058year = {2008},
     1059volume = {1},
     1060pages = {119--137},
     1061}
     1062
     1063@Article{wessel98,
     1064author = {Wessel, P. and Smith, W.H.F.},
     1065title = {New, improved version of Generic Mapping Tools released},
     1066journal = {EOS trans. AGU},
     1067year = {1998},
     1068OPTvolume = {79},
     1069OPTpages = {579},
     1070}
     1071
     1072@Article{wang03,
     1073author = {Wang, R. and  Martin, F. L.  and Roth, F.},
     1074title = {Computation of deformation induced by earthquakes in a multi-layered crust – FORTRAN programs EDGRN EDCMP},
     1075journal = {Comp. and Geosc.},
     1076year = {2006},
     1077OPTvolume = {2003},
     1078OPTpages = {195--207},
     1079}
     1080
     1081@Article{stein07,
     1082author = {Stein, S. and Okal, E.A.},
     1083title = {Ultralong Period Seismic Study of the December 2004 Indian Ocean Earthquake and Implications for Regional Tectonics and the Subduction Process},
     1084journal = {Bulletin of the Seismological Society of America},
     1085year = {2007},
     1086volume = {97},
     1087number = {1A},
     1088pages = {S279–S295},
     1089OPTdoi = {10.1785/0120050617.},
     1090}
     1091
Note: See TracChangeset for help on using the changeset viewer.