Scaling Properties of Arctic Sea Ice Deformation in a High-Resolution Viscous-Plastic Sea Ice Model and in Satellite Observations

Sea ice models with the traditional viscous-plastic (VP) rheology and very small horizontal grid spacing can resolve leads and deformation rates localized along Linear Kinematic Features (LKF). In a 1 km pan-Arctic sea ice-ocean simulation, the small-scale sea ice deformations are evaluated with a scaling analysis in relation to satellite observations of the Envisat Geophysical Processor System (EGPS) in the Central Arctic. A new coupled scaling analysis for data on Eulerian grids is used to determine the spatial and temporal scaling and the coupling between temporal and spatial scales. The spatial scaling of the modeled sea ice deformation implies multifractality. It is also coupled to temporal scales and varies realistically by region and season. The agreement of the spatial scaling with satellite observations challenges previous results with VP models at coarser resolution, which did not reproduce the observed scaling. The temporal scaling analysis shows that the VP model, as configured in this 1 km simulation, does not fully resolve the intermittency of sea ice deformation that is observed in satellite data

Adaptation of a Fast Optimal Interpolation Algorithm to the Mapping of Oceangraphic Data

A fast, recently developed, multiscale optimal interpolation algorithm has been adapted to the mapping of hydrographic and other oceanographic data. This algorithm produces solution and error estimates which are consistent with those obtained from exact least squares methods, but at a small fraction of the computational cost. Problems whose solution would be completely impractical using exact least squares, that is, problems with tens or hundreds of thousands of measurements and estimation grid points, can easily be solved on a small workstation using the multiscale algorithm. In contrast to methods previously proposed for solving large least squares problems, our approach provides estimation error statistics while permitting long-range correlations, using all measurements, and permitting arbitrary measurement locations. The multiscale algorithm itself, published elsewhere, is not the focus of this paper. However, the algorithm requires statistical models having a very particular multiscale structure; it is the development of a class of multiscale statistical models, appropriate for oceanographic mapping problems, with which we concern ourselves in this paper. The approach is illustrated by mapping temperature in the northeastern Pacific. The number of hydrographic stations is kept deliberately small to show that multiscale and exact least squares results are comparable. A portion of the data were not used in the analysis; these data serve to test the multiscale estimates. A major advantage of the present approach is the ability to repeat the estimation procedure a large number of times for sensitivity studies, parameter estimation, and model testing. We have made available by anonymous Ftp a set of MATLAB-callable routines which implement the multiscale algorithm and the statistical models developed in this paper

Scaling properties of Arctic sea ice deformation in high-resolution viscous-plastic sea-ice models

Many climate models use a rheology of the viscous-plastic type to simulate sea ice dynamics. With this rheology, large scale velocity and thickness fields can be realistically simulated, but the representation of small scale deformation rates and Linear Kinematic Features (LKF) is thought to be inadequate. However, at high resolution (< 5 km) the rheology starts to produce lines of localised deformation rates. In this study we use results from a 1-km Pan-Arctic model to investigate the influence of these deformation features on the scaling properties of sea ice deformation. For evaluation the EGPS satellite data set of small-scale sea ice kinematics for the Central Arctic (successor of RGPS) is used. The modelled sea ice deformation shows multi-fractal spatial scaling and, in this sense, agrees with the satellite data. In addition, the temporal coupling of the spatial scaling is reproduced as well. Furthermore, we examine the regional and seasonal variations of spatial scaling properties and its dependence on the ice condition, i.e. sea ice concentration and thickness, which are in agreement with previous RGPS studies

Numerical investigations of seasonal and interannual variability of North Pacific Subtropical Mode Water and its implications for Pacific climate variability

Author Posting. © American Meteorological Society, 2011. 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 24 (2011): 2648–2665, doi:10.1175/2010JCLI3435.1.North Pacific Subtropical Mode Water (NPSTMW) is an essential feature of the North Pacific subtropical gyre imparting significant influence on regional SST evolution on seasonal and longer time scales and, as such, is an important component of basin-scale North Pacific climate variability. This study examines the seasonal-to-interannual variability of NPSTMW, the physical processes responsible for this variability, and the connections between NPSTMW and basin-scale climate signals using an eddy-permitting 1979–2006 ocean simulation made available by the Estimating the Circulation and Climate of the Ocean, Phase II (ECCO2). The monthly mean seasonal cycle of NPSTMW in the simulation exhibits three distinct phases: (i) formation during November–March, (ii) isolation during March–June, and (iii) dissipation during June–November—each corresponding to significant changes in upper-ocean structure. An interannual signal is also evident in NPSTMW volume and other characteristic properties with volume minima occurring in 1979, 1988, and 1999. This volume variability is correlated with the Pacific decadal oscillation (PDO) with zero time lag. Further analyses demonstrate the connection of NPSTMW to the basin-scale ocean circulation. With this, modulations of upper-ocean structure driven by the varying strength and position of the westerlies as well as the regional air–sea heat flux pattern are seen to contribute to the variability of NPSTMW volume on interannual time scales.Support for this research was provided by the Partnership for Advancing Interdisciplinary Modeling (PARADIGM), a National Ocean Partnership Program and by a NASA Modeling, Analysis, and the Prediction (MAP) project called Estimating the Circulation and Climate of the Ocean, Phase II (ECCO2)

Altimetry-Based Diagnosis of Deep-Reaching Sub-Mesoscale Ocean Fronts

Recent studies demonstrate that energetic sub-mesoscale fronts (10–50 km width) extend in the ocean interior, driving large vertical velocities and associated fluxes. However, diagnosing the dynamics of these deep-reaching fronts from in situ observations remains challenging because of the lack of information on the 3-D structure of the horizontal velocity. Here, a realistic numerical simulation in the Antarctic Circumpolar Current (ACC) is used to study the dynamics of submesocale fronts in relation to velocity gradients, responsible for the formation of these fronts. Results highlight that the stirring properties of the flow at depth, which are related to the velocity gradients, can be inferred from finite-size Lyapunov exponent (FSLE) at the surface. Satellite altimetry observations of FSLE and velocity gradients are then used in combination with recent in situ observations collected by an elephant seal in the ACC to reconstruct frontal dynamics and their associated vertical velocities down to 500 m. The approach proposed here is well suited for the analysis of sub-mesoscale-resolving datasets and the design of future sub-mesoscale field campaigns

Characterizing the transition from balanced to unbalanced motions in the Southern California Current

Author Posting. © American Geophysical Union, 2019. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research-Oceans 124(3), (2019): 2088-2109, doi:10.1029/2018JC014583.As observations and models improve their resolution of oceanic motions at ever finer horizontal scales, interest has grown in characterizing the transition from the geostrophically balanced flows that dominate at large‐scale to submesoscale turbulence and waves that dominate at small scales. In this study we examine the mesoscale‐to‐submesoscale (100 to 10 km) transition in an eastern boundary current, the southern California Current System (CCS), using repeated acoustic Doppler current profiler transects, sea surface height from high‐resolution nadir altimetry and output from a (1/48)° global model simulation. In the CCS, the submesoscale is as energetic as in western boundary current regions, but the mesoscale is much weaker, and as a result the transition lacks the change in kinetic energy (KE) spectral slope observed for western boundary currents. Helmholtz and vortex‐wave decompositions of the KE spectra are used to identify balanced and unbalanced contributions. At horizontal scales greater than 70 km, we find that observed KE is dominated by balanced geostrophic motions. At scales from 40 to 10 km, unbalanced contributions such as inertia‐gravity waves contribute as much as balanced motions. The model KE transition occurs at longer scales, around 125 km. The altimeter spectra are consistent with acoustic Doppler current profiler/model spectra at scales longer than 70/125 km, respectively. Observed seasonality is weak. Taken together, our results suggest that geostrophic velocities can be diagnosed from sea surface height on scales larger than about 70 km in the southern CCS.This research was funded by NASA (NNX13AE44G, NNX13AE85G, NNX16AH67G, NNX16AO5OH, and NNX17AH53G). We thank Sung Yong Kim for providing the high‐frequency radar spectral estimates and the two anonymous reviewers for providing useful comments and suggestions that greatly improved the manuscript. High‐frequency ALES data for Jason‐1 and Jason‐2 altimeters are available upon request (https://openadb.dgfi.tum.de/en/contact/ALES). Both AltiKa and Sentinel‐3 altimeter products were produced and distributed by the Copernicus Marine and Environment Monitoring Service (CMEMS; http://www.marine.copernicus.eu). D. M. worked on the modeling component of this study at the Jet Propulsion Laboratory (JPL), California Institute of Technology, under a contract with the National Aeronautics and Space Administration (NASA). High‐end computing resources were provided by the NASA Advanced Supercomputing (NAS) Division of the Ames Research Center. The LLC output can be obtained from the ECCO project (ftp://ecco.jpl.nasa.gov/ECCO2/LLC4320/). The ADCP data are available at the Joint Archive for Shipboard ADCP data (JASADCP; http://ilikai.soest.hawaii.edu/sadcp).2019-08-2

Improved representation of river runoff in Estimating the Circulation and Climate of the Ocean Version 4 (ECCOv4) simulations: Implementation, evaluation, and impacts to coastal plume regions

In this study, we improve the representation of global river runoff in the Estimating the Circulation and Climate of the Ocean Version 4 (ECCOv4) framework, allowing for a more realistic treatment of coastal plume dynamics. We use a suite of experiments to explore the sensitivity of coastal plume regions to runoff forcing, model grid resolution, and grid type. The results show that simulated sea surface salinity (SSS) is reduced as the model grid resolution increases. Compared to Soil Moisture Active Passive (SMAP) observations, simulated SSS is closest to SMAP when using daily, point-source runoff (DPR) and the intermediate-resolution LLC270 grid. The Willmott skill score, which quantifies agreement between models and SMAP, yields up to 0.92 for large rivers such as the Amazon. There was no major difference in SSS for tropical and temperate coastal rivers when the model grid type was changed from the ECCO v4 latitude-longitude-polar-cap grid to the ECCO2 cube-sphere grid. We also found that using DPR forcing and increasing model resolution from the coarse-resolution LLC90 grid to the intermediate-resolution LLC270 grid elevated the river plume area, volume, stabilized the stratification and shoal the mixed layer depth (MLD). Additionally, we find that the impacts of increasing model resolution from the intermediate-resolution LLC270 grid to the high-resolution LLC540 grid are regionally dependent. The Mississippi River Plume is more sensitive than other regions, possibly because the wider and shallower Texas-Louisiana shelf drives a stronger baroclinic effect, as well as relatively weak sub-grid vertical mixing and adjustment in this region. Since rivers deliver large amounts of freshwater and anthropogenic materials to coastal regions, improving the representation of river runoff in global, high-resolution models will advance studies of coastal hypoxia, carbon cycling, and regional weather and climate and will ultimately help to predict land-ocean-atmospheric feedbacks seamlessly in the next generation of Earth system models