Engineering topics

Wirasaet, Damrongsak

Publications and source records attributed to Wirasaet, Damrongsak.

4 records

DOE OSTI2023

Barotropic tides in MPAS-Ocean (E3SM V2): impact of ice shelf cavities

Abstract. Oceanic tides are seldom represented in Earth system models (ESMs) owing to the need for high horizontal resolution to accurately represent the associated barotropic waves close to coasts. This paper presents results of tides implemented in the Model for Prediction Across Scales–Ocean or MPAS-Ocean, which is the ocean component within the U.S. Department of Energy developed Energy Exascale Earth System Model (E3SM). MPAS-Ocean circumvents the limitation of low resolution using unstructured global meshing. We are at this stage simulating the largest semidiurnal (M2, S2, N2) and diurnal (K1, O1) tidal constituents in a single-layer version of MPAS-O. First, we show that the tidal constituents calculated using MPAS-Ocean closely agree with the results of the global tidal prediction model TPXO8 when suitably tuned topographic wave drag and bottom drag coefficients are employed. Thereafter, we present the sensitivity of global tidal evolution due to the presence of Antarctic ice shelf cavities. The effect of ice shelves on the amplitude and phase of tidal constituents are presented. Lower values of complex errors (with respect to TPXO8 results) for the M2 tidal constituents are observed when the ice shelf is added in the simulations, with particularly strong improvement in the Southern Ocean. Our work points towards future research with varying Antarctic ice shelf geometries and sea ice coupling that might lead to better comparison and prediction of tides and thus better prediction of sea-level rise and also the future climate variability.

58 GEOSCIENCES↗

DOE OSTI2023

Establishing connectivity patterns of eastern oysters ( Crassostrea virginica ) on regional oceanographic scales

Localized depletion of marine and estuarine populations often results from large-scale natural and anthropogenic disturbances (e.g., hurricanes, oil spills) as well as overharvest of fisheries resources. Understanding how such localized depletions may affect populations at larger regional scales requires knowledge of connectivity among local populations within the larger regional landscape or metapopulation. Efforts to restore populations following such decreases require similar knowledge. During the 2010 Deepwater Horizon oil spill, dramatic declines of oysters occurred throughout estuaries west and east of the Mississippi River. We examined trajectories of particles, which were parameterized to mimic oyster larvae, using the ADvanced CIRCulation (ADCIRC; https://adcirc.org ) model to evaluate potential connectivity within and among embayments from Western Louisiana to Alabama. Patterns of larval settlement, which we defined as the intersection of a larval particle with known or expected oyster habitat at any point 13–21 days post-release, reflected much greater local contributions, with 50%–90% of settled particles originating within the same subdivision of embayments. Exchange among subdivision was much less (0%–40%) and settlement originating from outside the embayment of release was trivial under most scenarios (0%–14%). Connectivity between adjacent basins was greatest for larvae released in the southern portions of the embayment, whereas connectivity among nonadjacent basins was not predicted under the scenarios modeled. Because most local populations are relatively isolated on ecological timescales, areas suffering from extensive local depletion are likely to require extensive time to recover due to the lack of larval subsidy from the overall regional population. Restoration would require building stepping-stone populations or reefs within and among basin to restore a high degree of connectivity.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Scalable self attraction and loading calculations for unstructured ocean tide models

Self attraction and earth-loading effects are important for accurately modeling global tides. A common approach of handling this forcing is to expand mass anomalies into spherical harmonics, which are scaled by load Love numbers to account for elastic earth deformation. We investigate two different approaches to perform these calculations for ocean models that employ unstructured meshes and distributed memory parallelization. The first approach leverages a highly efficient spherical harmonics library, but requires all-to-one and one-to-all communications and interpolation operations between the unstructured and a structured mesh. This approach is compared to a parallel algorithm that computes the spherical harmonic transformations directly on the unstructured mesh with an all-reduce communication. Here, our results show that although the unstructured mesh calculations are more expensive, the scalability of the unstructured mesh approach allows for more efficient spherical harmonics transforms for high-resolution meshes and large processor counts. This methodology enables the efficient inclusion of tidal dynamics large-scale Earth system model simulations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Global Barotropic Tide Modeling Using Inline Self‐Attraction and Loading in MPAS‐Ocean

Abstract We examine ocean tides in the barotropic version of the Model for Prediction Across Scales (MPAS‐Ocean), the ocean component of the Department of Energy Earth system model. We focus on four factors that affect tidal accuracy: self‐attraction and loading (SAL), model resolution, details of the underlying bathymetry, and parameterized topographic wave drag. The SAL term accounts for the tidal loading of Earth's crust and the self‐gravitation of the ocean and the load‐deformed Earth. A common method for calculating SAL is to decompose mass anomalies into their spherical harmonic constituents. Here, we compare a scalar SAL approximation versus an inline SAL using a fast spherical harmonic transform package. Wave drag accounts for energy lost by breaking internal tides that are produced by barotropic tidal flow over topographic features. We compare a series of successively finer quasi‐uniform resolution meshes (62.9, 31.5, 15.7, and 7.87 km) to a variable resolution (45 to 5 km) configuration. We ran MPAS‐Ocean in a single‐layer barotropic mode forced by five tidal constituents. The 45 to 5 km variable resolution mesh obtained the best total root‐mean‐square error (5.4 cm) for the deep ocean ( 1,000 m) tide compared to TPXO8 and ran twice as fast as the quasi‐uniform 8 km mesh, which had an error of 5.8 cm. This error is comparable to those found in other forward (non‐assimilative) ocean tide models. In future work, we plan to use MPAS‐Ocean to study tidal interactions with other Earth system components, and the tidal response to climate change.

54 ENVIRONMENTAL SCIENCES↗