Ice-water Interface Research Articles

Ion migration at the ice–water interface during the freezing process of salt solution – Experimental investigation

Within the scope of this work, systematic investigations of the ion migration process at the ice–water interface during the freezing of salt solutions are performed via experimental investigations (part I) and numerical investigations (part II). In this study, the ion migration trend under different freezing temperatures was revealed, and the formation process of the brine pockets during phase transition and the mechanism of ion release in ice were explored. The results show that the migration capabilities of Na+ and K+ are greater than those of Ca2+ and Mg2+, whereas the migration capability of Cl− is greater than those of HCO3− and SO42−. At −2 °C, the ion concentration in the surface ice decreased by 47.64 %, indicating the occurrence of ion release within the ice. When the moving speed of the phase transition interface is greater than the ion diffusion ability, the brine pockets form in the ice. Under a temperature gradient, the brine pockets in ice become thermodynamically unstable, resulting in the inevitable ion release toward the high-temperature side. The ion release depends on the temperature gradient and the brine pocket concentration. Furthermore, the unsynchronous migration of ions is related to the freezing temperature, diffusion coefficient, hydration radius, and hydration free energy.

Reconsidering freeze-induced protein aggregation: Air bubbles as the root cause of ice-water interface stress

Freeze-induced stress causing aggregation of proteins has typically been primarily attributed to the ice-water interface. However, we hypothesize that the underlying observed and perceived detrimental effect of ice is, to some extent, attributed to air bubbles expelled from ice crystal lattices or to nanobubbles existing prior to freezing. The reduction of dissolved air was achieved via a deaeration process by placing samples in a reduced pressure chamber, while the reduction of nanobubbles was achieved by filtering samples via a syringe filter. The results showed that the reduction of both dissolved air molecules and stable colloidal nanobubbles in a bovine IgG solution prior to freezing led to a significant decrease in aggregation after thawing compared to untreated samples (∼6,000 vs. ∼ 40,000 particles/mL at a freezing rate of 100 K/s, respectively). The deaeration-filtration treatment works additively with cryoprotectants such as trehalose, further reducing the freeze-induced aggregation of IgG. The results also demonstrated that air–water interfacial aggregation of IgG in bulk liquid samples is a time-dependent process. The number of IgG subvisible particles increased with time and temperature, suggesting that random collisions of denatured molecules promoted the formation of aggregates with spherical morphology. In contrast, the IgG subvisible count after freeze-thawing had already reached its nominal value, suggesting a time-independent process where denatured protein molecules were compressed between ice crystals into filament-like aggregates. In summary, the findings shift the current paradigm from ice crystals being the main destabilizing factor during freezing to air bubbles, although the two are intertwined. From a translational aspect, this study underscores the value of deaeration-filtration as an essential supplemental process that can be applied in addition to formulation approaches such as the use of cryoprotectants to further reduce freezing stress on proteins and increase their stability.

Meso-scale investigation on the permeability of frozen soils with the lattice Boltzmann method

Complex composition and intricate pore-scale structure of frozen soils poses significant challenges in reliably and efficiently obtaining their permeability. In this study, we propose a modified quartet structure generation set (QSGS) numerical tool for generating frozen soils and present the development of a computational simulation code based on the multiple-relaxation-time lattice Boltzmann method (LBM). In the modified QSGS, the arc-shaped water-ice interface is depicted, and the influence of pore-scale geometry on freezing temperature is considered. The validity of combining the proposed QSGS model and the LBM code is proved by comparing calculated results to analytical and experimental results of porous media. Our objective was to investigate the effects of soil features, including porosity, grain diameter, shape anisotropy of soil particles, and ice content on the intrinsic permeability of frozen soil. Additionally, we examined the relationship between these features and the specific surface area and tortuosity. Numerical results show that the intrinsic permeability of frozen soils increases with increasing porosity, larger granular diameter, and anisotropy, which is identical with the pressure gradient. The presence of ice led to clogging flow pathways and drastically decreased the intrinsic permeability, which is significantly less than unfrozen soil with same effective porosity. This study provides a useful tool to investigate the intricate interplay between the pore-scale structure and the intrinsic permeability of frozen soils.

Underwater Acoustic Scattering from Multiple Ice Balls at the Ice–Water Interface

Underwater Acoustic Scattering from Multiple Ice Balls at the Ice–Water Interface

Stable isotopes of landfast sea ice as a record of La Grande River under-ice plume dispersal

Vertical profiles of salinity, and isotopic abundance ratios of hydrogen (δ2H) and oxygen (δ18O) of 18 landfast ice cores, collected along the northeast coast of James Bay in March 2019, and one ice core collected in Belcher Islands, were used to obtain the winter timeseries of the spatiotemporal evolution of the under-ice plume of La Grande River (LGR), the dominant river in the area. Variability in the isotopic composition and salinity of the ice cores indicated changes to the water source composition at the ice–water interface when the ice layers formed. The increased presence of river water beneath the ice during January–March was marked by more negative isotopic ratios in the lower portion of the ice as river discharge was increased for hydroelectric production. River water was the source of ∼43% of the ice in our ice core samples ( n = 320) with the interquartile range from 16% to 71%. The river water fractions incorporated into the ice indicate that LGR under-ice plume extended more than 75 km north and at least 30 km south of the river mouth for ∼3 months. These findings correspond well with more challenging to obtain hydrographic observations. End-of-winter ice core sampling and analysis for isotopic abundance has potential as a tool to monitor dispersal of LGR discharge into the ice-covered coastal environment.

Electric Field-Enhanced Ion Rejection Rate in Freeze Desalination.

Freeze desalination is an appealing method for seawater desalination through freezing seawater. The percentage of ions in the liquid phase, which is termed ion rejection rate, is a critical factor affecting the performance of freeze desalination. Improving the ion rejection rate is an important topic for freeze desalination. In this work, we investigate the effects of electric fields on the ion rejection rate during the freezing of seawater through molecular dynamics simulations. It is found that the ion rejection rate increases with increasing electric field strength. The enhanced ion rejection rate is due to the reduction of the energy barrier at the ice-water interface caused by the electric field, which affects the orientation of water molecules and ion-water interactions. However, the electric field hinders the ice growth rate, which affects the productivity of freeze desalination. Nevertheless, the finding in this work offers a new idea to improve the ion rejection rate. Practically, a trade-off needs to be found to optimize the overall performance of freeze desalination.

Molecular-Resolution Electron Imaging of Defects and Dynamics at the Ice-Water Interface

Molecular-Resolution Electron Imaging of Defects and Dynamics at the Ice-Water Interface

Rough ice cover modified hyporheic exchange: A numerical study on the mechanical effect propagating from the ice-water interface to the sediment–water interface

Many rivers are subject to ice-covered flow conditions in cold climate regions or during the winter months. The presence of ice cover can change stream flow hydrodynamics and thus affect the hyporheic exchange involved in many biogeochemical processes in aquatic environments. Among the various features of ice cover, its underside roughness is a main factor modifying the surface–subsurface flow hydrodynamics. To date, there has been little research on how ice cover affects the hyporheic exchange, even less on the effect of the ice cover roughness. This leaves open questions on the mechanical changes initiated at the ice-water interface and transferred to the water–sediment interface. In this study, a coupled surface–subsurface computational approach was used to model conjunctively the channel flow over dune bed and the underlying porous flow under different ice cover conditions. The results show that ice cover can enhance the hyporheic flux and reduce the hyporheic zone depth when compared to open channel flow at equivalent discharge. As the river flow depth increases to 4 times the height of the dune, the effect of the ice cover on the hyporheic flow becomes negligible. Within this effective range, the hyporheic flux via the channel water depth follows a power-law function. Ice cover roughness greatly augments the exchange rate while compresses the exchange space, therefore reduces flow residence time in the hyporheic zone. Based on the modeling results, prediction models were proposed for evaluating the impacts of ice cover roughness on the hyporheic exchange flux and depth. Our study indicates that hyporheic flow exchange is significantly impacted by the roughness of the ice cover which in turn the hyporheic exchange is likely to affect the river ice processes as well as the water quality and ecological functions throughout the river corridors.

Eutrophication in cold-arid lakes: molecular characteristics and transformation mechanism of DOM under microbial action at the ice-water interface

During freezing periods, nutrients (carbon and organic matter, etc.) are enriched in the water and sediment of lakes in cold-arid regions, leading to potential algal bloom outbreaks and other health risks to the ecosystem. Particularly, dissolved organic matter (DOM) is a critical component of the nutrients and plays an important role in the global carbon cycle. However, the mechanisms of DOM transfer between ice and water remain unclear. This study analyzed the influence of microbial community on DOM composition using 16 s RNA, 3DEEM, and FT-ICR MS in Daihai Lake and Wuliangsuhai Lake in the Yellow River Basin, China. According to the spectral analysis, the content of endogenous organic matter, such as humus, accounted for 40% of the total DOM in water, while the content of tryptophan and tyrosine accounted for 80% of the total DOM in ice. The results of mass spectrometry showed that lignin was the main component, and the content of organic matter in the ice was less than that in the water. Molecular structures of seven DOM coexisting in the lake ice and water were elucidated with adapted Kendrick-analogous network visualization, which clearly illustrates that long-chain DOM molecules are derived from small molecules, while other heteroatoms are complexed with the side groups. The positive correlations between CHO, CHNO, CHOS, CHOS and Actinomyces indicate that DOM actively interacted with the microbial community. 44% of CHO compounds have the same molecular formula in water, the content of CHOS in the water of the two lakes was closed to 7% higher than that in the ice. Meanwhile, DOM dynamically migrate between ice and water via interstitial water because of the solubility changes under microbial transformation, which has been proved by the decrease in the contents of the humus and tryptophan-like substances in the ice from the bottom to the surface and lower contents of carbohydrate and unsaturated aromatic hydrocarbon in the water than the ice. This study helps to predict the composition and structure of DOM during the migration in lakes and provides a scientific basis for environmental remediation with high concentration of carbon.

Bioenergetics of iron snow fueling life on Europa

The main sources of redox gradients supporting high-productivity life in the Europan and other icy ocean world oceans were proposed to be photolytically derived oxidants, such as reactive oxygen species (ROS) from the icy shell, and reductants (Fe(II), S(-II), CH4, H2) from bottom waters reacting with a (ultra)mafic seafloor. Important roadblocks to maintaining life, however, are that the degree of ocean mixing to combine redox species is unknown, and ROS damage biomolecules. Here, we envisage a unique solution using an acid mine drainage (AMD)-filled pit lakes analog system for the Europan ocean, which previous models predicted to be acidic. We hypothesize that surface-generated ROS oxidize dissolved Fe(II) resulting in Fe(III) (hydr)oxide precipitates, that settle to the seafloor as "iron snow." The iron snow provides a respiratory substrate for anaerobic microorganisms ("breathing iron"), and limits harmful ROS exposure since they are now neutralized at the ice-water interface. Based on this scenario, we calculated Gibbs energies and maximal biomass productivities of various anaerobic metabolisms for a range of pH, temperatures, and H2 fluxes. Productivity by iron reducers was greater for most environmental conditions considered, whereas sulfate reducers and methanogens were more favored at high pH. Participation of Fe in the metabolic redox processes is largely neglected in most models of Europan biogeochemistry. Our model overcomes important conceptual roadblocks to life in icy ocean worlds and broadens the potential metabolic diversity, thus increasing total primary productivity, the diversity and volume of habitable environmental niches and, ultimately, the probability of biosignature detection.

Acoustic observations of individual bubble release events from melting glacier ice in an arctic fjord

Melting below the waterline of marine-terminating glaciers and ice sheets has the potential to remove ice that is holding back land-based ice upstream, and increase the rate of global mean sea level rise. It is difficult to observe this melting directly, and the physical processes involved remain an active area of study. As glacier ice melts, the air bubbles trapped in the ice are released into the water. The air bubbles are often at greater pressure than the surrounding hydrostatic pressure in the water, causing them to be released explosively and generate sound through resulting monopole oscillations. Studying this sound at the level of individual bubbles provides an opportunity to learn more about the role of the bubbles in affecting heat flux across the ice-water boundary layer. It is also essential for proposed efforts to use acoustic emissions as a tool to quantify glacier-scale submarine melting. We will present and discuss observations made at the ice-water interface of floating pieces of glacier ice in an Arctic fjord. These observations included measurements of the water thermal structure, velocity, and salinity, as well as ablation of the ice face, made concurrently with acoustic measurements of bubble release pulses from a 2-element hydrophone array.

Forward modelling of synthetic-aperture radar (SAR) backscatter during lake ice melt conditions using the Snow Microwave Radiative Transfer (SMRT) model

Abstract. Monitoring of lake ice is important to maintain transportation routes, but in recent decades the number of in situ observations have declined. Remote sensing has worked to fill this gap in observations, with active microwave sensors, particularly synthetic-aperture radar (SAR), being a crucial technology. However, the impact of wet conditions on radar and how interactions change under these conditions have been largely ignored. It is important to understand these interactions as warming conditions are likely to lead to an increase in the occurrence of slush layers. This study works to address this gap using the Snow Microwave Radiative Transfer (SMRT) model to conduct forward-modelling experiments of backscatter for Lake Oulujärvi in Finland. Experiments were conducted under dry conditions, under moderate wet conditions, and under saturated conditions. These experiments reflected field observations during the 2020–2021 ice season. Results of the dry-snow experiments support the dominance of surface scattering from the ice–water interface. However, conditions where layers of wet snow are introduced show that the primary scattering interface changes depending on the location of the wet layer. The addition of a saturated layer at the ice surface results in the highest backscatter values due to the larger dielectric contrast created between the overlying dry snow and the slush layer. Improving the representation of these conditions in SMRT can also aid in more accurate retrievals of lake ice properties such as roughness, which is key for inversion modelling of other properties such as ice thickness.

Wavy Ice Patterns as a Result of Morphological Instability of an Ice–Water Interface with Allowance for the Convective–Conductive Heat Transfer Mechanism

In this research, the wavy ice patterns that form due to the evolution of morphological perturbations on the water–ice phase transition interface in the presence of a fluid flow are studied. The mathematical model of heat transport from a relatively warm fluid to a cold wall includes the mechanism of convective–conductive heat transfer in liquid and small sinusoidal perturbations of the water–ice interface. The analytical solutions describing the main state with a flat phase interface as well as its small morphological perturbations are derived. Namely, the migration velocity of perturbations and the dispersion relation are found. We show that the amplification rate of morphological perturbations changes its sign with variation of the wavenumber. This confirms the existence of two different crystallization regimes with (i) a stable (flat) interfacial boundary and (ii) a wavy interfacial boundary. The maximum of the amplification rate representing the most dangerous (quickly growing) perturbations is found. The theory is in agreement with experimental data.

Ice Recrystallization Unveils the Binding Mechanism Operating at a Diffused Interface.

Recrystallization of ice is a natural phenomenon that causes adverse effects in cryopreservation, agriculture, and in frozen food industry. It has long been recognized that ice recrystallization occurs through the Ostwald ripening and accretion processes. However, neither of these processes has been explored in microscopic detail by state-of-the-art experimental techniques. We carried out atomistic molecular dynamics (MD) simulations to explore ice recrystallization through the accretion process. Attempts have been made to elucidate the binding mechanism that is operating at the diffused ice-water interface. It is demonstrated that two ice crystals spontaneously recognize each other and bind together to form a large crystal in liquid water, resulting in ice recrystallization by accretion. Interestingly, the study reveals that the binding occurs due to the freezing of the interfacial water layer present between the two ice planes, even at a temperature above the melting point of the ice crystal. The synergistically enhanced ordering effect of two ice surfaces on the interfacial water leads to such freezing occurring during the binding process. However, proper crystallographic alignment is not necessarily required for the binding of the two crystals. Simulations have also been carried out to study the binding between an ice crystal and the model ice-binding surface (IBS) of an antifreeze protein above the melting point of the ice crystal. It is found that such binding at the IBS is accompanied by freezing of the interfacial water. This establishes that the synergetic ordering-driven freezing of interfacial water is a common binding mechanism at the diffused surfaces of ice crystals. We believe that this mechanism will provide a microscopic understanding of the process of recrystallization inhibition and thus help in designing suitable materials for potent applications in recrystallization inhibition.

An effect of a snow cover on solar heating and melting of lake or sea ice

Solar radiative heating and melting of lake and sea ice is a geophysical problem that has attracted the attention of researchers for many years. This problem is important in connection with the current global change of the climate. Physical and computational models of the process are suggested in the paper. Analytical solutions for the transfer of solar radiation in light-scattering snow cover and ice are combined with numerical calculations of heat transfer in a multilayer system. The thermal boundary conditions take into account convective heat losses to the ambient air and radiative cooling in the mid-infrared window of transparency of the cloudless atmosphere. The study begins with an anomalous spring melting of ice on the large high-mountain lakes of Tibet. It was found that a thick ice layer not covered with snow starts to melt at the ice-water interface due to volumetric solar heating of ice. The results of the calculations are in good agreement with the field observations. The computational analysis showed a dramatic change in the process when the ice is covered with snow. A qualitative change in the physical picture of the process occurs when the snow cover thickness increases to 20–30 cm. In this case, the snow melting precedes ice melting and water ponds are formed on the ice surface. This is typical for the Arctic Sea in polar summer. Known experimental data are used to estimate the melting of sea ice under the melt pond. Positive or negative feedback related to the specific optical and thermal properties of snow, ice, and water are discussed.

Brief communication: A technique for making in situ measurements at the ice–water boundary of small pieces of floating glacier ice

Abstract. This paper presents an apparatus and associated methods for making direct in situ measurements of the ice–water boundary of small pieces of floating glacier ice. The method involves approaching ice pieces in a small boat and attaching a frame with instruments on it to them using ice screws. These types of measurements provide an opportunity to study small-scale processes at the ice–water interface which control heat flux across the boundary. Recent studies have suggested that current parameterizations of these processes may be performing poorly. Improving understanding of these processes may allow for more accurate theoretical and model descriptions of submarine melting.

Biomolecular profiles of Arctic sea-ice diatoms highlight the role of under-ice light in cellular energy allocation.

Arctic sea-ice diatoms fuel polar marine food webs as they emerge from winter darkness into spring. Through their photosynthetic activity they manufacture the nutrients and energy that underpin secondary production. Sea-ice diatom abundance and biomolecular composition vary in space and time. With climate change causing short-term extremes and long-term shifts in environmental conditions, understanding how and in what way diatoms adjust biomolecular stores with environmental perturbation is important to gain insight into future ecosystem energy production and nutrient transfer. Using synchrotron-based Fourier transform infrared microspectroscopy, we examined the biomolecular composition of five dominant sea-ice diatom taxa from landfast ice communities covering a range of under-ice light conditions during spring, in Svalbard, Norway. In all five taxa, we saw a doubling of lipid and fatty acid content when light transmitted to the ice-water interface was >5% but <15% (85%-95% attenuation through snow and ice). We determined a threshold around 15% light transmittance after which biomolecular synthesis plateaued, likely because of photoinhibitory effects, except for Navicula spp., which continued to accumulate lipids. Increasing under-ice light availability led to increased energy allocation towards carbohydrates, but this was secondary to lipid synthesis, whereas protein content remained stable. It is predicted that under-ice light availability will change in the Arctic, increasing because of sea-ice thinning and potentially decreasing with higher snowfall. Our findings show that the nutritional content of sea-ice diatoms is taxon-specific and linked to these changes, highlighting potential implications for future energy and nutrient supply for the polar marine food web.

Ultrafast Axial Freezing in a Liquid-Filled Capillary Tube.

Liquid-filled capillary tubes are a kind of standard component in life science (e.g., blood vessels, interstitial pores, and plant vessels) and engineering (e.g., MEMS microchannel resonators, heat pipe wicks, and water-saturated soils). Under sufficiently low temperatures, the liquid in a capillary tube undergoes phase transition, forming an ice nucleus randomly on its inner wall. However, how an ice layer forms from the nucleus and then expands, either axially or radially to the tube inner wall, remains obscure. We demonstrated, both experimentally and theoretically, that axial freezing along the inner wall of a water-filled capillary tube occurs way ahead of radial freezing, at a nearly constant velocity 3 orders in magnitude faster than the latter. Rapid release of latent heat during axial freezing was identified as the determining factor for the short duration of recalescence, resulting in an exponential rise of the supercooling temperature from ice nucleation via axial freezing to radial freezing. The profile of the ice-water interface is strongly dependent upon the length-to-radius ratio of the capillary tube and the supercooling degree at ice nucleation. The results obtained in this study bridge the knowledge gap between the classical nucleation theory and the Stefan solution of phase transition.

Trajectories of Liquid Particles in a Dark Soliton Field in a Fluid Beneath an Ice Cover

A fluid layer of finite depth described by Euler equations is considered. The ice cover is modeled by a geometrically non-linear elastic Kirchhoff–Love plate. The trajectories of liquid particles under the ice cover are found in the field of a nonlinear surface traveling wave tending to a periodic wave at infinity: the so-called “dark soliton” (which is a nonlinear product of a step-wave and a periodic wave) of small but finite amplitude. The presence of dark solitons in the system is an indicator of the modulation stability of the carrier periodic wave (defocusing). The analysis uses explicit asymptotic expressions for solutions describing wave structures on the water–ice interface such as dark solitons, as well as asymptotic solutions for the velocity field in the liquid column generated by these waves.

Research of thermoelectric null-thermostat

Objective. The purpose of the study is to develop and study the design of a thermoelectric null-thermostat, characterized by low power consumption, in which high accuracy of temperature stabilization of the reference junctions of differential thermocouples is achieved by placing them directly near the water-ice interface.Method. The study was carried out on an experimental model of a thermoelectric null-thermostat. To study heat transfer processes, visual observation of the movement of the ice-water interface in the null-thermostat chamber was carried out. Temperature values were recorded at the lower and upper bases of the chamber, as well as at the hot and cold junctions of the thermoelectric module (TEM) using thermocouples.Result. The dependences of the change in temperature of the upper and lower base of the working chamber, as well as the phase boundary over time, the duration of time for complete ice penetration on the TEM supply current, and the temperature difference, respectively, between the surface of the cold junction of the TEM and the lower base of the working chamber were obtained , the surface of the hot junction of the thermal module and the upper base of the working chamber.Conclusion. It has been established that the speed of movement of the phase boundary strongly depends on the value of the supply current of the thermoelectric battery, which corresponds to the heat flow at the upper and lower base of the working chamber. As follows from the data presented, the change in temperature of the upper base is more noticeable than the lower one. At a current of 2, 4 and 6 A, the speed of movement of the phase boundary is 0.007 m/h, 0.01 m/h and 0.013 m/h, respectively. In this case, the duration of complete melting of ice, corresponding to the duration of operation of the null thermostat, when the TEM supply current changes from 2 to 7 A, is reduced from 1,91⋅104 s to 1,38⋅104.