Water Ice Research Articles

The Pre-perihelion Evolution of the Activity of Comet C/2017 K2 (Pan-STARRS) during the Water Ice-line Crossover

Comets, relics from the early solar system, consist of dust and ice. The ice sublimates as comets approach the Sun, ejecting dust from their nuclei seen as activity. Different volatiles sublimate at different Sun–comet distances and eject dust of unique sizes, structures, and compositions. In this study, we present new polarimetric observations of Oort cloud comet C/2017 K2 (Pan-STARRS) in R- and I-filter domains before, during, and after its crossover of the water-ice sublimation regime at phase angles of 15.°9, 10.°5, and 20.°0, respectively. Combining multiband optical imaging data covering a wide range of heliocentric distances (∼14−2.3 au), we aim to characterize the pre-perihelion evolution of cometary activity as well as the properties of its coma dust. Two discontinuous brightening events were observed: at ∼6 au presumably associated with changes in CO-like supervolatile ice activity, and at ∼2.9 au when water ice took over. Particularly, the latter activation is accompanied by changes in coma morphology and color whose trends differ between the inner (∼103 km) and outer (∼104 km) parts of the coma. No polarimetric discontinuities on the comet were observed over the inner coma region, all epochs showing phase-angle and wavelength dependencies compatible with those of active comets observed in similar observing geometry. During this period, the underlying dust continuum overwhelmed Hα emission at around 656.3 nm, suggesting less water ice on the comet’s surface than expected. We discuss K2's coma environment by combining numerical simulations of light scattered by dust and place the observations within the context of the comet’s evolution.

Lunar Exploration Based on Ground-Based Radar: Current Research Progress and Future Prospects

Lunar exploration is of significant importance in the development and utilization of in situ lunar resources, water ice exploration, and astronomical science. In recent years, ground-based radar (GBR) has gained increasing attention in the field of lunar exploration due to its flexibility, low cost, and penetrating capabilities. This paper reviews the scientific research on lunar exploration using GBR, outlining the basic principles of GBR and the progress made in lunar exploration studies. Our paper introduces the fundamental principles of lunar imaging using GBR and systematically reviews studies on lunar surface/subsurface detection, the dielectric properties inversion of the lunar regolith, and polar water ice detection using GBR. In particular, the paper summarizes the current development status of the Chinese GBR and forecasts future development trends in China. This review will enhance the understanding of lunar exploration results using GBR radar, systematically demonstrate the main applications and scientific achievements of GBR in lunar exploration, and provide a reference for GBR radar in future lunar exploration missions.

The Ensemble Tangent Linear Model Applied to a Cloud Physics Parameterization

Abstract An Ensemble Tangent Linear Model (ETLM) is applied to a cloud physics scheme used in the Navy Global Environmental Model (NAVGEM). The ensemble is created using 3-hour forecasts from the Ensemble Transform method used in the NAVGEM data assimilation system. The model states are saved before and after applying the cloud physics parameterization (which includes condensation/evaporation of cloud ice and cloud liquid water and stratiform precipitation), and these states are used to construct linearized model tendencies for temperature, specific humidity, cloud liquid water, and cloud ice water. We examine separately the application of the ETLM to cloud physics components that are explicitly local versus non-local. For the local components, an ETLM is built using a single grid point. ETLMs from 50 to 1000 members are tested, and skillful forecasts can be obtained for both local and non-local physics even with a moderate sized ensemble (e.g., 100 members). At 1000 members, the globally-averaged forecast error reductions (relative to persistence errors) are ∼40% for temperature, water vapor, and cloud liquid water and ∼30% for cloud ice. When initial perturbations are reduced by a factor of 0.1, the error reductions increase to ∼65% for all variables. For physics with non-local components (stratiform precipitation) the covariances that comprise the ETLM are localized with a Schur product matrix using a Gaussian localization shape with tunable length. The optimal lengths increased with ensemble size from ∼2-3 km for 50 members to ∼10 km for 1000 members. ETLMs for “all cloud physics” are also constructed and evaluated.

Towards the perfect soft-boiled chicken eggs: defining cooking conditions, quality criteria, and safety assessments

In recent years, ready-to-eat soft-boiled chicken eggs, with coagulated whites and semi-solid yolks have become popular among Chinese consumers due to their convenience, tender texture, and nutritional benefits. However, the standards for these products are currently inconsistent, and quality evaluation parameters and food safety issues remain unclear. Softness ratio, representing the semi-liquid yolk proportion in a cross-section, was defined through a survey of different brands of ready-to-eat soft-boiled chicken eggs. By accounting for egg weight and cooking time, optimal softness under specific conditions was determined. A quality evaluation method was established based on the softness ratio and shelling score. Finally, the safety of soft-boiled chicken eggs was assessed by measuring yolk center temperature. The optimal softness ratio was 0.46 to 0.64. Optimal cooking times in boiling water (100°C) with a 5:1 water-to-egg ratio were 300 s for 43 to 48 g eggs, 330 seconds for 48 to 53g, 350 s for 53 to 58 g, and 370 s for 58 to 63g. After cooking, eggs were cooled for 6 min in an ice water mixture with a 3:1 water-to-egg ratio. Shelling scores (0-5) depended on the egg white surface exposed post-shelling, peaking at an air cell diameter of 21.55 ± 2.26 mm. Egg weight and shelling score had a correlation of -0.41, while egg white springiness and shelling score had a correlation of 0.86. Ensuring core yolk pasteurization standards was difficult, thus stricter management or whole egg irradiation was suggested for safety.

Requirements for the Development and Operation of a Freeze-Up Ice-Jam Flood Forecasting System

This article provides a comprehensive overview of ice-jam flood forecasting methodologies applicable to rivers during freezing. It emphasizes the importance of understanding river ice processes and fluvial geomorphology for developing a freeze-up ice-jam flood forecasting system. The article showcases a stochastic modelling approach, which involves simulating a deterministic river ice model multiple times with varying parameters and boundary conditions. This approach has been applied to the Exploits River at Badger in Newfoundland, Canada, a river that has experienced several freeze-up ice-jam floods. The forecasting involves two approaches: predicting the extent of the ice cover during river freezing and using an ensemble method to determine backwater flood level elevations. Other examples of current ice-jam flood forecasting systems for the Kokemäenjoki River (Pori, Finland), Saint John River (Edmundston, NB, Canada), and Churchill River (Mud Lake, NL, Canada) that are operational are also presented. The text provides a detailed explanation of the processes involved in river freeze-up and ice-jam formation, as well as the methodologies used for freeze-up ice-jam flood forecasting. Ice-jam flood forecasting systems used for freeze-up were compared to those employed for spring breakup. Spring breakup and freeze-up ice-jam flood forecasting systems differ in their driving factors and methodologies. Spring breakup, driven by snowmelt runoff, typically relies on deterministic and probabilistic approaches to predict peak flows. Freeze-up, driven by cold temperatures, focuses on the complex interactions between atmospheric conditions, river flow, and ice dynamics. Both systems require air temperature forecasts, but snowpack data are more crucial for spring breakup forecasting. To account for uncertainty, both approaches may employ ensemble forecasting techniques, generating multiple forecasts using slightly different initial conditions or model parameters. The objective of this review is to provide an overview of the current state-of-the-art in ice-jam flood forecasting systems and to identify gaps and areas for improvement in existing ice-jam flood forecasting approaches, with a focus on enhancing their accuracy, reliability, and decision-making potential. In conclusion, an effective freeze-up ice-jam flood forecasting system requires real-time data collection and analysis, historical data analysis, ice jam modeling, user interface design, alert systems, and integration with other relevant systems. This combination allows operators to better understand ice jam behavior and make informed decisions about potential risks or mitigation measures to protect people and property along rivers. The key findings of this review are as follows: (i) Ice-jam flood forecasting systems are often based on simple, empirical models that rely heavily on historical data and limited real-time monitoring information. (ii) There is a need for more sophisticated modeling techniques that can better capture the complex interactions between ice cover, water levels, and channel geometry. (iii) Combining data from multiple sources such as satellite imagery, ground-based sensors, numerical models, and machine learning algorithms can significantly improve the accuracy and reliability of ice-jam flood forecasts. (iv) Effective decision-support tools are crucial for integrating ice-jam flood forecasts into emergency response and mitigation strategies.

The LCROSS Impact Crater as Seen by ShadowCam and Mini‐RF: Size, Context, and Excavation of Copernican Volatiles

AbstractThe Lunar CRater Observations and Sensing Satellite (LCROSS) impacted a Centaur rocket stage into a permanently shadowed region (PSR) in Cabeus crater, excavating water ice and other volatiles. We used the Miniature Radio Frequency (Mini‐RF) instrument on the Lunar Reconnaissance Orbiter and the ShadowCam instrument on the Korean Pathfinder Lunar Orbiter to detect the probable 22‐m diameter crater that resulted from the LCROSS impact. The crater formed superposed upon a dense small crater population along a crater ray from a larger pre‐existing crater. From its geologic context, the ice and regolith excavated by LCROSS were likely modified within the last 0.1–0.5 Gyr. An upper limit for the excavated volatiles is ~0.9 Gyr, as the location was not a PSR prior to that time. A young age for the LCROSS‐detected volatiles supports the idea that they were mostly emplaced by an exogenic mechanism, such as from comets or the solar wind.

Late Pleistocene polar bear genomes reveal the timing of allele fixation in key genes associated with Arctic adaptation

The polar bear (Ursus maritimus) occupies a relatively narrow ecological niche, with many traits adapted for cold temperatures, movement across snow, ice and open water, and for consuming highly lipid-dense prey species. The divergence of polar bears from brown bears (Ursus arctos) and their adaptation to their Arctic lifestyle is a well-known example of rapid evolution. Previous research investigating whole genomes uncovered twelve key genes that are highly differentiated between polar and brown bears, show signatures of selection in the polar bear lineage, and are associated with polar bear adaptation to the Arctic environment. Further research suggested fixed derived alleles in these genes arose from selection on both standing variation and de novo mutations in the evolution of polar bears. Here, we reevaluate these findings based on a larger and geographically more representative dataset of 119 polar bears and 135 brown bears, and assess the timing of derived allele fixation in polar bears by incorporating the genomes of two Late Pleistocene individuals (aged 130–100,000 years old and 100–70,000 years old). In contrast with previous results, we found no evidence of derived alleles fixed in present-day polar bears within the key genes arising from de novo mutation. Most derived alleles fixed in present-day polar bears were also fixed in the Late Pleistocene polar bears, suggesting selection occurred prior to 70,000 years ago. However, some derived alleles fixed in present-day polar bears were not fixed in the two Late Pleistocene polar bears, including at sites within APOB, LYST, and TTN. These three genes are associated with cardiovascular function, metabolism, and pigmentation, suggesting selection may have acted on different loci at different times.

Evaluation of the Surface Downward Longwave Radiation Estimation Models over Land Surface

Surface downward longwave radiation (SDLR) is crucial for maintaining the global radiative budget balance. Due to their ease of practicality, SDLR parameterization models are widely used, making their objective evaluation essential. In this study, against comprehensive ground measurements collected from more than 300 globally distributed sites, four SDLR parameterization models, including three popular existing ones and a newly proposed model, were evaluated under clear- and cloudy-sky conditions at hourly (daytime and nighttime) and daily scales, respectively. The validation results indicated that the new model, namely the Peng model, originally proposed for SDLR estimation at the sea surface and applied for the first time to the land surface, outperformed all three existing models in nearly all cases, especially under cloudy-sky conditions. Moreover, the Peng model demonstrated robustness across various land cover types, elevation zones, and seasons. All four SDLR models outperformed the Global Land Surface Satellite product from Advanced Very High-Resolution Radiometer Data (GLASS-AVHRR), ERA5, and CERES_SYN1de-g_Ed4A products. The Peng model achieved the highest accuracy, with validated RMSE values of 13.552 and 14.055 W/m2 and biases of −0.25 and −0.025 W/m2 under clear- and cloudy-sky conditions at daily scale, respectively. Its superior performance can be attributed to the inclusion of two cloud parameters, total column cloud liquid water and ice water, besides the cloud fraction. However, the optimal combination of these three parameters may vary depending on specific cases. In addition, all SDLR models require improvements for wetlands, bare soil, ice-covered surfaces, and high-elevation regions. Overall, the Peng model demonstrates significant potential for widespread use in SDLR estimation for both land and sea surfaces.

Vertical structures of typhoon cloud microphysical and radiative features associated with the precipitation type over the western North Pacific

Intensity of different precipitation types (convective, stratiform and shallow) and associated cloud vertical microphysical and radiative heating features are analyzed considering typhoon development, maturity, and decaying stages over the western North Pacific using the CloudSat Tropical Cyclone and China Meteorological Administration tropical cyclone best-track datasets from 2 June 2006 to 31 December 2015. At all three stages, the convective precipitation intensity, approximately twice that of stratiform precipitation, is the highest and peaks at development stage. The strongest stratiform precipitation occurs at typhoon maturity stage. Shallow precipitation is the weakest throughout the typhoon lifespan. Although the cloud microphysical parameters (radar reflectivity, cloud ice particle number concentration and effective radius) of both convective and stratiform precipitation tend to increase with precipitation intensity, convective precipitation contains more ice water of larger sizes in upper layers than stratiform precipitation. Unlike convective and stratiform precipitation, dominated by cold clouds, shallow precipitation is dominated by warm clouds with weak vertical contrast in the radiative distribution but strong radiation nearby 5 km. Our results show that cloud ice particle number concentration is important not only in precipitation intensity enhancement but also in determining the shortwave radiative heating center vertical location. More and larger ice particles in convective precipitation profiles result in stronger or comparable shortwave radiative heating than those in stratiform precipitation profiles, while the longwave radiative cooling rates in convective and stratiform precipitation profiles exhibit very similar features, likely attributable to similar infrared radiation levels due to comparable temperatures in these profiles.

The diurnal variation of dust and water ice aerosol optical depth at Jezero crater observed by MEDA/TIRS over a full Martian year

The Thermal InfraRed Sensor (TIRS) on the Perseverance rover has provided nearly two full Mars years of systematic monitoring of the total aerosol optical depth above Jezero Crater. These observations span a wide range of timescales, capturing seasonal patterns, diurnal variations, and minute-to-minute fluctuations in aerosol loading. By combining TIRS retrievals with orbital observations, the relative contributions of dust and water ice aerosols can be estimated, revealing their different seasonal and diurnal behaviors. The TIRS record shows distinct periods of dust storm activity, including strong regional storms during the perihelion season as well as short-lived but intense dust events outside the typical dust storm season. Water ice clouds exhibit pronounced seasonal and diurnal variability, with peak activity occurring during the aphelion season but with a presence throughout the year. The diurnal variation of clouds differs significantly between the aphelion and perihelion seasons, with clouds persisting throughout the night during the aphelion season, while largely absent outside of specific periods after sunrise and sunset during the perihelion season. These results provide new insights into the complex behavior of aerosols at Jezero Crater and their connections to atmospheric dynamics and the Martian dust and water cycles.

Thermal modeling of the lunar South Pole: Application to the PROSPECT landing site

Water ice is distributed on the surface and in the subsurface of the Moon, as confirmed by observational data, and predicted by several numerical models. In this respect, the direct search for lunar water is the main objective of the ESA’s PROSPECT package, that aims to analyze a region of interest at the lunar South Pole. PROSPECT, originally on the Russian Luna 27, is now on the CLPS (Commercial Lunar Provider Service) “CP” 22 mission. In this work, we applied our 3-D FEM thermophysical model to investigate the landing site selected for the CP 22 mission, which is centred at −84.496°S, 31.588°E, and located on the Leibnitz Plateau and within an area of high elevation. The purpose of our model is to investigate regions of interest (ROI) on the lunar surface by working with the real topography at the scale of 5 m, by using the DEM (Digital Elevation Model) of the region. Since the lunar surface is characterized by topographic variations such as craters or boulders, a 3-D model is preferable over a 1-D numerical model. We produced temperature maps of the surface and 1-D temperature vs depth, as well as we produced illumination maps, computing also the indirect contribution. These simulations will provide a complete thermophysical vision of the landing site, offering a theoretical support to the researchers and engineers of the CP 22 mission, and of future lunar missions. In addition, this model can be applied to every site of the Moon surface and subsurface and, in general, to any airless body of the Solar System.

Coastal Sea Ice Concentration Derived from Marine Radar Images: A Case Study from Utqiaġvik, Alaska

We apply the Canny edge algorithm to imagery from the Utqiaġvik coastal sea ice radar system (CSIRS) to identify regions of open water and sea ice and quantify ice concentration. The radar-derived sea ice concentration (SIC) is compared against the (closest to the radar field of view) 25 km resolution NSIDC Climate Data Record (CDR) and the 1 km merged MODIS-AMSR2 sea ice concentrations within the ∼11 km field of view for the year 2022–2023, when improved image contrast was first implemented. The algorithm was first optimized using sea ice concentration from 14 different images and 10 ice analysts (140 analyses in total) covering a range of ice conditions with landfast ice, drifting ice, and open water. The algorithm is also validated quantitatively against high-resolution MODIS-Terra in the visible range. Results show a correlation coefficient and mean bias error between the optimized algorithm, the CDR and MODIS-AMSR2 daily SIC of 0.18 and 0.54, and ∼−1.0 and 0.7%, respectively, with an averaged inter-analyst error of ±3%. In general, the CDR captures the melt period correctly and overestimates the SIC during the winter and freeze-up period, while the merged MODIS-AMSR2 better captures the punctual break-out events in winter, including those during the freeze-up events (reduction in SIC). Remnant issues with the detection algorithm include the false detection of sea ice in the presence of fog or precipitation (up to 20%), quantified from the summer reconstruction with known open water conditions. The proposed technique allows for the derivation of the SIC from CSIRS data at spatial and temporal scales that coincide with those at which coastal communities members interact with sea ice. Moreover, by measuring the SIC in nearshore waters adjacent to the shoreline, we can quantify the effect of land contamination that detracts from the usefulness of satellite-derived SIC for coastal communities.

Potential cloud precipitation capacity in typical regions over China

Based on the CLoud, Albedo and RAdiation dataset, AVHRR-based, version 2 (CLARA-A2), Tropical Rainfall Measuring Mission 3B43 (TRMM-3B43), and European Centre for Medium-Range Weather Forecasts Reanalysis v5 (ERA5) reanalysis data, the potential cloud precipitation capacity (PCPA) of typical regions in China is compared, and the relationship between impact factors and PCPA is discussed. Results have suggested that the Tarim Basin (TB) has scarce cloud water resources, while cloud water path (CWP) values are higher in South China (SC) and Sichuan Basin (SB) under the influence of the East Asian monsoon. Moreover, different typical regions of China exhibit varying dependencies on the ice water path (IWP) and liquid water path (LWP). There is a strong correlation between the IWP and precipitation in the Tibet Plateau (TP), Northeast China (NE), SC, and SB. The precipitation in TB demonstrates a more pronounced correlation with the LWP. Through a comparison of the correlation between PCPA and influencing factors in different typical regions of China, it is found that convective available potential energy (CAPE), surface latent heat flux (SLHF), surface sensible heat flux (SSHF), and 0–3 km relative humidity (RH) exhibit stronger correlation with PCPA than 2 m temperature (T2m) and 2–5 km vertical wind shear (SHEAR). Further investigation revealed that the joint effect of CAPE, RH, and SLHF has a pronounced effect on PCPA, particularly during spring and autumn. Additionally, the PCPA of TP exhibits significant dependency on the joint effect of these three influential factors. Furthermore, the ratio of LWP to IWP (RLI) also affects PCPA. In spring and autumn, the PCPA of TB and NC exhibits a positive correlation with RLI, whereas the PCPA of TP, SC, NE, and SB shows a negative correlation with RLI. In summer, the PCPA of TB and SC exhibits a notably negative correlation with RLI. This study deepens the understanding of the formation mechanism of cloud precipitation in typical regions of China, provides the basis for climate forecast and improves the accuracy of weather forecast.

Interstellar Glycolaldehyde, Methyl Formate, and Acetic Acid. II. Chemical Modeling of the Bimodal Abundance Pattern in NGC 6334I

Gas-phase abundance ratios between C2H4O2 isomers methyl formate (MF), glycolaldehyde (GA), and acetic acid (AA) are typically on the order of 100:10:1 in star-forming regions. However, an unexplained divergence from this neat relationship was recently observed toward a collection of sources in the massive protocluster NGC 6334I; some sources exhibited extreme MF:GA ratios, producing a bimodal behavior between different sources, while the MF:AA ratio remained stable. Here, we use a three-phase gas-grain hot-core chemical model to study the effects of a large parameter space on the simulated C2H4O2 abundances. A combination of high gas densities and long timescales during ice-mantle desorption (∼125–160 K) appears to be the physical cause of the high MF:GA ratios. The main chemical mechanism for GA destruction occurring under these conditions is the rapid adsorption and reaction of atomic H with GA on the ice surfaces before it has time to desorb. The different binding energies of MF and GA on water ice are crucial to the selectivity of the surface destruction mechanism; individual MF molecules rapidly escape the surface when exposed by water loss, while GA lingers and is destroyed by H. Moderately elevated cosmic-ray ionization rates can increase absolute levels of “complex organic molecule” (COM) production in the ices and increase the MF:GA ratio, but extreme values are destructive for gas-phase COMs. We speculate that the high densities required for extreme MF:GA ratios could be evidence of COM emission dominated by COMs desorbing within a circumstellar disk.

Competitive Entrapment of Hypervolatiles in Interstellar and Cometary Water Ice Analogs

The distribution of chemical species in protoplanetary disks around young stars, especially their division between gas and solid phases, fundamentally shapes the composition of future planets and planetesimals. This distribution is likely affected by entrapment, a mechanism whereby volatile species are mechanically or chemically bound within a less volatile ice. In this study, we investigate the entrapment efficiencies of four hypervolatiles (CH4, CO, N2, and Ar) in multicomponent water ice mixtures deposited at different temperatures and mixture ratios. At low ice deposition temperatures, we observe small differences in entrapment efficiency (CH4>CO>N2∼Ar) up to a factor of two across species. The differences in entrapment between species increase by up to an order of magnitude with increasing deposition temperature. The relative entrapment efficiencies are also impacted by changes in the overall hypervolatile concentration of the ice mixtures. Collectively, these experiments suggest that relative entrapment efficiencies are mainly regulated by small differences in binding energies to the ice matrix, though competition for the best sites also influences entrapment in more concentrated ices. We use these results to better inform interpretations of hypervolatile observations in comets and related objects.

Sufficiency of near-surface water ice as a driver of dust activity on comets

Context. Nearly all contemporary theoretical research on cometary dust activity relies on models depicting heat transfer and sublimation products within the near-surface porous layer. Gas flow exerts a pressure drag to the crust agglomerates, counteracting weak gravity and the tensile strength of that layer. Our interpretation of data from the Rosetta mission, and our broader comprehension of cometary activity, hinges significantly on the study of this process. Aims. We investigate the role played by the structure of the near-surface porous layer and its associated resistance to gas flow, tensile strength, pressure distribution, and other characteristics in the scenario of the potential release of dust agglomerates and the resulting dust activity. Methods. We employ a thermophysical model that factors in the microstructure of this layer and radiative heat conductivity. We consider gas flow in both the Knudsen and transition regimes. To accomplish this, we use methods such as test-particles Monte Carlo, direct-simulation Monte Carlo, and transmission probability. Our study encompasses a broad spectrum of dust-particle sizes. Results. We evaluated the permeability of a dust layer composed of porous aggregates in the submillimetre and millimetre ranges. We carried out comparisons among various models that describe gas diffusion in a porous dust layer. For both the transition and Knudsen regimes, we obtained pressure profiles within a non-isothermal layer. We discuss how the gaps in our understanding of the structure and composition could impact tensile strength estimates. We demonstrate that for particles in the millimetre range, the lifting force of the sublimation products of water ice is adequate to remove the layer. This scenario remains feasible even for particles on the scale of hundreds of microns. This finding is crucial as the sublimation of water ice continues to be the most probable mechanism for dust removal. Conclusions. This study partially overturns the previously held, pessimistic view regarding the possibility of dust removal via water sublimation. We demonstrate that a more precise consideration of various physical processes allows elevation of the matter of dust activity to a practical plane, necessitating a fresh quantitative analysis.

Turbulent heat fluxes in the North Water Polynya and ice estimated based on ASRv2 data and their impact on cloud

The presence or absence of sea ice introduces a substantial perturbation to surface‒atmosphere energy exchanges. Comprehending the effect of varying sea ice cover on surface‒atmosphere interactions is an important consideration for understanding the Arctic climate system. The recurring North Water Polynya (NOW) serves as a natural laboratory for isolating cloud responses to a rapid, near-step perturbation in sea ice. In this study, we employed high-resolution Arctic System Reanalysis version 2 (ASRv2) data to estimate turbulent heat fluxes over the NOW and nearby sea ice (NSI) area between 2005/2006 and 2015/2016. The results indicate that the average turbulent heat fluxes in the polynya are about 87% and 86% higher than in the NSI area over the 10 years during the entire duration of the polynya and during polar night, respectively. Enhanced turbulent heat fluxes from the polynya tend to produce more low-level clouds. The relationship between the polynya and low cloud in winter was examined based on Cloud‒Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO). The low-cloud fraction (0–2 km) was about 7–34% larger over the polynya than the NSI area, and the ice water content below 200 m was about 250%–413% higher over the former than the latter. The correlation between cloud fraction and turbulent heat fluxes in the polynya peaks around the altitude of 200–300 m. These results suggest that the NOW affects the Arctic boundary layer cloudiness and structure in wintertime. Furthermore, higher horizontal resolution reanalysis data can advance our understanding of the cloud-polynya response.

Understanding the Impact of H2 Diffusion Energy on the Formation Efficiency of H2 on the Interstellar Dust Grain Surface

We use microscopic Monte Carlo simulation techniques to investigate the impact of H2 diffusion energy on the recombination efficiency of H2 on interstellar dust grain surfaces under diffuse and translucent cloud conditions. We constructed five models representing different possible conditions encountered by adsorbed H and H2 on interstellar dust grains. We implemented adsorption sites with multiple binding energies for surface species; the Encounter-Desorption mechanism was also included. The study focused on silicate surfaces in diffuse clouds and water ice surfaces in translucent clouds. The results show that the recombination efficiency of H2 on dust surfaces decreases as H2 diffusion energy increases. An interesting finding of this work is that considering different binding sites for H and H2 gives rise to multiple steady phases, during which the recombination efficiency remains constant with a change in H2 diffusion energy.

Numerical modelling and measurement of the E-I characteristics of MgB2 wire in sub-cooled water ice

This work presents a comprehensive 3D numerical model of MgB2 multi-filamentary superconducting wires using the Finite Element Method (FEM) software, COMSOL Multiphysics® 6.0. The study aims to investigate the electro-thermal behavior of MgB2 composite wires during standard transport measurements at various initial temperatures under subcooled water ice conditions. By solving a series of partial differential equations governing heat transfer and dynamic current transport, the model provides detailed insights into the wire’s performance. The simulation results are rigorously compared with experimental E-I characteristics measured for 6-filament MgB2 wires with internal copper stabilization. This comparison validates the model and highlights its capability to predict the behavior of superconducting wires under cryogenic conditions. The findings offer valuable data on the current distribution, ohmic losses, and overall thermal stability of the composite wires, contributing to the advancement of cryogen-free superconducting technologies. This study bridges the gap in the literature regarding the electrothermal dynamics of MgB2 wires cooled by subcooled water ice, providing a foundation for further research and practical applications in high-field generation devices.

Mathematical Modeling of Reverse Processes with Consideration of the Water–Ice Phase Transition

Mathematical Modeling of Reverse Processes with Consideration of the Water–Ice Phase Transition