Search NASA⌕ Search

Engineering topics

Roberts, Andrew F.

Publications and source records attributed to Roberts, Andrew F..

Southern Ocean polynyas and dense water formation in a high-resolution, coupled Earth system model

Antarctic coastal polynyas produce dense shelf water, a primary source of Antarctic Bottom Water that contributes to the global overturning circulation. This paper investigates Antarctic dense water formation in the high-resolution version of the Energy Exascale Earth System Model (E3SM-HR). The model is able to reproduce the main Antarctic coastal polynyas, although the polynyas are smaller in area compared to observations. E3SM-HR also simulates several occurrences of open-ocean polynyas (OOPs) in the Weddell Sea at a higher rate than what the last 50 years of the satellite sea ice observational record suggests, but similarly to other high-resolution Earth system model simulations. Furthermore, the densest water masses in the model are formed within the OOPs rather than on the continental shelf as is typically observed. Biases related to the lack of dense water formation on the continental shelf are associated with overly strong atmospheric polar easterlies, which lead to a strong Antarctic Slope Front and too little exchange between on- and off-continental shelf water masses. Strong polar easterlies also produce excessive southward Ekman transport, causing a build-up of sea ice over the continental shelf and enhanced ice melting in the summer season. This, in turn, produces water masses on the continental shelf that are overly fresh and less dense relative to observations. Our results indicate that high resolution alone is insufficient for models to properly reproduce Antarctic dense water; the large-scale polar atmospheric circulation around Antarctica must also be accurately simulated.

54 ENVIRONMENTAL SCIENCES↗

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↗

MPAS - Ocean Simulation Quality for Variable-Resolution North American Coastal Meshes

Climate model components utilizing unstructured meshes enable variableresolution, regionally enhanced simulations within global domains. Here we investigate the relationship between mesh quality and simulation statistics using the JIGSAW unstructured meshing library and the Model for Prediction Across ScalesOcean (MPASOcean) with a focus on Gulf Stream dynamics. In the base configuration, the refined region employs 8 km cells that extend 400 km from the coast of North America. This coastal refined region is embedded within a lowresolution global domain, with cell size varying latitudinally between 30 and 60 km. The resolution transition region between the refined region and background mesh is 600 km wide. Three sensitivity tests are conducted: 1) the quality of meshes is intentionally degraded so that horizontal cells are progressively more distorted; 2) the transition region from high to low resolution is steepened; and 3) resolution of the coastal refinement region is varied from 30 km to 8 km. Overall, the ocean simulations are shown to be robust to mesh resolution and quality alterations. Meshes that are substantially degraded still produce realistic currents, with Southern Ocean transports within 0.4% and Gulf Stream transports within 12% of highquality mesh results. The narrowest transition case of 100 km did not produce any spurious effects. Refined regions with high resolution produce eddy kinetic energy and sea surface height variability that are similar to the highresolution reference simulation. These results provide heuristics for the design criteria of variableresolution climate model domains.

Hoch, Kristin E.↗