Ice Core Research Articles

Holocene climate variability in Slovenia: A review

The Slovenian climate has undergone significant fluctuations, and an understanding of the past climate is necessary to improve models and recognise long-term patterns. The cryosphere environment, such as ice core samples, provides valuable palaeoclimate data. Palynology and dendroclimatology are also effective ways to study long-term changes in vegetation and reconstruct past climates using pollen and tree proxies. Sediment cores from various locations in Slovenia have been studied to understand past environmental changes. Borehole temperature profiles as well as historical records were also used to reconstruct past climate conditions. Studies have shown specific periods when climatic changes likely played a major role, but a complete timeline of the Slovenian climate throughout the Holocene has not yet been fully developed.

Giant viral signatures on the Greenland ice sheet.

Dark pigmented snow and glacier ice algae on glaciers and ice sheets contribute to accelerating melt. The biological controls on these algae, particularly the role of viruses, remain poorly understood. Giant viruses, classified under the nucleocytoplasmic large DNA viruses (NCLDV) supergroup (phylum Nucleocytoviricota), are diverse and globally distributed. NCLDVs are known to infect eukaryotic cells in marine and freshwater environments, providing a biological control on the algal population in these ecosystems. However, there is very limited information on the diversity and ecosystem function of NCLDVs in terrestrial icy habitats. In this study, we investigate for the first time giant viruses and their host connections on ice and snow habitats, such as cryoconite, dark ice, ice core, red and green snow, and genomic assemblies of five cultivated Chlorophyta snow algae. Giant virus marker genes were present in almost all samples; the highest abundances were recovered from red snow and the snow algae genomic assemblies, followed by green snow and dark ice. The variety of active algae and protists in these GrIS habitats containing NCLDV marker genes suggests that infection can occur on a range of eukaryotic hosts. Metagenomic data from red and green snow contained evidence of giant virus metagenome-assembled genomes from the orders Imitervirales, Asfuvirales, and Algavirales. Our study highlights NCLDV family signatures in snow and ice samples from the Greenland ice sheet. Giant virus metagenome-assembled genomes (GVMAGs) were found in red snow samples, and related NCLDV marker genes were identified for the first time in snow algal culture genomic assemblies; implying a relationship between the NCLDVs and snow algae. Metatranscriptomic viral genes also aligned with metagenomic sequences, suggesting that NCLDVs are an active component of the microbial community and are potential "top-down" controls of the eukaryotic algal and protistan members. This study reveals the unprecedented presence of a diverse community of NCLDVs in a variety of glacial habitats dominated by algae.

Mapping the vertical heterogeneity of Greenland's firn from 2011–2019 using airborne radar and laser altimetry

Abstract. The firn layer on the Greenland Ice Sheet (GrIS) plays a crucial role in buffering surface meltwater runoff, which is constrained by the available firn pore space and impermeable ice layers that limit deeper meltwater percolation. Understanding these firn properties is essential for predicting current and future meltwater runoff and its contribution to global sea-level rise. While very-high-frequency (VHF) radars have been extensively used for surveying the GrIS, their lower bandwidth restricts direct firn stratigraphy extraction. In this study, we use concurrent VHF airborne radar and laser altimetry data collected as part of Operation IceBridge over the 2011–2019 period to investigate our hypothesis that vertical heterogeneities in firn (i.e. ice layers) cause vertical offsets in the radar surface reflection (dz). Our results, corroborated by modelling and firn core analyses, show that a dz larger than 1 m is strongly related to the vertical heterogeneity of a firn profile and effectively delineates between vertically homogeneous and vertically heterogeneous firn profiles over a depth range of ∼ 4 m. Temporal variations in dz align with climatic events and reveal an expansion of heterogeneous firn between 2011–2013 covering an area of ∼ 350 815 km2, followed by firn replenishment over the years 2014–2019 spanning an area of ∼ 667 725 km2. Our approach reveals the firn evolution of key regions on the Greenland Ice Sheet, providing valuable insights for detecting potential alterations in meltwater runoff patterns.

Impact of Aerosols on the Macrophysical and Microphysical Characteristics of Ice-Phase and Mixed-Phase Clouds over the Tibetan Plateau

Using CloudSat/CALIPSO satellite data and ERA5 reanalysis data from 2006 to 2010, the effects of aerosols on ice- and mixed-phase, single-layer, non-precipitating clouds over the Tibetan Plateau during nighttime in the MAM (March to May), JJA (June to August), SON (September to November), and DJF (December to February) seasons were examined. The results indicated the following: (1) The macrophysical and microphysical characteristics of ice- and mixed-phase clouds exhibit a nonlinear trend with increasing aerosol optical depth (AOD). When the logarithm of AOD (lnAOD) was ≤−4.0, with increasing AOD during MAM and JJA nights, the cloud thickness and ice particle effective radius of ice-phase clouds and mixed-phase clouds, the ice water path and ice particle number concentration of ice-phase clouds, and the liquid water path and cloud fraction of mixed-phase clouds all decreased; during SON and DJF nights, the cloud thickness of ice-phase clouds, cloud top height, liquid droplet number concentration, and liquid water path of mixed-phase clouds all decreased. When the lnAOD was >−4.0, with increasing AOD during MAM and JJA nights, the cloud top height, cloud base height, cloud fraction, and ice particle number concentration of ice-phase clouds, and the ice water path of mixed-phase clouds all increased; during SON and DJF nights, the cloud fraction of mixed-phase clouds and the ice water path of ice-phase clouds all increased. (2) Under the condition of excluding meteorological factors, including the U-component of wind, V-component of wind, pressure vertical velocity, temperature, and relative humidity at the atmospheric pressure heights near the average cloud top height, within the cloud, and the average cloud base height, as well as precipitable water vapor, convective available potential energy, and surface pressure. During MAM and JJA nights. When the lnAOD was ≤−4.0, an increase in aerosols may have led to a decrease in the thickness of ice and mixed-phase cloud layers, as well as a reduction in cloud water path values. In contrast, when the lnAOD was >−4.0, an increase in aerosols may contribute to elevated cloud base and cloud top heights for ice-phase clouds. During SON and DJF nights, changes in various cloud characteristics may be influenced by both aerosols and meteorological factors.

Estimation of the suspended sediment concentration in ice-covered channels based on the gravitational theory

Rivers in cold high latitude northern regions usually freeze during winter, thereby forming an ice cover on the water surface. The additional ice cover significantly changes the flow characteristics of the rivers and hence affects the process of water and sediment transport. Understanding the vertical profile of suspended sediment concentration (SSC) in ice-covered channel flows is helpful for investigating the sediment deposition and the bed morphology of rivers located in cold regions. To predict the SSC distribution in ice-covered channels, this study is based on the gravitational theory and puts forward a new SSC model by dividing the flow underneath the ice cover into the upper ice layer and the lower bed layer with the maximum streamwise velocity and the zero Reynolds shear stress as the interface. The calculated results of the new SSC model agree well with the collected experimental data, and the accuracy of the prediction results is evaluated by error analysis. Besides, the differences in the governing equations and the prediction accuracy between the gravitational theory and the diffusion theory are discussed, and thus the pros and cons of the proposed SSC model are illustrated. Although the gravitational theory has the controversial point, the validation results and discussions demonstrate that adopting the gravitational theory to estimate the SSC distribution in ice-covered channels is effective. The theoretical approach used herein can provide some guidance for studying the sediment-laden flow of freezing rivers in high-latitude northern regions.

Southern Ocean drives multidecadal atmospheric CO2 rise during Heinrich Stadials

The last glacial period was punctuated by cold intervals in the North Atlantic region that culminated in extensive iceberg discharge events. These cold intervals, known as Heinrich Stadials, are associated with abrupt climate shifts worldwide. Here, we present CO2 measurements from the West Antarctic Ice Sheet Divide ice core across Heinrich Stadials 2 to 5 at decadal-scale resolution. Our results reveal multi-decadal-scale jumps in atmospheric CO2 concentrations within each Heinrich Stadial. The largest magnitude of change (14.0 ± 0.8 ppm within 55 ± 10 y) occurred during Heinrich Stadial 4. Abrupt rises in atmospheric CO2 are concurrent with jumps in atmospheric CH4 and abrupt changes in the water isotopologs in multiple Antarctic ice cores, the latter of which suggest rapid warming of both Antarctica and Southern Ocean vapor source regions. The synchroneity of these rapid shifts points to wind-driven upwelling of relatively warm, carbon-rich waters in the Southern Ocean, likely linked to a poleward intensification of the Southern Hemisphere westerly winds. Using an isotope-enabled atmospheric circulation model, we show that observed changes in Antarctic water isotopologs can be explained by abrupt and widespread Southern Ocean warming. Our work presents evidence for a multi-decadal- to century-scale response of the Southern Ocean to changes in atmospheric circulation, demonstrating the potential for dynamic changes in Southern Ocean biogeochemistry and circulation on human timescales. Furthermore, it suggests that anthropogenic CO2 uptake in the Southern Ocean may weaken with poleward strengthening westerlies today and into the future.

The Australian ice core programme: history, context, and bibliometric analysis

ABSTRACT East Antarctic ice cores are key to Australia’s Antarctic science objectives. Australia is internationally recognised for its climate contributions relating to impurities, greenhouse gases, and water stable isotope studies. These factors make the Australian ice core programme an effective case study to demonstrate the usefulness of bibliometric studies in understanding the broader impacts of scientific research. This work presents a bibliometric analysis of Australian-authored and/or funded ice core publications from within the Australian Antarctic Territory over the programme’s 50 years (1972–2022). This work also discusses the history and chronology of Australia’s ice core programme to contextualise the results of the bibliometric analysis. These are illustrated by examining the key drilling locations, such as Law Dome, and the resultant research which establishes Australia as a global leader in producing high-resolution ice core records from coastal East Antarctic sites. This article demonstrates the applicability of bibliometric studies to identify themes of research that are intended to inform international and domestic policies for the betterment of society.

Robust and Superhydrophobic Polydimethylsiloxane/Ni@Ti3C2Tx Nanocomposite Coatings with Assembled Eyelash-Like Microstructure Array: A New Approach for Effective Passive Anti-Icing and Active Photothermal Deicing.

To solve the problem of ice condensation and adhesion, it is urgent to develop new anti-icing and deicing technologies. This study presented the development of a highly efficient photothermal-enhanced superhydrophobic PDMS/Ni@Ti3C2Tx composite film (m-NMPA) fabricated cost-effectively and straightforwardly. This film was fabricated utilizing PDMS as a hydrophobic agent, adhesive, and surface protector, while Ni@Ti3C2Tx as a magnetic photothermal filler innovatively. Through a simple spraying method, the filler is guided by a strong magnetic field to self-assemble into an eyelash-like microstructure array. The unique structure not only imparts superhydrophobic properties to the surface but also constructs an efficient "light-capturing" architecture. Remarkably, the m-NMPA film demonstrates outstanding superhydrophobic passive anti-icing and efficient photothermal active deicing performance without the use of fluorinated chemicals. The micro-/nanostructure of the film forms a gas layer, significantly delaying the freezing time of water. Particularly under extreme cold conditions (-30 °C), the freezing time is extended by a factor of 7.3 compared to the bare substrate. Furthermore, under sunlight exposure, surface droplets do not freeze. The excellent photothermal performance is attributed to the firm anchoring of nickel particles on the MXene surface, facilitating effective "point-to-face" photothermal synergy. The eyelash-like microarray structure enhances light-capturing capability, resulting in a high light absorption rate of 98%. Furthermore, the microstructure aids in maintaining heat at the uppermost layer of the surface, maximizing the utilization of thermal energy for ice melting and frost thawing. Under solar irradiation, the m-NMPA film can rapidly melt approximately a 4 mm thick ice layer within 558 s and expel the melted water promptly, reducing the risk of secondary icing. Additionally, the ice adhesion force on the surface of the m-NMPA film is remarkably low, with an adhesion strength of approximately 4.7 kPa for a 1 × 1 cm2 ice column. After undergoing rigorous durability tests, including xenon lamp weathering test, pressure resistance test, repeated adhesive tape testing, xenon lamp irradiation, water drop impact testing, and repeated brushing with hydrochloric acid and particles, the film's surface structure and superhydrophobic performance have remained exceptional. The photothermal superhydrophobic passive anti-icing and active deicing technology in this work rely on sustainable solar energy for efficient heat generation. It presents broad prospects for practical applications with advantages such as simple processing method, environmental friendliness, outstanding anti-icing effects, and exceptional durability.

Biophysical characterization of summer Arctic sea-ice habitats using a remotely operated vehicle-mounted underwater hyperspectral imager

The impact of a rapidly shifting sea-ice cover on climate, ecosystem processes and biophysical habitat properties is not yet fully understood, particularly in the central Arctic Ocean, due to a lack of spatially representative observations. From June to July 2020 during the year-long Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC, leg 4) in the Transpolar Drift, we deployed an underwater hyperspectral imager (UHI) mounted on a remotely operated vehicle (ROV) to characterize the biophysical properties of different sea-ice habitats. We conducted UHI surveys along two transects: i) under level first-year sea ice (FYI), which had a mean sea-ice draft of 1.4 m and was composed of primarily level FYI but also had a relatively shallow ridge (keel depth ∼2.6 m); and ii) under the flank of a ridge, named Jaridge, with a mean ice draft of 1.7 m, which was composed of a mix of level ice and thicker ridge blocks with over 3 m draft. We present a new unsupervised bio-optical quantification algorithm for hyperspectral surveys, the relative ice algal biomass index (RBI), using spectral mixture analysis (SMA). We compare this method to the supervised machine learning habitat classification algorithm, Support Vector Machine (SVM). The RBI showed good agreement to literature-based normalized difference indices (NDI) and PCA analyses, which confirm the RBI as a reliable unsupervised index for ice algal biomass. Our biophysical characterization of the two surveyed regions showed an association of sea-ice algal biomass with sea-ice ridge features. Our surveys also indicate that ice algal spatial distribution may be influenced by ice melt rates, and the formation of under ice meltwater layers and false bottoms. With high spatial coverage (>100 m) at microscale resolution (∼cm) we documented large spatial variability of summer Arctic sea-ice algal biomass and different patterns between adjacent ice habitats. We further demonstrate the need for improved understanding of sea-ice algal spatial variability as a complementary tool for sea-ice biogeochemical sampling using destructive ice core sampling.

Global Ocean Cooling of 2.3°C During the Last Glacial Maximum

AbstractQuantitative constraints on past mean ocean temperature (MOT) critically inform our historical understanding of Earth's energy balance. A recently developed MOT proxy based on paleoatmospheric Xe, Kr, and N2 ratios in ice core air bubbles is a promising tool rooted in the temperature dependences of gas solubilities. However, these inert gases are systematically undersaturated in the modern ocean interior, and it remains unclear how air‐sea disequilibrium may have changed in the past. Here, we carry out 30 tracer‐enabled model simulations under varying circulation, sea ice cover, and wind stress regimes to evaluate air‐sea disequilibrium in the Last Glacial Maximum (LGM) ocean. We find that undersaturation of all three gases was likely reduced, primarily due to strengthened high‐latitude winds, biasing reconstructed MOT by −0.38 ± 0.37°C (1σ). Accounting for air‐sea disequilibrium, paleoatmospheric inert gases indicate that LGM MOT was 2.27 ± 0.46°C (1σ) colder than the pre‐industrial era.

Characterization of Freezing Processes in Drug Substance Bottles by Ice Core Sampling

Freezing of biological drug substance (DS) is a critical unit operation that may impact product quality, potentially leading to protein aggregation and sub-visible particle formation. Cryo-concentration has been identified as a critical parameter to impact protein stability during freezing and should therefore be minimized. The macroscopic cryo-concentration, in the following only referred to as cryo-concentration, is majorly influenced by the freezing rate, which is in turn impacted by product independent process parameters such as the DS container, its size and fill level, and the freezing equipment. (At-scale) process characterization studies are crucial to understand and optimize freezing processes. However, evaluating cryo-concentration requires sampling of the frozen bulk, which is typically performed by cutting the ice block into pieces for subsequent analysis. Also, the large amount of product requirement for these studies is a major limitation. In this study, we report the development of a simple methodology for experimental characterization of frozen DS in bottles at relevant scale using a surrogate solution. The novel ice core sampling technique identifies the axial ice core in the center to be indicative for cryo-concentration, which was measured by osmolality, and concentrations of histidine and polysorbate 80 (PS80), whereas osmolality revealed to be a sensitive read-out. Finally, we exemplify the suitability of the method to study cryo-concentration in DS bottles by comparing cryo-concentrations from different freezing protocols (-80°C vs -40°C). Prolonged stress times during freezing correlated to a higher extent of cryo-concentration quantified by osmolality in the axial center of a 2 L DS bottle.Graphical

20th‐Century Antarctic Sea Level Mitigation Driven by Uncertain East Antarctic Accumulation History

AbstractIncreasing snow accumulation over the Antarctic Ice Sheet may mitigate future sea level rise. However, current estimates of mitigation potential are poorly constrained due to limited records of past variability. We present an annually resolved reconstruction of Antarctic snow accumulation from 1801 to 2000 CE, employing a paleoclimate data assimilation methodology to integrate ice core records with a multi‐model ensemble of climate simulations. Our reconstruction correlates well with instrumental reanalysis, and we find that Antarctic accumulation rates increased over the 20th‐century, resulting in a modest amount (∼1 mm) of sea level mitigation. Mitigation is primarily driven by an accelerating trend since around 1970. Our results contrast with previous mitigation estimates of ∼10–12 mm; this discrepancy is due to poorly constrained baseline estimates of 19th‐century accumulation in East Antarctica. Our reconstruction suggests that the uncertainty of future sea level mitigation from increasing Antarctic accumulation has been underestimated.

Aeolian dynamics at the northern edge of Deliblato (Banat) Sand Sea, Vojvodina, Serbia, at the time of the last deglaciation

Abstract The Deliblato (Banat) Sand Sea, which is one of the largest areas of аeolian sand in Europe, is located near the Iron Gate, which marks the crossing of the Danube River through the biggest gorge of this river. Here, Danubian alluvium has served as the sand source for the Banat Sand Sea, which was formed primarily through southeasterly (Košava) winds. Utilizing a multi-proxy approach, the objective of this study is to gain a better understanding of the environmental dynamics of the Banat Sand Sea. To achieve this goal, we conducted an analysis of an archive representing an approximately 20-m-thick dune formation on the northern edge of this dune field. Using optically stimulated luminescence (OSL) dating, we calculated aeolian sedimentation rates and dune ages. Sand was deposited here approximately between 17 ka and 13 ka. Magnetic susceptibility, grain size, and colorimetric analyses were interpreted in terms of local paleoenvironmental conditions. Calculated sedimentation rates (SR) indicate intensive aeolian deposition during the study period that range from 483 cm/ka to 502 cm/ka. We compared our data with regional and other European archives, as well as with climatic variations recorded in the Greenland ice core North Greenland Ice Core Project (NGRIP).

Quasi-equilibrium chemical evolution in starless cores

ABSTRACT The chemistry of H2O, CO, and other small molecular species in an isolated pre-stellar core, L1544, has been assessed in the context of a comprehensive gas-grain chemical model, coupled to an empirically constrained physical/dynamical model. Our main findings are (i) that the chemical network remains in near equilibrium as the core evolves towards star formation and the molecular abundances change in response to the evolving physical conditions. The gas-phase abundances at any time can be calculated accurately with equilibrium chemistry, and the concept of chemical clocks is meaningless in molecular clouds with similar conditions and dynamical time-scales, and (ii) A comparison of the results of complex and simple chemical networks indicates that the abundances of the dominant oxygen and carbon species, H2O, CO, C, and C+ are reasonably approximated by simple networks. In chemical equilibrium, the time-dependent differential terms vanish, and a simple network reduces to a few algebraic equations. This allows rapid calculation of the abundances most responsible for spectral line radiative cooling in molecular clouds with long dynamical time-scales. The dust ice mantles are highly structured and the ice layers retain a memory of the gas-phase abundances at the time of their deposition. A complex (gas-phase and gas-grain) chemical structure therefore exists, with cosmic-ray induced processes dominating in the inner regions. The inferred H2O abundance profiles for L1544 require that the outer parts of the core and also any medium exterior to the core are essentially transparent to the interstellar radiation field.

Impact of subsurface crevassing on the depth–age relationship of high-Alpine ice cores extracted at Col du Dôme between 1994 and 2012

Abstract. Three seasonally resolved ice core records covering the 20th century were extracted in 1994, 2004, and 2012 at a nearly identical location from the Col du Dôme (4250 m above sea level, m a.s.l.; Mont Blanc, French Alps) drill site. Here, we complete and combine chemical records of major ions and radiometric measurements of 3H and 210Pb obtained from these three cores with a 3D ice flow model of the Col du Dôme glacier to investigate in detail the origin of discontinuities observed in the depth–age relation of the ice cores drilled in 2004 and 2012. Taking advantage of the granitic bedrock at Col du Dôme, which makes the ice core 210Pb records sensitive to the presence of upstream crevasses, and the fact that the depth–age disturbances are observed at depths for which absolute time markers are available, we draw an overall picture of a dynamic crevasse formation. This can explain the non-disturbed depth–age relation of the ice core drilled in 1994 and the perturbations observed in those drilled in 2004 and 2012. Since crevasses are common at high-Alpine glacier sites, our study points to the important need for rigorous investigations of the depth–age scale and glaciological conditions upstream of drill sites before interpreting high-alpine ice core records in terms of atmospheric changes.

Theoretical and numerical study on transient acoustic wave propagation across ice layers in the Arctic Ocean.

The study of transient acoustic wave propagation across the Arctic Ocean ice layer provides theoretical guidance for the design of trans-ice acoustic communication systems. In this study, the Arctic Ocean was modeled as an ice-water composite structure, where the ice and water are regarded as an elastic solid and liquid, respectively. An analytical transient solution for acoustic wave propagation in this structure was derived using the eigenfunction expansion method. Further, the numerical procedures were presented and used to analyze the acoustic wave propagation characteristics across the ice layer. The results show that waveforms corresponding to the radial displacements are more severely distorted than the axial displacements. The amplitudes of the radial and axial displacements decreased rapidly with increasing propagation distance. The ice thickness had a greater impact on the radial displacement than axial displacement; the thicker the ice, the greater the distortion for both radial and axial displacements.

Himalayas is barrier for transportation of light-absorbing particles from South Asia to inner Third Pole

Through a comprehensive investigation into the historical profiles of black carbon derived from ice cores, the spatial distributions of light-absorbing impurities in snowpit samples, and carbon isotopic compositions of black carbon in snowpit samples of the Third Pole, we have identified that due to barriers of the Himalayas and remove of wet deposition, local sources rather than those from seriously the polluted South Asia are main contributors of light-absorbing impurities in the inner part of the Third Pole. Therefore, reducing emissions from residents of the Third Pole themselves is a more effective way of protecting the glaciers of the inner Third Pole in terms of reducing concentrations of light-absorbing particles in the atmosphere and on glaciers.

Probing the Rock Mass Fraction and Transport Efficiency inside Uranus Using 40Ar Measurements

The bulk of Uranus consists of a rock–ice core, but the relative proportions of rock and ice are unknown. Radioactive decay of potassium in the silicates produces 40Ar. If transport of argon from the core to the gaseous envelope is efficient, a measurement of 40Ar in the envelope will provide a direct constraint on the rock mass present (assuming a chondritic rock composition). The expected 40Ar concentrations in this case would be readily detectable by a mass spectrometer carried by a future atmospheric probe. For a given envelope concentration there is a trade-off between the rock mass present and the transport efficiency; this degeneracy could be overcome by making independent determinations of the rock mass (e.g., by gravity and seismology). Primordial 40Ar is a potential confounding factor, especially if Ar/H2 is significantly enhanced above solar or if degassing of radiogenic 40Ar were inefficient. Unfortunately, the primordial 40Ar/36Ar ratio is very uncertain; better constraints on this ratio through measurement or theory would be very helpful. Pollution of the envelope by silicates is another confounding factor but can be overcome by a measurement of the alkali metals in the envelope.

Developing a Hybrid Wavelet-Artificial Neural Network model for simulating high-resolution Antarctic ice core CO2 concentration records during 9–120 thousand years ago

Developing a Hybrid Wavelet-Artificial Neural Network model for simulating high-resolution Antarctic ice core CO₂ concentration records during 9–120 thousand years ago

Size of icy grains after a heating event in a protoplanetary disk

A heating event in a protoplanetary disk could sublimate H2O ice from icy aggregates. After the heating event, the temperature decreased, and the H2O vapor recondensed to form pure H2O grains through homogeneous nucleation or to form H2O ice layer covering the silicate aggregates left behind during the heating event. As a result, the size of the pure H2O ice grains could change from the commonly assumed value of 0.1μm of the aggregate composing ice-covered grains. Here, we numerically determined the final radii of the grains and the thickness of the H2O ice layer on the aggregates after the cooling. We also estimated the maximum radius of the aggregate consisting of the processed grains. If the size of the silicate aggregate left behind in the heating event was smaller than 0.073μm, continuous collisional growth of the aggregates of the processed grains was possible. If the size was larger than 0.073μm, a short cooling time depending on the heliocentric distance was required to achieve high mechanical strengths of the aggregates for the collisional growth.