O atual estudo visa cooperar na compreensão dos efeitos da passagem de tempestades extremas sobre os sistemas litorâneos através do estudo de caso na Desembocadura de Cananéia, localizada na Ilha Comprida - SP. O trabalho contou com campanhas de campo para levantamento de dados topográficos e coletas de sedimentos, análises das condições ambientais atuantes durante o período dos campos, determinação da evolução da linha de costa através da análise de fotos aéreas e imagens de satélite e simulações numéricas para caracterizar o regime de ondas atuante na região. Foram produzidos mapas temáticos para representar as cotas de inundação associadas a ocorrência de marés meteorológicas .Os resultados indicaram que a ação das ondas de tempestade se da principalmente por meio do transporte onshore-offshore nas escalas temporais diária e mensal. A evolução da linha de costa obtida apontou para maior eficiência das correntes de deriva litorânea na escala decadal. A comparação entre os resultados topográficos/volumétricos e da evolução da linha de costa evidenciam uma tendência erosiva nas adjacências do Pontal de Fora associada à incidência de eventos atmosféricos extremos para o período entre 2001 e início de 2012. Os resultados das determinações das cotas de inundação apontam que os maiores impactos locais associados às mudanças climáticas serão provenientes das alterações na frequência de ocorrência e intensidade das tempestades / The current study aims to cooperate in the knowledge about the effects of the occurrence of extreme storms on coastal systems through the case study of the Desembocadura de Cananéia, located on Ilha Comprida - SP. The work included campaigns of field survey to collect topographic data and sediment samples, analysis of environmental conditions during the period of surveys, determination of the shoreline evolution through analysis of aerial photos and satellite images and numerical simulations to characterize the wave regime present in the region. Thematic maps were produced to represent the flood hazard associated with the occurrence of storm surges. Results indicated that the action of storm waves occurs primarily by means of the onshore-offshore transport on daily and monthly time scales. The evolution of the shoreline obtained pointed to greater efficiency of longshore currents on the decadal scale. The comparison between the topographic/volumetric results and the evolution of the coastline shows an erosive trend in the vicinity of the Pontal de Fora associated with the incidence of extreme weather events for the period between 2001 and early 2012. The results of the determinations of flood hazard areas points out that the major local impacts related to the global climate changes will come from variations in the frequency of occurrence and intensity of storms

O atual estudo visa cooperar na compreensão dos efeitos da passagem de tempestades extremas sobre os sistemas litorâneos através do estudo de caso na Desembocadura de Cananéia, localizada na Ilha Comprida - SP. O trabalho contou com campanhas de campo para levantamento de dados topográficos e coletas de sedimentos, análises das condições ambientais atuantes durante o período dos campos, determinação da evolução da linha de costa através da análise de fotos aéreas e imagens de satélite e simulações numéricas para caracterizar o regime de ondas atuante na região. Foram produzidos mapas temáticos para representar as cotas de inundação associadas a ocorrência de marés meteorológicas .Os resultados indicaram que a ação das ondas de tempestade se da principalmente por meio do transporte onshore-offshore nas escalas temporais diária e mensal. A evolução da linha de costa obtida apontou para maior eficiência das correntes de deriva litorânea na escala decadal. A comparação entre os resultados topográficos/volumétricos e da evolução da linha de costa evidenciam uma tendência erosiva nas adjacências do Pontal de Fora associada à incidência de eventos atmosféricos extremos para o período entre 2001 e início de 2012. Os resultados das determinações das cotas de inundação apontam que os maiores impactos locais associados às mudanças climáticas serão provenientes das alterações na frequência de ocorrência e intensidade das tempestades / The current study aims to cooperate in the knowledge about the effects of the occurrence of extreme storms on coastal systems through the case study of the Desembocadura de Cananéia, located on Ilha Comprida - SP. The work included campaigns of field survey to collect topographic data and sediment samples, analysis of environmental conditions during the period of surveys, determination of the shoreline evolution through analysis of aerial photos and satellite images and numerical simulations to characterize the wave regime present in the region. Thematic maps were produced to represent the flood hazard associated with the occurrence of storm surges. Results indicated that the action of storm waves occurs primarily by means of the onshore-offshore transport on daily and monthly time scales. The evolution of the shoreline obtained pointed to greater efficiency of longshore currents on the decadal scale. The comparison between the topographic/volumetric results and the evolution of the coastline shows an erosive trend in the vicinity of the Pontal de Fora associated with the incidence of extreme weather events for the period between 2001 and early 2012. The results of the determinations of flood hazard areas points out that the major local impacts related to the global climate changes will come from variations in the frequency of occurrence and intensity of storms

Sea levels have risen approximately 20 cm since the beginning of the 20th century and more than 3 cm in the past 20 years, suggesting that global sea level rise is accelerating. As sea levels continue to rise and storms become more intense, coastal property and populations become more susceptible to damage. Florida is especially vulnerable to hurricane-induced storm surge (HSS) and the onset of accelerated sea-level rise (ASLR) due to its extensive coastline and high population density along the coast. The main purpose of this research is to assess the potential economic impacts of ASLR and HSS for two of western Pasco County’s municipalities, Port Richey and New Port Richey. A Geographic Information System is used to determine the spatial extent at a high-resolution of coastal inundation, the economic loss based on property value and road expenditure due to this inundation, and its impact on critical infrastructure. The results from this study showed coastal flooding generated by 0.5m SLR amounted to 48.8% land loss and $217,108,692 of property loss. Monetary losses from inundated properties shifted dramatically from 1.0m to 1.5m SLR, from $295 million to $417 million, suggesting that the tipping point could only be a half-meter SLR. Based on the 2.0m SLR results, most of major highway US-19 was completely flooded, property tax losses amounted to approximately $7.1 million, and road expenditure was approximately $158 million. Data provided in this study can be useful for coastal management and planning in Port Richey and New Port Richey.

Conventional coastal ocean models solve the shallow water equations, which describe the conservation of mass and momentum when the horizontal length scale is much greater than the vertical length scale. In this case vertical pressure gradients in the momentum equations are nearly hydrostatic. The outputs of coastal ocean models are thus sensitive to the bottom stress terms defined through the formulation of Manning’s n coefficients. This thesis considers the Bayesian inference problem of the Manning’s n coefficient in the context of storm surge based on the coastal ocean ADCIRC model. In the first part if the thesis, we apply an ensemble-based Kalman filter, the singular evolutive interpolated Kalman (SEIK) filter to estimate both a constant Manning’s n coefficient and a 2-D parameterized Manning’s coefficient on one ideal and one of more realistic domain using observation system simulation experiments (OSSEs). We study the sensitivity of the system to the ensemble size. we also access the benefits from using an inflation factor on the filter performance. To study the limitation of the Guassian restricted assumption on the SEIK filter, we also implemented in the second part of this thesis a Markov Chain Monte Carlo (MCMC) method based on a Generalized Polynomial chaos (gPc) approach for the estimation of the 1-D and 2-D Mannning’s n coefficient. The gPc is used to build a surrogate model that imitate the ADCIRC model in order to make the computational cost of implementing the MCMC with the ADCIRC model reasonable. We evaluate the performance of the MCMC-gPc approach and study its robustness to different OSSEs scenario. we also compare its estimates with those resulting from SEIK in term of parameter estimates and full distributions. we present a full analysis of the solution of these two methods, of the contexts of their algorithms, and make recommendation for fully realistic application.

Home to the world’s longest coastline, Canada has experienced devastating economic and social from coastal flooding events. While there have been a variety of mitigation methods employed over the years to increase a community’s resistance to coastal hazards, it is unrealistic to think that there exists a solution to guarantee a community’s safety under all possible flood hazards. Instead, the community’s efforts to raise their resistance to flood hazards should be augmented with careful planning and management to increase a community’s resilience to flood hazards, allowing them to recover quickly after a natural disaster. The first step in elevating a community’s resilience is to better understand the expected hazards that it may experience. This thesis presents two unique case studies to better understand the flooding hazards present on the Canadian coastline. A large-scale numerical model that accounts for the presence of ice was developed to investigate storm surges in Canada’s western Arctic. It was found that the quality of the climatic forcing data used, ERA5, was poor in capturing peak wind speeds, but could be compensated for by using elevated wind drag coefficients. The use of non-traditional high-water marks such as driftwood lines were validated and were shown to significantly alter expected flood return periods compared to the return periods estimated from only the incomplete tide gauge measurements present on the Arctic coastline. The second case study extends the results of a tsunami hydrodynamic simulation on Canada’s Pacific coastline through a life safety assessment. The performance between an agent-based and GIS-based approach to modelling tsunami evacuation were directly compared and were shown to yield different magnitudes in fatality rate and facility demand, but similar trends. Both models agreed on a mitigation option that can significantly reduce the loss of life during a tsunami.

京都大学 / 0048 / 新制・課程博士 / 博士(工学) / 甲第20680号 / 工博第4377号 / 新制||工||1680(附属図書館) / 京都大学大学院工学研究科社会基盤工学専攻 / (主査)教授 平石 哲也, 教授 中北 英一, 准教授 森 信人 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DFAM

Recently, a highly resolved, finite element mesh was developed for the purpose of performing hydrodynamic calculations in the Western North Atlantic Tidal (WNAT) model domain. The WNAT model domain consists of the Gulf of Mexico, the Caribbean Sea, and the entire portion of the North Atlantic Ocean found west of the 60° W meridian. This high resolution mesh (333K) employs 332,582 computational nodes and 647,018 triangular elements to provide approximately 1.0 to 25 km node spacing. In the previous work, the 333K mesh was applied in a Localized Truncation Error Analysis (LTEA) to produce nodal density requirements for the WNAT model domain. The goal of the work herein is to use these LTEA-based element sizing guidelines in order to obtain a more optimal finite element mesh for the WNAT model domain, where optimal refers to minimizing nodes (to enhance computational efficiency) while maintaining model accuracy, through an automated procedure. Initially, three finite element meshes are constructed: 95K, 60K, and 53K. The 95K mesh consists of 95,062 computational nodes and 182,941 triangular elements providing about 0.5 to 120 km node spacing. The 60K mesh contains 60,487 computational nodes and 108,987 triangular elements. It has roughly 0.5 to 185 km node spacing. The 53K mesh includes 52,774 computational nodes and 98,365 triangular elements. This is a particularly coarse mesh, consisting of approximately 0.5 to 160 km node spacing. It is important to note that these three finite element meshes were produced automatically, with each employing the bathymetry and coastline (of various levels of resolution) of the 333K mesh, thereby enabling progress towards an optimal finite element mesh. Tidal simulations are then performed for the WNAT model domain by solving the shallow water equations in a time marching manner for the deviation from mean sea level and depth-integrated velocities at each computational node of the different finite element meshes. In order to verify the model output and compare the performance of the various finite element mesh applications, historical tidal constituent data from 150 tidal stations located within the WNAT model domain are collected and examined. These historical harmonic data are applied in two types of comparative analyses to evaluate the accuracy of the simulation results. First, qualitative comparisons are based on visual sense by utilizing plots of resynthesized model output and historical tidal constituents. Second, quantitative comparisons are performed via a statistical analysis of the errors between model response and historical data. The latter method elicits average phase errors and goodness of average amplitude fits in terms of numerical values, thus providing a quantifiable way to present model error. The error analysis establishes the 53K finite element mesh as optimal when compared to the 333K, 95K, and 60K meshes. However, its required time step of less than ten seconds constrains its application. Therefore, the 53K mesh is manually edited to uphold accurate simulation results and to produce a more computationally efficient mesh, by increasing its time step, so that it can be applied to forecast tide and storm surge in the Western North Atlantic Ocean on a real-time basis.

Florida's Department of Transportation requires design storm tide hydrographs for coastal waters surrounding tidal inlets along the coast of Florida. These hydrographs are used as open ocean boundary conditions for local bridge scour models. At present, very little information is available on the effect that tidal inlets have on these open coast storm tide hydrographs. Furthermore, current modeling practice enforces a single design hydrograph along the open coast boundary for bridge scour models. This thesis expands on these concepts and provides a more fundamental understanding on both of these modeling areas. A numerical parameter study is undertaken to elucidate the influence of tidal inlets on open coast storm tide hydrographs. Four different inlet-bay configurations are developed based on a statistical analysis of existing tidal inlets along the Florida coast. The length and depth of the inlet are held constant in each configuration, but the widths are modified to include the following four inlet profiles: 1) average Florida inlet width; 2) 100 meter inlet width; 3) 500 meter inlet width; and 4) 1000 meter inlet width. In addition, two unique continental shelf profiles are used to design the ocean bathymetry in the model domains: a bathymetry profile consistent with the west/northeast coast of Florida (wide continental shelf width), and a bathymetry profile similar to the southeast coast of Florida (narrow continental shelf width). The four inlet-bay configurations are paired with each of the bathymetry profiles to arrive at eight model domains employed in this study. Results from these domains are compared to control cases that do not include any inlet-bay system in the computational domain. The ADCIRC-2DDI numerical code is used to obtain water surface elevations for all studies performed herein. The code is driven by astronomic tides at the open ocean boundary, and wind velocities and atmospheric pressure profiles over the surface of the computational domains. Model results clearly indicate that the four inlet-bay configurations do not have a significant impact on the open coast storm tide hydrographs. Furthermore, a spatial variance amongst the storm tide hydrographs is recognized for open coast boundary locations extending seaward from the mouth of the inlet. The results and conclusions presented herein have implications toward future bridge scour modeling efforts. In addition, a hindcast study of Hurricane Ivan in the vicinity of Escambia Bay along the Panhandle of Florida is performed to assess the findings of the numerical parameter study in a real-life scenario. Initially, emphasis is placed on domain scale by comparing model results with historical data for three computational domains: an ocean-based domain, a shelf-based domain, and an inlet-based domain. Results indicate that the ocean-based domain favorably simulates storm surge levels within the bay compared to the other model domains. Furthermore, the main conclusions from the numerical parameter study are verified in the hindcast study: 1) the Pensacola Pass-Escambia Bay system has a minimal effect on the open coast storm tide hydrographs; and 2) the open coast storm tide hydrographs exhibit spatial dependence along typical open coast boundary locations.

The city of Pascagoula and its coastal areas along the United States Gulf Coast have experienced many catastrophic hurricanes and were devastated by high storm surges caused by Hurricane Katrina (August 23 to 30, 2005). The National Hurricane Center reported high water marks exceeding 6 meters near the port of Pascagoula with a near 10-meter high water mark recorded near the Hurricane Katrina landfall location in Waveland, MS. Although the Pascagoula River is located 105 km east of the landfall location of Hurricane Katrina, the area was devastated by storm surge-induced inundation because of its low elevation. Building on a preliminary finite element mesh for the Pascagoula River, the work presented herein is aimed at incorporating the marsh areas lying adjacent to the Lower Pascagoula and Escatawpa Rivers for the purpose of simulating the inland inundation which occurred during Hurricane Katrina. ADCIRC-2DDI (ADvanced CIRCulation Model for Shelves, Coasts and Estuaries, Two-Dimensional Depth Integrated) is employed as the hydrodynamic circulation code. The simulations performed in this study apply high-resolution winds and pressures over the 7-day period associated with Hurricane Katrina. The high resolution of the meteorological inputs to the problem coupled with the highly detailed description of the adjacent inundation areas will provide an appropriate modeling tool for studying storm surge dynamics within the Pascagoula River. All simulation results discussed herein are directed towards providing for a full accounting of the hydrodynamics within the Pascagoula River in support of ongoing flood/river forecasting efforts. In order to better understand the hydrodynamics within the Pascagoula River when driven by an extreme storm surge event, the following tasks were completed as a part of this study: 1) Develop an inlet-based floodplain DEM (Digital Elevation Model) for the Pascagoula River. The model employs topography up to the 1.5-meter contour extracted from the Southern Louisiana Gulf Coast Mesh (SL15 Mesh) developed by the Federal Emergency Management Agency (FEMA). 2) Incorporate the inlet-based floodplain model into the Western North Atlantic Tidal (WNAT) model domain, which consists of the Gulf of Mexico, the Caribbean Sea, and the entire portion of the North Atlantic Ocean found west of the 60 degree West meridian, in order to more fully account for the storm surge dynamics occurring within the Pascagoula River. This large-scale modeling approach will utilize high-resolution wind and pressure fields associated with Hurricane Katrina, so that storm surge hydrographs (elevation variance) at the open-ocean boundary locations associated with the localized domain can be adequately obtained. 3) Understand the importance of the various meteorological forcings that are attributable to the storm surge dynamics that are setup within the Pascagoula River. Different implementations of the two model domains (large-scale, including the WNAT model domain; localized, with its focus concentrated solely on the Pascagoula River) will involve the application of tides, storm surge hydrographs and meteorological forcing (winds and pressures) in isolation (i.e., as the single forcing mechanism) and collectively (i.e., together in combination). The following conclusions are drawn from the research presented in this thesis: 1) Incorporating the marsh areas into the preliminary in-bank mesh provides for significant improvement in the astronomic tide simulation; 2) the large-scale modeling approach (i.e., the localized floodplain mesh incorporated into the WNAT model domain) is shown to be most adequate towards simulating storm surge dynamics within the Pascagoula River. Further, we demonstrate the utility of the large-scale model domain towards providing storm surge hydrographs for the open-ocean boundary of the localized domain. Only when the localized domain is forced with the storm surge hydrograph (generated by the large-scale model domain) does it most adequately capture the full behavior of the storm surge. Finally, we discover that while the floodplain description up to the 1.5-m contour greatly improves the model response by allowing for the overtopping of the river banks, a true recreation of the water levels caused by Hurricane Katrina will require a floodplain description up to the 5-m contour.

Storm surge flooding caused by tropical cyclones is a devastating threat to coastal regions, and this threat is growing due to sea-level rise (SLR). Therefore, accurate and rapid projection of the storm surge hazard is critical for coastal communities. This study focuses on developing a new framework that can rapidly predict storm surges under SLR scenarios for any random synthetic storms of interest and assign a probability to its likelihood. The framework leverages the Joint Probability Method with Response Surfaces (JPM-RS) for probabilistic hazard characterization, a storm surge machine learning model, and a SLR model. The JPM probabilities are based on historical tropical cyclone track observations. The storm surge machine learning model was trained based on high-fidelity storm surge simulations provided by the U.S. Army Corps of Engineers (USACE). The SLR was considered by adding the product of the normalized nonlinearity, arising from surge-SLR interaction, and the sea-level change from 1992 to the target year, where nonlinearities are based on high-fidelity storm surge simulations and subsequent analysis by USACE. In this study, this framework was applied to the Chesapeake Bay region of the U.S. and used to estimate the SLR-adjusted probabilistic tropical cyclone flood hazard in two areas: one is an urban Virginia site, and the other is a rural Maryland site. This new framework has the potential to aid in reducing future coastal storm risks in coastal communities by providing robust and rapid hazard assessment that accounts for future sea-level rise. / Master of Science / Storm surge flooding, which is the rise in sea level caused by tropical cyclones and other storms, is a devastating threat to coastal regions, and its impact is increasing due to sea-level rise (SLR). This poses a considerable risk to communities living near the coast. Therefore, it is crucial to accurately and quickly predict the potential for storm surge flooding. This study aimed to develop a new way that can rapidly estimate peak storm surges under different sea-level rise scenarios for any random synthetic storms of interest and assess the likelihood of their occurrence. The approach is based on historical tropical cyclone datasets and a machine learning model trained on high-quality simulations provided by the US Army Corps of Engineers (USACE). The study focused on the Chesapeake Bay area of the US and estimated the probabilistic tropical cyclone flood hazard in two locations, an urban site in Virginia and a rural site in Maryland. This new approach has the potential to assist in reducing coastal storm risks in vulnerable communities by providing a quick and reliable assessment of the hazard that takes into account the effects of future sea-level rise.