Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sea ice simulator”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A Hierarchy of Global Ocean Models Coupled to CESM1

Abstract We develop a hierarchy of simplified ocean models for coupled ocean, atmosphere, and sea ice climate simulations using the Community Earth System Model version 1 (CESM1). The hierarchy has four members: a slab ocean model, a mixed‐layer model (MLM) with entrainment and detrainment, an Ekman MLM, and an ocean general circulation model (OGCM). Flux corrections of heat and salt are applied to the simplified models ensuring that all hierarchy members have the same climatology. We diagnose the needed flux corrections from auxiliary simulations in which we restore the temperature and salinity to the daily climatology obtained from a target CESM1 simulation. The resulting three‐dimensional corrections contain the interannual variability fluxes that maintain the correct vertical gradients of temperature and salinity in the tropics. We find that the inclusion of mixed‐layer entrainment and Ekman flow produces sea surface temperature and surface air temperature fields whose means and variances are progressively more similar to those produced by the target CESM1 simulation. We illustrate the application of the hierarchy to the problem of understanding the response of the climate system to the loss of Arctic sea ice. We find that the shifts in the positions of the mid‐latitude westerly jet and of the Inter‐tropical Convergence Zone (ITCZ) in response to sea‐ice loss depend critically on upper ocean processes. Specifically, heat uptake associated with the mixed‐layer entrainment influences the shift in the westerly jet and ITCZ. Moreover, the shift of ITCZ is sensitive to the form of Ekman flow parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A hierarchy of global ocean models coupled to CESM1

We develop a hierarchy of simplified ocean models for coupled ocean, atmosphere, and sea ice climate simulations using the Community Earth System Model version 1 (CESM1). The hierarchy has four members: a slab ocean model, a mixed-layer model (MLM) with entrainment and detrainment, an Ekman MLM, and an ocean general circulation model (OGCM). Flux corrections of heat and salt are applied to the simplified models ensuring that all hierarchy members have the same climatology. We diagnose the needed flux corrections from auxiliary simulations in which we restore the temperature and salinity to the daily climatology obtained from a target CESM1 simulation. The resulting three-dimensional corrections contain the interannual variability fluxes that maintain the correct vertical gradients of temperature and salinity in the tropics. We find that the inclusion of mixed-layer entrainment and Ekman flow produces sea surface temperature and surface air temperature fields whose means and variances are progressively more similar to those produced by the target CESM1 simulation. We illustrate the application of the hierarchy to the problem of understanding the response of the climate system to the loss of Arctic sea ice. We find that the shifts in the positions of the mid-latitude westerly jet and of the Inter-tropical Convergence Zone (ITCZ) in response to sea-ice loss depend critically on upper ocean processes. Specifically, heat uptake associated with the mixed-layer entrainment influences the shift in the westerly jet and ITCZ. Moreover, the shift of ITCZ is sensitive to the form of Ekman flow parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A modified thermodynamic sea ice model and its application

This paper presents a modified thermodynamic sea ice model based on the Winton’s three-layer model framework. Several improvements focused on the vertical thermodynamics have been made in the model. Results from a series of one-dimensional experiments show that equilibrium ice thickness in our modified model is increased by about 45 cm when compared with the original Winton’s model. Sensitivity tests indicate that all modifications mentioned above contribute to this change of ice thickness, among which the increase of ice layer has the most significant effect followed by the sea ice conductivity parameterization. The modified thermodynamic sea ice model is then coupled into the sea ice component SIS in GFDL MOM4 to exert long-term integration experiments forced by CORE2.0 data. Results indicate that the modified model shows significant improvement in simulating the Arctic sea ice, including an increase in both the sea ice thickness and volume over the whole Arctic region. This confirms the above founding from 1-D simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Large Eddy Simulations of the Atmospheric Boundary Layer Over Satellite‐Sensed Sea Ice Maps

Surface heterogeneity in the marginal ice zone (MIZ) causes multiscale secondary atmospheric circulations that are challenging to model or observe. The absence or inadequate representation of these circulations in ocean‐atmosphere exchange schemes in climate models is partially responsible for the underestimation of Arctic sea ice loss. Observationally, such circulations obfuscate the interpretation of polar atmospheric chemistry measurements, among others. To address this open challenge, large‐eddy simulations are conducted over real‐world satellite‐sensed sea ice maps with an ice‐sea temperature contrast, as well as over idealizations of these maps that alter the ice pattern but conserve its fraction, showing that the ice fraction of a surface is not sufficient to predict the surface heat flux. In a second suite of simulations, three other heterogeneity metrics (representing the surface fragmentation, patch edge tortuosity, and patch size variability) are introduced to complement the ice fraction in describing the surface. Simulations varying these parameters suggest that they also significantly impact surface‐air interactions. A roughness contrast is then added to the surface temperature contrast, indicating that the contribution of roughness changes to the resulting atmospheric circulations is less pronounced than that of thermal heterogeneity. Based on these results, we illustrate, using a multi‐linear regression on these map features, that MIZ surface parameterizations in Earth Systems Models can be improved if they account for these various characteristics of the sea‐ice patterns.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Simple Hybrid Sea Ice Nudging Method for Improving Control Over Partitioning of Sea Ice Concentration and Thickness

Abstract To assess the effect of ocean‐atmosphere coupling in the climate response to forced sea ice loss, the Polar Amplification Model Intercomparison Project protocol includes centennial coupled atmosphere‐ocean general circulation model simulations with imposed sea ice loss. The protocol, which specifies sea ice concentration and thickness distribution targets, does not prescribe a method for achieving them. Although different methods for imposing sea ice loss (or growth) in models have been documented, testing of the method‐dependence of the resulting climate responses has been limited. Achieving the targeted sea ice state has proven to be challenging using the so‐called ghost‐flux nudging method, which induces ice melt from below, as this method does not constrain the partitioning between thickness and concentration. We propose, describe and test a simple method that combines the advantages of direct sea ice nudging and ghost‐flux nudging. The hybrid nudging method better captures the partitioning between thickness and concentration while conserving total water content. We document how this novel sea ice constraining method reaches specific targets, enhances surface turbulent heat flux responses to sea ice loss, and induces tropospheric warming for both polar regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Separating the Influences of Low-Latitude Warming and Sea Ice Loss on Northern Hemisphere Climate Change

Abstract Analyzing a multimodel ensemble of coupled climate model simulations forced with Arctic sea ice loss using a two-parameter pattern-scaling technique to remove the cross-coupling between low- and high-latitude responses, the sensitivity to high-latitude sea ice loss is isolated and contrasted to the sensitivity to low-latitude warming. Despite some differences in experimental design, the Northern Hemisphere near-surface atmospheric sensitivity to sea ice loss is found to be robust across models in the cold season; however, a larger intermodel spread is found at the surface in boreal summer, and in the free tropospheric circulation. In contrast, the sensitivity to low-latitude warming is most robust in the free troposphere and in the warm season, with more intermodel spread in the surface ocean and surface heat flux over the Northern Hemisphere. The robust signals associated with sea ice loss include upward turbulent and longwave heat fluxes where sea ice is lost, warming and freshening of the Arctic Ocean, warming of the eastern North Pacific Ocean relative to the western North Pacific with upward turbulent heat fluxes in the Kuroshio Extension, and salinification of the shallow shelf seas of the Arctic Ocean alongside freshening in the subpolar North Atlantic Ocean. In contrast, the robust signals associated with low-latitude warming include intensified ocean warming and upward latent heat fluxes near the western boundary currents, freshening of the Pacific Ocean, salinification of the North Atlantic, and downward sensible and longwave fluxes over the ocean.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Arctic ice-ocean interactions in an 8-to-2 kilometer resolution global model

In the last decades, the Arctic climate has changed dramatically, with the loss of multiyear sea ice one of the clearest consequences. These changes have occurred on relatively rapid timescales, and both accurate short-term Arctic prediction (e.g., 10 days to three months) and climate projection of future Arctic scenarios present ongoing challenges. Here we describe a representation of the Arctic ocean and sea ice in a ultrahigh resolution simulation in which the horizontal grid mesh reduces from 8 km at the equator to 2 km at the poles (UH8to2) for the years 2017-2020. We find the simulation reproduces observed distributions of seasonal sea-ice thickness and concentration realistically, although concentration is biased low in the spring and summer and low biases in thickness are found in the central and eastern basins in the fall. Volume, fresh water, and heat transports through key passages are realistic, lying within observationally determined ranges. Climatological comparisons reveal that the UH8to2 Atlantic Water is shallower, warmer, and saltier than the World Ocean Atlas 2018 climatology for 2005-2017 in the eastern basin. Our analysis suggests that these biases, combined with a lack of stratification in the upper 100 m of the simulated ocean, contribute to the winter biases in modeled sea ice thickness. This relationship between biases in the sea ice and ocean points to a potential positive feedback within the model, illuminating challenges for long term model predictive power in a changing Arctic climate.

58 GEOSCIENCES↗
DOE OSTI2023

A Machine Learning Correction Model of the Winter Clear-Sky Temperature Bias over the Arctic Sea Ice in Atmospheric Reanalyses

Atmospheric reanalyses are widely used to estimate the past atmospheric near-surface state over sea ice. They provide boundary conditions for sea ice and ocean numerical simulations and relevant information for studying polar variability and anthropogenic climate change. Previous research revealed the existence of large near-surface temperature biases (mostly warm) over the Arctic sea ice in the current generation of atmospheric reanalyses, which is linked to a poor representation of the snow over the sea ice and the stably stratified boundary layer in the forecast models used to produce the reanalyses. These errors can compromise the employment of reanalysis products in support of polar research. Here, we train a fully connected neural network that learns from remote sensing infrared temperature observations to correct the existing generation of uncoupled atmospheric reanalyses (ERA5, JRA-55) based on a set of sea ice and atmospheric predictors, which are themselves reanalysis products. The advantages of the proposed correction scheme over previous calibration attempts are the consideration of the synoptic weather and cloud state, compatibility of the predictors with the mechanism responsible for the bias, and a self-emerging seasonality and multidecadal trend consistent with the declining sea ice state in the Arctic. The correction leads on average to a 27% temperature bias reduction for ERA5 and 7% for JRA-55 if compared to independent in situ observations from the MOSAiC campaign (respectively, 32% and 10% under clear-sky conditions). These improvements can be beneficial for forced sea ice and ocean simulations, which rely on reanalyses surface fields as boundary conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Role of Sea‐Surface Salinity in Simulating Historical Decadal Variations of Atlantic Meridional Overturning Circulation in a Coupled Climate Model

Abstract Extending climate prediction from seasonal to decadal timescales requires realistic initialization of not only upper ocean heat content but also the Atlantic meridional overturning circulation (AMOC). However, it remains a major challenge to realistically initialize AMOC in a coupled system while also maintaining a balanced atmosphere‐ocean initial state. This study demonstrates the feasibility of generating fully coupled historical states with realistic AMOC variability. Employing a forced ocean—sea‐ice (FOSI) model simulation as the “truth,” we show reproducibility of key features of historical AMOC decadal variability in a fully coupled model by restoring sea‐surface salinity, in addition to sea‐surface temperature restoring widely used in seasonal prediction. The atmospheric state of the restored coupled model solution is much closer to that of the free coupled simulation than to the observations used in FOSI, pointing to potential advantages of using this approach for initializing decadal predictions with reduced ocean initialization shock.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Impact of a New Sea Ice Thermodynamic Formulation in the CESM2 Sea Ice Component

Abstract The sea ice component of the Community Earth System Model version 2 (CESM2) contains new “mushy‐layer” physics that simulates prognostic salinity in the sea ice, with consequent modifications to sea ice thermodynamics and the treatment of melt ponds. The changes to the sea ice model and their influence on coupled model simulations are described here. Two simulations were performed to assess the changes in the vertical thermodynamics formulation with prognostic salinity compared to a constant salinity profile. Inclusion of the mushy layer thermodynamics of Turner et al. (2013, https://doi.org/10.1002/jgrc.20171 ) in a fully coupled Earth system model produces thicker and more extensive sea ice in the Arctic, with relatively unchanged sea ice in the Antarctic compared to simulations using a constant salinity profile. While this is consistent with the findings of uncoupled ice‐ocean model studies, the role of the frazil and congelation growth is more important in fully coupled simulations. Melt pond drainage is also an important contribution to simulated ice thickness differences as also found in the uncoupled simulations of Turner and Hunke (2015; https://doi.org/10.1002/2014JC010358 ). However, it is an interaction of the ponds and the snow fraction that impacts the surface albedo and hence the top melt. The changes in the thermodynamics and resulting ice state modify the ice‐ocean‐atmosphere fluxes with impacts on the atmosphere and ocean states, particularly temperature.

Bailey, David A.↗
DOE OSTI2022

An evaluation of the E3SMv1 Arctic ocean and sea-ice regionally refined model

The Energy Exascale Earth System Model (E3SM) is a state-of-the-science Earth system model (ESM) with the ability to focus horizontal resolution of its multiple components in specific areas. Regionally refined global ESMs are motivated by the need to explicitly resolve, rather than parameterize, relevant physics within the regions of refined resolution, while offering significant computational cost savings relative to the respective cost of configurations with high-resolution (HR) everywhere on the globe. In this paper, we document results from the first Arctic regionally refined E3SM configuration for the ocean and sea-ice components (E3SM-Arctic-OSI), while employing data-based atmosphere, land, and hydrology components. Our aim is an improved representation of the Arctic coupled ocean and sea-ice state, its variability and trends, and the exchanges of mass and property fluxes between the Arctic and the sub-Arctic. We find that E3SM-Arctic-OSI increases the realism of simulated Arctic ocean and sea-ice conditions compared to a similar low-resolution E3SM simulation without the Arctic regional refinement in ocean and sea-ice components (E3SM-LR-OSI). In particular, exchanges through the main Arctic gateways are greatly improved with respect to E3SM-LR-OSI. Other aspects, such as the Arctic freshwater content variability and sea-ice trends, are also satisfactorily simulated. Yet, other features, such as the upper-ocean stratification and the sea-ice thickness distribution, need further improvements, involving either more advanced parameterizations, model tuning, or additional grid refinements. Overall, E3SM-Arctic-OSI offers an improved representation of the Arctic system relative to E3SM-LR-OSI, at a fraction (15 %) of the computational cost of comparable global high-resolution configurations, while permitting exchanges with the lower-latitude oceans that cannot be directly accounted for in Arctic regional models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Mid-Holocene Antarctic sea-ice increase driven by marine ice sheet retreat

Abstract. Over recent decades Antarctic sea-ice extent has increased, alongside widespread ice shelf thinning and freshening of waters along the Antarctic margin. In contrast, Earth system models generally simulate a decrease in sea ice. Circulation of water masses beneath large-cavity ice shelves is not included in current Earth System models and may be a driver of this phenomena. We examine a Holocene sediment core off East Antarctica that records the Neoglacial transition, the last major baseline shift of Antarctic sea ice, and part of a late-Holocene global cooling trend. We provide a multi-proxy record of Holocene glacial meltwater input, sediment transport, and sea-ice variability. Our record, supported by high-resolution ocean modelling, shows that a rapid Antarctic sea-ice increase during the mid-Holocene (∼ 4.5 ka) occurred against a backdrop of increasing glacial meltwater input and gradual climate warming. We suggest that mid-Holocene ice shelf cavity expansion led to cooling of surface waters and sea-ice growth that slowed basal ice shelf melting. Incorporating this feedback mechanism into global climate models will be important for future projections of Antarctic changes.

58 GEOSCIENCES↗
DOE OSTI2022

On the variability of the Bering Sea Cold Pool and implications for the biophysical environment

The Bering Sea experiences a seasonal sea ice cover, which is important to the biophysical environment found there. A pool of cold bottom water (<2°C) is formed on the shelf each winter as a result of cooling and vertical mixing due to brine rejection during the predominately local sea ice growth. The extent and distribution of this Cold Pool (CP) is largely controlled by the winter extent of sea ice in the Bering Sea, which can vary considerably and recently has been much lower than average. The cold bottom water of the CP is important for food security because it delineates the boundary between arctic and subarctic demersal fish species. A northward retreat of the CP will likely be associated with migration of subarctic species toward the Chukchi Sea. We use the fully-coupled Regional Arctic System Model (RASM) to examine variability of the extent and distribution of the CP and its relation to change in the sea ice cover in the Bering Sea during the period 1980–2018. RASM results confirm the direct correlation between the extent of sea ice and the CP and show a smaller CP as a consequence of realistically simulated recent declines of the sea ice cover in the Bering Sea. In fact, the area of the CP was found to be only 31% of the long-term mean in July of 2018. In addition, we also find that a low ice year is followed by a later diatom bloom, while a heavy ice year is followed by an early diatom bloom. Finally, the RASM probabilistic intra-annual forecast capability is reviewed, based on 31-member ensembles for 2019–2021, for its potential use for prediction of the winter sea ice cover and the subsequent summer CP area in the Bering Sea.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Reducing uncertainty of polar to midlatitude linkages using DOE’s E3SM in a coordinated model-experiment setting

This project brought DOE’s climate modeling effort with the Energy Exascale Earth System Model into the Polar Amplification Model Intercomparison Project (PAMIP), which is part of the sixth and latest Coupled Model Intercomparison Project, CMIP6. PAMIP examines the causes and consequences of polar amplification, when external forcing results in a larger temperature increase in high latitudes than the global average, in a coordinated set of model experiments with a common modeling protocol. Our teams from UC Irvine and the University of Toronto have designed, carried out, analyzed, and disseminated PAMIP output from the Energy Exascale Earth System Model (E3SM) and the Community Earth System Model (CESM). PAMIP’s ongoing stream of significant new results have advanced progress in the community’s understanding and led to new outstanding research questions that have motivated further work. PAMIP has led to improved consensus on the atmospheric response to sea ice loss. The important finding is that for a similar sea ice anomaly forcing, the simulated atmospheric response in the troposphere is remarkably consistent among the 16 models’ runs analyzed. The zonal-mean tropospheric response consists of a very robust equatorward shift of the westerly flow in mid-latitudes. However, while the multi-model mean response is robust, it has a weak amplitude relative to internal variability. We identified a weakness in the models (including E3SM) in terms of their response to sea-ice forcing that is related to eddy forcing (or nonlinear dynamical effects) at mid- to high latitudes. In fact, E3SM is an outlier in terms of the models participating, and in that sense, it turned out to be a vital participant model. We found that reductions in energy transport due sea-ice loss and involving dry air only are compensated by increases in moist energy transport from warmer sea surface temperature in midlatitudes. This leads to a large spread in energy transport into the Arctic and is a potential source of spread in Arctic amplification. We identified an important role that climate modes, including tropical modes of variability (El Nino and the Southern Oscillation (ENSO); the Quasi-biennial Oscillation (QBO)) play in the response to sea-ice anomalies, including in ocean coupled experiments. Similarly, we identified and quantified the contribution of sea-ice thickness to the atmospheric response compared to the response to sea-ice extent only. We found that it is important to run large ensembles and even with an ensemble size of 100 simulations the response is largely influenced by internal variability. We demonstrated convincingly that Ural blocking, not sea-ice loss, provides the weakening of the stratospheric polar vortex in fall/early winter and a negative phase of the North Atlantic Oscillation that can last for up to two months. However, sea-ice anomalies can influence the background flow so that the response to Ural blocking is more persistent under low sea-ice conditions in the Barents/Kara Sea than high sea ice Atmospheric model hierarchies that progressively add individual processes have a long history in providing dynamical insight for modeling the atmosphere. Similarly, coupled model hierarchies that progressively add individual ocean processes can provide insights into the workings of the coupled climate system, however such hierarchies have not been available except for a non-dynamical slab ocean model. Because of the missing processes, surface flux corrections must be added to produce a target climate. In this project, we managed to overcome this problem and develop a globally coupled ocean model hierarchy in CESM that can turn on and off the processes of mixed-layer entrainment and Ekman flow. We used the hierarchy to study the impact of Arctic sea-ice loss on the climate system. We find that the effect of mixed-layer entrainment on ocean heat uptake influences the atmospheric circulation by shifting the latitudinal positions of the mid-latitude westerly jet and the Intertropical Convergence Zone (ITCZ). In quadrupled CO 2 experiments, we studied how air-sea coupling affects the response of tropical rainfall under global warming. In order to identify the importance of individual ocean processes, we used the hierarchy of ocean models to separate the effects of seasonal mixed-layer entrainment, wind-driven Ekman flows, and frictional flows. We showed that including Ekman and frictional flows allows our simulation to produce the Pacific Ocean's enhanced equatorial warming pattern and equatorward ITCZ contraction noted in previous climate simulations. We also showed that the frictional flow, which has yet to receive much attention, is as important as the Ekman flow in generating equatorial heat convergence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Investigating controls on sea ice algal production using E3SMv1.1-BGC

We present the analysis of global sympagic primary production (PP) from 300 years of pre-industrial and historical simulations of the E3SMv1.1-BGC model. The model includes a novel, eight-element sea ice biogeochemical component, MPAS-Seaice zbgc, which is resolved in three spatial dimensions and uses a vertical transport scheme based on internal brine dynamics. Modeled ice algal chlorophyll-a concentrations and column-integrated values are broadly consistent with observations, though chl-a profile fractions indicate that upper ice communities of the Southern Ocean are underestimated. Simulations of polar integrated sea ice PP support the lower bound in published estimates for both polar regions with mean Arctic values of 7.5 and 15.5 TgC/a in the Southern Ocean. However, comparisons of the polar climate state with observations, using a maximal bound for ice algal growth rates, suggest that the Arctic lower bound is a significant underestimation driven by biases in ocean surface nitrate, and that correction of these biases supports as much as 60.7 TgC/a of net Arctic PP. Simulated Southern Ocean sympagic PP is predominantly light-limited, and regional patterns, particularly in the coastal high production band, are found to be negatively correlated with snow thickness.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Sea Ice Rheology Experiment (SIREx): 1. Scaling and Statistical Properties of Sea-Ice Deformation Fields

As the sea-ice modeling community is shifting to advanced numerical frameworks, developing new sea-ice rheologies, and increasing model spatial resolution, ubiquitous deformation features in the Arctic sea ice are now being resolved by sea-ice models. Initiated at the Forum for Arctic Modeling and Observational Synthesis, the Sea Ice Rheology Experiment (SIREx) aims at evaluating state-of-the-art sea-ice models using existing and new metrics to understand how the simulated deformation fields are affected by different representations of sea-ice physics (rheology) and by model configuration. Part 1 of the SIREx analysis is concerned with evaluation of the statistical distribution and scaling properties of sea-ice deformation fields from 35 different simulations against those from the RADARSAT Geophysical Processor System (RGPS). For the first time, the viscous-plastic (and the elastic-viscous-plastic variant), elastic-anisotropic-plastic, and Maxwell-elasto-brittle rheologies are compared in a single study. We find that both plastic and brittle sea-ice rheologies have the potential to reproduce the observed RGPS deformation statistics, including multi-fractality. Model configuration (e.g., numerical convergence, atmospheric representation, spatial resolution) and physical parameterizations (e.g., ice strength parameters and ice thickness distribution) both have effects as important as the choice of sea-ice rheology on the deformation statistics. It is therefore not straightforward to attribute model performance to a specific rheological framework using current deformation metrics. In light of these results, we further evaluate the statistical properties of simulated Linear Kinematic Features in a SIREx Part 2 companion paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Influence of Antarctic and Greenland Continental Shelf Circulation on High‐Latitude Oceans in E3SM

The science objectives of this project are to simulate and understand the impacts of both deep-basin warm-water intrusions and land-ice melt on the continental shelf circulations and sea-ice distributions around the margins of Greenland and Antarctica. As well, the role of subsurface ocean heat from the Atlantic on declining sea-ice cover in the Arctic is explored. Mesoscale processes and fine bathymetry are implicated in cross-shelf property transports around both Greenland and Antarctica. Therefore, we configured and ran an atmospheric reanalysis-forced global ocean/sea-ice simulation on a grid that reduces from 8 km at the Equator to 2 km at the poles (UH8to2) with 60 vertical levels. It was produced using the Energy Exascale Earth System Model ‘‘HiLAT’’ code (E3SMv0-HiLAT) that uses the Parallel Ocean Program (POP) and CICE5 as its ocean and sea-ice components, respectively. Two main UH8to2 simulations were carried out: one for 1975-2009 and the other for July 2016-2020 after it was initialized from a 1/25° data-assimilative ocean/sea-ice prediction system ocean/sea-ice state. The UH8to2 is not coupled to an active land-ice model. Rather, land-ice melt is represented by observationally informed freshwater fluxes (FWFs). Short (multi-year) UH8to2 simulations were conducted to understand sensitivities when Greenland ice sheet (GrIS) melt is released only at the ocean surface or when it is distributed over the upper water column in accordance with fjord melt plume behavior; these cases were compared with a no GrIS melt case. West Greenland continental shelf currents were fastest in the vertical distribution case and an increase in baroclinic conversion at the shelf break associated with increased eddy kinetic energy was found relative to the surface release case. Further, salinity is lower and meltwater volume greater in the eastern Labrador Sea in the vertical distribution case. For the Arctic, the veracity of the UH8to2 was evaluated for 2017-2020 using available observations. Simulated seasonal sea-ice thickness and concentration are realistic, but the ice is unrealistically thin in the central and eastern Arctic in the fall. Comparisons of vertical sections of ocean temperature, salinity, and buoyancy collected from Ice-Tethered Profilers (ITPs) in the eastern Arctic in the fall and winter of 2019/2020 and co-located/concurrent UH8to2 fields show the stratification over the top 100 m of the water column is too low in the model, the simulated mixed layer too deep, and the simulated subsurface Atlantic Water (AW) too warm; these biases may contribute to the sea-ice biases. A model intercomparison study using the UH8to2 and a forced 1/25° regional Arctic ocean/sea-ice (uses the HYbrid Coordinate Ocean Model and CICE5) simulation further investigates the relationship between AW and sea-ice in the eastern Arctic. The models show a mesoscale-rich pulse of Atlantic Water extending into the eastern basin that reaches maximum intensity in late winter of 2018, after which it decreases in strength. Concurrent and co-located sea-ice melt or the inhibition of sea-ice growth is seen and is attributed to halocline mesoscale eddies doming into the mixed layer with convection bringing this heat into the vicinity of the sea-ice.

58 GEOSCIENCES↗
DOE OSTI2020

Expanding the design space of stratospheric aerosol geoengineering to include precipitation-based objectives and explore trade-offs

Previous climate modeling studies demonstrate the ability of feedback-regulated, stratospheric aerosol geoengineering with injection at multiple independent latitudes to meet multiple simultaneous temperature-based objectives in the presence of anthropogenic climate change. However, the impacts of climate change are not limited to rising temperatures but also include changes in precipitation, loss of sea ice, and many more; knowing how a given geoengineering strategy will affect each of these climate metrics is vital to understanding the limits and trade-offs of geoengineering. Here, we first introduce a new method of visualizing the design space in which desired climate outcomes are represented by 2-D surfaces on a 3-D graph. Surface orientations represent how different injection choices influence that objective, and intersecting surfaces represent objectives which can be met simultaneously. Using this representation as a guide, we present simulations of two new strategies for feedback-regulated aerosol injection, using the Community Earth System Model with the Whole Atmosphere Community Climate Model – CESM1(WACCM). The first simultaneously manages global mean temperature, tropical precipitation centroid, and Arctic sea ice extent, while the second manages global mean precipitation, tropical precipitation centroid, and Arctic sea ice extent. Both simulations control the tropical precipitation centroid to within 5 % of the goal, and the latter controls global mean precipitation to within 1% of the goal. Additionally, the first simulation overcompensates sea ice, while the second undercompensates sea ice; all of these results are consistent with the expectations of our design space model. In addition to showing that precipitation-based climate metrics can be managed using feedback alongside other goals, our simulations validate the utility of our design space visualization in predicting our climate model behavior under a given geoengineering strategy, and together they help illustrate the fundamental limits and trade-offs of stratospheric aerosol geoengineering.

58 GEOSCIENCES↗