Search NASA⌕ Search

SEARCH · Search NASA

Results for “sea ice concentration”

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

Ocean Variability Beneath the Filchner-Ronne Ice Shelf Inferred From Basal Melt Rate Time Series

Fourteen phase-sensitive radars (ApRES) were deployed on the Filchner-Ronne Ice Shelf (FRIS) to measure variability in its basal melt rate. Melt rates from sites along the Ronne Depression vary seasonally, consistent with the dynamics of the propagation of seasonal dense water from the western ice front into the cavity. Several sites at the back of the FRIS cavity feature a signal with two seasonal maxima. Sub-ice shelf oceanographic data available from one of the sites indicate that this signal is caused by two different pathways followed by the same source water. Inter-annual variability is strongest along a direct flow pathway between western Ronne Ice Front and western Berkner Island. Highest melting occurred in 1999 and 2018, following anomalously low summer sea-ice concentrations in front of the ice shelf. Inter-annual melt rate variability at the back of the FRIS cavity is limited. If present, it is expressed as a suppression or delay in the arrival of the seasonal melt rate minimum, which can be understood in terms of inter-annual stratification changes and variable inflow pathways toward the sites. Long term mean ApRES melt rates agree with estimates from satellite data over eastern FRIS. However, the satellite estimates overstate the area of active basal freezing in the western part of the ice shelf. Furthermore, the temporal melt rate variability from the satellite estimates exaggerates the range of variability at both seasonal and inter-annual time scales with any correspondence between the in-situ and remotely derived inter-annual variability being limited to a single site.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of the Planetary Boundary Layer Height From ERA5 Reanalysis With MOSAiC Observations Over the Arctic Ocean

The planetary boundary layer height (PBLH) is a crucial indicator reflecting the region of the atmosphere characterized by continuous turbulence. Here, we use radiosonde and surface meteorological observations (4–7 times per day, year-round measurements) during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition to derive the PBLH (PBLH MOSAiC ), and further evaluate the PBLH from the ERA5 reanalysis (PBLH ERA5 ). Comparisons between PBLH MOSAiC and PBLH ERA5 from different perspectives reveal that: (a) The overestimation of PBLH ERA5 when the sea ice concentration is >90% is significant with the centered root mean squared error reaching up to 201 m; (b) The difference between the two products is notably pronounced in cold seasons, while it is comparatively diminished in warm seasons; (c) In neutral boundary layers, differences in PBLH ERA5 are larger compared with stable and convective boundary layers. In addition, the analysis of error sources indicates that the bias of PBLH ERA5 is sensitive to the bias of vertical thermal structure and wind speed profiles in ERA5 data sets in all conditions. Finally, we find a Random Forest model effectively reduces the bias of PBLH ERA5 with the index of agreement reaching up to 0.71 in the test data set, while a multiple linear regression demonstrates comparable performance to the Random Forest model.

54 ENVIRONMENTAL SCIENCES↗

Model-observation discrepancies in Arctic moisture intrusions: causes and pathways for improved simulation

Arctic moisture intrusions (MIs), narrow filaments of strong moisture transport, are key drivers of poleward moisture flux and Arctic weather extremes, yet their representation in climate models is poorly understood. Using a new Arctic MI detection algorithm, we document persistent biases across three CMIP generations (CMIP3–CMIP6): models overestimate MI occurrence over the Pacific sector and underestimate it over the Atlantic sector. These errors stem from misrepresented midlatitude westerly jets, with an equatorward North Atlantic jet associated with too few Atlantic MIs, and a poleward, weakened North Pacific jet linked to too many Pacific MIs. Experiments that correct sea surface temperature and sea ice concentration biases and increase atmospheric resolution improve jet structure and MI statistics, while a cloud-locking simulation indicates that better high-frequency cloud–radiation–circulation interactions can yield further gains. Our results clarify pathways to reducing long-standing MI and jet biases, providing guidance for improving simulations of Arctic and midlatitude climate.

54 ENVIRONMENTAL SCIENCES↗

Northern Hemisphere Winter Air Temperature Patterns and Their Associated Atmospheric and Ocean Conditions

The Northern Hemisphere (NH) has experienced winter Arctic warming and continental cooling in recent decades, but the dominant patterns in winter surface air temperature (SAT) are not well understood. Here, a self-organizing map (SOM) analysis is performed to identify the leading patterns in winter daily SAT fields from 1979 to 2018, and their associated atmospheric and ocean conditions are also examined. Three distinct winter SAT patterns with two phases of nearly opposite signs and a time scale of 7–12 days are found: one pattern exhibits concurrent SAT anomalies of the same sign over North America (NA) and northern Eurasia, while the other two patterns show SAT anomalies of opposite signs between, respectively, NA and the Bering Sea, and the Kara Sea and East Asia (EA). Winter SAT variations may arise from changes in the SOM frequencies. Specifically, the observed increasing trends of winter cold extremes over NA, central Eurasia, and EA during 1998–2013 can be understood as a result of the increasing occurrences of some specific SAT patterns. These SOMs are closely related to poleward advection of midlatitude warm air and equatorward movements of polar cold airmass. These meridional displacements of cold and warm airmasses cause concurrent anomalies over different regions not only in SAT but also in water vapor and surface downward longwave radiation. Anomalous sea surface temperatures in the tropical Pacific, midlatitude North Pacific, and North Atlantic and anomalous Arctic sea ice concentrations also concur to support and maintain the anomalous atmospheric circulation that causes the SAT anomalies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Collaborative Research: Advancing Arctic Climate Projection Capability at Seasonal to Decadal Scales (Final Technical Report)

The Regional Arctic System Model (RASM) at process resolving configurations has been used to (i) advance understanding of physical processes and feedbacks involved in Arctic amplification and (ii) understand and potentially reduce uncertainty in prediction of arctic climate change at seasonal to decadal scales. RASM consists the atmosphere (Weather and Research Forecasting model, WRF), ocean (Parallel Ocean Program, POP), sea ice (CICE), land hydrology (Variable Infiltration Capacity model, VIC), river routing scheme (RVIC), marine biogeochemistry components and the coupling framework (CPL7). Its domain is pan-Arctic, with the atmosphere and land components configured on a 50-km or 25-km grid and four configurations of the ocean and sea ice components: 1/12°(~9.3km) or 1/48°(~2.4km) and 45 or 60 vertical layers. These RASM configurations have been motivated by the emerging exascale capability for high performance computing to improve model fidelity. The dynamical downscaling of reanalysis allows comparison of RASM results with observations in place and time to: (i) advance system level understanding of physical processes and coupling involved in an event, (ii) optimize model parameter space, (iii) diagnose and reduce model biases and (iv) produce realistic and consistent across all the components initial conditions for predictions and predictability studies, which are all unique capabilities not available in global Earth System Models (ESMs). An evaluation of RASM 1.0 (Cassano et al. 2017) revealed that it had a cold bias over the oceans and a warm bias over land areas due largely to cloud and radiation biases in the model, with too little cloud cover simulated over land and too much cloud cover simulated over sub-polar oceans. This study has motivated an upgrade to WRF version 3.7.1 in RASM and allowed for the inclusion of the radiative impact of convective clouds. A variety of atmospheric physics parameterizations were evaluated against observations (e.g. data from the Arctic Clouds in Summer Experiment (ACSE); Sedlar et al. 2020) to identify an optimal suite of WRF physics options in RASM. The RASM with the optimized WRF physics were used to study the impact of strong mesoscale winds over the ocean around the southern tip of Greenland (DuVivier and Cassano 2016) and their impact on oceanic convection (DuVivier et al. 2017a). Data from the PolarWinds field campaign were used to evaluate WRF boundary layer physics and resolution impacts on the simulation of a Greenland barrier wind event (DuVivier et al. 2017b). The RVIC streamflow routing model has been implemented in RASM to realistically represent high-resolution streamflow processes (Hamman et al. 2017) and to couple the land buoyancy fluxes to the ocean. The RASM-RVIC high-resolution data set of all coastal freshwater fluxes for the Arctic drainage basin and surrounding areas for 1979-2014 was published as a separate product (https://doi.org/10.5281/zenodo.293037). The fidelity of atmospheric momentum transfer to and the response of polar marine Ekman layer in RASM and Community Earth System Model (CESM) was investigated by Roberts et al. (2015). The increased frequency of oceanic flux exchange in CESM, following the RASM guidance, caused a considerable increase in the median inertial ice speed across the Southern Ocean and parts of the Arctic. A comprehensive evaluation of the RASM1.0 atmosphere-ocean-sea ice-land interface was completed by Brunke et al. (2018). RASM was also demonstrated for its capability to simulate extreme events in agreement with observations in space and time (Lee et al. submitted). In particular, the development of three open water events, known as polynyas, have been simulated north of Greenland in February of 2011, 2017 and 2018, in agreement with satellite observations for the past four decades. The optimized RASM sea ice results have been favorably evaluated against satellite observations and a subset of eleven CMIP6 models (Watts et al. submitted). In a complementary project, Jin et al. (2018) have shown that RASM with higher-resolution and new sea-ice processes contributed to lower model errors in sea-ice conditions, concentrations of nutrients and ice algae, in comparison to results from the coarse-resolution (1°) CESM. In two other complementary studies, RASM results were used (i) to explain areas of concentrated use by bowhead whales, the seasonal progression in the use, and the physical environment within those areas (Citta et al. 2015) and (ii) for a synthesis of fall bowhead whales distribution and migration in the Bering-Chukchi-Beaufort (BCB) Sea to investigate whale movements and feeding to the local ocean hydrography and currents (Citta et al. 2018). However, the multi-decadal output from the CESM Large Ensemble yielded unrealistic forcing. Instead, the shorter NCEP CFSv2 9-month forecasts were successfully tested and afforded an increased ensemble size (~30) to demonstrate gains of dynamical downscaling at sub-seasonal to intra-annual time scales. The improved model physics and coupling among RASM model components have yielded more realistic representation of the sea ice cover and consistent across all model components initial conditions. Consequently, RASM demonstrates significant gains compared to simulation of sea ice in the NCEP reanalysis. In addition, RASM 6-month ensemble forecasts yield very realistic sea ice distribution, which demonstrates both significant gains of dynamical downscaling and the continued impact of the initial conditions on forecasts out to 6 months (https://nps.edu/web/rasm/predictions). A follow up study using RASM for dynamical downscaling of the more realistic CESM initialized Decadal Prediction Large Ensemble output is currently ongoing as part of the DOE RGMA HiLAT-RASM project.

54 ENVIRONMENTAL SCIENCES↗

Expanded Understanding of the Western Antarctic Peninsula Sea‐Ice Environment Through Local and Regional Observations at Palmer Station

Abstract The Western Antarctic Peninsula (WAP) has been experiencing rapid regional warming since at least the 1950s, however, the impacts of this warming at the local scale are variable and nuanced. Previous studies that have linked sea‐ice variability to biogeochemical cycles and food web dynamics often combine local‐scale biogeochemical data with coarse‐resolution regional satellite sea‐ice data, which may not adequately capture local sea‐ice conditions. In this study, we analyzed local‐scale in situ sea‐ice observations collected as part of a 28‐year record (1992–2020) from the Palmer Long‐Term Ecological Research site at Anvers Island, mid‐WAP, in conjunction with isotopically‐derived sea‐ice meltwater (SIM) fractions and satellite‐derived sea‐ice motion and concentration, to quantify the variability and long‐term trends in local sea‐ice behavior. In situ sea ice observations at Palmer Station displayed higher variability than satellite observations and showed no significant declines over this time, despite region‐wide declines identified in prior studies. Higher spring SIM fractions were attributed to strong northward sea‐ice motion throughout the winter. Applying these local‐scale sea‐ice insights to similarly scaled stratification and chlorophyll‐ a measurements, we found that a longer‐lasting, more consistent sea‐ice pack led to greater water column stratification following the spring sea‐ice retreat. Greater sea‐ice persistence and stronger stratification led to larger peaks in chlorophyll‐ a , though sea‐ice metrics did not explain the positive temporal trends in either stratification strength or chlorophyll‐ a . Through this study, we identify how local sea‐ice observations and meltwater data can enhance satellite data to build an understanding of the intricate connections between ice, water column dynamics, and phytoplankton.

Goodell, E.↗

Sea Ice Collision Risk Assessment for Tidal Turbine Siting in Cook Inlet, Alaska

Cook Inlet has great potential for tidal stream energy development. However, the presence of drifting sea ice could create hazardous collision risks for tidal turbine farms. Before turbines can be installed in Cook Inlet, sites must be surveyed to determine how often sea ice is present, how fast it will be moved by the current, and where the trajectories of drifting sea ice will be concentrated. In this study, we use remotely sensed data to characterize the seasonal sea ice conditions in Cook Inlet, a hydrodynamic model to map the water velocities, and a particle trajectory tracking model to calculate collision risks. The resulting sea ice coverage and collision risk maps will enable tidal energy developers to choose the best locations for the deployment of tidal turbines and other offshore platforms in Cook Inlet.

16 TIDAL AND WAVE POWER↗

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↗

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↗

Hidden Production: On the Importance of Pelagic Phytoplankton Blooms Beneath Arctic Sea Ice

Abstract Recent observations suggest that substantial phytoplankton blooms occur under sea ice on Arctic continental shelves during June and July. This is opposed to the traditional view that no significant biomass is produced in sea‐ice covered waters. However, no observational estimates are available on the Arctic‐wide primary production beneath sea ice. Here, using a fully coupled Arctic system model, we estimate that 63%/41% of the total primary production in the central Arctic occurs in waters covered by sea ice that is ≥50%/≥85% concentration. The total primary production there is increasing at a rate of 5.2% per decade during 1980–2018. Increased light transmission, due to the removal of sea ice, more extensive melt ponds, and thinner sea ice, is implicated as the main cause of increasing trends in primary production.

54 ENVIRONMENTAL SCIENCES↗

The annual cycle and sources of relevant aerosol precursor vapors in the central Arctic during the MOSAiC expedition

Abstract. In this study, we present and analyze the first continuous time series of relevant aerosol precursor vapors from the central Arctic (north of 80° N) during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition. These precursor vapors include sulfuric acid (SA), methanesulfonic acid (MSA), and iodic acid (IA). We use FLEXPART simulations, inverse modeling, sulfur dioxide (SO2) mixing ratios, and chlorophyll a (chl a) observations to interpret the seasonal variability in the vapor concentrations and identify dominant sources. Our results show that both natural and anthropogenic sources are relevant for the concentrations of SA in the Arctic, but anthropogenic sources associated with Arctic haze are the most prevalent. MSA concentrations are an order of magnitude higher during polar day than during polar night due to seasonal changes in biological activity. Peak MSA concentrations were observed in May, which corresponds with the timing of the annual peak in chl a concentrations north of 75° N. IA concentrations exhibit two distinct peaks during the year, namely a dominant peak in spring and a secondary peak in autumn, suggesting that seasonal IA concentrations depend on both solar radiation and sea ice conditions. In general, the seasonal cycles of SA, MSA, and IA in the central Arctic Ocean are related to sea ice conditions, and we expect that changes in the Arctic environment will affect the concentrations of these vapors in the future. The magnitude of these changes and the subsequent influence on aerosol processes remains uncertain, highlighting the need for continued observations of these precursor vapors in the Arctic.

Boyer, Matthew↗

A numerical framework for simulating the atmospheric variability of supermicron marine biogenic ice nucleating particles

Abstract. We present a framework for estimating concentrations of episodically elevated high-temperature marine ice nucleating particles (INPs) in the sea surface microlayer and their subsequent emission into the atmospheric boundary layer. These episodic INPs have been observed in multiple ship-based and coastal field campaigns, but the processes controlling their ocean concentrations and transfer to the atmosphere are not yet fully understood. We use a combination of empirical constraints and simulation outputs from an Earth system model to explore different hypotheses for explaining the variability of INP concentrations, and the occurrence of episodic INPs, in the marine atmosphere. In our calculations, we examine the following two proposed oceanic sources of high-temperature INPs: heterotrophic bacteria and marine biopolymer aggregates (MBPAs). Furthermore, we assume that the emission of these INPs is determined by the production of supermicron sea spray aerosol formed from jet drops, with an entrainment probability that is described by Poisson statistics. The concentration of jet drops is derived from the number concentration of supermicron sea spray aerosol calculated from model runs. We then derive the resulting number concentrations of marine high-temperature INPs (at 253 K) in the atmospheric boundary layer and compare their variability to atmospheric observations of INP variability. Specifically, we compare against concentrations of episodically occurring high-temperature INPs observed during field campaigns in the Southern Ocean, the Equatorial Pacific, and the North Atlantic. In this case study, we evaluate our framework at 253 K because reliable observational data at this temperature are available across three different ocean regions, but suitable data are sparse at higher temperatures. We find that heterotrophic bacteria and MBPAs acting as INPs provide only a partial explanation for the observed high INP concentrations. We note, however, that there are still substantial knowledge gaps, particularly concerning the identity of the oceanic INPs contributing most frequently to episodic high-temperature INPs, their specific ice nucleation activity, and the enrichment of their concentrations during the sea–air transfer process. Therefore, targeted measurements investigating the composition of these marine INPs and drivers for their emissions are needed, ideally in combination with modeling studies focused on the potential cloud impacts of these high-temperature INPs.

54 ENVIRONMENTAL SCIENCES↗

Reversible ice sheet thinning in the Amundsen Sea Embayment during the Late Holocene

Cosmogenic-nuclide concentrations in subglacial bedrock cores show that the West Antarctic Ice Sheet (WAIS) at a site between Thwaites and Pope glaciers was at least 35 m thinner than present in the past several thousand years and then subsequently thickened. This is important because of concern that present thinning and grounding line retreat at these and nearby glaciers in the Amundsen Sea Embayment may irreversibly lead to deglaciation of significant portions of the WAIS, with decimeter- to meter-scale sea level rise within decades to centuries. A past episode of ice sheet thinning that took place in a similar, although not identical, climate was not irreversible. We propose that the past thinning–thickening cycle was due to a glacioisostatic rebound feedback, similar to that invoked as a possible stabilizing mechanism for current grounding line retreat, in which isostatic uplift caused by Early Holocene thinning led to relative sea level fall favoring grounding line advance.

54 ENVIRONMENTAL SCIENCES↗

Environmental Factors Associated With Fall Phytoplankton Blooms in the Northern Bering and Chukchi Seas

This study investigates environmental drivers of fall phytoplankton blooms in the Arctic, focusing on the northern Bering and Chukchi seas. Random Forests models were used to analyze covariates of fall phytoplankton blooms from 2013 to 2018, incorporating shipboard, remote sensing, and modeled environmental properties. Four regional models and one comprehensive all-station model considered fall as well as midsummer conditions. Midsummer properties included suspended particulate matter, chlorophyll-a, and the proportion of degraded pheophytin to chlorophyll-a used as a proxy for bloom stage. Open water duration was one of the highest ranked factors in predicting fall blooms. Open water duration also influences the stage of midsummer (July) blooms as indicated by pheophytin proportions, which in turn were the highest-ranked factor for predicting fall bloom events in the Chirikov Basin (northern Bering Sea between St. Lawrence Island and the Bering Strait) and the Chukchi Sea. Wind direction, specifically easterly winds, was an important predictor in the northern Bering Sea. Maximum wind speed ranked highly at stations located within the nutrient-poor Alaska Coastal Current in the Chukchi Sea. However, stormy days, average and maximum wind speeds generally ranked low in importance as a predictor of fall bloom events. Other parameters, including photosynthetic active radiation, modeled nutrient concentrations, mixed layer depth, and time since sea ice breakup date showed strong but regionally varying relationships with fall blooms. Altogether, results from these Random Forests models suggest that high wind events and storms in the absence of sea ice provide an incomplete narrative for initiating fall bloom events.

Gaffey, C. B. [Clark University, Worcester, MA (Un↗

Geometric remapping of particle distributions in the Discrete Element Model for Sea Ice (DEMSI v0.0)

Abstract. A new sea ice dynamical core, the Discrete Element Model for Sea Ice (DEMSI), is under development for use in coupled Earth system models. DEMSI is based on the discrete element method, which models collections of ice floes as interacting Lagrangian particles. In basin-scale sea ice simulations the Lagrangian motion results in significant convergence and ridging, which requires periodic remapping of sea ice variables from a deformed particle configuration back to an undeformed initial distribution. At the resolution required for Earth system models we cannot resolve individual sea ice floes, so we adopt the sub-grid-scale thickness distribution used in continuum sea ice models. This choice leads to a series of hierarchical tracers depending on ice fractional area or concentration that must be remapped consistently. The circular discrete elements employed in DEMSI help improve the computational efficiency at the cost of increased complexity in the effective element area definitions for sea ice cover that are required for the accurate enforcement of conservation. An additional challenge is the accurate remapping of element values along the ice edge, the location of which varies due to the Lagrangian motion of the particles. In this paper we describe a particle-to-particle remapping approach based on well-established geometric remapping ideas that enforces conservation, bounds preservation, and compatibility between associated tracer quantities, while also robustly managing remapping at the ice edge. One element of the remapping algorithm is a novel optimization-based flux correction that enforces concentration bounds in the case of nonuniform motion. We demonstrate the accuracy and utility of the algorithm in a series of numerical test cases.

58 GEOSCIENCES↗

Sea ice evolution along the Northern Sea Route and implications for trans-Arctic shipping from 2021 through 2060

Arctic surface temperatures warmed at twice the global average in the second half of the 20 th century due to Arctic amplification (AA), a phenomenon predominantly caused by regional polar changes, like the melting of perennial sea ice and reduced sea ice extent (leading to more solar radiation being absorbed by the ocean surface as opposed to being reflected back to space by the ice surface). AA is projected to reach a factor of three even if the climate is stabilized by the mid 21 st century by reduced greenhouse gas emissions. In all emissions scenarios, AA is projected to lead to temperature changes at least 2.4 times larger than global mean surface temperature changes occurring between 2070 and 2100. The ice-albedo feedback, which occurs when the polar-marine surface absorbs more radiation as highly reflective sea ice melts, is reversible such that the premise of a runaway process is no longer accepted as a realistic possibility. No matter what actions are taken to reduce CO 2 concentrations in the atmosphere from now on, two different methods of predicting an ice-free Arctic suggest that perennial sea ice will mostly disappear in September by the year 2050. If there is no reduction in anthropogenic CO 2 and methane emissions, that scenario could occur sooner than 2030. Defining Arctic navigability as safe and economic passage of Polar Class 7 cargo ships without need of an escorting icebreaker, no single trans Arctic ship route will be navigable year-round in the first half of the 21 st century, including in the strong emission scenarios. However, seasonal trans-Arctic navigability will increase this century. An estimate on the number of days per year that the Northern See Route (NSR) will be navigable in the future is beyond the scope of this report. Along Northern Sea Routes 5 and 6, which run close to the Russian coast and Yamal LNG plant, an ARC 7 ice class LNG tanker, the equivalent of a Polar Class 3 (PC3) vessel, will be at low risk in December through April at some point during the current decade. May will continue to entail some risk (more than April) along relatively short segments of these routes through the end of the next decade (2030-2039). Come June, snow rapidly melts away, and thereafter the underlying sea ice begins becomes thinner and less concentrated, greatly reducing risk. However, neither path is desirable for ARC 7 tankers due to shallow bathymetry in Sannikov Strait, and a more desirable path for these ships passes to the north of the New Siberian Islands (route 20, discussed further below). Conventional LNG tankers (i.e., non-ice-strengthened vessels according to the IMO classification), without icebreaker escorts, will continue to encounter dangerous or impassable conditions along many sections of the NSR through the end of this decade for many months of the year. Through 2039, these conventional tankers will be able to operate safely along NSR 6 from August through October. By 2040-2049, the span of safe operation increases to August through November, and by 2050-2059 it increases to July through November, assuming a northward deviation to avoid the Sannikov Strait.

54 ENVIRONMENTAL SCIENCES↗

Low ozone dry deposition rates to sea ice during the MOSAiC field campaign: Implications for the Arctic boundary layer ozone budget

Dry deposition to the surface is one of the main removal pathways of tropospheric ozone (O 3 ). We quantified for the first time the impact of O 3 deposition to the Arctic sea ice on the planetary boundary layer (PBL) O 3 concentration and budget using year-round flux and concentration observations from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) campaign and simulations with a single-column atmospheric chemistry and meteorological model (SCM). Based on eddy-covariance O 3 surface flux observations, we find a median surface resistance on the order of 20,000 s m –1 , resulting in a dry deposition velocity of approximately 0.005 cm s –1 . This surface resistance is up to an order of magnitude larger than traditionally used values in many atmospheric chemistry and transport models. The SCM is able to accurately represent the yearly cycle, with maxima above 40 ppb in the winter and minima around 15 ppb at the end of summer. However, the observed springtime ozone depletion events are not captured by the SCM. In winter, the modelled PBL O 3 budget is governed by dry deposition at the surface mostly compensated by downward turbulent transport of O 3 towards the surface. Advection, which is accounted for implicitly by nudging to reanalysis data, poses a substantial, mostly negative, contribution to the simulated PBL O 3 budget in summer. During episodes with low wind speed (<5 m s –1 ) and shallow PBL (<50 m), the 7-day mean dry deposition removal rate can reach up to 1.0 ppb h –1 . Our study highlights the importance of an accurate description of dry deposition to Arctic sea ice in models to quantify the current and future O 3 sink in the Arctic, impacting the tropospheric O 3 budget, which has been modified in the last century largely due to anthropogenic activities.

54 ENVIRONMENTAL SCIENCES↗

A review of black carbon in snow and ice and its impact on the cryosphere

Black carbon (BC) has emerged as an important short-lived climate forcer. Due to its light absorption properties, BC can darken the snow/ice surface, affect the energy balance, and further lead to acceleration of the cryospheric melting. By reviewing the recent published literatures, we present an overview of the historical changes, spatial distribution of BC in snow/ice, and how these changes are related to the melting of the cryosphere (e.g., glaciers, snow cover, and sea ice). Ice core records show an abrupt increase of BC concentrations that began in the 1850s and continued throughout the 20th century, which is consistent with an increase in BC emissions owing to industrialization. A decrease in the amount of BC since 1970s in Arctic and European ice cores has been partially attributed to the Clean Air Act. However, in the Himalayas, BC records show a continuous increase during this period. BC concentrations in snow and ice in the mid-latitude regions are one to two orders of magnitude higher than those in the polar regions. In particular, BC concentrations in aged snow and granular ice in ablation areas of mountain glaciers are one to three orders of magnitude higher than those in fresh snow or snowpits in the glacier accumulation areas due to BC accumulation during melting season. BC in the surface snow/ice is responsible for about 20% of the albedo reduction in the Tibetan Plateau during glacier melt season. Globally, observations and modeling results indicate that radiative forcing (RF) induced by BC in snow and ice is highest in the mid-latitudes, ranging from several W m–2 in fresh snow to hundreds of W m–2 in aged snow and granular ice in glacier ablation areas. The large BC-in-snow RF and associated snow albedo feedback lead to an acceleration in the total glacier melt (approximately 20%) and/or a reduction in the duration of the snow cover by several days, resulting in an increase of glacier discharge. Given our limited understanding of the role of BC in cryospheric melting, it is important to synthesize the existing research on the multi-scale processes related to BC in snow and ice to identify the gaps in our understanding of these processes and to propose a path forward to improve the quality of our observations of the aforementioned phenomena to fill these gaps.

Kang, Shichang↗