Changeset 7522 for anuga_work/publications/boxing_day_validation_2008
- Timestamp:
- Sep 22, 2009, 6:09:32 PM (15 years ago)
- Location:
- anuga_work/publications/boxing_day_validation_2008
- Files:
-
- 9 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/publications/boxing_day_validation_2008/appendix.tex
r7480 r7522 16 16 provided by the \textsc{ursga} model. The left image shows the maximum 17 17 modelled depth while the right hand image shows the maximum modelled 18 flow velocities.}18 flow speeds.} 19 19 \label{fig:reference_model} 20 20 \end{figure} … … 34 34 friction value while the higher slows the flow and decreases the 35 35 inundation extent. Ideally, friction should vary across the entire 36 domain depending on terrain and vegetation , but thisis beyond the36 domain depending on terrain and vegetation. This, however, is beyond the 37 37 scope of this study.} 38 38 \label{fig:sensitivity_friction} … … 44 44 \end{center} 45 45 46 \caption{The maximal flow speeds for the same model parameterisations46 \caption{The maximal modelled flow speeds for the same model parameterisations 47 47 found in Figure \protect \ref{fig:sensitivity_friction}. 48 48 The reference flow speeds for a -
anuga_work/publications/boxing_day_validation_2008/conclusion.tex
r7521 r7522 1 1 \section{Conclusion} 2 2 This paper proposes a new field data benchmark for the 3 v erification of tsunami inundation models. Currently, there is a3 validation of tsunami inundation models. Currently, there is a 4 4 scarcity of appropriate validation datasets due to a lack of well-documented 5 5 historical tsunami impacts. The benchmark proposed here … … 20 20 and able to predict detailed inundation extents and dynamics with reasonable accuracy. 21 21 Model predictions matched well a detailed inundation survey 22 of Patong Bay, Thailand as well as altimetry data from the \textsc{jason} satellite,22 of Patong City, Thailand as well as altimetry data from the \textsc{jason} satellite, 23 23 eye-witness accounts of wave front arrival times and onshore flow speeds. 24 24 -
anuga_work/publications/boxing_day_validation_2008/data.tex
r7521 r7522 7 7 subsequent coastal field surveys of run-up and flooding, and 8 8 measurements of coseismic displacements as well as bathymetry from ship-based 9 expeditions , have now been made9 expeditions and high quality topographic data, have now been made 10 10 available. %~\cite{vigny05,amnon05,kawata05,liu05}. 11 11 … … 22 22 cause of any discrepancies between modelled and observed inundation. 23 23 Consequently, in this section we present data not only to facilitate 24 validation of inundation but to also aid the assessment of tsunami24 validation of inundation extent but to also aid the assessment of tsunami 25 25 generation and propagation. 26 26 … … 69 69 it encounters the shoreline bordering coastal regions. This period 70 70 of the tsunami evolution is referred to as the propagation stage. The 71 height and velocityof the tsunami is dependent on the local71 height and speed of the tsunami is dependent on the local 72 72 bathymetry in the regions through which the wave travels and the size 73 73 of the initial wave. This section details the bathymetry data needed … … 79 79 sources: 80 80 \begin{itemize} 81 \item a two arc minute grid data setcovering the Bay of Bengal,81 \item a two arc minute data grid covering the Bay of Bengal, 82 82 DBDB2, obtained from US Naval Research Labs 83 83 (\url{http://www7320.nrlssc.navy.mil/DBDB2_WWW}); 84 \item a 3 second arcgrid obtained directly from NOAA covering the84 \item a three arc second data grid obtained directly from NOAA covering the 85 85 whole of the Andaman Sea based on the 86 Smith \& Sandwell 2-minute86 Smith \& Sandwell two minute 87 87 dataset (\url{http://topex.ucsd.edu/WWW_html/srtm30_plus.html}), 88 88 coastline constrained using SRTM data (\url{http://srtm.csi.cgiar.org}) 89 89 as well as Thai Navy charts no.\ 45 and no.\ 362; and 90 \item a one second grid created from the digitised Thai Navy 91 bathymetry chart, no. 358, which covers Patong Bay and the 92 immediately adjacent regions. The digitised points and contour lines 93 from this chart are shown in Figure~\ref{fig:patong_bathymetry}. 94 The gridding of data was performed using \textsc{Intrepid}, a commercial 95 geophysical processing package developed by Intrepid Geophysics. The 96 gridding scheme employed the nearest neighbour algorithm followed by 97 an application of minimum curvature akima spline smoothing. 98 See \url{http://www.intrepid-geophysics.com/ig/manuals/english/gridding.pdf} 99 for details on the Intrepid gridding scheme. 90 \item Thai Navy chart no.\ 358 providing water depths in Patong Bay. 100 91 \end{itemize} 101 92 102 These sets were combined via 103 interpolation and resampling to produce four nested grids 104 which are relatively coarse in the deeper water and 105 progressively finer as the distance to 106 Patong Beach decreases as shown in Figure~\ref{fig:nested_grids}. 107 108 The coarsest 109 bathymetry was obtained by interpolating the DBDB2 grid to a 27 second 110 arc grid. A subsection of this region was then replaced by nine second 111 data which was generated by sub-sampling the three second of arc grid from 112 NOAA. It is an artificially generated data set which is a subset of the original data. 113 114 A subset of the nine second grid was replaced by the three second 115 data. Finally, the one second grid was used to approximate the 116 bathymetry in Patong Bay. Any 117 points that deviated from the general trend near the boundary were 118 deleted as a quality check. 119 120 A one second grid was used to approximate the bathymetry in Patong 121 Bay. This bathymetry data was created from the digitised Thai 122 Navy bathymetry chart, no 358. 123 124 93 These data sets were combined via gridding, interpolation and resampling to produce 94 four nested grids which are relatively coarse in the deeper water and 95 progressively finer as the distance to shore Patong Beach decreases as 96 shown in Figure~\ref{fig:nested_grids}. This progression was chosen 97 to match model resolution requirements according to the principle that 98 shallow water flows are more sensitive to variations in elevation data 99 than deep water flows. Consequently, the elevation data in shallow 100 waters and on-shore need to be resolved better than elevation data 101 further off-shore. 102 103 The coarsest bathymetry was obtained by interpolating the DBDB2 grid 104 to a 27~second arc grid. A subsection of this region was then replaced 105 by nine second data which was generated by sub-sampling the three 106 second of arc grid from NOAA. It is an artificially generated data set 107 which is a subset of the original data. A subset of the nine second 108 grid was replaced by the three second data. Finally, a one arc second 109 grid approximating the bathymetry in Patong Bay and the immediately 110 adjacent regions was created by digitising Thai Navy bathymetry chart, 111 no.\ 358. The digitised points and contour lines from this chart are 112 shown in Figure~\ref{fig:patong_bathymetry}. The gridding was 113 performed using \textsc{Intrepid}, a commercial geophysical processing 114 package developed by Intrepid Geophysics\footnote{ 115 See 116 \url{http://www.intrepid-geophysics.com/ig/manuals/english/gridding.pdf} 117 for details on the Intrepid gridding scheme.}. 118 Any points that deviated from the general trend near the boundary were 119 deleted through a quality control process. 125 120 The sub-sampling of larger grids was performed by using \textsc{resample}, 126 121 a Generic Mapping Tools (\textsc{GMT}) program \cite{wessel98}. … … 132 127 \end{center} 133 128 134 \caption{Nested bathymetry grids.} 129 \caption{Nested elevation grids of the Andaman Sea with 130 highest resolution at and around Patong Bay.} 135 131 \label{fig:nested_grids} 136 132 \end{figure} … … 138 134 \subsubsection{JASON Satellite Altimetry}\label{sec:data_jason} 139 135 During the 26 December 2004 event, the \textsc{jason} satellite tracked from 140 north to south and over the equator at 02:55 136 north to south and over the equator at 02:55~UTC nearly two hours 141 137 after the earthquake \cite{gower05}. The satellite recorded the sea 142 138 level anomaly compared to the average sea level from its previous five … … 162 158 Coordinating Committee Co-ordinating Committee for Geoscience 163 159 Programmes in East and Southeast Asia (CCOP) \cite{szczucinski06} 164 was obtained to validate model inundation. See also acknowledgements 165 at the end of this paper. In this section we also present eye-witness 160 was obtained to validate model inundation. In this section we also present eye-witness 166 161 accounts which can be used to qualitatively validate tsunami 167 162 inundation. … … 172 167 (described in Section \ref{sec:bathymetry data}) and from 1~m and 10~m 173 168 elevation contours provided by the CCOP. The 1~second terrain model 174 for the and community as shown in Figure~\ref{fig:patong_bathymetry}. 175 176 Two 1/3~second grids were created: One for the saddle point covering 177 Merlin and Tri Trang Beaches and one for Patong City and its immediate 178 shore area. These grids were based on the same data used for 169 for the community is shown in Figure~\ref{fig:patong_bathymetry}. 170 171 To provide increased resolution for the surveyed area, 172 two 1/3~second grids were created: One for the saddle point covering 173 Merlin and Tri Trang Beaches (separate survey patch to the left in 174 Figure~\ref{fig:patongescapemap}) 175 and one for Patong City and its immediate 176 shore area (main surveyed area in Figure~\ref{fig:patongescapemap}). 177 These grids were based on the same data used for 179 178 the 1~second data grid. The Patong city grid was further modified based on 180 179 satellite imagery to include the river and lakes towards the south of … … 188 187 \end{center} 189 188 190 \caption{3D vi sualisation of the elevation data set used for the nearshore propagation and inundation in Patong Bay showing189 \caption{3D view of the elevation data set used for the nearshore propagation and inundation in Patong City showing 191 190 digitised data points and contours as well as rivers and roads 192 191 draped over the data model.} … … 197 196 \subsubsection{Buildings and Other Structures} 198 197 Human-made buildings and structures can significantly affect tsunami 199 inundation. The footprint and number of floors of the 200 buildings in Patong Bay were extracted from the data provided by CCOP. 201 The heights of these 202 buildings were estimated assuming that each floor has a height of 3 m and they 203 were added to the topographic dataset. 198 inundation. The footprint and number of floors of the buildings in 199 Patong Bay were extracted from the data provided by CCOP. The heights 200 of these buildings were estimated assuming that each floor has a 201 height of 3~m and the resulting profiles were added to the topographic 202 dataset. The resulting elevation model and its interaction with one of 203 the tsunami waves can be seen in Figure~\ref{fig:anuga screenshot} in 204 Section~\ref{sec:inundation results}. 204 205 205 206 … … 225 226 226 227 \caption{Tsunami survey mapping the maximum observed inundation at 227 Patong beachcourtesy of the CCOP \protect \cite{szczucinski06}.}228 Patong City courtesy of the CCOP \protect \cite{szczucinski06}.} 228 229 \label{fig:patongescapemap} 229 230 \end{figure} … … 246 247 \end{center} 247 248 248 \caption{Location of time series extracted from the model output.}249 \caption{Location of time series extracted from the model output.} 249 250 \label{fig:gauge_locations} 250 251 \end{figure} … … 255 256 (Comfort Resort) and 256 257 \url{http://www.archive.org/download/tsunami_patong_beach/tsunami_patong_beach.wmv} 257 (Novotel) }258 (Novotel).} 258 259 %http://wizbangblog.com/content/2005/01/01/wizbang-tsunami.php 259 which include footage of the tsunami in Patong Bay on the day260 of the 2004 Indian Ocean Tsunami. Both videos show an already inundated260 which include footage of the tsunami in Patong City on the day 261 of the 2004 Indian Ocean tsunami. Both videos show an already inundated 261 262 street. They also show what is to be assumed as the second 262 263 and third waves approaching and further flooding of the buildings and 263 street. 264 street. The first video is in the very north, filmed from what is 264 265 believed to be the roof of the Novotel Hotel marked ``north'' in Figure 265 266 \ref{fig:gauge_locations}. The second video is in the very south, … … 304 305 should reproduce the following behaviour: 305 306 \begin{itemize} 306 \item reproducethe inundation survey map in Patong city307 \item the inundation survey map in Patong city 307 308 (Figure~\ref{fig:patongescapemap}), 308 \item simulatea leading depression followed by two distinct crests309 of decreasing magnitude at the beach, and309 \item a leading depression followed by two distinct crests 310 of decreasing magnitude at the beach, 310 311 \item predict the water depths and flow speeds, at the locations of 311 312 the eye-witness videos, that fall within the bounds obtained from 312 the videos .313 \item reproducethe \textsc{jason} satellite altimetry sea surface314 anomalies (see Section~\ref{sec:data_jason}), 315 \item reproducethe vertical deformation observed in north-western313 the videos, 314 \item the \textsc{jason} satellite altimetry sea surface 315 anomalies (see Section~\ref{sec:data_jason}), and 316 \item the vertical deformation observed in north-western 316 317 Sumatra and along the Nicobar--Andaman islands (see 317 Section~\ref{sec:gen_data}) ,318 Section~\ref{sec:gen_data}). 318 319 \end{itemize} 319 320 320 Ideally, the model should also be compared to measured time series of321 wave heights and velocities but the authors are not aware of the321 Ideally, the model should also be compared to measured time series of 322 wave heights and flow speeds but the authors are not aware of the 322 323 availability of such data near Patong Bay. -
anuga_work/publications/boxing_day_validation_2008/introduction.tex
r7521 r7522 20 20 These models are typically used to predict quantities such as arrival 21 21 times, wave speeds and heights, as well as inundation extents 22 whichcan be used to develop efficient hazard mitigation plans. Physics based23 models combine observed seismic, geodetic and sometimes t sunamidata to22 that can be used to develop efficient hazard mitigation plans. Physics based 23 models combine observed seismic, geodetic and sometimes tide gauge data to 24 24 provide estimates of initial sea floor and ocean surface 25 25 deformation. The shallow water wave equations~\cite{george06}, … … 29 29 30 30 Inaccuracies in model prediction can result in inappropriate 31 evacuation plans and town zoning, which may result in loss of life and 32 large financial losses. Consequently tsunami models must undergo 31 evacuation plans, town zoning and land use planning, 32 which ultimately may result in loss of life and infrastructure. 33 Consequently tsunami models must undergo 33 34 sufficient end-to-end testing to increase scientific and community 34 35 confidence in the model predictions. … … 63 64 64 65 Currently, the extent of tsunami-related field data is limited. The 65 cost of tsunami monitoring programs a s well as66 cost of tsunami monitoring programs and the rarity of events as well as 66 67 bathymetry and topography surveys 67 68 prohibits the collection of data in many of the regions in which … … 96 97 al~\cite{synolakis08} to validate and verify tsunami models. 97 98 The benchmark proposed here allows evaluation of 98 model structure during all three distinct stages tsunami evolution. 99 model components during three distinct stages tsunami 100 evolution, namely generation, propagation and inundation. 99 101 It consists of geodetic measurements of the 100 102 Sumatra--Andaman earthquake that are used to validate the description … … 115 117 tailored accordingly. 116 118 117 Unlike the existing field benchmarks the proposed test facilitates119 Unlike the existing field benchmarks, the proposed test facilitates 118 120 localised and highly detailed spatially distributed assessment of 119 121 modelled inundation. To the authors knowledge it is also the first benchmark to … … 130 132 The numerical models used to simulate tsunami impact 131 133 are computationally intensive and high resolution models of the entire 132 evolution process will often take a number of days to133 run. Consequently, the uncertainty in model predictions is difficult to134 quantify as it would require a very large number of runs.134 evolution process will often require a number of days to 135 complete. Consequently, the uncertainty in model predictions is difficult to 136 quantify as it would require a very large number of simulations. 135 137 However, model uncertainty should not be ignored. Section 136 138 ~\ref{sec:sensitivity} provides a simple analysis that can 137 139 be used to investigate the sensitivity of model predictions to a number 138 of model parameters.140 of key model parameters. -
anuga_work/publications/boxing_day_validation_2008/method.tex
r7521 r7522 1 1 \section{Modelling the Event}\label{sec:models} 2 2 Numerous models are currently used to model and predict tsunami 3 generation, propagation and run-up. These range in solving different3 generation, propagation and inundation. These range in solving different 4 4 equations and employing different methodologies with some examples 5 5 being~\cite{titov97a,satake95,zhang08}. Here we introduce the … … 65 65 ~\cite{chlieh07,asavanant08,arcas06,grilli07,ioualalen07}. Some are 66 66 determined from various geological surveys of the site. Others solve 67 an inverse problem whichcalibrates the source based upon the tsunami67 an inverse problem that calibrates the source based upon the tsunami 68 68 wave signal, the seismic signal and/or even the run-up. 69 69 The source … … 73 73 range of geodetic and seismic data. The slip model consists 74 74 of 686~20~km~x~20~km subsegments each with a different slip, strike and dip 75 angle. The dip subfaults gofrom $17.5^\circ$ in the north and $12^\circ$ in75 angle. The dip subfaults range from $17.5^\circ$ in the north and $12^\circ$ in 76 76 the south. Refer to Chlieh et al~\cite{chlieh07} for a detailed 77 77 discussion of this model and its derivation. %Note that the geodetic … … 84 84 %accurately. 85 85 86 \subsection{ Deepwater propagation}\label{sec:modelPropagation}86 \subsection{Open water propagation}\label{sec:modelPropagation} 87 87 The \textsc{ursga} model described below was used to simulate the 88 88 propagation of the 2004 Indian Ocean tsunami across the open ocean, based on a … … 107 107 grid system described in Section \ref{sec:propagation data} where 108 108 the coarse grids were used in the open ocean and the finest 109 resolution grid was employed in the region closest to Patong bay.109 resolution grid was employed in the region closest to Patong City. 110 110 \textsc{Ursga} is not publicly available. 111 111 … … 115 115 unless an intricate sequence of nested grids is employed. In 116 116 comparison \textsc{anuga}, described below, is designed to produce 117 robust and accurate predictions of on-shoreinundation, but is less117 robust and accurate predictions of inundation, but is less 118 118 suitable for earthquake source modelling and large study areas because 119 119 it is based on projected spatial coordinates. Consequently, the -
anuga_work/publications/boxing_day_validation_2008/paper.tex
r7521 r7522 14 14 %----------title-------------% 15 15 \title{Benchmarking Tsunami Models using the December 2004 Indian 16 Ocean Tsunami and its Impact at Patong Bay}16 Ocean Tsunami and its Impact at Patong City, Thailand} 17 17 \titlerunning{A tsunami model benchmark} 18 18 … … 28 28 \email{john.jakeman@anu.edu.au} 29 29 \and 30 O. Nielsen \and R. Mleczko \and D. Burbidge \and K. Van Putten \and N. Horspool \at 30 O. Nielsen \and K. Van Putten \and R. Mleczko 31 \and D. Burbidge \and N. Horspool \at 31 32 Geoscience Australia, Canberra, \textsc{Australia} 32 33 } … … 53 54 propagation and a detailed inundation survey of Patong city, Thailand 54 55 to compare model and observed inundation. Furthermore we utilise this 55 benchmark to further validate the \textsc{ursga--anuga}modelling methodology56 benchmark to further validate the modelling methodology 56 57 used by Geoscience Australia to simulate the tsunami 57 58 inundation. Important buildings and other structures were incorporated 58 into the underlying computational mesh and shown to have a large59 into the underlying computational mesh and are shown to have a large 59 60 influence on inundation extent. 60 61 -
anuga_work/publications/boxing_day_validation_2008/results.tex
r7521 r7522 118 118 After propagating the tsunami in the open ocean using \textsc{ursga}, 119 119 the approximated ocean and surface elevation and horizontal flow 120 velocities were extracted and used to construct a boundary condition 120 speeds were extracted and used to construct a boundary condition 121 121 for the \textsc{anuga} model. The interface between the \textsc{ursga} 122 122 and \textsc{anuga} models was chosen to roughly follow the 100~m depth 123 contour along the west coast of Phuket Island. Data from the 124 3 second grid was decimated to match the resolution chosen in ANUGA. 125 The computational123 contour along the west coast of Phuket Island. Data from the 124 three~second grid which is approximately 30~m apart was decimated to 125 match the resolution chosen in \textsc{Anuga}. The computational 126 126 domain is shown in Figure~\ref{fig:computational_domain}. 127 127 … … 164 164 The boundary condition at each side of the domain towards the south 165 165 and the north where no information about the incident wave or 166 its velocity field is available 166 its velocity field is available from the \textsc{Ursga} model 167 167 was chosen as a transmissive 168 168 boundary condition, effectively replicating the time dependent wave … … 170 170 The velocity field on these boundaries was kept at 171 171 to zero during the simulation. Other choices include applying the mean tide value as a 172 Dirichlet boundary condition. But experiments as well as the172 Dirichlet boundary condition. Experiments as well as the 173 173 result of the verification reported here showed that this approach 174 174 tends to underestimate the tsunami impact due to the tempering of the … … 176 176 condition robustly preserves the wave. 177 177 178 During the \textsc{anuga} simulation the tide was kept constant at178 During the \textsc{anuga} simulation the tide was kept constant in the offshore region at 179 179 $0.80$ m. This value was chosen to correspond to the tidal height 180 180 specified by the Thai Navy tide charts … … 184 184 wave propagated through the \textsc{anuga} domain is much 185 185 smaller of the order of 2 hours. Consequently the assumption of constant tide height is 186 reasonable. 186 reasonable. The initial water level for the river was set to 0. 187 187 188 188 \subsection{Inundation}\label{sec:inundation results} 189 189 The \textsc{anuga} simulation described in the previous section and used to 190 190 model shallow water propagation also predicts 191 inundation. Maximum onshore inundation depth was computed from the model 192 throughout the entire Patong Bay region and used to generate 193 a measure of the inundated area. 191 inundation. Maximum onshore inundation depth was computed from the inundation model 192 and used to generate a measure of the inundated area. 194 193 Figure~\ref{fig:inundationcomparison1cm} (left) shows very good 195 194 agreement between the measured and simulated inundation. However, … … 252 251 \end{equation} 253 252 These values for the two aforementioned simulations are given in 254 Table~\ref{table:inundationAreas}. High value of both $\rho_{in}$ and $\rho_{out}$ indicate 253 Table~\ref{table:inundationAreas} along with results from the sensitivity analysis in 254 Section~\ref{sec:sensitivity}. High values of both $\rho_{in}$ and $\rho_{out}$ indicate 255 255 that the model overestimates the impact whereas low values of both quantities would indicate 256 256 an underestimation. A high value of $\rho_{in}$ combined with a low value of $\rho_{out}$ … … 267 267 flow, as well as uncertainties in the elevation data including effects of 268 268 erosion and deposition by the tsunami event. 269 The impacts of some of the model uncertainties are is investigated in269 The impacts of some of the model uncertainties are as investigated in 270 270 Section~\ref{sec:sensitivity}. 271 271 … … 274 274 As one aim of this paper is to provide a new benchmark for tsunami 275 275 inundation modelling we have made the datasets available 276 available on \textsc{Sourceforge} in \textsc{anuga}276 available on \textsc{Sourceforge} in the \textsc{anuga} 277 277 project (\url{http://sourceforge.net/projects/anuga}) under the directory 278 278 \url{validation\_data/patong-1.0}. … … 284 284 will need to run the validation scripts (\url{anuga\_validation/automated\_validation\_tests/patong\_beach\_validation}) which are part of the 285 285 \textsc{anuga} distribution also available from 286 Sourceforge\url{http://sourceforge.net/projects/anuga}.286 \textsc{Sourceforge} \url{http://sourceforge.net/projects/anuga}. 287 287 288 288 … … 291 291 \subsubsection{Arrival time} 292 292 The arrival time of the first wave took place between 9:55 and 10:05 as described in 293 Section~\ref{sec:eyewitness data}. The modelled arrival time at the beach is 10:02293 Section~\ref{sec:eyewitness data}. The modelled arrival time at the beach is around 10:02 294 294 as can be verified from the animation provided in 295 Section \ref{sec:inundation results} .295 Section \ref{sec:inundation results} or from Figure~\ref{fig:onshore_timeseries} below. 296 296 Subsequent waves of variable magnitude appear over the next two hours 297 297 approximately 20-30 minutes apart. … … 304 304 series have been extracted from the model. These are the locations where video footage from the event is 305 305 available as described in Section \ref{sec:eyewitness data}. 306 The corresponding are shown in Figure \ref{fig:onshore_timeseries}.306 The corresponding time series are shown in Figure \ref{fig:onshore_timeseries}. 307 307 308 308 … … 325 325 326 326 \caption{Time series obtained from the two onshore locations, North and South, 327 shown in Figure \protect \ref{fig:gauge_locations}. }327 shown in Figure \protect \ref{fig:gauge_locations}. Time is given in hours since the earthquake event (7:59am).} 328 328 \label{fig:onshore_timeseries} 329 329 \end{figure} … … 335 335 Table \ref{tab:depth and flow comparisons}. 336 336 The predicted maximum depths and speeds are all of the same order 337 of what was observed . However, unlike the real event,337 of what was observed as is the approximate arrival time at the two locations. However, unlike the real event, 338 338 the model estimates complete withdrawal of the water between waves at the 339 339 chosen locations and shows that the model must be used with caution at this 340 340 level of detail. 341 Nonetheless, this comparison serves to check that depths and speeds341 Nonetheless, this comparison serves to check that the peak depths and speeds 342 342 predicted are within the range of what is expected. 343 343 -
anuga_work/publications/boxing_day_validation_2008/sensitivity.tex
r7480 r7522 5 5 are computationally intensive and high resolution models of the entire 6 6 evolution process will often take a number of days to 7 run. Consequently, the uncertainty in model predictions is difficult to7 compute. Consequently, the uncertainty in model predictions is difficult to 8 8 quantify as it would require a very large number of runs. 9 9 However, model uncertainty should not be ignored. The aim of this section is … … 28 28 appropriate values of Manning's coefficient range from 0.007 to 0.03 29 29 for tsunami propagation over a sandy sea floor and the reference model 30 uses a value of 0.01. 30 uses a value of 0.01. To investigate sensitivity to this parameter, 31 31 we simulated the maximum onshore inundation using a Manning's 32 32 coefficient of 0.0003 and 0.03. The resulting inundation maps are … … 34 34 and the maximum flow speeds in Figure~\ref{fig:sensitivity_friction_speed}. 35 35 The figure, along with Table~\ref{table:inundationAreas}, 36 shows that the on-shore inundation extent decreases with increasing36 shows, as expected, that the on-shore inundation extent decreases with increasing 37 37 friction and that small perturbations in the friction cause bounded 38 38 changes in the output. This is consistent with the conclusions of … … 46 46 generally well predicted by the generation and propagation models such as 47 47 \textsc{ursga} as demonstrated in Section \ref{sec:resultsPropagation} 48 and also in \cite{thomas2009}. Nevertheless, the effect of errors in 48 and also in \cite{thomas2009} assuming that the source parameters chosen appropriately. 49 Nevertheless, the effect of errors in 49 50 the wave height used as input to the inundation model \textsc{anuga} 50 51 was investigated by perturbing the 51 amplitude of the input wave by $\pm$10 cm. This value was chosen to be larger 52 than the expected error in the amplitude predicted by the propagation model. 52 amplitude of the input wave by $\pm$10 cm. This value was chosen to be consistent 53 with the expected error in the amplitude predicted by the propagation model 54 and amounts to about $\pm$5\% of the maximal waveheight at the boundary. 53 55 54 Figure~\ref{fig:sensitivity_boundary}, Figure~\ref{fig:sensitivity_boundary_speed} ,55 and 56 Figure~\ref{fig:sensitivity_boundary}, Figure~\ref{fig:sensitivity_boundary_speed} 57 and Table~\ref{table:inundationAreas} 56 58 indicate that the inundation severity is directly proportional to the 57 59 boundary waveheight but small -
anuga_work/publications/boxing_day_validation_2008/tsunami07.bib
r7521 r7522 1048 1048 author = {Burbidge, D. and Cummins, P.R. and Mleczko, R. and Thio, H.K.}, 1049 1049 title = "{A Probabilistic Tsunami Hazard Assessment for {W}estern {A}ustralia}", 1050 journal = {Pure appl. geophys.},1050 journal = "{Pure Appl. Geophys.}", 1051 1051 year = {2008}, 1052 1052 volume = {165},
Note: See TracChangeset
for help on using the changeset viewer.