Changeset 7463


Ignore:
Timestamp:
Sep 2, 2009, 2:15:02 PM (15 years ago)
Author:
ole
Message:

Added info from conversations with Richard.

Location:
anuga_work/publications/boxing_day_validation_2008
Files:
2 edited

Legend:

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

    r7451 r7463  
    7575used here to validate open ocean tsunami models.
    7676
    77 \subsubsection{Bathymetry Data}
     77\subsubsection{Bathymetry Data}\label{sec:bathymetry data}
    7878The bathymetry data used in this study was derived from the following
    7979sources:
    8080\begin{itemize}
    8181\item a two arc minute grid data set covering the Bay of Bengal,
    82   DBDB2, obtained from US Naval Research Labs;
     82  DBDB2, obtained from US Naval Research Labs
     83  (\url{http://www7320.nrlssc.navy.mil/DBDB2_WWW});
    8384\item a 3 second arc grid covering the whole of the Andaman Sea based
    84   on Thai Navy charts no. 45 and no. 362; and
     85  on Thai Navy charts no. 45 and no. 362; and 
     86  (FIXME (OLE): wait for DB's reply)
    8587\item a one second grid created from the digitised Thai Navy
    8688  bathymetry chart, no. 358, which covers Patong Bay and the
    87   immediately adjacent regions.
    88   (FIXME (Ole): How was the grid created from these digitised points?) RICHARD
     89  immediately adjacent regions.
     90  The gridding of data was performed using {\bf Intrepid}, a commercial
     91  geophysical processing package developed by Intrepid Geophysics. The
     92  gridding scheme employed the nearest neighbour algorithm followed by
     93  an application of minimum curvature akima spline smoothing.
     94  See \url{http://www.intrepid-geophysics.com/ig/manuals/english/gridding.pdf}
     95  for details on the Intrepid model. 
    8996\end{itemize}
    90 FIXME (Jane): Refs for all these. RICHARD
    9197
    9298These sets were combined via
     
    100106arc grid. A subsection of this region was then replaced by nine second
    101107data which was generated by sub-sampling the three second of arc grid from
    102 NOAA (FIXME (Jane): This was not mentioned in the dots above. RICHARD). It is an artificially generated data set which is a subset of the original data.
     108NOAA. It is an artificially generated data set which is a subset of the original data.
    103109
    104110A subset of the nine second grid was replaced by the three second
     
    110116A one second grid was used to approximate the bathymetry in Patong
    111117Bay. This elevation data was created from the digitised Thai
    112 Navy bathymetry chart, no 358. The digitised points are shown in
    113 Figure~\ref{fig:patong_bathymetry}. The contours were obtained from...(FIXME? RICHARD)
     118Navy bathymetry chart, no 358. The digitised points and contour lines
     119from this chart are shown in Figure~\ref{fig:patong_bathymetry}.
     120
    114121
    115122The sub-sampling of larger grids was performed by using {\bf resample},
    116 a Generic Mapping Tools (\textsc{GMT}) program (\cite{wessel98}). The
    117 gridding of data was performed using {\bf Intrepid}, a commercial
    118 geophysical processing package developed by Intrepid Geophysics. The
    119 gridding scheme employed the nearest neighbour algorithm followed by
    120 an application of minimum curvature akima spline smoothing.
    121 See \url{http://www.intrepid-geophysics.com/ig/manuals/english/gridding.pdf}
    122 for details on the Intrepid model.
     123a Generic Mapping Tools (\textsc{GMT}) program (\cite{wessel98}).
    123124
    124125
     
    153154data which is often not available. In the case of model validation
    154155high quality field measurements are also required. For the proposed
    155 benchmark a high resolution topography data set and a high quality inundation
    156  survey map from the
     156benchmark a high resolution (1 second) topography data set and a
     157tsunami inundation survey map from the
    157158Coordinating Committee Co-ordinating Committee for Geoscience Programmes
    158159in East and Southeast Asia (CCOP) (\cite{szczucinski06}) was obtained
     
    161162
    162163\subsubsection{Topography Data}
    163 A one second grid was used to approximate the topography in Patong
    164 Bay. This elevation data was again created from the digitised Thai
    165 Navy bathymetry chart, no 358.
    166 FIXME (Ole): I don't think so. The Navy chart is only offshore. RICHARD
    167 
    168  A visualisation of the elevation data
    169 set used in Patong Bay is shown in
    170 Figure~\ref{fig:patong_bathymetry}.
     164The 1 second onshore topography for Patong Beach provided by the CCOP was
     165merged with the nearshore 1 second bathymetry described in Section
     166\ref{sec:bathymetry data} to provide a seamless terrain model for the
     167bay and community as shown in Figure~\ref{fig:patong_bathymetry}.
     168
    171169
    172170\begin{figure}[ht]
     
    272270Flow rates were estimated using landmarks found in both videos and
    273271were found to be in the range of 5 to 7 metres per second (+/- 2 m/s)
    274 in the north and 0.5 to 2 metres per second (+/- 1 m/s) in the south.
    275 FIXME (Jane): How were these error bounds derived? RICHARD
     272in the north and 0.5 to 2 metres per second (+/- 1 m/s) in the
     273south\footnote{These error bounds were estimated from uncertainty in aligning the debris with building boundaries in the videos.}.
    276274Water depths could also
    277275be estimated from the videos by the level at which water rose up the
    278276sides of buildings such as shops. Our estimates are in the order of
    279 1.5 to 2.0 metres (+/- 0.5 m).
     2771.5 to 2.0 metres (+/- 0.5 m estimated error bounds).
    280278Fritz ~\cite{fritz06} performed a detailed
    281279analysis of video frames taken around Banda Aceh and arrived at flow
  • anuga_work/publications/boxing_day_validation_2008/results.tex

    r7451 r7463  
    8787shown in Figure~\ref{fig:computational_domain}.
    8888
    89 FIXME (Ole): I know that a nested ursga model was trialled for the
    90 end-to-end modelling. However, for the study done here, where models
    91 were coupled, I didn't think nested grids were used with URSGA -
    92 and certainly not down to 1 arc second. Can someone shed some light
    93 on this please? RICHARD
    94 
    9589\begin{figure}[ht]
    9690\begin{center}
     
    152146in the small regions surrounding the inundation
    153147region in Patong Bay. The coarse resolution was chosen to be
    154 commensurate with the model output from the \textsc{ursga} model
     148commensurate with the model output from the \textsc{ursga} model (FIXME (Ole): Richard says that the ursga model used all four grids which would mean that the
     149resolution at the ANUGA boundary was 1 second or about 30m.
     150This is not consistent with my memory and certainly not with us choosing a
     151resolution of 440m. John, do you remember what the spacing was between the
     152URSGA points? Did we weed them out or did we take them as they were?)
     153
    155154while the latter was chosen to match the available resolution of topographic
    156155data and building data in Patong city.
Note: See TracChangeset for help on using the changeset viewer.