skip to main content

Title: Is Computational Oceanography Coming of Age?
Abstract Computational oceanography is the study of ocean phenomena by numerical simulation, especially dynamical and physical phenomena. Progress in information technology has driven exponential growth in the number of global ocean observations and the fidelity of numerical simulations of the ocean in the past few decades. The growth has been exponentially faster for ocean simulations, however. We argue that this faster growth is shifting the importance of field measurements and numerical simulations for oceanographic research. It is leading to the maturation of computational oceanography as a branch of marine science on par with observational oceanography. One implication is that ultraresolved ocean simulations are only loosely constrained by observations. Another implication is that barriers to analyzing the output of such simulations should be removed. Although some specific limits and challenges exist, many opportunities are identified for the future of computational oceanography. Most important is the prospect of hybrid computational and observational approaches to advance understanding of the ocean.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Bulletin of the American Meteorological Society
Page Range / eLocation ID:
E1481 to E1493
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Surface boundaries in supercells have been suspected of being important in the arrangement and concentration of vorticity for the development and intensification of tornadoes, but there has been little attention given to the effects of the underlying surface roughness on their behavior. This study investigates the impact of surface drag on the structure and evolution of these boundaries, their associated distribution of near-surface vorticity, and tornadogenesis and maintenance. Comparisons between idealized simulations without and with drag introduced in the mature stage of the storm prior to tornadogenesis reveal that the inclusion of surface drag substantially alters the low-level structure, particularly with respect to the number, location, and intensity of surface convergence boundaries. Substantial drag-generated horizontal vorticity induces rotor structures near the surface associated with the convergence boundaries in both the forward and rear flanks of the storm. Stretching of horizontal vorticity and subsequent tilting into the vertical along the convergence boundaries lead to elongated positive vertical vorticity sheets on the ascending branch of the rotors and the opposite on the descending branch. The larger near-surface pressure deficit associated with the faster development of the near-surface cyclone when drag is active creates a downward dynamic vertical pressure gradient force that suppresses vertical growth, leading to a weaker and wider tornado detached from the surrounding convergence boundaries. A conceptual model of the low-level structure of the tornadic supercell is presented that focuses on the contribution of surface drag, with the aim of adding more insight and complexity to previous conceptual models.

    Significance Statement

    Tornado development is sensitive to near-surface processes, including those associated with front-like boundaries between regions of airflow within the parent storm. However, observations and theory are insufficient to understand these phenomena, and numerical simulation remains vital. In our simulations, we find that a change in a parameter that controls how much the near-surface winds are reduced by friction (or drag) can substantially alter the storm behavior and tornado potential. We investigate how surface drag affects the low-level storm structure, the distribution of regions of near-surface rotation, and the development of tornadoes within the simulation. Our results provide insight into the role of surface drag and lead to an improved conceptual model of the near-surface structure of a tornadic storm.

    more » « less
  2. null (Ed.)
    Abstract Antarctic sea ice extent (SIE) has slightly increased over the satellite observational period (1979 to the present) despite global warming. Several mechanisms have been invoked to explain this trend, such as changes in winds, precipitation, or ocean stratification, yet there is no widespread consensus. Additionally, fully coupled Earth system models run under historic and anthropogenic forcing generally fail to simulate positive SIE trends over this time period. In this work, we quantify the role of winds and Southern Ocean SSTs on sea ice trends and variability with an Earth system model run under historic and anthropogenic forcing that nudges winds over the polar regions and Southern Ocean SSTs north of the sea ice to observations from 1979 to 2018. Simulations with nudged winds alone capture the observed interannual variability in SIE and the observed long-term trends from the early 1990s onward, yet for the longer 1979–2018 period they simulate a negative SIE trend, in part due to faster-than-observed warming at the global and hemispheric scale in the model. Simulations with both nudged winds and SSTs show no significant SIE trends over 1979–2018, in agreement with observations. At the regional scale, simulated sea ice shows higher skill compared to the pan-Antarctic scale both in capturing trends and interannual variability in all nudged simulations. We additionally find negligible impact of the initial conditions in 1979 on long-term trends. 
    more » « less
  3. Arctic Ocean properties and processes are highly relevant to the regional and global coupled climate system, yet still scarcely observed, especially in winter. Team OCEAN conducted a full year of physical oceanography observations as part of the Multidisciplinary drifting Observatory for the Study of the Arctic Climate (MOSAiC), a drift with the Arctic sea ice from October 2019 to September 2020. An international team designed and implemented the program to characterize the Arctic Ocean system in unprecedented detail, from the seafloor to the air-sea ice-ocean interface, from sub-mesoscales to pan-Arctic. The oceanographic measurements were coordinated with the other teams to explore the ocean physics and linkages to the climate and ecosystem. This paper introduces the major components of the physical oceanography program and complements the other team overviews of the MOSAiC observational program. Team OCEAN’s sampling strategy was designed around hydrographic ship-, ice- and autonomous platform-based measurements to improve the understanding of regional circulation and mixing processes. Measurements were carried out both routinely, with a regular schedule, and in response to storms or opening leads. Here we present along-drift time series of hydrographic properties, allowing insights into the seasonal and regional evolution of the water column from winter in the Laptev Sea to early summer in Fram Strait: freshening of the surface, deepening of the mixed layer, increase in temperature and salinity of the Atlantic Water. We also highlight the presence of Canada Basin deep water intrusions and a surface meltwater layer in leads. MOSAiC most likely was the most comprehensive program ever conducted over the ice-covered Arctic Ocean. While data analysis and interpretation are ongoing, the acquired datasets will support a wide range of physical oceanography and multi-disciplinary research. They will provide a significant foundation for assessing and advancing modeling capabilities in the Arctic Ocean. 
    more » « less
  4. null (Ed.)
    Abstract. We present a newly developed upper-thermocline, open-ocean biogeochemical flux model that is complex and flexible enough to capture open-ocean ecosystem dynamics but reduced enough to incorporate into highly resolved numerical simulations and parameter optimization studies with limited additional computational cost. The model, which is derived from the full 56-state-variable Biogeochemical Flux Model (BFM56; Vichi et al., 2007), follows a biological and chemical functional group approach and allows for the development of critical non-Redfield nutrient ratios. Matter is expressed in units of carbon, nitrogen, and phosphate, following techniques used in more complex models. To reduce the overall computational cost and to focus on upper-thermocline, open-ocean, and non-iron-limited or non-silicate-limited conditions, the reduced model eliminates certain processes, such as benthic, silicate, and iron influences, and parameterizes others, such as the bacterial loop. The model explicitly tracks 17 state variables, divided into phytoplankton, zooplankton, dissolved organic matter, particulate organic matter, and nutrient groups. It is correspondingly called the Biogeochemical Flux Model 17 (BFM17). After describing BFM17, we couple it with the one-dimensional Princeton Ocean Model for validation using observational data from the Sargasso Sea. The results agree closely with observational data, giving correlations above 0.85, except for chlorophyll (0.63) and oxygen (0.37), as well as with corresponding results from BFM56, with correlations above 0.85, except for oxygen (0.56), including the ability to capture the subsurface chlorophyll maximum and bloom intensity. In comparison to previous models of similar size, BFM17 provides improved correlations between several model output fields and observational data, indicating that reproduction of in situ data can be achieved with a low number of variables, while maintaining the functional group approach. Notable additions to BFM17 over similar complexity models are the explicit tracking of dissolved oxygen, allowance for non-Redfield nutrient ratios, and both dissolved and particulate organic matter, all within the functional group framework. 
    more » « less
  5. Abstract Quasi-decadal climate of the Kuroshio Extension (KE) is pivotal to understanding the North Pacific coupled ocean–atmosphere dynamics and their predictability. Recent observational studies suggest that extratropical-tropical coupling between the KE and the central tropical Pacific El Niño Southern Oscillation (CP-ENSO) leads to the observed preferred decadal time-scale of Pacific climate variability. By combining reanalysis data with numerical simulations from a high-resolution climate model and a linear inverse model (LIM), we confirm that KE and CP-ENSO dynamics are linked through extratropical-tropical teleconnections. Specifically, the atmospheric response to the KE excites Meridional Modes that energize the CP-ENSO (extratropicstropics), and in turn, CP-ENSO teleconnections energize the extratropical atmospheric forcing of the KE (tropicsextratropics). However, both observations and the model show that the KE/CP-ENSO coupling is non-stationary and has intensified in recent decades after the mid-1980. Given the short length of the observational and climate model record, it is difficult to attribute this shift to anthropogenic forcing. However, using a large-ensemble of the LIM we show that the intensification in the KE/CP-ENSO coupling after the mid-1980 is significant and linked to changes in the KE atmospheric downstream response, which exhibit a stronger imprint on the subtropical winds that excite the Pacific Meridional modes and CP-ENSO. 
    more » « less