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 145 records · Page 8

The Aqua-planet Experiment (APE): Response to Changed Meridional SST Profile

This paper explores the sensitivity of Atmospheric General Circulation Model (AGCM) simulations to changes in the meridional distribution of sea surface temperature (SST). The simulations are for an aqua-planet, a water covered Earth with no land, orography or sea- ice and with specified zonally symmetric SST. Simulations from 14 AGCMs developed for Numerical Weather Prediction and climate applications are compared. Four experiments are performed to study the sensitivity to the meridional SST profile. These profiles range from one in which the SST gradient continues to the equator to one which is flat approaching the equator, all with the same maximum SST at the equator. The zonal mean circulation of all models shows strong sensitivity to latitudinal distribution of SST. The Hadley circulation weakens and shifts poleward as the SST profile flattens in the tropics. One question of interest is the formation of a double versus a single ITCZ. There is a large variation between models of the strength of the ITCZ and where in the SST experiment sequence they transition from a single to double ITCZ. The SST profiles are defined such that as the equatorial SST gradient flattens, the maximum gradient increases and moves poleward. This leads to a weakening of the mid-latitude jet accompanied by a poleward shift of the jet core. Also considered are tropical wave activity and tropical precipitation frequency distributions. The details of each vary greatly between models, both with a given SST and in the response to the change in SST. One additional experiment is included to examine the sensitivity to an off-equatorial SST maximum. The upward branch of the Hadley circulation follows the SST maximum off the equator. The models that form a single precipitation maximum when the maximum SST is on the equator shift the precipitation maximum off equator and keep it centered over the SST maximum. Those that form a double with minimum on the equatorial maximum SST shift the double structure off the equator, keeping the minimum over the maximum SST. In both situations only modest changes appear in the shifted profile of zonal average precipitation. When the upward branch of the Hadley circulation moves into the hemisphere with SST maximum, the zonal average zonal, meridional and vertical winds all indicate that the Hadley cell in the other hemisphere dominates.

AGCM↗

ISMIP6-based Antarctic projections to 2100: simulations with the BISICLES ice sheet model

The contribution of the Antarctic Ice Sheet is one of the most uncertain components of sea level rise to 2100. Ice sheet models are the primary tool for projecting future sea level contribution from continental ice sheets. The Ice Sheet Model Intercomparison for the Coupled Model Intercomparison Phase 6 (ISMIP6) provided projections of the ice sheet contribution to sea level over the 21st century, quantifying uncertainty due to ice sheet model, climate model, emission scenario, and uncertain parameters. We present simulations following the ISMIP6 framework with the BISICLES ice sheet model and new experiments extending the ISMIP6 protocol to more comprehensively sample uncertainties in future climate, ice shelf sensitivity to ocean melting, and their interactions. These results contributed to the land ice projections of Edwards et al. (2021), which formed the basis of sea level projections for the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (AR6). Our experiments show the important interplay between surface mass balance processes and ocean-driven melt in determining Antarctic sea level contribution. Under higher-warming scenarios, high accumulation offsets more ocean-driven mass loss when sensitivity to ocean-driven melt is low. Conversely, we show that when sensitivity to ocean warming is high, ocean melting drives increased mass loss despite high accumulation. Overall, we simulate a sea level contribution range across our experiments from 2 to 178 mm. Finally, we show that collapse of ice shelves due to surface warming increases sea level contribution by 25 mm relative to the no-collapse experiments, for both moderate and high sensitivity of ice shelf melting to ocean forcing.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of Arctic Warming in Mid-Pliocene Climate Simulations

Palaeoclimate simulations improve our understanding of the climate, inform us about the performance of climate models in a different climate scenario, and help to identify robust features of the climate system. Here, we analyse Arctic warming in an ensemble of 16 simulations of the mid-Pliocene Warm Period (mPWP), derived from the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2). The PlioMIP2 ensemble simulates Arctic (60–90∘ N) annual mean surface air temperature (SAT) increases of 3.7 to 11.6 ∘C compared to the pre-industrial period, with a multi-model mean (MMM) increase of 7.2 ∘C. The Arctic warming amplification ratio relative to global SAT anomalies in the ensemble ranges from 1.8 to 3.1 (MMM is 2.3). Sea ice extent anomalies range from −3.0 to −10.4×106 km2, with a MMM anomaly of −5.6×106 km2, which constitutes a decrease of 53 % compared to the pre-industrial period. The majority (11 out of 16) of models simulate summer sea-ice-free conditions (≤1×106 km2) in their mPWP simulation. The ensemble tends to underestimate SAT in the Arctic when compared to available reconstructions, although the degree of underestimation varies strongly between the simulations. The simulations with the highest Arctic SAT anomalies tend to match the proxy dataset in its current form better. The ensemble shows some agreement with reconstructions of sea ice, particularly with regard to seasonal sea ice. Large uncertainties limit the confidence that can be placed in the findings and the compatibility of the different proxy datasets. We show that while reducing uncertainties in the reconstructions could decrease the SAT data–model discord substantially, further improvements are likely to be found in enhanced boundary conditions or model physics. Lastly, we compare the Arctic warming in the mPWP to projections of future Arctic warming and find that the PlioMIP2 ensemble simulates greater Arctic amplification than CMIP5 future climate simulations and an increase instead of a decrease in Atlantic Meridional Overturning Circulation (AMOC) strength compared to pre-industrial period. The results highlight the importance of slow feedbacks in equilibrium climate simulations, and that caution must be taken when using simulations of the mPWP as an analogue for future climate change.

Palaeoclimate simulations↗

Multilayer sea ice in the Arctic - Model- and satellite-derived

The behavior of sea ice is a strong function of ice strength which, in the central Arctic, depends on the relative fractions of multiyear and first-year ice coverage. A dynamic-thermodynamic sea ice model has been configured to include both multiyear and first-year ice. The model has been run for a 7-year period, 1979-1985, for which a satellite microwave record of multiyear sea ice coverage is available. The simulated concentrations of multiyear ice are slightly higher than the satellite-derived concentrations, although the total coverage of simulated multiyear ice agrees well with the satellite-derived coverage. Both the simulated and the satellite-derived coverages in the early 1980s show large interannual variations that represent potentially large perturbations in the mass balance of Arctic pack ice. These fluctuations are generally consistent in the model and satellite data, and they are generally consistent with year-to-year differences in the large-scale wind forcing. In at least some cases the wintertime distribution of multiyear ice is influenced by the pattern of ice drift as far back as the previous summer. The late December coverage of multiyear ice in the Alaskan region appears to foreshadow summer ice conditions in the same region. The short-term response to an intense synoptic system is considerably larger in the satellite data than in the model results.

Walsh, John E.↗

Radiative Properties of Smoke and Aerosol Over Land Surfaces

This talk discusses smoke and aerosol's radiative properties with particular attention to distinguishing the measurement over clear sky from clouds over land, sea, snow, etc. surfaces, using MODIS Airborne Simulator data from (Brazil, arctic sea ice and tundra and southern Africa, west Africa, and other ecosystems. This talk also discusses the surface bidirectional reflectance using Cloud Absorption Radiometer, BRDF measurements of Saudi Arabian desert, Persian Gulf, cerrado and rain forests in Brazil, sea ice, tundra, Atlantic Ocean, Great Dismal Swamp, Kuwait oil fire smoke. Recent upgrades to instrument (new TOMS UVA channels at 340 and 380 planned use in Africa (SAFARI 2000) and possibly for MEIDEX will also be discussed. This talk also plans to discuss the spectral variation of surface reflectance over land and the sensitivity of off-nadir view angles to correlation between visible near-infrared reflectance for use in remote sensing of aerosol over land.

King, Michael D.↗

Status and Plans for Reanalysis at NASA/GMAO

Reanalysis plays a critical role in GMAOs goal to enhance NASA's program of Earth observations, providing vital data sets for climate research and the development of future missions. As the breadth of NASAs observations expands to include multiple components of the Earth system, so does the need to assimilate observations from currently uncoupled components of the system in a more physically consistent manner. GMAOs most recent reanalysis of the satellite era, MERRA-2, has completed the period 1980-present, and is now running as a continuing global climate analysis with two- to three-week latency. MERRA-2 assimilates meteorological and aerosol observations as a weakly coupled assimilation system as a first step toward GMAOs longer term goal of developing an integrated Earth system analysis (IESA) capability that will couple assimilation systems for the atmosphere, ocean, land and chemistry. The GMAO strategy is to progress incrementally toward an IESA through an evolving combination of coupled systems and offline component reanalyses driven by, for example, MERRA-2 atmospheric forcing. Most recently, the GMAO has implemented a weakly coupled assimilation scheme for analyzing ocean skin temperature within the existing atmospheric analysis. The scheme uses background fields from a near-surface ocean diurnal layer model to assimilate surface-sensitive radiances plus in-situ observations along with all other observations in the atmospheric assimilation system. In addition, MERRA-2-driven simulations of the ocean (plus sea ice) and atmospheric chemistry (for the EOS period) are currently underway, as is the development of a coupled atmosphere-ocean assimilation system. This talk will describe the status of these ongoing efforts and the planned steps toward an IESA capability for climate research.

MERRA-2↗

Projected Impact of Climate Change on the Energy Budget of the Arctic Ocean by a Global Climate Model

The annual energy budget of the Arctic Ocean is characterized by a net heat loss at the air-sea interface that is balanced by oceanic heat transport into the Arctic. The energy loss at the air-sea interface is due to the combined effects of radiative, sensible, and latent heat fluxes. The inflow of heat by the ocean can be divided into two components: the transport of water masses of different temperatures between the Arctic and the Atlantic and Pacific Oceans and the export of sea ice, primarily through Fram Strait. Two 150-year simulations (1950-2099) of a global climate model are used to examine how this balance might change if atmospheric greenhouse gases (GHGs) increase. One is a control simulation for the present climate with constant 1950 atmospheric composition, and the other is a transient experiment with observed GHGs from 1950 to 1990 and 0.5% annual compounded increases of CO2 after 1990. For the present climate the model agrees well with observations of radiative fluxes at the top of the atmosphere, atmospheric advective energy transport into the Arctic, and surface air temperature. It also simulates the seasonal cycle and summer increase of cloud cover and the seasonal cycle of sea-ice cover. In addition, the changes in high-latitude surface air temperature and sea-ice cover in the GHG experiment are consistent with observed changes during the last 40 and 20 years, respectively. Relative to the control, the last 50-year period of the GHG experiment indicates that even though the net annual incident solar radiation at the surface decreases by 4.6 W(per square meters) (because of greater cloud cover and increased cloud optical depth), the absorbed solar radiation increases by 2.8 W(per square meters) (because of less sea ice). Increased cloud cover and warmer air also cause increased downward thermal radiation at the surface so that the net radiation into the ocean increases by 5.0 Wm-2. The annual increase in radiation into the ocean, however, is compensated by larger increases in sensible and latent heat fluxes out of the ocean. Although the net energy loss from the ocean surface increases by 0.8 W (per square meters), this is less than the interannual variability, and the increase may not indicate a long-term trend. The seasonal cycle of heat fluxes is significantly enhanced. The downward surface heat flux increases in summer (maximum 2 of 19 W per square meters or 23% in June) while the upward heat flux increases in winter (maximum of 16 W per square meters or 28% in November). The increased downward flux in summer is due to a combination of increases in absorbed solar and thermal radiation and smaller losses of sensible and latent heat. The increased heat loss in winter is due to increased sensible and latent heat fluxes, which in turn are due to reduced sea-ice cover. On the other hand, the seasonal cycle of surface air temperature is damped, as there is a large increase in winter temperature but little change in summer.

Miller, James R.↗

LongRunMIP - Motivation and Design for a Large Collection of Millennial-Length AO-GCM Simulations

LongRunMIP is the first collection of millennial-length simulations of complex coupled climate models and enables investigations of how these models equilibrate in response to radiative perturbations.We present a model intercomparison project, LongRunMIP, the first collection of millennial-length (1000+ year) simulations of complex coupled climate models with a representation of ocean, atmosphere, sea ice, and land surface, and their interactions. Standard model simulations are generally only a few hundred years long. However, modeling the long-term equilibration in response to radiative forcing perturbation is important for understanding many climate phenomena, such as the evolution of ocean circulation, time-and temperature-dependent feedbacks, and the differentiation of forced signal and internal variability. The aim of LongRunMIP is to facilitate research into these questions by serving as an archive for simulations that capture as much of this equilibration as possible. The only requirement to participate in LongRunMIP is to contribute a simulation with elevated, constant CO2 forcing that lasts at least 1000 years. LongRunMIP is a MIP of opportunity in that the simulations were mostly performed prior to the conception of the archive without an agreed-upon set of experiments. For most models, the archive contains a preindustrial control simulation and simulations with an idealized (typically abrupt) CO2 forcing. We collect 2D surface and top-of-atmosphere fields, and 3D ocean temperature and salinity fields. Here, we document the collection of simulations and discuss initial results, including the evolution of surface and deep ocean temperature and cloud radiative effects. As of summer 2019, the collection includes 50 simulations of 15 models by 10 modeling centers. The data of LongRunMIP are publicly available. We encourage submission of more simulations in the future.

Maria Rugenstein↗

Spatiotemporal evolution of melt ponds on Arctic sea ice

Melt ponds on sea ice play an important role in the Arctic climate system. Their presence alters the partitioning of solar radiation: decreasing reflection, increasing absorption and transmission to the ice and ocean, and enhancing melt. The spatiotemporal properties of melt ponds thus modify ice albedo feedbacks and the mass balance of Arctic sea ice. The Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition presented a valuable opportunity to investigate the seasonal evolution of melt ponds through a rich array of atmosphere-ice-ocean measurements across spatial and temporal scales. In this study, we characterize the seasonal behavior and variability in the snow, surface scattering layer, and melt ponds from spring melt to autumn freeze-up using in situ surveys and auxiliary observations. We compare the results to satellite retrievals and output from two models: the Community Earth System Model (CESM2) and the Marginal Ice Zone Modeling and Assimilation System (MIZMAS). During the melt season, the maximum pond coverage and depth were 21% and 22 ± 13 cm, respectively, with distribution and depth corresponding to surface roughness and ice thickness. Compared to observations, both models overestimate melt pond coverage in summer, with maximum values of approximately 41% (MIZMAS) and 51% (CESM2). This overestimation has important implications for accurately simulating albedo feedbacks. During the observed freeze-up, weather events, including rain on snow, caused high-frequency variability in snow depth, while pond coverage and depth remained relatively constant until continuous freezing ensued. Both models accurately simulate the abrupt cessation of melt ponds during freeze-up, but the dates of freeze-up differ. MIZMAS accurately simulates the observed date of freeze-up, while CESM2 simulates freeze-up one-to-two weeks earlier. This work demonstrates areas that warrant future observation-model synthesis for improving the representation of sea-ice processes and properties, which can aid accurate simulations of albedo feedbacks in a warming climate.

54 ENVIRONMENTAL SCIENCES↗

Should Sea-Ice Modeling Tools Designed for Climate Research Be Used for Short-Term Forecasting?

In theory, the same sea-ice models could be used for both research and operations, but in practice, differences in scientific and software requirements and computational and human resources complicate the matter. Although sea-ice modeling tools developed for climate studies and other research applications produce output of interest to operational forecast users, such as ice motion, convergence, and internal ice pressure, the relevant spatial and temporal scales may not be sufficiently resolved. For instance, sea-ice research codes are typically run with horizontal resolution of more than 3 km, while mariners need information on scales less than 300 m. Certain sea-ice processes and coupled feedbacks that are critical to simulating the Earth system may not be relevant on these scales; and therefore, the most important model upgrades for improving sea-ice predictions might be made in the atmosphere and ocean components of coupled models or in their coupling mechanisms, rather than in the sea-ice model itself. This paper discusses some of the challenges in applying sea-ice modeling tools developed for research purposes for operational forecasting on short time scales, and highlights promising new directions in sea-ice modeling.

54 ENVIRONMENTAL SCIENCES↗

Benchmark estimate of the effect of anthropogenic emissions on the ocean surface

Abstract Investigations into the role of anthropogenic emissions in the occurrence of extreme weather often use a method that compares simulations of atmospheric climate models run under a factual scenario of historical boundary conditions observed during the period of the event against simulations run under a counterfactual scenario of what those boundary conditions might naturally have been over that same period in the absence of anthropogenic emissions. A particular requirement for this experiment design is an accurate estimation of ocean surface boundary conditions for use by the counterfactual natural simulations. Here we use output from the CMIP5 multi‐climate‐model archive to develop a robust estimate of sea surface temperatures and sea ice conditions for use in counterfactual natural simulations, intended as a benchmark estimate to facilitate comparison across climate models and across studies. This development includes tests to ensure that the final estimate is stable from year‐to‐year and stable against other perturbations to the methodology, as well as consideration of the strengths and weaknesses in comparison to other available attributable warming estimates. While this estimate is tailored specifically for the International CLIVAR C20C+ Detection and Attribution Project, it can be used by related projects as well.

Stone, Dáithí A.↗

The sensitivity of the general circulation to Arctic Sea ice boundaries - A numerical experiment

Results are presented for a set of numerical experiments conducted with the Goddard (formerly GISS) general circulation model. The experiments were designed to test the model atmospheric response to a single fixed and specified parameter, the total ice cover in the Davis Strait, Barents Sea, East Greenland Sea, Sea of Okhotsk and Bering Sea. Margin variations are considered that are substantially smaller than those involved in ice age or ice-free Arctic simulations. Anomaly is defined as the mean of two runs corresponding to climatological maximum sea ice conditions. Model results indicate that the ice margin anomalies are capable of altering local climates in certain regions of high and middle latitudes. Possible interactions between high latitudes and subtropical regions are suggested.

Herman, G. F.↗

Seasonal-to-Interannual Variability in Antarctic Sea-Ice Dynamics, and Its Impact on Surface Fluxes and Water Mass Production

Strong seasonal and interannual signals in Antarctic bottom-water outflow remain unexplained yet are highly correlated with anomalies in net sea-ice growth in coastal polynyas. The mechanisms responsible for driving salination and replenishment and rejuvenation of the dense shelf "source" waters likely also generate pulses of bottom water outflow. The objective of this research is to investigate time-scales of variability in the dynamics of sea-ice in the Southern Ocean in order to determine the primary sites for production of dense shelf waters. We are using a merged satellite/buoy sea-ice motion data set for the period 1978-present day to compute the dynamics of opening and closing of coastal polynyas over the continental shelf. The Ocean Circulation and Climate Advanced Model (OCCAM) ocean general circulation model with coupled sea-ice dynamics is presently forced using National Center for Environmental Prediction (NCEP) data to simulate fluxes and the salination impact of the ocean shelf regions. This work is relevant in the context of measuring the influence of polar sea-ice dynamics upon polar ocean characteristics, and thereby upon global thermohaline ocean circulation. Interannual variability in simulated net freezing rate in the Southern Weddell Sea is shown for the period 1986-1993. There is a pronounced maximum of ice production in 1988 and minimum in 1991 in response to anomalies in equatorward meridional wind velocity. This follows a similar approximate 8-year interannual cycle in Sea Surface Temperature (SST) and satellite-derived ice-edge anomalies reported elsewhere as the "Antarctic Circumpolar Wave." The amplitude of interannual fluctuations in annual net ice production are about 40% of the mean value, implying significant interannual variance in brine rejection and upper ocean heat loss. Southward anomalies in wind stress induce negative anomalies in open water production, which are observed in passive microwave satellite images. Thus, cycles of enhanced poleward wind stress reduce ice growth by compacting the ice along the coastline and closing open water in leads and polynyas. Model simulations confirm that years of low ice production, such as 1991, coincide with years of lower than normal bottom water outflow. Future plans include the assimilation of satellite ice concentrations and ice drift dynamics to more accurately constrain boundary conditions in the model.

Drinkwater, Mark R.↗

Sea Ice and Cloud Processes Mediating Compensation between Atmospheric and Oceanic Meridional Heat Transports across the CMIP6 Preindustrial Control Experiment

Abstract Bjerknes compensation (BJC) refers to the anticorrelation observed between atmospheric and oceanic heat transport (AHT/OHT) variability, particularly on decadal to longer time scales that may be important to the predictability of the climate system. This study investigates the spread in BJC across fully coupled simulations of phase 6 of the Coupled Model Intercomparison Project (CMIP6) and critical processes (particularly related to sea ice and clouds) that may contribute to that spread. BJC on decadal to longer time scales is confirmed across all the simulations evaluated, and it is strongest in the Northern Hemisphere (NH) between 60° and 70°N. At these latitudes, BJC appears to be primarily driven by the exchange of turbulent fluxes (sensible and latent) in the Greenland, Iceland, and Barents Seas. Metrics to break down how sea ice and clouds uniquely modify the radiative balance of the polar atmosphere during anomalous OHT events are presented. These metrics quantify the impacts of sea ice and clouds on surface and top of atmosphere (latent, sensible, longwave, and shortwave radiative) energy fluxes. Cloud responses tend to counter the clear sky impacts over the Marginal Ice Zone (MIZ). It is further shown that the degree of BJC present in a simulation at high latitudes is heavily influenced by the sensitivity of the sea ice to OHT, which is most influential over the MIZ. These results are qualitatively robust across models and explain the intermodel spread in NH BJC in the preindustrial control experiment.

54 ENVIRONMENTAL SCIENCES↗

Advancing Sea Ice Predictability in E3SM with Machine Learning

Focal area(s): To improve predictions of sea ice in E3SM we propose to develop a hierarchy of data-driven models using observational and simulation data to investigate the most important Earth system drivers of sea ice variability and loss, develop surrogates that build on the reduced parameter space of important drivers, and, where appropriate, couple machine learning models with standard PDE models to capture important physical behavior at different scales. This work falls under Focal Area 2. Predictive modeling through the use of AI techniques.

54 ENVIRONMENTAL SCIENCES↗

Earth System Model Parameter Adjustment Using a Green's Functions Approach

We demonstrate the practicality and effectiveness of using a Green's functions estimation approach for adjusting uncertain parameters in an Earth system model (ESM). This estimation approach has previously been applied to an intermediate-complexity climate model and to individual ESM components, e.g., ocean, sea ice, or carbon cycle components. Here, the Green's functions approach is applied to a state-of-the-art ESM that comprises a global atmosphere/land configuration of the Goddard Earth Observing System (GEOS) coupled to an ocean and sea ice configuration of the Massachusetts Institute of Technology general circulation model (MITgcm). Horizontal grid spacing is approximately 110 km for GEOS and 37–110 km for MITgcm. In addition to the reference GEOS-MITgcm simulation, we carried out a series of model sensitivity experiments, in which 20 uncertain parameters are perturbed. These “control” parameters can be used to adjust sea ice, microphysics, turbulence, radiation, and surface schemes in the coupled simulation. We defined eight observational targets: sea ice fraction, net surface shortwave radiation, downward longwave radiation, near-surface temperature, sea surface temperature, sea surface salinity, and ocean temperature and salinity at 300 m. We applied the Green's functions approach to optimize the values of the 20 control parameters so as to minimize a weighted least-squares distance between the model and the eight observational targets. The new experiment with the optimized parameters resulted in a total cost reduction of 9 % relative to a simulation that had already been adjusted using other methods. The optimized experiment attained a balanced cost reduction over most of the observational targets. We also report on results from a set of sensitivity experiments that are not used in the final optimized simulation but helped explore options and guided the optimization process. These experiments include an assessment of sensitivity to the number of control parameters and to the selection of observational targets and weights in the cost function. Based on these sensitivity experiments, we selected a specific definition for the cost function. The sensitivity experiments also revealed a decreasing overall cost as the number of control variables was increased. In summary, we recommend using the Green's functions estimation approach as an additional fine-tuning step in the model development process. The method is not a replacement for modelers' experience in choosing and adjusting sensitive model parameters. Instead, it is an additional practical and effective tool for carrying out final adjustments of uncertain ESM parameters.

Green's Function↗

A simulated Antarctic fast ice ecosystem

A 2D numerical ecosystem model of Antarctic land fast ice is developed to elucidate the primary production with the Antarctic sea ice zone. The physical component employs atmospheric data to simulate congelation ice growth, initial brine entrapment, desalination, and nutrient flux. The biological component is based on the concept of a maximum temperature-dependent algal growth rate which is reduced by limitations imposed from insufficient light or nutrients, as well as suboptimal salinity. Preliminary simulations indicate that, during a bloom, microalgae are able to maintain their vertical position relative to the lower congelation ice margin and are not incorporated into the crystal matrix as the ice sheet thickens. It is inferred that land fast sea ice contains numerous microhabitats that are functionally distinct based upon the unique set of processes that control microalgal growth and accumulation within each.

Arrigo, Kevin R.↗

Modelling and simulation of brinicle formation

Below the Arctic sea ice, under the right conditions, a flux of icy brine flows down into the sea. The icy brine has a much lower fusion point and is denser than normal seawater. As a result, it sinks while freezing everything around it, forming an ice channel called a brinicle (also known as ice stalactite). In this paper, we develop a mathematical model for this phenomenon, assuming cylindrical symmetry. The fluid is considered to be viscous and quasi-stationary. The heat and salt transport are weakly coupled to the fluid motion and are modelled with the corresponding conservation equations, accounting for diffusive and convective effects. Finite-element discretization is employed to solve the coupled system of partial differential equations. We find that the model can capture the general behaviour of the physical system and generate brinicle-like structures while also recovering dendrite composition, which is a physically expected feature aligned with previous experimental results. This represents, to our knowledge, the first complete model proposed that captures the global structure of the physical phenomenon even though it has some discrepancies, such as brine accumulation.

97 MATHEMATICS AND COMPUTING↗