Richard A. Luettich

A Basin- to Channel-Scale Unstructured Grid Hurricane Storm Surge Model Applied to Southern Louisiana

Abstract Southern Louisiana is characterized by low-lying topography and an extensive network of sounds, bays, marshes, lakes, rivers, and inlets that permit widespread inundation during hurricanes. A basin- to channel-scale implementation of the Advanced Circulation (ADCIRC) unstructured grid hydrodynamic model has been developed that accurately simulates hurricane storm surge, tides, and river flow in this complex region. This is accomplished by defining a domain and computational resolution appropriate for the relevant processes, specifying realistic boundary conditions, and implementing accurate, robust, and highly parallel unstructured grid numerical algorithms. The model domain incorporates the western North Atlantic, the Gulf of Mexico, and the Caribbean Sea so that interactions between basins and the shelf are explicitly modeled and the boundary condition specification of tidal and hurricane processes can be readily defined at the deep water open boundary. The unstructured grid enables highly refi...

A High-Resolution Coupled Riverine Flow, Tide, Wind, Wind Wave, and Storm Surge Model for Southern Louisiana and Mississippi. Part I: Model Development and Validation

Abstract A coupled system of wind, wind wave, and coastal circulation models has been implemented for southern Louisiana and Mississippi to simulate riverine flows, tides, wind waves, and hurricane storm surge in the region. The system combines the NOAA Hurricane Research Division Wind Analysis System (H*WIND) and the Interactive Objective Kinematic Analysis (IOKA) kinematic wind analyses, the Wave Model (WAM) offshore and Steady-State Irregular Wave (STWAVE) nearshore wind wave models, and the Advanced Circulation (ADCIRC) basin to channel-scale unstructured grid circulation model. The system emphasizes a high-resolution (down to 50 m) representation of the geometry, bathymetry, and topography; nonlinear coupling of all processes including wind wave radiation stress-induced set up; and objective specification of frictional parameters based on land-cover databases and commonly used parameters. Riverine flows and tides are validated for no storm conditions, while winds, wind waves, hydrographs, and high wa...

A High-Resolution Coupled Riverine Flow, Tide, Wind, Wind Wave, and Storm Surge Model for Southern Louisiana and Mississippi. Part II: Synoptic Description and Analysis of Hurricanes Katrina and Rita

Abstract Hurricanes Katrina and Rita were powerful storms that impacted southern Louisiana and Mississippi during the 2005 hurricane season. In Part I, the authors describe and validate a high-resolution coupled riverine flow, tide, wind, wave, and storm surge model for this region. Herein, the model is used to examine the evolution of these hurricanes in more detail. Synoptic histories show how storm tracks, winds, and waves interacted with the topography, the protruding Mississippi River delta, east–west shorelines, manmade structures, and low-lying marshes to develop and propagate storm surge. Perturbations of the model, in which the waves are not included, show the proportional importance of the wave radiation stress gradient induced setup.

The influence of domain size on the response characteristics of a hurricane storm surge model

The influence of domain size on boundary condition specification and on computed storm surge response is investigated. Storm surge response along the Florida shelf in the Gulf of Mexico due to Hurricane Kate is examined over three domains using two different open ocean boundary forcing functions, a still water (or zero elevation) condition and an inverted barometer condition which accounts for the atmospheric pressure component of the meteorological forcing. The first domain is relatively small and is situated primarily on the continental shelf in the region of intense storm surge generation. A second domain includes the entire Gulf of Mexico basin. The final domain covers the Gulf of Mexico, contiguous basins, and extends out into the deep Atlantic Ocean. The computed storm surge response indicates that the small domain is inadequate, since cross-shelf boundaries are in regions of significant storm surge generation where surge and therefore boundary conditions are not known a priori. Also, the behavior of resonant modes that are physically excited within the Gulf of Mexico due to the passage of the hurricane is unknown at the boundaries of this small domain. The domain that includes the entire Gulf of Mexico captures the primary storm surge well but may not correctly model resonant modes. In general, these resonant modes are difficult to accurately set up by boundary condition specification, since they may be dependent on interactions between the Gulf and contiguous basins. The primary storm surge response as well as resonant modes excited by the storm are best represented using a domain which encompasses the western North Atlantic Ocean, the Caribbean Sea, and the Gulf of Mexico. This domain with deep Atlantic Ocean boundaries facilitates simple boundary condition specification and minimizes the influence of boundary conditions on storm surge generation in coastal regions. Basin resonant modes and basin to basin interactions are also captured.

Hurricane Gustav (2008) Waves and Storm Surge: Hindcast, Synoptic Analysis, and Validation in Southern Louisiana

AbstractHurricane Gustav (2008) made landfall in southern Louisiana on 1 September 2008 with its eye never closer than 75 km to New Orleans, but its waves and storm surge threatened to flood the city. Easterly tropical-storm-strength winds impacted the region east of the Mississippi River for 12–15 h, allowing for early surge to develop up to 3.5 m there and enter the river and the city’s navigation canals. During landfall, winds shifted from easterly to southerly, resulting in late surge development and propagation over more than 70 km of marshes on the river’s west bank, over more than 40 km of Caernarvon marsh on the east bank, and into Lake Pontchartrain to the north. Wind waves with estimated significant heights of 15 m developed in the deep Gulf of Mexico but were reduced in size once they reached the continental shelf. The barrier islands further dissipated the waves, and locally generated seas existed behind these effective breaking zones.The hardening and innovative deployment of gauges since Hur...

Abundance of estuarine crab larvae is associated with tidal hydrologic variables

Abundances of brachyuran megalopae and juveniles were measured throughout consecutive tidal cycles during six 2 to 3 d sampling periods in summer 1992, and associated with rates of change of tidal hydrologic variables in the Newport River Estuary. Current speeds and rates of pressure change fitted sinusoidal (tidal) models well; however, rates of salinity and temperature change did not. Analysis of plankton samples taken during spring and neap tides showed peak abundances during nighttime rising tides for all taxonomic groups: Callinectes sapidus, Uca spp., Xanthidae, and Pinnixa spp. megalopae, and Pinnotheres spp. juveniles. Megalopal and juvenile abundances from time-intensive sampling were related to rates of changes in the hydrologic variables using stepwise logistic regression. No hydrologic variable accounted well for the presence of Uca spp. megalopae. Megalopal presence was best predicted by current speed for Pinnixa spp. megalopae, and rates of changes in pressure for xanthid megalopae and Pinnotheres spp. juveniles, and salinity for C. sapidus megalopae. These variables might act as cues causing megalopae to ascend into the water column at a particular point in the flooding tide, and subsequently descend to or near the bottom prior to ebb flow. In this way, larvae which develop on the continental shelf or lower estuary undergo transport up the estuary by behaviorally altering their swimming activity and depth concurrent with tidal changes.

Grid convergence studies for the prediction of hurricane storm surge

SUMMARY The focus of this paper is a systematic determination of the relationship between grid resolution and errors associated with computations of hurricane storm surge. A grid structure is sought that provides the spatial resolution necessary to capture pertinent storm surge physics and does not overdiscretize. A set of numerical experiments simulating storm surge generation over 14 grid discretizations of idealized domains examines the influence of grid spacing, shoreline detail, coastline resolution and characteristics of the meteorological forcing on storm surge computations. Errors associated with a given grid are estimated using a Richardson-based error estimator. Analysis of the magnitude and location of estimated errors indicates that underresolution on the continental shelf leads to significant overprediction of the primary storm surge. In deeper waters, underresolution causes smearing or damping of the inverted barometer forcing function, which in turn results in underprediction of the surge elevation. In order to maintain a specified error level throughout the duration of the storm, the highest grid resolution is required on the continental shelf and particularly in nearshore areas. The disparity of discretization requirements between deep waters and coastal regions is best met using a graded grid. Application of the graded gridding strategy to the hindcast of Hurricane Camille reinforces the necessity of using a grid that has high levels of resolution in nearshore regions and areas of complex coastal geometry. # 1998 John Wiley & Sons, Ltd. Int. J. Numer. Meth. Fluids, 26: 369‐401 (1998).

Ecological Response to Hurricane Events in the Pamlico Sound System, North Carolina, and Implications for Assessment and Management in a Regime of Increased Frequency

Since the mid 1990s, the Atlantic and Gulf Coast regions have experienced a dramatic increase in the number of hurricane landfalls. In eastern North Carolina alone, eight hurricances have affected the coast in the past 9 years. These storms have exhibited individualistic hydrologic, nutrient, and sediment loading effects and represent a formidable challenge to nutrient management aimed at reducing eutrophication in the Pamlico Sound and its estuarine tributaries. Different rainfall amounts among hurricanes lead to variable freshwater and nutrient discharge and variable nutrient, organic matter, and sediment enrichment. These enrichments differentially affected physical and chemical properties (salinity, water residence time, transparency, stratification, dissolved oxygen), phytoplankton primary production, and phytoplankton community composition. Contrasting ecological responses were accompanied, by changes in nutrient and oxygen cycling, habitat, and higher trophic levels, including different direct effects on fish populations. Floodwaters from the two largest hurricances, Fran (1996) and Floyd (1999), exerted, multi-month to multi-annual effects on hydrology, nutrient loads, productivity, and biotic composition. Relatively low rainfall coastal hurricanes like Isabel (2003) and Ophelia (2005) caused strong vertical mixing and storm surges, but relatively minor hydrologic and nutrient effects. Both hydrologic loading and wind forcing are important drivers and must be integrated with nutrient loading in assessing short-term and long-term ecological effects of these storms. These climatic forcings cannot be managed but should be considered in the development of water quality management strategies for these and other large estuarine ecosystems faced with increasing frequencies and intensities of hurricane activity.

A three-dimensional model of tidal currents in the mouth of the Tagus estuary

Abstract Three-dimensional simulations of barotropic tides are used to investigate the complex circulation at the mouth of the Tagus estuary. Comparison with field data shows that elevations are well represented in the main area of interest, but velocities are slightly over-predicted due to a simplified treatment of intertidal flats. Model results show the existence of strong velocity phase lags (up to 1.75 hours for the semi-diurnal constituents) between the shallow areas and the deep channel. These phase lags are partially responsible for the generation by advective accelerations of a strong residual velocity field (velocities reach 0.5 m s −1 ) with well defined eddies. The interaction between the residual and tidal velocity fields in turn generates strong chaotic stirring. Localized sigma coordinates (LSC), a recently proposed method which allows the number of nodes per vertical to vary horizontally, are used for the first time in a three-dimensional application. A previously proposed criterion for the horizontal distribution of the total number of vertical nodes is shown to be inadequate when advective accelerations are important. However, with an alternative criterion, the use of LSC reduces maximum errors by a factor of two relative to traditional sigma coordinates.

A comparison of the nonlinear frictional characteristics of two-dimensional and three-dimensional models of a shallow tidal embayment

The nonlinear frictional behavior of two-dimensional (2-D) and three-dimensional (3-D) models are compared in this study of tides in the Bight of Abaco. The shallow depths and the existence of an extensive set of tidal elevation data (five astronomical and two overtide constituents at 25 stations) from Filloux and Snyder (1979) offer an excellent opportunity to compare the effects of different frictional formulations. In addition, previous modeling efforts in the bight have consistently overpredicted the M6 and generally overdamped the O1, K1, and S2 tides. The results indicate that although the 2-D and 3-D models may be calibrated to produce nearly identical responses for the dominant M2 tide, there are systematic differences in the responses of the primary overrides. These differences are explained using analytical expansions of the friction terms and are shown to be due to differences in the terms that are nonlinear in velocity and in water level. The investigation concludes that the overgeneration of M6 and the overdamping of secondary astronomical tides will occur in 3-D models as well as 2-D models. Although several causes for these problems were considered, improvement in these constituents could be achieved only by modifying the standard quadratic friction or flow-dependent eddy viscosity relations to reduce the nonlinear frictional effect relative to the linear frictional effect. The required modifications suggest the presence of a constant background velocity, residual turbulence field, or possibly the need for a more advanced frictional closure.