High-Resolution In Situ Oxygen-Argon Studies of Surface Biological and Physical Processes in the Polar Oceans

The Arctic Ocean and Western Antarctic Peninsula (WAP) are the fastest warming regions on the planet and are undergoing rapid climate and ecosystem changes. Until we can fully resolve the coupling between biological and physical processes we cannot predict how warming will influence carbon cycling and ecosystem function and structure in these sensitive and climactically important regions. My dissertation centers on the use of high-resolution measurements of surface dissolved gases, primarily O2 and Ar, as tracers or physical and biological functioning that we measure underway using an optode and Equilibrator Inlet Mass Spectrometry (EIMS). Total O2 measurements are common throughout the historical and autonomous record but are influenced by biological (net metabolic balance) and physical (temperature, salinity, pressure changes, ice melt/freeze, mixing, bubbles and diffusive gas exchange) processes. We use Ar, an inert gas with similar solubility properties to O2, to devolve distinct records of biological (O2/Ar) and physical (Ar) oxygen. These high-resolution measurements that expose intersystem coupling and submesoscale variability were central to studies in the Arctic Ocean, WAP and open Southern Ocean that make up this dissertation.

Key findings of this work include the documentation of under ice and ice-edge blooms and basin scale net sea ice freeze/melt processes in the Arctic Ocean. In the WAP O2 and pCO2 are both biologically driven and net community production (NCP) variability is controlled by Fe and light availability tied to glacial and sea ice meltwater input. Further, we present a feasibility study that shows the ability to use modeled Ar to derive NCP from total O2 records. This approach has the potential to unlock critical carbon flux estimates from historical and autonomous O2 measurements in the global oceans.

A New Method for Modeling Free Surface Flows and Fluid-structure Interaction with Ocean Applications

The computational modeling of ocean waves and ocean-faring devices poses numerous challenges. Among these are the need to stably and accurately represent both the fluid-fluid interface between water and air as well as the fluid-structure interfaces arising between solid devices and one or more fluids. As techniques are developed to stably and accurately balance the interactions between fluid and structural solvers at these boundaries, a similarly pressing challenge is the development of algorithms that are massively scalable and capable of performing large-scale three-dimensional simulations on reasonable time scales. This dissertation introduces two separate methods for approaching this problem, with the first focusing on the development of sophisticated fluid-fluid interface representations and the second focusing primarily on scalability and extensibility to higher-order methods.

We begin by introducing the narrow-band gradient-augmented level set method (GALSM) for incompressible multiphase Navier-Stokes flow. This is the first use of the high-order GALSM for a fluid flow application, and its reliability and accuracy in modeling ocean environments is tested extensively. The method demonstrates numerous advantages over the traditional level set method, among these a heightened conservation of fluid volume and the representation of subgrid structures.

Next, we present a finite-volume algorithm for solving the incompressible Euler equations in two and three dimensions in the presence of a flow-driven free surface and a dynamic rigid body. In this development, the chief concerns are efficiency, scalability, and extensibility (to higher-order and truly conservative methods). These priorities informed a number of important choices: The air phase is substituted by a pressure boundary condition in order to greatly reduce the size of the computational domain, a cut-cell finite-volume approach is chosen in order to minimize fluid volume loss and open the door to higher-order methods, and adaptive mesh refinement (AMR) is employed to focus computational effort and make large-scale 3D simulations possible. This algorithm is shown to produce robust and accurate results that are well-suited for the study of ocean waves and the development of wave energy conversion (WEC) devices.