Development of the CSOMIO Coupled Ocean-Oil-Sediment- Biology Model

The fate and dispersal of oil in the ocean is dependent upon ocean dynamics, as well as transformations resulting from the interaction with the microbial community and suspended particles. These interaction processes are parameterized in many models limiting their ability to accurately simulate the fate and dispersal of oil for subsurface oil spill events. This paper presents a coupled ocean-oil-biology-sediment modeling system developed by the Consortium for Simulation of Oil-Microbial Interactions in the Ocean (CSOMIO) project. A key objective of the CSOMIO project was to develop and evaluate a modeling framework for simulating oil in the marine environment, including its interaction with microbial food webs and sediments. The modeling system developed is based on the Coupled Ocean-Atmosphere-Wave-Sediment Transport model (COAWST). Central to CSOMIO’s coupled modeling system is an oil plume model coupled to the hydrodynamic model (Regional Ocean Modeling System, ROMS). The oil plume model is based on a Lagrangian approach that describes the oil plume dynamics including advection and diffusion of individual Lagrangian elements, each representing a cluster of oil droplets. The chemical composition of oil is described in terms of three classes of compounds: saturates, aromatics, and heavy oil (resins and asphaltenes). The oil plume model simulates the rise of oil droplets based on ambient ocean flow and density fields, as well as the density and size of the oil droplets. The oil model also includes surface evaporation and surface wind drift. A novel component of the CSOMIO model is two-way Lagrangian-Eulerian mapping of the oil characteristics. This mapping is necessary for implementing interactions between the ocean-oil module and the Eulerian sediment and biogeochemical modules. The sediment module is a modification of the Community Sediment Transport Modeling System. The module simulates formation of oil-particle aggregates in the water column. The biogeochemical module simulates microbial communities adapted to the local environment and to elevated concentrations of oil components in the water column. The sediment and biogeochemical modules both reduce water column oil components. This paper provides an overview of the CSOMIO coupled modeling system components and demonstrates the capabilities of the modeling system in the test experiments

Skill metrics for evaluation and comparison of sea ice models

© The Author(s), 2015. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Journal of Geophysical Research: Oceans 120 (2015): 5910–5931, doi:10.1002/2015JC010989.Five quantitative methodologies (metrics) that may be used to assess the skill of sea ice models against a control field are analyzed. The methodologies are Absolute Deviation, Root-Mean-Square Deviation, Mean Displacement, Hausdorff Distance, and Modified Hausdorff Distance. The methodologies are employed to quantify similarity between spatial distribution of the simulated and control scalar fields providing measures of model performance. To analyze their response to dissimilarities in two-dimensional fields (contours), the metrics undergo sensitivity tests (scale, rotation, translation, and noise). Furthermore, in order to assess their ability to quantify resemblance of three-dimensional fields, the metrics are subjected to sensitivity tests where tested fields have continuous random spatial patterns inside the contours. The Modified Hausdorff Distance approach demonstrates the best response to tested differences, with the other methods limited by weak responses to scale and translation. Both Hausdorff Distance and Modified Hausdorff Distance metrics are robust to noise, as opposed to the other methods. The metrics are then employed in realistic cases that validate sea ice concentration fields from numerical models and sea ice mean outlook against control data and observations. The Modified Hausdorff Distance method again exhibits high skill in quantifying similarity between both two-dimensional (ice contour) and three-dimensional (ice concentration) sea ice fields. The study demonstrates that the Modified Hausdorff Distance is a mathematically tractable and efficient method for model skill assessment and comparison providing effective and objective evaluation of both two-dimensional and three-dimensional sea ice characteristics across data sets.U.S. National Science Foundation (NSF) Grant Number: PLR-0804017, NASA JPL OVWST, Bureau of Ocean Energy Management (BOEM), FSU Grant Number: M12PC00003, NSF Grant Numbers: projects PLR-0804010 , PLR-1313614 , PLR-1203720, BP/The Gulf of Mexico Research Initiative Grant Number: SA12-12, GoMRI-008, DoD High Performance Computing Modernization Progra

Time scales of the Greenland freshwater anomaly in the subpolar North Atlantic

Author Posting. © American Meteorological Society, 2021. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 34(22), (2021): 8971–8987, https://doi.org/10.1175/JCLI-D-20-0610.1.The impact of increasing Greenland freshwater discharge on the subpolar North Atlantic (SPNA) remains unknown as there are uncertainties associated with the time scales of the Greenland freshwater anomaly (GFWA) in the SPNA. Results from numerical simulations tracking GFWA and an analytical approach are employed to estimate the response time, suggesting that a decadal time scale (13 years) is required for the SPNA to adjust for increasing GFWA. Analytical solutions obtained for a long-lasting increase of freshwater discharge show a non-steady-state response of the SPNA with increasing content of the GFWA. In contrast, solutions for a short-lived pulse of freshwater demonstrate different responses of the SPNA with a rapid increase of freshwater in the domain followed by an exponential decay after the pulse has passed. The derived theoretical relation between time scales shows that residence time scales are time dependent for a non-steady-state case and asymptote the response time scale with time. The residence time of the GFWA deduced from Lagrangian experiments is close to and smaller than the response time, in agreement with the theory. The Lagrangian analysis shows dependence of the residence time on the entrance route of the GFWA and on the depth. The fraction of the GFWA exported through Davis Strait has limited impact on the interior basins, whereas the fraction entering the SPNA from the southwest Greenland shelf spreads into the interior regions. In both cases, the residence time of the GFWA increases with depth demonstrating long persistence of the freshwater anomaly in the subsurface layers.D. S. Dukhovskoy and E. P. Chassignet were funded by the DOE (Award DE-SC0014378) and HYCOM NOPP (Award N00014-19-1-2674). The HYCOM-CICE simulations were supported by a grant of computer time from the DoD High-Performance Computing Modernization Program at NRL SSC. G. Platov was funded by the RSF N19-17-00154. P. G. Myers was funded by an NSERC Discovery Grant (Grant RGPIN 04357). A. Proshutinsky was funded by FAMOS project (NSF Grant NSF 14-584)

Arctic decadal variability : an auto-oscillatory system of heat and fresh water exchange

Author Posting. © American Geophysical Union, 2004. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 31 (2004): L03302, doi:10.1029/2003GL019023.This paper presents a mechanism of decadal variability in the Artic Ocean–GIN Sea (Greenland, Iceland and Norwegian Seas) atmosphere-ice-ocean system. We hypothesize that Arctic variability is regulated by heat and freshwater exchange between the Arctic Ocean and the GIN Sea. The interaction between basins is weak during anticyclonic circulation regimes (low AO/NAO) and strong during cyclonic circulation regimes (high AO/NAO). Regime shifts are controlled by the system itself through oceanic and atmospheric gradients (dynamic height and surface air temperature) that increase during the anticyclonic regime and decrease during the cyclonic regime. This conceptual mechanism for Arctic decadal variability has been reproduced in a model experiment. Both model results and observational data support the suggested mechanism.This research has been supported by the National Science Foundation and by the International Arctic Research Center, University of Alaska Fairbanks, under auspices of the United States National Science Foundation and from the Alaska Sea Grant through the Center for Global Change, University of Alaska Fairbanks

Assessment of numerical simulations of deep circulation and variability in the Gulf of Mexico using recent observations

Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 50(4), (2020): 1045-1064, doi:10.1175/JPO-D-19-0137.1.Three simulations of the circulation in the Gulf of Mexico (the “Gulf”) using different numerical general circulation models are compared with results of recent large-scale observational campaigns conducted throughout the deep (>1500 m) Gulf. Analyses of these observations have provided new understanding of large-scale mean circulation features and variability throughout the deep Gulf. Important features include cyclonic flow along the continental slope, deep cyclonic circulation in the western Gulf, a counterrotating pair of cells under the Loop Current region, and a cyclonic cell to the south of this pair. These dominant circulation features are represented in each of the ocean model simulations, although with some obvious differences. A striking difference between all the models and the observations is that the simulated deep eddy kinetic energy under the Loop Current region is generally less than one-half of that computed from observations. A multidecadal integration of one of these numerical simulations is used to evaluate the uncertainty of estimates of velocity statistics in the deep Gulf computed from limited-length (4 years) observational or model records. This analysis shows that the main deep circulation features identified from the observational studies appear to be robust and are not substantially impacted by variability on time scales longer than the observational records. Differences in strengths and structures of the circulation features are identified, however, and quantified through standard error analysis of the statistical estimates using the model solutions.This work was supported by the Gulf Research Program of the National Academy of Sciences under Awards 2000006422 and 2000009966. The content is solely the responsibility of the authors and does not necessarily represent the official views of the Gulf Research Program or the National Academy of Sciences. The authors acknowledge the GLORYS project for providing the ocean reanalysis data used in the ROMS simulation. GLORYS is jointly conducted by MERCATOR OCEAN, CORIOLIS, and CNRS/INSU. Installation, recovery, data acquisition, and processing of the CANEK group current-meter moorings were possible because of CICESE-PetróleosMexicanos Grant PEP-CICESE 428229851 and the dedicated work of the crew of the B/O Justo Sierra and scientists of the CANEK group. The authors thank Dr. Aljaz Maslo, CICESE, for assistance with analysis of model data. The Bureau of Ocean Energy Management (BOEM), U.S. Dept. of the Interior, provided funding for the Lagrangian Study of the Deep Circulation in the Gulf of Mexico and the Observations and Dynamics of the Loop Current study. HYCOM simulation data are available from the HYCOM data server (https://www.hycom.org/data/goml0pt04/expt-02pt2), MITgcm data are available from the ECCO data server (http://ecco.ucsd.edu/gom_results2.html), and the ROMS simulation data are available from GRIIDC (NA.x837.000:0001)

Impact of upstream variability on the Loop Current dynamics in numerical simulations of the Gulf of Mexico

The Loop Current (LC), which is the main mesoscale dynamic feature of the Gulf of Mexico (GoM), has a major impact on the circulation and its variability in the interior Gulf. The LC is a highly variable and dynamic feature. It changes shape from a short jet connecting the two openings of the GoM in an almost straight line ("retracted phase") to a long loop invading most of the eastern part of the GoM ("extended phase"). When it is in the extended phase, it can shed large anticyclonic eddies, called Loop Current Eddies, which then migrate to the western GoM. In this study, the processes controlling the LC dynamics are investigated using two multi-decadal simulations of the Gulf of Mexico HYbrid Coordinate Ocean Model differing in their open boundary conditions (BCs) and altimetry-derived gridded fields. The LC in the simulation with BCs derived from monthly climatology state variables frequently remains in its retracted phase significantly longer than observed. In contrast, the duration of the retracted phase is notably shorter in the simulation in which the BCs have realistic daily variability. By examining the flow properties through the Yucatan Channel from which the LC originates, we find that increased intensity of this current and a westward shift of the mean core is associated with the LC transitions from the retracted to the extended phase. This transition is accompanied by an increase of both cyclonicity of the flow in the west and anticyclonicity in the east of the core of this jet. Moreover, the number of anticyclonic eddies entering in the GoM through the Yucatan Channel is significantly higher when the LC extends in the GoM. Consequently, this study demonstrates the importance of realistic flow variability at the lateral boundaries for accurate simulation of the LC system in a model, and reveals characteristics of the upstream flow associated with different LC behavior that can potentially aid in forecasting the LC system

Measurement Characteristics of Near-Surface Currents from Ultra-Thin Drifters, Drogued Drifters, and HF Radar

Concurrent measurements by satellite tracked drifters of different hull and drogue configurations and coastal high-frequency radar reveal substantial differences in estimates of the near-surface velocity. These measurements are important for understanding and predicting material transport on the ocean surface as well as the vertical structure of the near-surface currents. These near-surface current observations were obtained during a field experiment in the northern Gulf of Mexico intended to test a new ultra-thin drifter design. During the experiment, thirty small cylindrical drifters with 5 cm height, twenty-eight similar drifters with 10 cm hull height, and fourteen drifters with 91 cm tall drogues centered at 100 cm depth were deployed within the footprint of coastal High-Frequency (HF) radar. Comparison of collocated velocity measurements reveals systematic differences in surface velocity estimates obtained from the different measurement techniques, as well as provides information on properties of the drifter behavior and near-surface shear. Results show that the HF radar velocity estimates had magnitudes significantly lower than the 5 cm and 10 cm drifter velocity of approximately 45% and 35%, respectively. The HF radar velocity magnitudes were similar to the drogued drifter velocity. Analysis of wave directional spectra measurements reveals that surface Stokes drift accounts for much of the velocity difference between the drogued drifters and the thin surface drifters except during times of wave breaking

Analysis of the Beaufort Gyre freshwater content in 2003-2018

© The Author(s), 2019. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Proshutinsky, A., Krishfield, R., Toole, J. M., Timmermans, M-L., Williams, W. J., Zimmermann, S., Yamamoto-Kawai, M., Armitage, T. W. K., Dukhovskoy, D., Golubeva, E., Manucharyan, G. E., Platov, G., Watanabe, E., Kikuchi, T., Nishino, S., Itoh, M., Kang, S-H., Cho, K-H., Tateyama, K., & Zhao, J. Analysis of the Beaufort Gyre freshwater content in 2003-2018. Journal of Geophysical Research-Oceans, 124(12), (2019): 9658-9689, doi:10.1029/2019JC015281.Hydrographic data collected from research cruises, bottom‐anchored moorings, drifting Ice‐Tethered Profilers, and satellite altimetry in the Beaufort Gyre region of the Arctic Ocean document an increase of more than 6,400 km3 of liquid freshwater content from 2003 to 2018: a 40% growth relative to the climatology of the 1970s. This fresh water accumulation is shown to result from persistent anticyclonic atmospheric wind forcing (1997–2018) accompanied by sea ice melt, a wind‐forced redirection of Mackenzie River discharge from predominantly eastward to westward flow, and a contribution of low salinity waters of Pacific Ocean origin via Bering Strait. Despite significant uncertainties in the different observations, this study has demonstrated the synergistic value of having multiple diverse datasets to obtain a more comprehensive understanding of Beaufort Gyre freshwater content variability. For example, Beaufort Gyre Observational System (BGOS) surveys clearly show the interannual increase in freshwater content, but without satellite or Ice‐Tethered Profiler measurements, it is not possible to resolve the seasonal cycle of freshwater content, which in fact is larger than the year‐to‐year variability, or the more subtle interannual variations.National Science Foundation. Grant Numbers: PLR‐1302884,OPP‐1719280, and OPP‐1845877, PLR‐1303644 and OPP‐1756100, OPP‐1756100, PLR‐1303644, OPP‐1845877, OPP‐1719280, PLR‐1302884 Key Program of National Natural Science Foundation of China. Grant Number: 41330960 Global Change Research Program of China. Grant Number: 2015CB953900 Ministry of Education, Korea Japan Aerospace Exploration Agency (JAXA) /Earth Observation Research Center (EORC) Ministry of Education, Culture, Sports, Science and Technology of Japan (MEXT) Stanback Postdoctoral Fellowship Russian Foundation for Basic Research. Grant Number: 17‐05‐00382 Presidium of Russian Academy of Sciences HYCOM NOPP. Grant Number: N00014‐15‐1‐2594 DOE. Grant Number: DE‐SC0014378 National Aeronautics and Space Administration Tokyo University of Marine Science and Technology Department of Fisheries and Oceans Canada Woods Hole Oceanographic Institutio