Publications
- The GeoClaw software for depth-averaged flows with adaptive refinement
(with Berger, Marsha J. George, David L. and LeVeque, Randall J. and )
Advances in Water Resources ()
@article{Berger.2011, year = {2011}, rating = {0}, keywords = {#math.NA,#physics.geo-ph,math.NA,physics.geo-ph}, title = {{The GeoClaw software for depth-averaged flows with adaptive refinement}}, author = {Berger, Marsha J. and George, David L. and LeVeque, Randall J. and Mandli, Kyle T.}, journal = {Advances in Water Resources}, issn = {0309-1708}, doi = {10.1016/j.advwatres.2011.02.016}, eprint = {1008.0455}, url = {http://www.sciencedirect.com/science/article/pii/S0309170811000480}, pages = {1195--1206}, number = {9}, volume = {34}, note = {\#math.NA \#physics.geo-ph 18 pages, 11 figures, Animations and source code for some examples at http://www.clawpack.org/links/awr10 18 pages, 11 figures, Animations and source code for some examples at http://www.clawpack.org/links/awr10 18 pages, 11 figures, Animations and source code for some examples at http://www.clawpack.org/links/awr10 18 pages, 11 figures, Animations and source code for some examples at http://www.clawpack.org/links/awr10} }[1]M.J. Berger, D.L. George, R.J. LeVeque, K.T. Mandli, The GeoClaw software for depth-averaged flows with adaptive refinement, Advances in Water Resources. 34 (2011) 1195–1206. https://doi.org/10.1016/j.advwatres.2011.02.016.
Many geophysical flow or wave propagation problems can be modeled with two-dimensional depth-averaged equations, of which the shallow water equations are the simplest example. We describe the GeoClaw software that has been designed to solve problems of this nature, consisting of open source Fortran programs together with Python tools for the user interface and flow visualization. This software uses high-resolution shock-capturing finite volume methods on logically rectangular grids, including latitude–longitude grids on the sphere. Dry states are handled automatically to model inundation. The code incorporates adaptive mesh refinement to allow the efficient solution of large-scale geophysical problems. Examples are given illustrating its use for modeling tsunamis and dam-break flooding problems. Documentation and download information is available at www.clawpack.org/geoclaw. - PyClaw: Accessible, Extensible, Scalable Tools for Wave Propagation Problems
(with Ketcheson, David I. Kyle and Ahmadia, Aron J. and Alghamdi, Amal and Luna, Manuel Quezada de and Parsani, Matteo and Knepley, Matthew G. and Emmett, Matthew)
SIAM Journal on Scientific Computing ()
@article{Ketcheson.2012, year = {2012}, rating = {0}, title = {{PyClaw: Accessible, Extensible, Scalable Tools for Wave Propagation Problems}}, author = {Ketcheson, David I. and Mandli, Kyle and Ahmadia, Aron J. and Alghamdi, Amal and Luna, Manuel Quezada de and Parsani, Matteo and Knepley, Matthew G. and Emmett, Matthew}, journal = {SIAM Journal on Scientific Computing}, issn = {1064-8275}, doi = {10.1137/110856976}, eprint = {1111.6583}, pages = {C210--C231}, number = {4}, volume = {34}, language = {English}, month = nov }[2]D.I. Ketcheson, K. Mandli, A.J. Ahmadia, A. Alghamdi, M.Q.de Luna, M. Parsani, M.G. Knepley, M. Emmett, PyClaw: Accessible, Extensible, Scalable Tools for Wave Propagation Problems, SIAM Journal on Scientific Computing. 34 (2012) C210–C231. https://doi.org/10.1137/110856976.
Development of scientific software involves tradeoffs between ease of use, generality, and performance. We describe the design of a general hyperbolic PDE solver that can be operated with the convenience of MATLAB yet achieves efficiency near that of hand-coded Fortran and scales to the largest supercomputers. This is achieved by using Python for most of the code while employing automatically wrapped Fortran kernels for computationally intensive routines, and using Python bindings to interface with a parallel computing library and other numerical packages. The software described here is PyClaw, a Python-based structured grid solver for general systems of hyperbolic PDEs [K. T. Mandli et al., PyClaw Software, Version 1.0, http://numerics.kaust.edu.sa/pyclaw/ (2011)]. PyClaw provides a powerful and intuitive interface to the algorithms of the existing Fortran codes Clawpack and SharpClaw, simplifying code development and use while providing massive parallelism and scalable solvers via the PETSc library. The package is further augmented by use of PyWENO for generation of efficient high-order weighted essentially nonoscillatory reconstruction code. The simplicity, capability, and performance of this approach are demonstrated through application to example problems in shallow water flow, compressible flow, and elasticity.
(with Terrel, Andy R Kyle T)
()
@inproceedings{Terrel.2012, year = {2012}, rating = {0}, author = {Terrel, Andy R and Mandli, Kyle T}, title = {{ManyClaw: Slicing and dicing Riemann solvers for next generation highly parallel architectures}}, urldate = {0}, pages = {1 -- 6}, series = {TACC-Intel Symposium on Highly Parallel Architectures 2012}, keywords = {}, month = feb }[3]A.R. Terrel, K.T. Mandli, ManyClaw: Slicing and dicing Riemann solvers for next generation highly parallel architectures, in: 2012: pp. 1–6.
Next generation computer architectures will include order of magnitude more intra-node parallelism; however, many application programmers have a difficult time keeping their codes current with the state- of-the-art machines. In this context, we analyze Hyperbolic PDE solvers, which are used in the solution of many important applications in science and engineering. We present ManyClaw , a project intended to explore the exploitation of intra-node parallelism in hyperbolic PDE solvers via the Clawpack software package for solving hyperbolic PDEs. Our goal is to separate the low level parallelism and the physical equations thus providing users the capability to leverage intra-node parallelism without explicitly writing code to take advantage of newer architectures.
(with Burstedde, Carsten Calhoun, Donna A and Kyle and Terrel, Andy R)
()
@inproceedings{Burstedde.2013, year = {2013}, rating = {0}, keywords = {Adaptive mesh refinement,CLAWPACK,HPC,hyperbolic conservation laws,manycore}, author = {Burstedde, Carsten and Calhoun, Donna A and Mandli, Kyle and Terrel, Andy R}, title = {{ForestClaw: Hybrid forest-of-octrees AMR for hyperbolic conservation laws}}, doi = {10.3233/978-1-61499-381-0-253}, url = {http://www.ebooks.iospress.nl/volumearticle/35888}, urldate = {0}, pages = {253 -- 262}, volume = {25}, series = {ParCo 2013}, note = {\#Adaptive mesh refinement \#CLAWPACK \#HPC \#hyperbolic conservation laws \#manycore} }[4]C. Burstedde, D.A. Calhoun, K. Mandli, A.R. Terrel, ForestClaw: Hybrid forest-of-octrees AMR for hyperbolic conservation laws, in: 2013: pp. 253–262. https://doi.org/10.3233/978-1-61499-381-0-253.
We present a new hybrid paradigm for parallel adaptive mesh refinement (AMR) that combines the scalability and lightweight architecture of tree-based AMR with the computational efficiency of patch-based solvers for hyperbolic conservation laws. The key idea is to interpret each leaf of the AMR hierarchy as one uniform compute patch in Rd with md degrees of freedom, where m is customarily between 8 and 32. Thus, computation on each patch can be optimized for speed, while we inherit the flexibility of adaptive meshes. In our work we choose to integrate with the p4est AMR library since it allows us to compose the mesh from multiple mapped octrees and enables the cubed sphere and other nontrivial multiblock geometries. We describe aspects of the parallel implementation and close with scalings for both MPI-only and OpenMP/MPI hybrid runs, where the largest MPI run executes on 16,384 CPU cores.- A numerical method for the two layer shallow water equations with dry states
(with )
Ocean Modelling ()
@article{Mandli.2013, year = {2013}, rating = {0}, title = {{A numerical method for the two layer shallow water equations with dry states}}, author = {Mandli, Kyle T.}, journal = {Ocean Modelling}, issn = {1463-5003}, doi = {10.1016/j.ocemod.2013.08.001}, eprint = {1308.1905}, url = {http://www.sciencedirect.com/science/article/pii/S1463500313001212}, pages = {80--91}, volume = {72}, language = {English}, local-url = {file://localhost/Users/kmandli/Documents/Papers%20Library/Mandli/Mandli_2013.pdf} }[5]K.T. Mandli, A numerical method for the two layer shallow water equations with dry states, Ocean Modelling. 72 (2013) 80–91. https://doi.org/10.1016/j.ocemod.2013.08.001.
A numerical method is proposed for solving the two layer shallow water equations with variable bathymetry in one dimension based on high-resolution f-wave-propagation finite volume methods. The method splits the jump in the fluxes and source terms into waves propagating away from each grid cell interface. It addresses the required determination of the system’s eigenstructure and a scheme for evaluating the flux and source terms. It also handles dry states in the system where the bottom layer depth becomes zero, utilizing existing methods for the single layer solution and handling single layer dry states that can exist independently. Sample results are shown illustrating the method and its handling of dry states including an idealized ocean setting. - Uncertainty quantification and inference of Manning’s friction coefficients using DART buoy data during the Tōhoku tsunami
(with Sraj, Ihab and Knio, Omar M. and Dawson, Clint N. and Hoteit, Ibrahim)
Ocean Modelling ()
@article{Sraj.2014, year = {2014}, rating = {0}, title = {{Uncertainty quantification and inference of Manning’s friction coefficients using DART buoy data during the Tōhoku tsunami}}, author = {Sraj, Ihab and Mandli, Kyle T. and Knio, Omar M. and Dawson, Clint N. and Hoteit, Ibrahim}, journal = {Ocean Modelling}, issn = {1463-5003}, doi = {10.1016/j.ocemod.2014.09.001}, url = {http://www.sciencedirect.com/science/article/pii/S1463500314001322}, pages = {82--97}, volume = {83}, language = {English} }[6]I. Sraj, K.T. Mandli, O.M. Knio, C.N. Dawson, I. Hoteit, Uncertainty quantification and inference of Manning’s friction coefficients using DART buoy data during the Tōhoku tsunami, Ocean Modelling. 83 (2014) 82–97. https://doi.org/10.1016/j.ocemod.2014.09.001.
Tsunami computational models are employed to explore multiple flooding scenarios and to predict water elevations. However, accurate estimation of water elevations requires accurate estimation of many model parameters including the Manning’s n friction parameterization. Our objective is to develop an efficient approach for the uncertainty quantification and inference of the Manning’s n coefficient which we characterize here by three different parameters set to be constant in the on-shore, near-shore and deep-water regions as defined using iso-baths. We use Polynomial Chaos (PC) to build an inexpensive surrogate for the GeoClaw model and employ Bayesian inference to estimate and quantify uncertainties related to relevant parameters using the DART buoy data collected during the Tōhoku tsunami. The surrogate model significantly reduces the computational burden of the Markov Chain Monte-Carlo (MCMC) sampling of the Bayesian inference. The PC surrogate is also used to perform a sensitivity analysis. - Adaptive mesh refinement for storm surge
(with Dawson, Clint N.)
Ocean Modelling ()
@article{Mandli.2014, year = {2014}, rating = {0}, title = {{Adaptive mesh refinement for storm surge}}, author = {Mandli, Kyle T. and Dawson, Clint N.}, journal = {Ocean Modelling}, issn = {1463-5003}, doi = {10.1016/j.ocemod.2014.01.002}, eprint = {1401.5744}, pages = {36--50}, volume = {75}, keywords = {} }[7]K.T. Mandli, C.N. Dawson, Adaptive mesh refinement for storm surge, Ocean Modelling. 75 (2014) 36–50. https://doi.org/10.1016/j.ocemod.2014.01.002.
An approach to utilizing adaptive mesh refinement algorithms for storm surge modeling is proposed. Currently numerical models exist that can resolve the details of coastal regions but are often too costly to be run in an ensemble forecasting framework without significant computing resources. The application of adaptive mesh refinement algorithms substantially lowers the computational cost of a storm surge model run while retaining much of the desired coastal resolution. The approach presented is implemented in the GeoClaw framework and compared to ADCIRC for Hurricane Ike along with observed tide gauge data and the computational cost of each model run. - Visualizing uncertainties in a storm surge ensemble data assimilation and forecasting system
(with Höllt, Thomas Altaf, M. Umer and and Hadwiger, Markus and Dawson, Clint N. and Hoteit, Ibrahim)
Natural Hazards ()
@article{Höllt.2015, year = {2015}, rating = {0}, title = {{Visualizing uncertainties in a storm surge ensemble data assimilation and forecasting system}}, author = {Höllt, Thomas and Altaf, M. Umer and Mandli, Kyle T. and Hadwiger, Markus and Dawson, Clint N. and Hoteit, Ibrahim}, journal = {Natural Hazards}, issn = {0921-030X}, doi = {10.1007/s11069-015-1596-y}, url = {http://link.springer.com/article/10.1007/s11069-015-1596-y}, pages = {317--336}, number = {1}, volume = {77}, language = {English}, month = jan }[8]T. Höllt, M.U. Altaf, K.T. Mandli, M. Hadwiger, C.N. Dawson, I. Hoteit, Visualizing uncertainties in a storm surge ensemble data assimilation and forecasting system, Natural Hazards. 77 (2015) 317–336. https://doi.org/10.1007/s11069-015-1596-y.
We present a novel integrated visualization system that enables the interactive visual analysis of ensemble simulations and estimates of the sea surface height and other model variables that are used for storm surge prediction. Coastal inundation, caused by hurricanes and tropical storms, poses large risks for today’s societies. High-fidelity numerical models of water levels driven by hurricane-force winds are required to predict these events, posing a challenging computational problem, and even though computational models continue to improve, uncertainties in storm surge forecasts are inevitable. Today, this uncertainty is often exposed to the user by running the simulation many times with different parameters or inputs following a Monte-Carlo framework in which uncertainties are represented as stochastic quantities. This results in multidimensional, multivariate and multivalued data, so-called ensemble data. While the resulting datasets are very comprehensive, they are also huge in size and thus hard to visualize and interpret. In this paper, we tackle this problem by means of an interactive and integrated visual analysis system. By harnessing the power of modern graphics processing units for visualization as well as computation, our system allows the user to browse through the simulation ensembles in real time, view specific parameter settings or simulation models and move between different spatial and temporal regions without delay. In addition, our system provides advanced visualizations to highlight the uncertainty or show the complete distribution of the simulations at user-defined positions over the complete time series of the prediction. We highlight the benefits of our system by presenting its application in a real-world scenario using a simulation of Hurricane Ike. - Clawpack: building an open source ecosystem for solving hyperbolic PDEs
(with Kyle T Ahmadia, Aron J and Berger, Marsha and Calhoun, Donna and George, David L and Hadjimichael, Yiannis and Ketcheson, David I and Lemoine, Grady I and LeVeque, Randall J)
PeerJ Computer Science ()
@article{Mandli.2016, year = {2016}, rating = {0}, title = {{Clawpack: building an open source ecosystem for solving hyperbolic PDEs}}, author = {Mandli, Kyle T and Ahmadia, Aron J and Berger, Marsha and Calhoun, Donna and George, David L and Hadjimichael, Yiannis and Ketcheson, David I and Lemoine, Grady I and LeVeque, Randall J}, journal = {PeerJ Computer Science}, doi = {10.7717/peerj-cs.68}, url = {https://peerj.com/articles/cs-68/?td=wk}, pages = {e68}, number = {3}, volume = {2}, language = {English}, month = aug }[9]K.T. Mandli, A.J. Ahmadia, M. Berger, D. Calhoun, D.L. George, Y. Hadjimichael, D.I. Ketcheson, G.I. Lemoine, R.J. LeVeque, Clawpack: building an open source ecosystem for solving hyperbolic PDEs, PeerJ Computer Science. 2 (2016) e68. https://doi.org/10.7717/peerj-cs.68.
Clawpack is a software package designed to solve nonlinear hyperbolic partial differential equations using high-resolution finite volume methods based on Riemann solvers and limiters. The package includes a number of variants aimed at different applications and user communities. Clawpack has been actively developed as an open source project for over 20 years. The latest major release, Clawpack 5, introduces a number of new features and changes to the code base and a new development model based on GitHub and Git submodules. This article provides a summary of the most significant changes, the rationale behind some of these changes, and a description of our current development model. - Quantifying uncertainties in fault slip distribution during the Tōhoku tsunami using polynomial chaos
(with Sraj, Ihab and Knio, Omar M. and Dawson, Clint N. and Hoteit, Ibrahim)
Ocean Dynamics ()
@article{Sraj.2017, year = {2017}, rating = {0}, title = {{Quantifying uncertainties in fault slip distribution during the Tōhoku tsunami using polynomial chaos}}, author = {Sraj, Ihab and Mandli, Kyle T. and Knio, Omar M. and Dawson, Clint N. and Hoteit, Ibrahim}, journal = {Ocean Dynamics}, issn = {1616-7341}, doi = {10.1007/s10236-017-1105-9}, eprint = {1607.07414}, pages = {1535--1551}, number = {12}, volume = {67}, keywords = {} }[10]I. Sraj, K.T. Mandli, O.M. Knio, C.N. Dawson, I. Hoteit, Quantifying uncertainties in fault slip distribution during the Tōhoku tsunami using polynomial chaos, Ocean Dynamics. 67 (2017) 1535–1551. https://doi.org/10.1007/s10236-017-1105-9.
An efficient method for inferring Manning’s n coefficients using water surface elevation data was presented in Sraj et al. (Ocean Modell 83:82–97 2014a) focusing on a test case based on data collected during the Tōhoku earthquake and tsunami. Polynomial chaos (PC) expansions were used to build an inexpensive surrogate for the numerical model GeoClaw, which were then used to perform a sensitivity analysis in addition to the inversion. In this paper, a new analysis is performed with the goal of inferring the fault slip distribution of the Tōhoku earthquake using a similar problem setup. The same approach to constructing the PC surrogate did not lead to a converging expansion; however, an alternative approach based on basis pursuit denoising was found to be suitable. Our result shows that the fault slip distribution can be inferred using water surface elevation data whereas the inferred values minimize the error between observations and the numerical model. The numerical approach and the resulting inversion are presented in this work.
(with Huang, Yipeng Guo, Ning and Seok, Mingoo and Tsividis, Yannis and Kyle and Sethumadhavan, Simha)
()
@inproceedings{Huang.2017, year = {2017}, rating = {0}, author = {Huang, Yipeng and Guo, Ning and Seok, Mingoo and Tsividis, Yannis and Mandli, Kyle and Sethumadhavan, Simha}, title = {{Hybrid analog-digital solution of nonlinear partial differential equations}}, isbn = {978-1-4503-4952-9}, doi = {10.1145/3123939.3124550}, urldate = {0}, pages = {665 -- 678}, series = {the 50th Annual IEEE/ACM International Symposium}, month = oct }[11]Y. Huang, N. Guo, M. Seok, Y. Tsividis, K. Mandli, S. Sethumadhavan, Hybrid analog-digital solution of nonlinear partial differential equations, in: 2017: pp. 665–678. https://doi.org/10.1145/3123939.3124550.
… ACM ISBN 978-1-4503-4952-9/17/10... 15.00 https://doi.org/ 10.1145 / 3123939.3124550 ACM Reference format: Yipeng Huang, Ning Guo, Mingoo Seok, Yannis Tsividis, Kyle Mandli, and Simha Sethumadhavan. 2017 … https://doi.org/ 10.1145 / 3123939.3124550 …- Bayesian inference of earthquake parameters from buoy data using a polynomial chaos-based surrogate
(with Giraldi, Loïc Maître, Olivier P. Le and and Dawson, Clint N. and Hoteit, Ibrahim and Knio, Omar M.)
Computational Geosciences ()
@article{Giraldi.2017, year = {2017}, rating = {0}, title = {{Bayesian inference of earthquake parameters from buoy data using a polynomial chaos-based surrogate}}, author = {Giraldi, Loïc and Maître, Olivier P. Le and Mandli, Kyle T. and Dawson, Clint N. and Hoteit, Ibrahim and Knio, Omar M.}, journal = {Computational Geosciences}, issn = {1420-0597}, doi = {10.1007/s10596-017-9646-z}, url = {https://link.springer.com/article/10.1007\%2Fs10596-017-9646-z}, pages = {683--699}, number = {4}, volume = {21}, language = {English}, month = apr }[12]L. Giraldi, O.P.L. Maître, K.T. Mandli, C.N. Dawson, I. Hoteit, O.M. Knio, Bayesian inference of earthquake parameters from buoy data using a polynomial chaos-based surrogate, Computational Geosciences. 21 (2017) 683–699. https://doi.org/10.1007/s10596-017-9646-z.
This work addresses the estimation of the parameters of an earthquake model by the consequent tsunami, with an application to the Chile 2010 event. We are particularly interested in the Bayesian inference of the location, the orientation, and the slip of an Okada-based model of the earthquake ocean floor displacement. The tsunami numerical model is based on the GeoClaw software while the observational data is provided by a single DARTⓇ buoy. We propose in this paper a methodology based on polynomial chaos expansion to construct a surrogate model of the wave height at the buoy location. A correlated noise model is first proposed in order to represent the discrepancy between the computational model and the data. This step is necessary, as a classical independent Gaussian noise is shown to be unsuitable for modeling the error, and to prevent convergence of the Markov Chain Monte Carlo sampler. Second, the polynomial chaos model is subsequently improved to handle the variability of the arrival time of the wave, using a preconditioned non-intrusive spectral method. Finally, the construction of a reduced model dedicated to Bayesian inference is proposed. Numerical results are presented and discussed. - Evolution and Controls of Large Glacial Lakes in the Nepal Himalaya
(with Haritashya, Umesh K. Kargel, Jeffrey S. and Shugar, Dan H. and Leonard, Gregory J. and Strattman, Katherine and Watson, C. Scott and Shean, David and Harrison, Stephan and and Regmi, Dhananjay)
Remote Sensing ()
@article{Haritashya.2018, year = {2018}, rating = {0}, title = {{Evolution and Controls of Large Glacial Lakes in the Nepal Himalaya}}, author = {Haritashya, Umesh K. and Kargel, Jeffrey S. and Shugar, Dan H. and Leonard, Gregory J. and Strattman, Katherine and Watson, C. Scott and Shean, David and Harrison, Stephan and Mandli, Kyle T. and Regmi, Dhananjay}, journal = {Remote Sensing}, doi = {10.3390/rs10050798}, pages = {798}, number = {5}, volume = {10}, language = {English} }[13]U.K. Haritashya, J.S. Kargel, D.H. Shugar, G.J. Leonard, K. Strattman, C.S. Watson, D. Shean, S. Harrison, K.T. Mandli, D. Regmi, Evolution and Controls of Large Glacial Lakes in the Nepal Himalaya, Remote Sensing. 10 (2018) 798. https://doi.org/10.3390/rs10050798.
Glacier recession driven by climate change produces glacial lakes, some of which are hazardous. Our study assesses the evolution of three of the most hazardous moraine-dammed proglacial lakes in the Nepal Himalaya—Imja, Lower Barun, and Thulagi. Imja Lake (up to 150 m deep; 78.4 × 106 m3 volume; surveyed in October 2014) and Lower Barun Lake (205 m maximum observed depth; 112.3 × 106 m3 volume; surveyed in October 2015) are much deeper than previously measured, and their readily drainable volumes are slowly growing. Their surface areas have been increasing at an accelerating pace from a few small supraglacial lakes in the 1950s/1960s to 1.33 km2 and 1.79 km2 in 2017, respectively. In contrast, the surface area (0.89 km2) and volume of Thulagi lake (76 m maximum observed depth; 36.1 × 106 m3; surveyed in October 2017) has remained almost stable for about two decades. Analyses of changes in the moraine dams of the three lakes using digital elevation models (DEMs) quantifies the degradation of the dams due to the melting of their ice cores and hence their natural lowering rates as well as the potential for glacial lake outburst floods (GLOFs). We examined the likely future evolution of lake growth and hazard processes associated with lake instability, which suggests faster growth and increased hazard potential at Lower Barun lake.
(with Ferreira, Chaulio R Kyle T and Bader, Michael)
()
@inproceedings{Ferreira.2018, year = {2018}, rating = {0}, author = {Ferreira, Chaulio R and Mandli, Kyle T and Bader, Michael}, title = {{Vectorization of Riemann solvers for the single- and multi-layer shallow water equations}}, doi = {10.1109/hpcs.2018.00073}, url = {https://ieeexplore.ieee.org/abstract/document/8514378}, urldate = {0}, series = { International Conference on High Performance Computing Simulation}, month = jul }[14]C.R. Ferreira, K.T. Mandli, M. Bader, Vectorization of Riemann solvers for the single- and multi-layer shallow water equations, in: 2018. https://doi.org/10.1109/hpcs.2018.00073.
We discuss vectorization of normal and transverse Riemann solvers for the single- and multi-layer shallow water equations. Our approach is simple and portable, as it is based on auto-vectorization by the compiler, aided by OpenMP 4.0 directives. Despite the high complexity of the solver routines, the Intel Fortran Compiler proved itself able to successfully vectorize loops containing calls to these solvers, after only a few small changes in their code. We evaluate the performance of the vectorized Riemann solvers within the context of GeoClaw, a software designed for simulation of geophysical flows with finite volume methods. Our performance studies consider two platforms with different sets of SIMD instructions: a dual-socket Haswell system with the AVX2 instruction set (256-bit) and an Intel Xeon Phi (Knights Landing) with AVX-512 instructions (512-bit). The experimental results indicate performance im- provements of up to 2.1x on the former platform and up to 6.5x on the latter (with double-precision arithmetic). We also show that these speedups can easily compensate for the overhead introduced by the rearrangement of the simulation data structures, which might be necessary to achieve efficient vectorization.- An Implementation of Adaptive Mesh Refinement for Shallow Water Equations
(with Schwarzschild, Avi Kyle T)
arXiv ()
@article{Schwarzschild.2018, year = {2018}, rating = {0}, title = {{An Implementation of Adaptive Mesh Refinement for Shallow Water Equations}}, author = {Schwarzschild, Avi and Mandli, Kyle T}, journal = {arXiv}, doi = {10.48550/arxiv.1803.01450}, eprint = {1803.01450}, url = {https://arxiv.org/abs/1803.01450}, keywords = {}, month = mar }[15]A. Schwarzschild, K.T. Mandli, An Implementation of Adaptive Mesh Refinement for Shallow Water Equations, ArXiv. (2018). https://doi.org/10.48550/arxiv.1803.01450.
An implementation of adaptive mesh refinement algorithms is presented for use with multilayer shallow water equations. Currently, adaptive mesh refinement is implemented with a single layer shallow water model in the GeoClaw framework. This implementation, also in the GeoClaw framework, is for multilayer models, which have been implemented in GeoClaw previously. Until now, however, these models were too computationally expensive to run on large domains while resolving detail in coastal regions. - Model reduction of a parametrized scalar hyperbolic conservation law using displacement interpolation
(with Rim, Donsub Kyle T)
arXiv ()
@article{Rim.2018um, year = {2018}, rating = {0}, title = {{Model reduction of a parametrized scalar hyperbolic conservation law using displacement interpolation}}, author = {Rim, Donsub and Mandli, Kyle T}, journal = {arXiv}, doi = {10.48550/arxiv.1805.05938}, eprint = {1805.05938}, url = {https://arxiv.org/abs/1805.05938}, keywords = {}, month = may }[16]D. Rim, K.T. Mandli, Model reduction of a parametrized scalar hyperbolic conservation law using displacement interpolation, ArXiv. (2018). https://doi.org/10.48550/arxiv.1805.05938.
We propose a model reduction technique for parametrized partial differential equations arising from scalar hyperbolic conservation laws. The key idea of the technique is to construct basis functions that are local in parameter and time space via displacement interpolation. The construction is motivated by the observation that the derivative of solutions to hyperbolic conservation laws satisfy a contractive property with respect to the Wasserstein metric [Bolley et al. J. Hyperbolic Differ. Equ. 02 (2005), pp. 91-107]. We will discuss the approximation properties of the displacement interpolation, and show that it can naturally complement linear interpolation. Numerical experiments illustrate that we can successfully achieve the model reduction of a parametrized Burgers’ equation, and that the reduced order model is suitable for performing typical tasks in uncertainty quantification. - hp discontinuous Galerkin methods for parametric, wind-driven water wave models
(with Conroy, Colton J. Kubatko, Ethan J. and Nappi, Angela and Sebian, Rachel and West, Dustin and )
Advances in Water Resources ()
@article{Conroy.2018, year = {2018}, rating = {0}, title = {{hp discontinuous Galerkin methods for parametric, wind-driven water wave models}}, author = {Conroy, Colton J. and Kubatko, Ethan J. and Nappi, Angela and Sebian, Rachel and West, Dustin and Mandli, Kyle T.}, journal = {Advances in Water Resources}, issn = {0309-1708}, doi = {10.1016/j.advwatres.2018.04.008}, url = {https://www.sciencedirect.com/science/article/pii/S0309170817306474}, pages = {70--83}, volume = {119}, language = {English}, month = sep }[17]C.J. Conroy, E.J. Kubatko, A. Nappi, R. Sebian, D. West, K.T. Mandli, hp discontinuous Galerkin methods for parametric, wind-driven water wave models, Advances in Water Resources. 119 (2018) 70–83. https://doi.org/10.1016/j.advwatres.2018.04.008.
We recast the parametric, wind-driven water wave modeling paradigm into a weak form that is advantageous to discontinuous Galerkin (DG) methods and demonstrate some advantages of polynomial refinement versus standard mesh refinement. Hindcast studies performed over Lake Erie indicate that this simplified parametric approach, when paired with advanced numerics, produces similar error measures to the well-established spectral wave model known as SWAN while executing significantly faster in terms of CPU time. - Surrogate-based parameter inference in debris flow model
(with Navarro, Maria Maître, Olivier P. Le and Hoteit, Ibrahim and George, David L. and and Knio, Omar M.)
Computational Geosciences ()
- Displacement Interpolation Using Monotone Rearrangement
(with Rim, Donsub )
SIAM/ASA Journal on Uncertainty Quantification ()
@article{Rim.2018, year = {2018}, rating = {0}, title = {{Displacement Interpolation Using Monotone Rearrangement}}, author = {Rim, Donsub and Mandli, Kyle T.}, journal = {SIAM/ASA Journal on Uncertainty Quantification}, doi = {10.1137/18m1168315}, eprint = {1712.04028}, url = {https://epubs.siam.org/doi/abs/10.1137/18M1168315}, pages = {1503--1531}, number = {4}, volume = {6}, language = {English}, month = nov }[19]D. Rim, K.T. Mandli, Displacement Interpolation Using Monotone Rearrangement, SIAM/ASA Journal on Uncertainty Quantification. 6 (2018) 1503–1531. https://doi.org/10.1137/18m1168315.
When approximating a function that depends on a parameter, one encounters many practical examples where linear interpolation or linear approximation with respect to the parameters proves ineffective. This is particularly true for responses from hyperbolic partial differential equations (PDEs) where linear, low-dimensional bases are difficult to construct. We propose the use of displacement interpolation, where the interpolation is done on the optimal transport map between the functions at nearby parameters, to achieve an effective dimensionality reduction of hyperbolic phenomena. We further propose a multidimensional extension by using the intertwining property of the Radon transform. This extension is a generalization of the classical translational representation of Lax and Phillips [P. D. Lax and R. S. Phillips, Bull. Amer. Math. Soc., 70 (1964), pp. 130–142]. - Dynamically adaptive data-driven simulation of extreme hydrological flows
(with Jain, Pushkar Kumar Kyle and Hoteit, Ibrahim and Knio, Omar and Dawson, Clint)
Ocean Modelling ()
@article{Jain.2018, year = {2018}, rating = {0}, title = {{Dynamically adaptive data-driven simulation of extreme hydrological flows}}, author = {Jain, Pushkar Kumar and Mandli, Kyle and Hoteit, Ibrahim and Knio, Omar and Dawson, Clint}, journal = {Ocean Modelling}, issn = {1463-5003}, doi = {10.1016/j.ocemod.2017.12.004}, pages = {85--103}, number = {Estuar. Coasts 38 1 2015}, volume = {122}, language = {English}, keywords = {}, month = jan }[20]P.K. Jain, K. Mandli, I. Hoteit, O. Knio, C. Dawson, Dynamically adaptive data-driven simulation of extreme hydrological flows, Ocean Modelling. 122 (2018) 85–103. https://doi.org/10.1016/j.ocemod.2017.12.004.
Hydrological hazards such as storm surges, tsunamis, and rainfall-induced flooding are physically complex events that are costly in loss of human life and economic productivity. Many such disasters could be mitigated through improved emergency evacuation in real-time and through the development of resilient infrastructure based on knowledge of how systems respond to extreme events. Data-driven computational modeling is a critical technology underpinning these efforts. This investigation focuses on the novel combination of methodologies in forward simulation and data assimilation. The forward geophysical model utilizes adaptive mesh refinement (AMR), a process by which a computational mesh can adapt in time and space based on the current state of a simulation. The forward solution is combined with ensemble based data assimilation methods, whereby observations from an event are assimilated into the forward simulation to improve the veracity of the solution, or used to invert for uncertain physical parameters. The novelty in our approach is the tight two-way coupling of AMR and ensemble filtering techniques. The technology is tested using actual data from the Chile tsunami event of February 27, 2010. These advances offer the promise of significantly transforming data-driven, real-time modeling of hydrological hazards, with potentially broader applications in other science domains.
(with Barba, Lorena A. Barker, Lecia J. and Blank, Douglas S. and Brown, Jed and Downey, Allen B. and George, Timothy and Heagy, Lindsey J. and and Moore, Jason K. and Lippert, David and Niemeyer, Kyle E. and Watkins, Ryan R. and West, Richard H. and Wickes, Elizabeth and Willing, Carol and Zingale, Michael)
()
@book{Barba.2019, year = {2019}, title = {{Teaching and Learning with Jupyter}}, author = {Barba, Lorena A. and Barker, Lecia J. and Blank, Douglas S. and Brown, Jed and Downey, Allen B. and George, Timothy and Heagy, Lindsey J. and Mandli, Kyle T. and Moore, Jason K. and Lippert, David and Niemeyer, Kyle E. and Watkins, Ryan R. and West, Richard H. and Wickes, Elizabeth and Willing, Carol and Zingale, Michael}, url = {https://jupyter4edu.github.io/jupyter-edu-book/}, month = dec }[21]L.A. Barba, L.J. Barker, D.S. Blank, J. Brown, A.B. Downey, T. George, L.J. Heagy, K.T. Mandli, J.K. Moore, D. Lippert, K.E. Niemeyer, R.R. Watkins, R.H. West, E. Wickes, C. Willing, M. Zingale, Teaching and Learning with Jupyter, 2019. https://jupyter4edu.github.io/jupyter-edu-book/.
- Tropical cyclone hazard to Mumbai in the recent historical climate Tropical cyclone hazard to Mumbai in the recent historical climate
(with Sobel, Adam H. Lee, Chia-Ying and Camargo, Suzana J. and and Emanuel, Kerry A. and Mukhopadhyay, Parthasarathi and Mahakur, M.)
Monthly Weather Review ()
@article{Sobel.2019, year = {2019}, title = {{Tropical cyclone hazard to Mumbai in the recent historical climate Tropical cyclone hazard to Mumbai in the recent historical climate}}, author = {Sobel, Adam H. and Lee, Chia-Ying and Camargo, Suzana J. and Mandli, Kyle T. and Emanuel, Kerry A. and Mukhopadhyay, Parthasarathi and Mahakur, M.}, journal = {Monthly Weather Review}, issn = {0027-0644}, doi = {10.1175/mwr-d-18-0419.1}, pages = {2355--2366}, number = {7}, volume = {147} }[22]A.H. Sobel, C.-Y. Lee, S.J. Camargo, K.T. Mandli, K.A. Emanuel, P. Mukhopadhyay, M. Mahakur, Tropical cyclone hazard to Mumbai in the recent historical climate Tropical cyclone hazard to Mumbai in the recent historical climate, Monthly Weather Review. 147 (2019) 2355–2366. https://doi.org/10.1175/mwr-d-18-0419.1.
The hazard to the city of Mumbai, India from a possible severe tropical cyclone under the recent historical climate is considered. The authors first determine, based on a review of primary sources, that the Bombay cyclone of 1882, documented in a number of print and internet sources and claimed to have caused 100,000 or more deaths, did not occur. Two different tropical cyclone hazard models, both of which generate large numbers of synthetic cyclones using environmental data - here taken from reanalyses in the satellite era - as input, are then used to quantify the hazard, in conjunction with historical observations. Both models indicate that a severe cyclone landfall at or near Mumbai is possible, though unlikely in any given year. Return periods for wind speeds exceeding 100 kt (the threshold for category 3 in the Saffir-Simpson scale) at Mumbai itself are estimated to be in the range of thousands to greater than ten thousand years, while the return period for a storm with maximum wind speed of 100 kt or greater passing within 150 km of Mumbai (possibly close enough to generate a substantial storm surge at the city) is estimated to be around 500 years. Return periods for winds exceeding 65 kt (hurricane intensity on the Saffir-Simpson scale) are estimated to be around 200 years at Mumbai itself, 50–90 years within 150 km. Climate change is not explicitly considered in this study, but the hazard to the city is likely to be increasing due to sea level rise as well as changes in storm climatology. - Modeling and Simulation of Tsunami Impact: A Short Review of Recent Advances and Future Challenges
(with Marras, Simone )
Geosciences ()
@article{Marras.2020, year = {2020}, title = {{Modeling and Simulation of Tsunami Impact: A Short Review of Recent Advances and Future Challenges}}, author = {Marras, Simone and Mandli, Kyle T.}, journal = {Geosciences}, doi = {10.3390/geosciences11010005}, pages = {5}, number = {1}, volume = {11} }[23]S. Marras, K.T. Mandli, Modeling and Simulation of Tsunami Impact: A Short Review of Recent Advances and Future Challenges, Geosciences. 11 (2020) 5. https://doi.org/10.3390/geosciences11010005.
Tsunami modeling and simulation has changed in the past few years more than it has in decades, especially with respect to coastal inundation. Among other things, this change is supported by the approaching era of exa-scale computing, whether via GPU or more likely forms of hybrid computing whose presence is growing across the geosciences. For reasons identified in this review, exa-scale computing efforts will impact the on-shore, highly turbulent régime to a higher degree than the 2D shallow water equations used to model tsunami propagation in the open ocean. This short review describes the different approaches to tsunami modeling from generation to impact and underlines the limits of each model based on the flow régime. Moreover, from the perspective of a future comprehensive multi-scale modeling infrastructure to simulate a full tsunami, we underline the current challenges associated with this approach and review the few efforts that are currently underway to achieve this goal. A table of existing tsunami software packages is provided along with an open Github repository to allow developers and model users to update the table with additional models as they are published and help with model discoverability. - Continental Scale Heterogeneous Channel Flow Routing Strategy for Operational Forecasting Models
(with Meselhe, Ehab Lamjiri, Maryam A. and Flint, Kelly and Matus, Sean and White, Eric D. and Kyle)
JAWRA Journal of the American Water Resources Association ()
@article{Meselhe.2021, year = {2021}, title = {{Continental Scale Heterogeneous Channel Flow Routing Strategy for Operational Forecasting Models}}, author = {Meselhe, Ehab and Lamjiri, Maryam A. and Flint, Kelly and Matus, Sean and White, Eric D. and Mandli, Kyle}, journal = {JAWRA Journal of the American Water Resources Association}, issn = {1093-474X}, doi = {10.1111/1752-1688.12847}, pages = {209--221}, number = {2}, volume = {57} }[24]E. Meselhe, M.A. Lamjiri, K. Flint, S. Matus, E.D. White, K. Mandli, Continental Scale Heterogeneous Channel Flow Routing Strategy for Operational Forecasting Models, JAWRA Journal of the American Water Resources Association. 57 (2021) 209–221. https://doi.org/10.1111/1752-1688.12847.
The benefits of operational forecasting models to the general public are numerous. They support water management decisions, provide the opportunity to mitigate the impacts of weather‐ and flood‐related disasters and potentially save lives and properties. Channel flow routing is a key component of these models and affects their ability to forecast flood depth, duration, and extent. Continental scale channel flow routing within the operational forecasting environments encounters a broad spectrum of hydraulic characteristics. Deploying computationally demanding approaches, such as the dynamic wave, should be limited in time and space to conditions where the inertia terms are significant (typically in low‐gradient environments and whenever backwater effects are prominent); otherwise, efficient and robust methods, e.g., Kinematic, Muskingum‐Cunge or diffusive waves should be the default. The heterogeneous routing approach presented here provides a framework to evaluate the balance between friction, inertia, and pressure and strategically triggers the appropriate wave approximation. The strategy recommended here is to activate the appropriate wave approximation based on the ambient hydraulic conditions, and smoothly transitions among these approximations. This strategy, if successfully implemented, would strike a balance among the performance metrics of operational forecasting models, namely, computational efficiency, accuracy, and minimization of computational instabilities. Research Impact Statement: A heterogeneous routing strategy, driven by the surrounding hydraulic conditions, activates terms in the momentum equation as needed to strike a balance between computational speed and accuracy. - An \h\-Box Method for Shallow Water Equations Including Barriers
(with Li, Jiao Kyle T)
SIAM Journal on Scientific Computing ()
@article{Li.2021, year = {2021}, title = {{An \$h\$-Box Method for Shallow Water Equations Including Barriers}}, author = {Li, Jiao and Mandli, Kyle T}, journal = {SIAM Journal on Scientific Computing}, issn = {1064-8275}, doi = {10.1137/19m128363x}, pages = {B431--B454}, number = {2}, volume = {43}, local-url = {file://localhost/Users/kmandli/Documents/Papers%20Library/Li/Li_2021.pdf} }[25]J. Li, K.T. Mandli, An \h\-Box Method for Shallow Water Equations Including Barriers, SIAM Journal on Scientific Computing. 43 (2021) B431–B454. https://doi.org/10.1137/19m128363x.
- Optimization of Coastal Protections in the Presence of Climate Change
(with Miura, Yuki Dinenis, Philip C. and and Deodatis, George and Bienstock, Daniel)
Frontiers in Climate ()
@article{Miura.2021ilf, year = {2021}, title = {{Optimization of Coastal Protections in the Presence of Climate Change}}, author = {Miura, Yuki and Dinenis, Philip C. and Mandli, Kyle T. and Deodatis, George and Bienstock, Daniel}, journal = {Frontiers in Climate}, doi = {10.3389/fclim.2021.613293}, pages = {613293}, volume = {3} }[26]Y. Miura, P.C. Dinenis, K.T. Mandli, G. Deodatis, D. Bienstock, Optimization of Coastal Protections in the Presence of Climate Change, Frontiers in Climate. 3 (2021) 613293. https://doi.org/10.3389/fclim.2021.613293.
It is generally acknowledged that interdependent critical infrastructure in coastal urban areas is constantly threatened by storm-induced flooding. Due to changing climate effects, such as sea level rise (SLR), the occurrence of catastrophic events will be more frequent and may trigger an increased likelihood of severe hazards. Planning a protective measure or mitigation strategy is a complex problem given the constraints that it must fit within a prescribed and limited fiscal budget and be beneficial to the community it protects both socially and economically. This article proposes a methodology for optimizing protective measures and mitigation strategies for interdependent infrastructures subjected to storm-induced flooding and climate change impacts such as SLR. Optimality is defined in this methodology as a maximum reduction in overall expected losses within a prescribed budget (compared to the expected losses in the case of doing nothing for protection/mitigation). Protective measures can include seawalls, barriers, artificial dunes, restoration of wetlands, raising individual buildings, sealing parts of the infrastructure, strategic retreat, insurance, and many more. The optimal protective strategy can be a combination of several protective measures implemented over space and time. The optimization process starts with parameterizing the protective measures. Storm-induced flooding and SLR, and their corresponding consequences, are estimated using a GIS-based subdivision-redistribution methodology (GISSR) developed by the authors for finding a rough solution in the first brute-force iterations of the optimization loop. A storm surge computational model called GeoClaw is subsequently used to simulate ensembles of synthetic storms in order to fine-tune and achieve the optimal solution. Damage loss, including economic impacts, is quantified based on calculated flood estimates. The suitability of the potential optimal solution is examined and assessed with input from stakeholders’ interviews. It should be mentioned that the results and conclusions provided in this work depend on the assumptions made about future sea level rise (SLR). The authors acknowledge that there are other, more severe predictions for sea level rise (SLR), than the one used in this paper. - Quantifying air–water turbulence with moment field equations
(with Conroy, Colton J. and Kubatko, Ethan J.)
Journal of Fluid Mechanics ()
@article{Conroy.2021rah, year = {2021}, title = {{Quantifying air–water turbulence with moment field equations}}, author = {Conroy, Colton J. and Mandli, Kyle T. and Kubatko, Ethan J.}, journal = {Journal of Fluid Mechanics}, issn = {0022-1120}, doi = {10.1017/jfm.2021.242}, pages = {A39}, volume = {917} }[27]C.J. Conroy, K.T. Mandli, E.J. Kubatko, Quantifying air–water turbulence with moment field equations, Journal of Fluid Mechanics. 917 (2021) A39. https://doi.org/10.1017/jfm.2021.242.
Energy transfer in turbulent fluids is non-Gaussian. We quantify non-Gaussian energy transfer between the atmosphere and bodies of water using a turbulent diffusion operator coupled with temporally self-affine velocity distributions and a recursive integration method that produce multifractal measures. The measures serve as input to a system of moment field equations (derived from Navier–Stokes) that generate and track high-frequency gravity waves that propagate through the water surface (as a result of the air–water interactions). The dimension of the support of the air–water turbulence produced by our methods falls within the range of theory and observation, and correspondingly, hindcast statistical measures of the water-wave surface such as significant water-wave height and wave period are well correlated to observational buoy data. Further, our recursive integration method can be used by spectral resolving phase-averaged models to interpolate temporal wind data to smaller scales to capture the non-Gaussian behaviour of the air–water interaction. - A methodological framework for determining an optimal coastal protection strategy against storm surges and sea level rise
(with Miura, Yuki Qureshi, Huda and Ryoo, Chanyang and Dinenis, Philip C. and Li, Jiao and and Deodatis, George and Bienstock, Daniel and Lazrus, Heather and Morss, Rebecca)
Natural Hazards ()
@article{Miura.2021, year = {2021}, title = {{A methodological framework for determining an optimal coastal protection strategy against storm surges and sea level rise}}, author = {Miura, Yuki and Qureshi, Huda and Ryoo, Chanyang and Dinenis, Philip C. and Li, Jiao and Mandli, Kyle T. and Deodatis, George and Bienstock, Daniel and Lazrus, Heather and Morss, Rebecca}, journal = {Natural Hazards}, issn = {0921-030X}, doi = {10.1007/s11069-021-04661-5}, pages = {1821--1843}, number = {2}, volume = {107} }[28]Y. Miura, H. Qureshi, C. Ryoo, P.C. Dinenis, J. Li, K.T. Mandli, G. Deodatis, D. Bienstock, H. Lazrus, R. Morss, A methodological framework for determining an optimal coastal protection strategy against storm surges and sea level rise, Natural Hazards. 107 (2021) 1821–1843. https://doi.org/10.1007/s11069-021-04661-5.
Interdependent critical infrastructures in coastal regions, including transportation, electrical grid, and emergency services, are continually threatened by storm-induced flooding. This has been demonstrated a number of times, most recently by hurricanes such as Harvey and Maria, as well as Sandy and Katrina. The need to protect these infrastructures with robust protection mechanisms is critical for our continued existence along the world’s coastlines. Planning these protections is non-trivial given the rare-event nature of strong storms and climate change manifested through sea level rise. This article proposes a framework for a methodology that combines multiple computational models, stakeholder interviews, and optimization to find an optimal protective strategy over time for critical coastal infrastructure while being constrained by budgetary considerations. - High-Speed GIS-Based Simulation of Storm Surge–Induced Flooding Accounting for Sea Level Rise
(with Miura, Yuki Kyle T and Deodatis, George)
Natural Hazards Review ()
@article{Miura.2021a0t, year = {2021}, title = {{High-Speed GIS-Based Simulation of Storm Surge–Induced Flooding Accounting for Sea Level Rise}}, author = {Miura, Yuki and Mandli, Kyle T and Deodatis, George}, journal = {Natural Hazards Review}, issn = {1527-6988}, doi = {10.1061/(asce)nh.1527-6996.0000465}, number = {3}, volume = {22} }[29]Y. Miura, K.T. Mandli, G. Deodatis, High-Speed GIS-Based Simulation of Storm Surge–Induced Flooding Accounting for Sea Level Rise, Natural Hazards Review. 22 (2021). https://doi.org/10.1061/(asce)nh.1527-6996.0000465.
- A new tropical cyclone surge index incorporating the effects of coastal geometry, bathymetry and storm information
(with Islam, Md. Rezuanul Lee, Chia-Ying and and Takagi, Hiroshi)
Scientific Reports ()
@article{Islam.2021, year = {2021}, title = {{A new tropical cyclone surge index incorporating the effects of coastal geometry, bathymetry and storm information}}, author = {Islam, Md. Rezuanul and Lee, Chia-Ying and Mandli, Kyle T. and Takagi, Hiroshi}, journal = {Scientific Reports}, doi = {10.1038/s41598-021-95825-7}, pmid = {34408207}, pmcid = {PMC8373937}, pages = {16747}, number = {1}, volume = {11} }[30]M.R. Islam, C.-Y. Lee, K.T. Mandli, H. Takagi, A new tropical cyclone surge index incorporating the effects of coastal geometry, bathymetry and storm information, Scientific Reports. 11 (2021) 16747. https://doi.org/10.1038/s41598-021-95825-7.
This study presents a new storm surge hazard potential index (SSHPI) for estimating tropical cyclone (TC) induced peak surge levels at a coast. The SSHPI incorporates parameters that are often readily available at real-time: intensity in 10-min maximum wind speed, radius of 50-kt wind, translation speed, coastal geometry, and bathymetry information. The inclusion of translation speed and coastal geometry information lead to improvements of the SSHPI to other existing surge indices. A retrospective analysis of SSHPI using data from 1978–2019 in Japan suggests that this index captures historical events reasonably well. In particular, it explains \textbackslashtextasciitilde 66% of the observed variance and \textbackslashtextasciitilde 74% for those induced by TCs whose landfall intensity was larger than 79-kt. The performance of SSHPI is not sensitive to the type of coastal geometry (open coasts or semi-enclosed bays). Such a prediction methodology can decrease numerical computation requirements, improve public awareness of surge hazards, and may also be useful for communicating surge risk. - Numerical Considerations for Quantifying Air–Water Turbulence with Moment Field Equations
(with Conroy, Colton J. and Kubatko, Ethan J.)
Water Waves ()
@article{Conroy.2021, year = {2021}, title = {{Numerical Considerations for Quantifying Air–Water Turbulence with Moment Field Equations}}, author = {Conroy, Colton J. and Mandli, Kyle T. and Kubatko, Ethan J.}, journal = {Water Waves}, issn = {2523-367X}, doi = {10.1007/s42286-021-00048-y}, pages = {319--354}, number = {2}, volume = {3}, local-url = {file://localhost/Users/kmandli/Documents/Papers%20Library/Conroy/Conroy_2021.pdf} }[31]C.J. Conroy, K.T. Mandli, E.J. Kubatko, Numerical Considerations for Quantifying Air–Water Turbulence with Moment Field Equations, Water Waves. 3 (2021) 319–354. https://doi.org/10.1007/s42286-021-00048-y.
We investigate energy transfer of air–water interactions and develop a numerical method that captures its temporal variability and generates and tracks the short waves that form in the water surface as a result of the air–water turbulence. We solve a novel system of balance equations derived from the Navier–Stokes equations known as moment field equations. The main advantage of our approach is that we do not assume a priori that the stochastic random variables that quantify the turbulent energy transfer between air and water are Gaussian. We generate non-conservative multifractal measures of turbulent energy transfer using a recursive integration process and a self-affine velocity kernel. The kernel exactly satisfies the (duration limited) kinetic equation for waves as well as invariant scaling properties of the Navier–Stokes equations. This allows us to derive source terms for the moment field equations using a turbulent diffusion operator. The operator quantifies energy transfer along a space time path associated with pressure instabilities in the air–sea interface and transfers the statistical shape (or fractal dimension) of the atmosphere to the wind-sea. Because we use observational data to begin the recursive integration process, the ocean–atmosphere interaction is inherently built into the model. Numerical results from application of our methods to air–sea turbulence off the coast of New Jersey and New York indicate that our methods produce measures of turbulent energy transfer that match theory and observation, and, correspondingly, significant wave heights and average wave periods predicted by our model qualitatively match buoy data. - A Novel Framework for Parametric Analysis of Coastal Transition Zone Modeling
(with Chegini, Taher Coelho, Gustavo de Almeida and Ratcliff, John and Ferreira, Celso M. and Kyle and Burke, Patrick and Li, Hong‐Yi)
JAWRA Journal of the American Water Resources Association ()
@article{Chegini.2022, year = {2022}, title = {{A Novel Framework for Parametric Analysis of Coastal Transition Zone Modeling}}, author = {Chegini, Taher and Coelho, Gustavo de Almeida and Ratcliff, John and Ferreira, Celso M. and Mandli, Kyle and Burke, Patrick and Li, Hong‐Yi}, journal = {JAWRA Journal of the American Water Resources Association}, issn = {1093-474X}, doi = {10.1111/1752-1688.12983}, pages = {86--103}, number = {1}, volume = {58} }[32]T. Chegini, G.de A. Coelho, J. Ratcliff, C.M. Ferreira, K. Mandli, P. Burke, H.Y. Li, A Novel Framework for Parametric Analysis of Coastal Transition Zone Modeling, JAWRA Journal of the American Water Resources Association. 58 (2022) 86–103. https://doi.org/10.1111/1752-1688.12983.
Vulnerability of coastal regions to extreme events motivates an operational coupled inland‐coastal modeling strategy focusing on the coastal transition zone (CTZ), an area between the coast and upland river. To tackle this challenge, we propose a top‐down framework for investigating the contribution of different processes to the hydrodynamics of CTZs with various geometrical shapes, different physical properties, and under several forcing conditions. We further propose a novel method, called tidal vanishing point (TVP), for delineating the extent of CTZs through the upland. We demonstrate the applicability of our framework over the United States East and Gulf coasts. We categorize CTZs in the region into three classes, namely, without estuary (direct river–coast connection), triangular‐, and trapezoidal‐shaped estuary. The results show that although semidiurnal tidal constituents are dominant in most cases, diurnal tidal constituents become more prevalent in the river segment as the discharge increases. Also, decreasing the bed roughness value promotes more significant changes in the results than increasing it by the same value. Additionally, the estuary promotes tidal energy attenuation and consequently decreases the reach of tidal signals through the upland. The proposed framework is generic and extensible to any coastal region. - Under the surface: Pressure-induced planetary-scale waves, volcanic lightning, and gaseous clouds caused by the submarine eruption of Hunga Tonga-Hunga Ha’apai volcano
(with Yuen, David A. Scruggs, Melissa A. and Spera, Frank J. and Zheng, Yingcai and Hu, Hao and McNutt, Stephen R. and Thompson, Glenn and Kyle and Keller, Barry R. and Wei, Songqiao Shawn and Peng, Zhigang and Zhou, Zili and Mulargia, Francesco and Tanioka, Yuichiro)
Earthquake Research Advances ()
@article{Yuen.2022, year = {2022}, title = {{Under the surface: Pressure-induced planetary-scale waves, volcanic lightning, and gaseous clouds caused by the submarine eruption of Hunga Tonga-Hunga Ha'apai volcano}}, author = {Yuen, David A. and Scruggs, Melissa A. and Spera, Frank J. and Zheng, Yingcai and Hu, Hao and McNutt, Stephen R. and Thompson, Glenn and Mandli, Kyle and Keller, Barry R. and Wei, Songqiao Shawn and Peng, Zhigang and Zhou, Zili and Mulargia, Francesco and Tanioka, Yuichiro}, journal = {Earthquake Research Advances}, issn = {2772-4670}, doi = {10.1016/j.eqrea.2022.100134}, pages = {100134}, number = {3}, volume = {2} }[33]D.A. Yuen, M.A. Scruggs, F.J. Spera, Y. Zheng, H. Hu, S.R. McNutt, G. Thompson, K. Mandli, B.R. Keller, S.S. Wei, Z. Peng, Z. Zhou, F. Mulargia, Y. Tanioka, Under the surface: Pressure-induced planetary-scale waves, volcanic lightning, and gaseous clouds caused by the submarine eruption of Hunga Tonga-Hunga Ha’apai volcano, Earthquake Research Advances. 2 (2022) 100134. https://doi.org/10.1016/j.eqrea.2022.100134.
We present a narrative of the eruptive events culminating in the cataclysmic January 15, 2022 eruption of Hunga Tonga-Hunga Ha’apai Volcano by synthesizing diverse preliminary seismic, volcanological, sound wave, and lightning data available within the first few weeks after the eruption occurred. The first hour of eruptive activity produced fast-propagating tsunami waves, long-period seismic waves, loud audible sound waves, infrasonic waves, exceptionally intense volcanic lightning and an unsteady volcanic plume that transiently reached—at 58 km—the Earth’s mesosphere. Energetic seismic signals were recorded worldwide and the globally stacked seismogram showed episodic seismic events within the most intense periods of phreatoplinian activity, and they correlated well with the infrasound pressure waveform recorded in Fiji. Gravity wave signals were strong enough to be observed over the entire planet in just the first few hours, with some circling the Earth multiple times subsequently. These large-amplitude, long-wavelength atmospheric disturbances come from the Earth’s atmosphere being forced by the magmatic mixture of tephra, melt and gasses emitted by the unsteady but quasi-continuous eruption from 0402±1–1800 UTC on January 15, 2022. Atmospheric forcing lasted much longer than rupturing from large earthquakes recorded on modern instruments, producing a type of shock wave that originated from the interaction between compressed air and ambient (wavy) sea surface. This scenario differs from conventional ideas of earthquake slip, landslides, or caldera collapse-generated tsunami waves because of the enormous (∼1000x) volumetric change due to the supercritical nature of volatiles associated with the hot, volatile-rich phreatoplinian plume. The time series of plume altitude can be translated to volumetric discharge and mass flow rate. For an eruption duration of ∼12 h, the eruptive volume and mass are estimated at 1.9 km3 and ∼2 900 Tg, respectively, corresponding to a VEI of 5–6 for this event. The high frequency and intensity of lightning was enhanced by the production of fine ash due to magma—seawater interaction with concomitant high charge per unit mass and the high pre-eruptive concentration of dissolved volatiles. Analysis of lightning flash frequencies provides a rapid metric for plume activity and eruption magnitude. Many aspects of this eruption await further investigation by multidisciplinary teams. It represents a unique opportunity for fundamental research regarding the complex, non-linear behavior of high energetic volcanic eruptions and attendant phenomena, with critical implications for hazard mitigation, volcano forecasting, and first-response efforts in future disasters. - Advancing Interdisciplinary and Convergent Science for Communities: Lessons Learned through the NCAR Early-Career Faculty Innovator Program
(with Bukvic, Anamaria Kyle and Finn, Donovan and Mayo, Talea and Wong-Parodi, Gabrielle and Merdjanoff, Alexis and Alland, Joshua and Davis, Christopher and Haacker, Rebecca and Morss, Rebecca and O’Lenick, Cassandra and Wilhelmi, Olga and Lombardozzi, Danica)
Bulletin of the American Meteorological Society ()
@article{Bukvic.2022, year = {2022}, title = {{Advancing Interdisciplinary and Convergent Science for Communities: Lessons Learned through the NCAR Early-Career Faculty Innovator Program}}, author = {Bukvic, Anamaria and Mandli, Kyle and Finn, Donovan and Mayo, Talea and Wong-Parodi, Gabrielle and Merdjanoff, Alexis and Alland, Joshua and Davis, Christopher and Haacker, Rebecca and Morss, Rebecca and O’Lenick, Cassandra and Wilhelmi, Olga and Lombardozzi, Danica}, journal = {Bulletin of the American Meteorological Society}, issn = {0003-0007}, doi = {10.1175/bams-d-21-0265.1}, pages = {E2513--E2532}, number = {11}, volume = {103} }[34]A. Bukvic, K. Mandli, D. Finn, T. Mayo, G. Wong-Parodi, A. Merdjanoff, J. Alland, C. Davis, R. Haacker, R. Morss, C. O’Lenick, O. Wilhelmi, D. Lombardozzi, Advancing Interdisciplinary and Convergent Science for Communities: Lessons Learned through the NCAR Early-Career Faculty Innovator Program, Bulletin of the American Meteorological Society. 103 (2022) E2513–E2532. https://doi.org/10.1175/bams-d-21-0265.1.
Abstract The authors introduce the National Center for Atmospheric Research’s Early-Career Faculty Innovator Program and present lessons learned about advancing interdisciplinary and convergent science with and for society. The Innovator Program brings together faculty and students from the social sciences with NCAR researchers to conduct interdisciplinary and convergent research on problems motivated by societal challenges in the face of climate change and environmental hazards. This article discusses aspects of program structure and the research being conducted. The article also emphasizes the challenges and successes of the research collaborations within the Innovator Program, along with lessons learned about engaging in highly interdisciplinary, potentially convergent work, particularly from the early-career perspective. Many projects involve faculty PIs from racially, ethnically, or otherwise minoritized groups, and minority serving institutions (MSIs), or those who engage with marginalized communities. Hence, the Innovator Program is contributing to the development of a growing research community pursuing science with and for society that also broadens participation in research related to the atmospheric sciences. - Extreme Water Level Simulation and Component Analysis in Delaware Estuary during Hurricane Isabel
(with Yin, Dongxiao Muñoz, David F. and Bakhtyar, Roham and Xue, Z. George and Moftakhari, Hamed and Ferreira, Celso and Kyle)
JAWRA Journal of the American Water Resources Association ()
@article{Yin.2022, year = {2022}, title = {{Extreme Water Level Simulation and Component Analysis in Delaware Estuary during Hurricane Isabel}}, author = {Yin, Dongxiao and Muñoz, David F. and Bakhtyar, Roham and Xue, Z. George and Moftakhari, Hamed and Ferreira, Celso and Mandli, Kyle}, journal = {JAWRA Journal of the American Water Resources Association}, issn = {1093-474X}, doi = {10.1111/1752-1688.12947}, pages = {19--33}, number = {1}, volume = {58} }[35]D. Yin, D.F. Muñoz, R. Bakhtyar, Z.G. Xue, H. Moftakhari, C. Ferreira, K. Mandli, Extreme Water Level Simulation and Component Analysis in Delaware Estuary during Hurricane Isabel, JAWRA Journal of the American Water Resources Association. 58 (2022) 19–33. https://doi.org/10.1111/1752-1688.12947.
Sea level rise and intense hurricane events make the East and Gulf Coasts of the United States increasingly vulnerable to flooding, which necessitates the development of computational models for accurate water level simulation in these areas to safeguard the coastal wellbeing. With this regard, a model framework for water level simulation over coastal transition zone during hurricane events is built in this study. The model takes advantage of the National Water Model’s strength in simulating rainfall–runoff process, and D‐Flow Flexible Mesh’s ability to support unstructured grid in hydrodynamic processes simulation with storm surges/tides information from the Advanced CIRCulation model. We apply the model on the Delaware Estuary to simulate extreme water level and to investigate the contribution of different physical components to it during Hurricane Isabel (2003). The model shows satisfactory performance with an average Willmott skill of 0.965. Model results suggest that storm surge is the most dominating component of extreme water level with an average contribution of 78.16%, second by the astronomical tide with 19.52%. While the contribution of rivers is mainly restricted to the upper part of the estuary upstream of Schuylkill River, local wind‐induced water level is more pronounced with values larger than 0.2 m over most part of the estuary. - Inter‐Model Comparison of Delft3D‐FM and 2D HEC‐RAS for Total Water Level Prediction in Coastal to Inland Transition Zones
(with Muñoz, David F. Yin, Dongxiao and Bakhtyar, Roham and Moftakhari, Hamed and Xue, Zuo and Kyle and Ferreira, Celso)
JAWRA Journal of the American Water Resources Association ()
@article{Muñoz.2022, year = {2022}, title = {{Inter‐Model Comparison of Delft3D‐FM and 2D HEC‐RAS for Total Water Level Prediction in Coastal to Inland Transition Zones}}, author = {Muñoz, David F. and Yin, Dongxiao and Bakhtyar, Roham and Moftakhari, Hamed and Xue, Zuo and Mandli, Kyle and Ferreira, Celso}, journal = {JAWRA Journal of the American Water Resources Association}, issn = {1093-474X}, doi = {10.1111/1752-1688.12952}, pages = {34--49}, number = {1}, volume = {58} }[36]D.F. Muñoz, D. Yin, R. Bakhtyar, H. Moftakhari, Z. Xue, K. Mandli, C. Ferreira, Inter‐Model Comparison of Delft3D‐FM and 2D HEC‐RAS for Total Water Level Prediction in Coastal to Inland Transition Zones, JAWRA Journal of the American Water Resources Association. 58 (2022) 34–49. https://doi.org/10.1111/1752-1688.12952.
Hydrodynamic models play a key role in simulating total water level (TWL), that is, a combination of river flow, tide, surge, wind and wave‐induced water level, and representing flood inundation dynamics in coastal areas. An appropriate selection of two‐dimensional (2D) models that integrate riverine and estuarine interactions with ocean dynamics is crucial to generate accurate TWL predictions and assist stakeholders and federal agencies in decision making and flood emergency responses. In this study, we compare the performance of two widely used hydrodynamic models (e.g., 2D HEC‐RAS and Delft3D‐Flexible Mesh [FM]) with respect to their ability of predicting TWL in Delaware Bay, United States. Based on a previously established model configuration, we simulate Hurricane Sandy and Isabel that affected the Bay and led to considerable damages and economic losses. We then evaluate model capabilities with tidal analysis, compare observed vs. simulated TWL and analyze spatiotemporal variations of TWL through scenario‐based simulations. Our results suggest that atmospheric forcing input in Delft3D‐FM significantly improves TWL predictions as compared to those of 2D HEC‐RAS. Furthermore, model simulations with Delft3D‐FM can be faster than 2D HEC‐RAS by a factor of 6–10. Despite these advantages, 2D HEC‐RAS (version 5.07) is a noncommercial software easier to implement and can be a simpler alternative for modeling extreme events when atmospheric forcing is not relevant in the model domain. - Moving from interdisciplinary to convergent research across geoscience and social sciences: challenges and strategies
(with Finn, Donovan Kyle and Bukvic, Anamaria and Davis, Christopher A and Haacker, Rebecca and Morss, Rebecca E and O’Lenick, Cassandra R and Wilhelmi, Olga and Wong-Parodi, Gabrielle and Merdjanoff, Alexis A and Mayo, Talea L)
Environmental Research Letters ()
@article{Finn.2022, year = {2022}, title = {{Moving from interdisciplinary to convergent research across geoscience and social sciences: challenges and strategies}}, author = {Finn, Donovan and Mandli, Kyle and Bukvic, Anamaria and Davis, Christopher A and Haacker, Rebecca and Morss, Rebecca E and O’Lenick, Cassandra R and Wilhelmi, Olga and Wong-Parodi, Gabrielle and Merdjanoff, Alexis A and Mayo, Talea L}, journal = {Environmental Research Letters}, doi = {10.1088/1748-9326/ac7409}, pages = {061002}, number = {6}, volume = {17} }[37]D. Finn, K. Mandli, A. Bukvic, C.A. Davis, R. Haacker, R.E. Morss, C.R. O’Lenick, O. Wilhelmi, G. Wong-Parodi, A.A. Merdjanoff, T.L. Mayo, Moving from interdisciplinary to convergent research across geoscience and social sciences: challenges and strategies, Environmental Research Letters. 17 (2022) 061002. https://doi.org/10.1088/1748-9326/ac7409.
- Episodic Magma Hammers for the 15 January 2022 Cataclysmic Eruption of Hunga Tonga‐Hunga Ha’apai
(with Zheng, Yingcai Hu, Hao and Spera, Frank J. and Scruggs, Melissa and Thompson, Glenn and Jin, Yuesu and Lapen, Tom and McNutt, Stephen R. and Kyle and Peng, Zhigang and Yuen, Dave A.)
Geophysical Research Letters ()
@article{Zheng.2023, year = {2023}, title = {{Episodic Magma Hammers for the 15 January 2022 Cataclysmic Eruption of Hunga Tonga‐Hunga Ha'apai}}, author = {Zheng, Yingcai and Hu, Hao and Spera, Frank J. and Scruggs, Melissa and Thompson, Glenn and Jin, Yuesu and Lapen, Tom and McNutt, Stephen R. and Mandli, Kyle and Peng, Zhigang and Yuen, Dave A.}, journal = {Geophysical Research Letters}, issn = {0094-8276}, doi = {10.1029/2023gl102763}, number = {8}, volume = {50} }[38]Y. Zheng, H. Hu, F.J. Spera, M. Scruggs, G. Thompson, Y. Jin, T. Lapen, S.R. McNutt, K. Mandli, Z. Peng, D.A. Yuen, Episodic Magma Hammers for the 15 January 2022 Cataclysmic Eruption of Hunga Tonga‐Hunga Ha’apai, Geophysical Research Letters. 50 (2023). https://doi.org/10.1029/2023gl102763.
Understanding the forces and magma system dynamics on timescales of seconds to minutes remains challenging. In the January 2022 phreatoplinian Hunga Tonga‐Hunga Ha’apai eruption, four remarkably similar seismic subevents within a 5‐min interval occurred during the intensifying early eruptive phase. The subevents are similar in waveforms and durations (∼25 s each). Each subevent begins with an unusual negative P‐wave polarity which is inferred, using full‐wave seismic modeling, to be caused by an upward single‐force mechanism at the volcano created by a magma hammer likely in response to magma flow blockage/constriction during the early part of the eruption as discharge rapidly increased over orders of magnitude with concomitant conduit geometry evolution and instability. Our proposed episodic magma hammer model is consistent with thermodynamic and phase properties of the magmatic mixture, and yields an estimate of conduit mass flow in agreement with vent discharge rates derived from satellite imagery of plume heights. The seismic record of the 15 January 2022 Hunga Tonga‐Hunga Ha’apai explosive eruption exhibited a remarkably regular pattern, recording repeating volcanic processes. Within an interval of ∼300 s during an early eruptive phase, four strong seismic subevents occurred and were recorded by global seismic stations. Detailed seismic analyses showed that each of these subevents is similar in waveform and duration and is characterized by a sequence of four forces: upward, downward, upward, and downward. We suggest that the first upward seismic force at the volcano was likely created by an ascending magma colliding with a blockage or conduit constriction that occurred during and because of the ongoing eruption. We attribute the other forces to upward magma backflow and piston motion in the conduit, owing to their similar time durations. The magma hammer mechanism allows us to estimate the magma flow rate in the subsurface conduit, which is consistent with the vent discharge rate observed by satellite imagery. The 15 January 2022 Hunga Tonga‐Hunga Ha’apai eruption had four episodic seismic subevents with similar waveforms within ∼300 s An unusual upward force jump‐started each subevent A magma hammer explains the force and estimates the subsurface magma mass flux which fits the vent discharge rate based on satellite data The 15 January 2022 Hunga Tonga‐Hunga Ha’apai eruption had four episodic seismic subevents with similar waveforms within ∼300 s An unusual upward force jump‐started each subevent A magma hammer explains the force and estimates the subsurface magma mass flux which fits the vent discharge rate based on satellite data - Coastal Sediments 2023, The Proceedings of the Coastal Sediments 2023
(with Wang, Ping Royer, Elizabeth and Rosati, Julie D and BEEVER, MEGAN A and IRISH, JENNIFER L and WEISS, ROBERT and MANDLI, KYLE T and JEFFRIES, CATHERINE R)
In Coastal Sediments 2023, WORLD SCIENTIFIC ()
@incollection{Wang.20237lc, year = {2023}, title = {{Coastal Sediments 2023, The Proceedings of the Coastal Sediments 2023}}, author = {Wang, Ping and Royer, Elizabeth and Rosati, Julie D and BEEVER, MEGAN A and IRISH, JENNIFER L and WEISS, ROBERT and MANDLI, KYLE T and JEFFRIES, CATHERINE R}, booktitle = {Coastal Sediments 2023}, isbn = {9789811279898}, urldate = {2024-05-17}, pages = {2462--2467}, series = {Coastal Sediments 2023}, publisher = {WORLD SCIENTIFIC}, copyright = {All rights reserved}, doi = {10.1142/9789811275135\_0224} }[39]P. Wang, E. Royer, J.D. Rosati, M.E.G.A.N.A. BEEVER, J.E.N.N.I.F.E.R.L. IRISH, R.O.B.E.R.T. WEISS, K.Y.L.E.T. MANDLI, C.A.T.H.E.R.I.N.E.R. JEFFRIES, Coastal Sediments 2023, The Proceedings of the Coastal Sediments 2023, in: Coastal Sediments 2023, WORLD SCIENTIFIC, 2023: pp. 2462–2467. https://doi.org/10.1142/9789811275135_0224.
This paper will discuss the beginnings of a sensitivity analysis of barrier island breaching. The study area of Mantoloking, New Jersey, USA is used as the barrier island breached significantly during Hurricane Sandy in 2012. The numerical model XBeach is used to conduct this study. The study investigates the affects that back-bay currents, water-level timing, and barrier-island configuration have on barrier island breaching. The results will lead to a better understanding of how certain parameters such as dune geometry and channelization of the back-bay contribute to the short-term morphological process of breaching during coastal storm events. - Manifold Approximations via Transported Subspaces: Model Reduction for Transport-Dominated Problems
(with Rim, Donsub Peherstorfer, Benjamin and )
SIAM Journal on Scientific Computing ()
@article{Rim.2023, year = {2023}, title = {{Manifold Approximations via Transported Subspaces: Model Reduction for Transport-Dominated Problems}}, author = {Rim, Donsub and Peherstorfer, Benjamin and Mandli, Kyle T.}, journal = {SIAM Journal on Scientific Computing}, issn = {1064-8275}, doi = {10.1137/20m1316998}, shorttitle = {Manifold Approximations via Transported Subspaces}, pages = {A170--A199}, number = {1}, volume = {45}, copyright = {All rights reserved}, local-url = {file://localhost/Users/kmandli/Documents/Papers%20Library/Rim/Rim_2023.pdf} }[40]D. Rim, B. Peherstorfer, K.T. Mandli, Manifold Approximations via Transported Subspaces: Model Reduction for Transport-Dominated Problems, SIAM Journal on Scientific Computing. 45 (2023) A170–A199. https://doi.org/10.1137/20m1316998.
This work presents a method for constructing online-efficient reduced models of large-scale systems governed by parametrized nonlinear scalar conservation laws. The solution manifolds induced by transport-dominated problems such as hyperbolic conservation laws typically exhibit nonlinear structures, which means that traditional model reduction methods based on linear approximations are inefficient when applied to these problems. In contrast, the approach introduced in this work derives reduced approximations that are nonlinear by explicitly composing global transport dynamics with locally linear approximations of the solution manifolds. A time-stepping scheme evolves the nonlinear reduced models by transporting local approximation spaces along the characteristic curves of the governing equations. The proposed computational procedure allows an offline/online decomposition and is online-efficient in the sense that the complexity of accurately time stepping the nonlinear reduced model is independent of that of the full model. Numerical experiments with transport through heterogeneous media and the Burgers equation show orders of magnitude speedups of the proposed nonlinear reduced models based on transported subspaces compared to traditional linear reduced models and full models. - Advances in Morphodynamic Modeling of Coastal Barriers: A Review
(with Hoagland, Steven W.H. Jeffries, Catherine R. and Irish, Jennifer L. and Weiss, Robert and Kyle and Vitousek, Sean and Johnson, Catherine M. and Cialone, Mary A.)
Journal of Waterway, Port, Coastal, and Ocean Engineering ()
@article{Hoagland.2023, year = {2023}, title = {{Advances in Morphodynamic Modeling of Coastal Barriers: A Review}}, author = {Hoagland, Steven W.H. and Jeffries, Catherine R. and Irish, Jennifer L. and Weiss, Robert and Mandli, Kyle and Vitousek, Sean and Johnson, Catherine M. and Cialone, Mary A.}, journal = {Journal of Waterway, Port, Coastal, and Ocean Engineering}, issn = {0733-950X}, doi = {10.1061/jwped5.wweng-1825}, shorttitle = {Advances in Morphodynamic Modeling of Coastal Barriers}, pages = {03123001}, number = {5}, volume = {149}, copyright = {All rights reserved}, local-url = {file://localhost/Users/kmandli/Documents/Papers%20Library/Hoagland/Hoagland_2023.pdf} }[41]S.W.H. Hoagland, C.R. Jeffries, J.L. Irish, R. Weiss, K. Mandli, S. Vitousek, C.M. Johnson, M.A. Cialone, Advances in Morphodynamic Modeling of Coastal Barriers: A Review, Journal of Waterway, Port, Coastal, and Ocean Engineering. 149 (2023) 03123001. https://doi.org/10.1061/jwped5.wweng-1825.
As scientific understanding of barrier morphodynamics has improved, so has the ability to reproduce observed phenomena and predict future barrier states using mathematical models. To use existing models effectively and improve them, it is important to understand the current state of morphodynamic modeling and the progress that has been made in the field. This manuscript offers a review of the literature regarding advancements in morphodynamic modeling of coastal barrier systems and summarizes current modeling abilities and limitations. Broadly, this review covers both event-scale and long-term morphodynamics. Each of these sections begins with an overview of commonly modeled phenomena and processes, followed by a review of modeling developments. After summarizing the advancements toward the stated modeling goals, we identify research gaps and suggestions for future research under the broad categories of improving our abilities to acquire and access data, furthering our scientific understanding of relevant processes, and advancing our modeling frameworks and approaches.
(with Miura, Yuki and Deodatis, George)
()
@unpublished{Miura.2023yv, year = {2023}, author = {Miura, Yuki and Mandli, Kyle T. and Deodatis, George}, title = {{Modeling Storm Surges with a Bounded Probability Distribution}} }[42]Y. Miura, K.T. Mandli, G. Deodatis, Modeling Storm Surges with a Bounded Probability Distribution, (2023).
(with Miura, Yuki Blackshaw, C. Y. and Zhang, M. S. and and Deodatis, George)
()
@unpublished{Miura.2023, year = {2023}, author = {Miura, Yuki and Blackshaw, C. Y. and Zhang, M. S. and Mandli, Kyle T. and Deodatis, George}, title = {{Flood Risk Assessment for Subway Systems in NYC with GIS-based Subdivision-Redistribution Methodology}} }[43]Y. Miura, C.Y. Blackshaw, M.S. Zhang, K.T. Mandli, G. Deodatis, Flood Risk Assessment for Subway Systems in NYC with GIS-based Subdivision-Redistribution Methodology, (2023).
- Climate Change Contributions to Increasing Compound Flooding Risk in New York City
(with Sarhadi, Ali Rousseau-Rizzi, Raphaël and Kyle and Neal, Jeffrey and Wiper, Michael P and Feldmann, Monika and Emanuel, Kerry)
Bulletin of the American Meteorological Society ()
@article{Sarhadi.2024, year = {2024}, title = {{Climate Change Contributions to Increasing Compound Flooding Risk in New York City}}, author = {Sarhadi, Ali and Rousseau-Rizzi, Raphaël and Mandli, Kyle and Neal, Jeffrey and Wiper, Michael P and Feldmann, Monika and Emanuel, Kerry}, journal = {Bulletin of the American Meteorological Society}, issn = {0003-0007}, doi = {10.1175/bams-d-23-0177.1}, pages = {E337--E356}, number = {2}, volume = {105}, copyright = {http://www.ametsoc.org/PUBSReuseLicenses}, local-url = {file://localhost/Users/kmandli/Documents/Papers%20Library/Sarhadi/Sarhadi_2024.pdf} }[44]A. Sarhadi, R. Rousseau-Rizzi, K. Mandli, J. Neal, M.P. Wiper, M. Feldmann, K. Emanuel, Climate Change Contributions to Increasing Compound Flooding Risk in New York City, Bulletin of the American Meteorological Society. 105 (2024) E337–E356. https://doi.org/10.1175/bams-d-23-0177.1.
Abstract Efforts to meaningfully quantify the changes in coastal compound surge- and rainfall-driven flooding hazard associated with tropical cyclones (TCs) and extratropical cyclones (ETCs) in a warming climate have increased in recent years. Despite substantial progress, however, obtaining actionable details such as the spatially and temporally varying distribution and proximal causes of changing flooding hazard in cities remains a persistent challenge. Here, for the first time, physics-based hydrodynamic flood models driven by rainfall and storm surge simultaneously are used to estimate the magnitude and frequency of compound flooding events. We apply this to the particular case of New York City. We find that sea level rise (SLR) alone will increase the TC and ETC compound flooding hazard more significantly than changes in storm climatology as the climate warms. We also project that the probability of destructive Sandy-like compound flooding will increase by up to 5 times by the end of the century. Our results have strong implications for climate change adaptation in coastal communities.
Submitted
(with Miura, Yuki and Deodatis, George)
()
@unpublished{Miura.2023yv, year = {2023}, author = {Miura, Yuki and Mandli, Kyle T. and Deodatis, George}, title = {{Modeling Storm Surges with a Bounded Probability Distribution}} }[1]Y. Miura, K.T. Mandli, G. Deodatis, Modeling Storm Surges with a Bounded Probability Distribution, (2023).
(with Miura, Yuki Blackshaw, C. Y. and Zhang, M. S. and and Deodatis, George)
()
@unpublished{Miura.2023, year = {2023}, author = {Miura, Yuki and Blackshaw, C. Y. and Zhang, M. S. and Mandli, Kyle T. and Deodatis, George}, title = {{Flood Risk Assessment for Subway Systems in NYC with GIS-based Subdivision-Redistribution Methodology}} }[2]Y. Miura, C.Y. Blackshaw, M.S. Zhang, K.T. Mandli, G. Deodatis, Flood Risk Assessment for Subway Systems in NYC with GIS-based Subdivision-Redistribution Methodology, (2023).