Nicole Jeffery

The Southern Ocean and Its Climate in CCSM4

AbstractThe new Community Climate System Model, version 4 (CCSM4), provides a powerful tool to understand and predict the earth’s climate system. Several aspects of the Southern Ocean in the CCSM4 are explored, including the surface climatology and interannual variability, simulation of key climate water masses (Antarctic Bottom Water, Subantarctic Mode Water, and Antarctic Intermediate Water), the transport and structure of the Antarctic Circumpolar Current, and interbasin exchange via the Agulhas and Tasman leakages and at the Brazil–Malvinas Confluence. It is found that the CCSM4 has varying degrees of accuracy in the simulation of the climate of the Southern Ocean when compared with observations. This study has identified aspects of the model that warrant further analysis that will result in a more comprehensive understanding of ocean–atmosphere–ice dynamics and interactions that control the earth’s climate and its variability.

Pan-Arctic simulation of coupled nutrient-sulfur cycling due to sea ice biology : Preliminary results

A dynamic model is constructed for interactive silicon, nitrogen, sulfur processing in and below Arctic sea ice, by ecosystems residing in the lower few centimeters of the distributed pack. A biogeochemically active bottom layer supporting sources/sinks for the pennate diatoms is appended to thickness categories of a global sea ice code. Nutrients transfer from the ocean mixed layer to drive algal growth, while sulfur metabolites are reinjected from the ice interface. Freeze, flux, flush and melt processes are linked to multielement geocycling for the entire high-latitude regime. Major element kinetics are optimized initially to reproduce chlorophyll observations, which extend across the seasons. Principal influences on biomass are solute exchange velocity at the solid interface, optical averaging in active ice and cell retention against ablation. The sulfur mechanism encompasses open water features such as accumulation of particulate dimethyl sulfoniopropionate, grazing and other disruptive releases, plus bacterial/enzymatic conversion to volatile dimethyl sulfide. For baseline settings, the mixed layer trace gas distribution matches sparging measurements where they are available. However, concentrations rise to well over 10 nM in remote, unsampled locations. Peak contributions are supported by ice grazing, mortality and fractional melting. The model bottom layer adds substantially to a ring maximum of reduced sulfur chemistry that may be dominant across the marginal Arctic environment. Sensitivity tests on this scenario include variation of cell sulfur composition and remineralization, routings/chemical time scales, and the physical dimension of water layers. An alternate possibility that peripheral additions are small cannot be excluded from the outcomes. It is concluded that seagoing dimethyl sulfide data are far too sparse at the present time to distinguish sulfur-ice production levels. Citation: Elliott, S., C. Deal, G. Humphries, E. Hunke, N. Jeffery, M. Jin, M. Levasseur, and J. Stefels (2012), Pan-Arctic simulation of coupled nutrient-sulfur cycling due to sea ice biology: Preliminary results, J. Geophys. Res., 117, G01016, doi:10.1029/2011JG001649.

The Effect of Tilted Rotation on Shear Instabilities at Low Stratifications

Abstract A linear stability analysis of the inviscid stratified Boussinesq equations is presented given a steady zonal flow with constant vertical shear in a tilted f plane. Full nonhydrostatic terms are included: 1) acceleration of vertical velocity and 2) Coriolis force terms arising from the meridional component of earth’s rotation vector. Calculations of growth rates, critical wavenumbers, and dominance regimes for baroclinic and symmetric instabilities are compared with results from the traditional nonhydrostatic equations, which include a strictly vertical rotation vector, as well as results from the hydrostatic equations. The authors find that for positive zonal z shear, tilted rotation enhances the dominance regime of symmetric instabilities at the expense of baroclinic instabilities and maintains symmetric instabilities at larger scales than previously indicated. Furthermore, in contrast to former studies, it is determined that hydrostatic growth rates for both instabilities are not maximal. Rath...

Sea ice thermohaline dynamics and biogeochemistry in the Arctic Ocean: Empirical and model results

Large changes in the sea ice regime of the Arctic Ocean have occurred over the last decades justifying the development of models to forecast sea ice physics and biogeochemistry. The main goal of this study is to evaluate the performance of the Los Alamos Sea Ice Model (CICE) to simulate physical and biogeochemical properties at time scales of a few weeks and to use the model to analyze ice algal bloom dynamics in different types of ice. Ocean and atmospheric forcing data and observations of the evolution of the sea ice properties collected from 18 April to 4 June 2015, during the Norwegian young sea ICE expedition, were used to test the CICE model. Our results show the following: (i) model performance is reasonable for sea ice thickness and bulk salinity; good for vertically resolved temperature, vertically averaged Chl a concentrations, and standing stocks; and poor for vertically resolved Chl a concentrations. (ii) Improving current knowledge about nutrient exchanges, ice algal recruitment, and motion is critical to improve sea ice biogeochemical modeling. (iii) Ice algae may bloom despite some degree of basal melting. (iv) Ice algal motility driven by gradients in limiting factors is a plausible mechanism to explain their vertical distribution. (v) Different ice algal bloom and net primary production (NPP) patterns were identified in the ice types studied, suggesting that ice algal maximal growth rates will increase, while sea ice vertically integrated NPP and biomass will decrease as a result of the predictable increase in the area covered by refrozen leads in the Arctic Ocean.

Modeling the winter‐spring transition of first‐year ice in the western Weddell Sea

A new halodynamic scheme is coupled with the Los Alamos sea ice model to simulate western Weddell Sea ice during the winter-spring transition. One-dimensional temperature and salinity profiles are consistent with the warming and melt stages exhibited in first-year ice cores from the 2004 Ice Station POLarstern (ISPOL) expedition. Results are highly sensitive to snowfall. Simulations which use reanalysis precipitation data do not retain a snow cover beyond mid-December, and the warming transition occurs too rapidly. Model performance is greatly improved by prescribing a snowfall rate based on reported snow thicknesses. During ice growth prior to ISPOL, simulations indicate a period of thick snow and upper ice salinity enrichment. Gravity drainage model parameters impact the simulation immediately, while effects from the flushing parameter (snow porosity at the ice top) appear as the freeboard becomes negative. Simulations using a snow porosity of 0.3, consistent with that of wet snow, agree with salinity observations. The model does not include lateral sources of sea-water flooding, but vertical transport processes account for the high upper-ice salinities observed in ice cores at the start of the expedition. As the ice warms, a fresh upper-ice layer forms, and the high salinity layer migrates downward. This pattern is consistent with the early spring development stages of high-porosity layers observed in Antarctic sea ice that are associated with rich biological production. Future extensions of the model may be valuable in Antarctic ice-biogeochemical applications.

Local Atmospheric Response to an Open-Ocean Polynya in a High-Resolution Climate Model

AbstractIn this paper the atmospheric response to an open-ocean polynya in the Southern Ocean is studied by analyzing the results from an atmospheric and oceanic synoptic-scale resolving Community Earth System Model (CESM) simulation. While coarser-resolution versions of CESM generally do not produce open-ocean polynyas in the Southern Ocean, they do emerge and disappear on interannual time scales in the synoptic-scale simulation. This provides an ideal opportunity to study the polynya’s impact on the overlying and surrounding atmosphere. This has been pursued here by investigating the seasonal cycle of differences of surface and air-column variables between polynya and nonpolynya years. The results indicate significant local impacts on turbulent heat fluxes, precipitation, cloud characteristics, and radiative fluxes. In particular, it is found that clouds over polynyas are optically thicker and higher than clouds over sea ice during nonpolynya years. Although the lower albedo of polynyas significantly in...

Effects of Model Resolution and Ocean Mixing on Forced Ice‐Ocean Physical and Biogeochemical Simulations Using Global and Regional System Models

The current coarse-resolution global Community Earth System Model (CESM) can reproduce major and large-scale patterns but is still missing some key biogeochemical features in the Arctic Ocean, e.g., low surface nutrients in the Canada Basin. We incorporated the CESM Version 1 ocean biogeochemical code into the Regional Arctic System Model (RASM) and coupled it with a sea-ice algal module to investigate model limitations. Four ice-ocean hindcast cases are compared with various observations: two in a global 18 (40 60 km in the Arctic) grid: G1deg and G1deg-OLD with/without new sea-ice processes incorporated; two on RASM’s 1/128 ( 9 km) grid R9km and R9km-NB with/without a subgrid scale brine rejection parameterization which improves ocean vertical mixing under sea ice. Higher-resolution and new sea-ice processes contributed to lower model errors in sea-ice extent, ice thickness, and ice algae. In the Bering Sea shelf, only higher resolution contributed to lower model errors in salinity, nitrate (NO3), and chlorophyll-a (Chl-a). In the Arctic Basin, model errors in mixed layer depth (MLD) were reduced 36% by brine rejection parameterization, 20% by new sea-ice processes, and 6% by higher resolution. The NO3 concentration biases were caused by both MLD bias and coarse resolution, because of excessive horizontal mixing of high NO3 from the Chukchi Sea into the Canada Basin in coarse resolution models. R9km showed improvements over G1deg on NO3, but not on Chl-a, likely due to light limitation under snow and ice cover in the Arctic Basin.

Validation of sea ice models using an uncertainty-based distance metric for multiple model variables: NEW METRIC FOR SEA ICE MODEL VALIDATION

We implement a variance-based distance metric (Dn) to objectively assess skill of sea ice models when multiple output variables or uncertainties in both model predictions and observations need to be considered. The metric compares observations and model data pairs on common spatial and temporal grids improving upon highly aggregated metrics (e.g., total sea ice extent or volume) by capturing the spatial character of model skill. The Dn metric is a gamma-distributed statistic that is more general than the χ2 statistic commonly used to assess model fit, which requires the assumption that the model is unbiased and can only incorporate observational error in the analysis. The Dn statistic does not assume that the model is unbiased, and allows the incorporation of multiple observational data sets for the same variable and simultaneously for different variables, along with different types of variances that can characterize uncertainties in both observations and the model. This approach represents a step to establish a systematic framework for probabilistic validation of sea ice models. The methodology is also useful for model tuning by using the Dn metric as a cost function and incorporating model parametric uncertainty as part of a scheme to optimize model functionality. We apply this approach to evaluate different configurations of the standalone Los Alamos sea ice model (CICE) encompassing the parametric uncertainty in the model, and to find new sets of model configurations that produce better agreement than previous configurations between model and observational estimates of sea ice concentration and thickness.

Biogeochemistry of CICE: the Los Alamos Sea Ice Model Documentation and Software User's Manual zbgc_colpkg modifications to Version 5

This is the software Users Manual for the zbgc_colpkg modifications in version 5. It is used to model the effect of aerosols on ice deposits.