VIS-NIR/SWIR Spectral Properties of H 2 O Ice Depending on Particle Size and Surface Temperature

: Laboratory measurements were performed to study the spectral signature of H 2 O ice be-tween 0.4 and 4.2 µ m depending on varying temperatures between 70 and 220 K. Spectral parameters of samples with particle sizes up to ~1360 µ m, particle size mixtures, and different particle shapes were analyzed. The band depth (BD) of the major H 2 O-ice absorptions at 1.04, 1.25, 1.5, and 2 µ m offers an excellent indicator for varying particle sizes in pure H 2 O ice. The spectral changes due to temperature rather, but not exclusively, affect the H 2 O-ice absorptions located at 1.31, 1.57, and 1.65 µ m and the Fresnel reﬂection peaks at 3.1 and 3.2 µ m, which strongly weaken with increasing temperature. As the BDs of the H 2 O-ice absorptions at 1.31, 1.57, and 1.65 µ m increase, the band centers (BCs) of the H 2 O-ice absorptions at 1.25 and 1.5 µ m slightly shift to shorter wavelengths. However, the BCs of the strong H 2 O-ice absorptions can also be affected by saturation in the case of large particles. The collected spectra provide a useful spectral library for future investigations of icy satellites such as Ganymede and Callisto, the major targets of ESA’s JUICE mission.


Introduction and Motivation
H 2 O ice is among the dominant chemical compounds on the surfaces of satellites in the outer solar system.Its varying abundance as well as its physical properties such as particle size and crystallinity are important indicators for understanding the geological and microphysical evolution of the icy satellites' surfaces, including how their environment influences the surface properties.Particularly, the surfaces of the Jovian satellites Ganymede and Callisto, the major targets of the upcoming JUICE mission [1], were found to be covered by a large range of H 2 O-ice particle sizes and different phases of crystallinity (amorphous and crystalline).The properties vary on both satellites depending on geographic latitude and point to a potential relationship with the local surface temperature [2][3][4] and/or sputtering-induced H 2 O-ice enrichment [5].However, major changes in the spectral signature of H 2 O ice only due to temperature have also been previously reported [6][7][8] and may add up to spectral effects related to variations in particle size and crystallinity.Therefore, an accurate characterization of how surface temperature affects the H 2 O-ice spectral signatures in association with varying particle sizes is essential to disentangle both effects and to avoid any misinterpretation of icy surface properties.This would in particular enable one to directly compare the particle sizes of geological structures characterized by different surface temperatures, on the same satellite and between different icy satellites.
So far, modelled H 2 O-ice spectra have mostly been used to estimate or quantify the physical H 2 O-ice parameters.However, these models are based on optical constants Minerals 2021, 11, 1328 2 of 27 derived from transmission spectra [6][7][8], and they need to assume the scattering properties of particles (single-scattering albedo and particle phase function) and their dependencies on particle size, shape, and wavelength.The single-scattering albedo is estimated from optical constants using various simplified models of the scattering of a single particle, but it is recognized that these different models lead to wide divergences in particle size and that one of the most-used versions of the Hapke 'equivalent slab' model [9] is known to be inapplicable, where the imaginary index is large.Even its improved version [10], although providing improvement in strong absorptions, still has difficulties to fit the shape of the Fresnel reflection peak of H 2 O ice at 3.1 µm in some particle size ranges and, on the other hand, leads to a factor of 2.5 overestimation of the particle size (e.g., see Figure 1 in the study of [11]).In addition, no simple model is able to correctly predict the wavelengthdependent particle phase function as function of particle size and shape, despite attempts to model it empirically using sets of exact Monte Carlo calculations [12].Therefore, it is generally assumed to be constant and fully independent of all these parameters.All these assumptions can be critical to correctly simulate spectra, particularly in the case of large H 2 O-ice particles (>500 µm in diameter), as expected in the equatorial regions of Ganymede and Callisto, for which indeed the modelled spectra often have difficulties to match the H 2 O-ice signature, particularly, where absorptions start to saturate [4].So far, modelled H2O-ice spectra have mostly been used to estimate or quantify the physical H2O-ice parameters.However, these models are based on optical constants derived from transmission spectra [6][7][8], and they need to assume the scattering properties of particles (single-scattering albedo and particle phase function) and their dependencies on particle size, shape, and wavelength.The single-scattering albedo is estimated from optical constants using various simplified models of the scattering of a single particle, but it is recognized that these different models lead to wide divergences in particle size and that one of the most-used versions of the Hapke 'equivalent slab' model [9] is known to be inapplicable, where the imaginary index is large.Even its improved version [10], although providing improvement in strong absorptions, still has difficulties to fit the shape of the Fresnel reflection peak of H2O ice at 3.1 µm in some particle size ranges and, on the other hand, leads to a factor of 2.5 overestimation of the particle size (e.g., see Figure 1 in the study of [11]).In addition, no simple model is able to correctly predict the wavelengthdependent particle phase function as function of particle size and shape, despite attempts to model it empirically using sets of exact Monte Carlo calculations [12].Therefore, it is generally assumed to be constant and fully independent of all these parameters.All these assumptions can be critical to correctly simulate spectra, particularly in the case of large H2O-ice particles (>500 µm in diameter), as expected in the equatorial regions of Ganymede and Callisto, for which indeed the modelled spectra often have difficulties to match the H2O-ice signature, particularly, where absorptions start to saturate [4].In order to evaluate in detail how variations in particle size and surface temperature together influence the spectral H 2 O-ice signature, we measured H 2 O-ice samples with particle sizes ranging between ~70 and ~1600 µm in the 0.4-4.2µm wavelength range at temperatures between 70 and 150 K (and 220 K, for selected cases)-the expected temperature range occurring on the icy Galilean satellites [13].Since these spectra were measured in reflectance, they provide the opportunity to directly compare the results to spectral remote sensing data sets of icy satellites.In particular, we investigated: (1) how the spectral signature of H 2 O-ice changes with particle size in the above mentioned size range; (2) how the spectrum of a specific particle size changes with temperature between 70 and 150-220 K, (3) including if there are recognizable differences when the sample is cooled or heated; further, (4) how the spectral signature changes with different particle shapes; and (5) how the relative abundances of small and large particles influence the spectral signature in particle size mixtures.

Sample Preparation
Samples of five different particle sizes (70 ± 30, 300 ± 100, 680 ± 120, 1060 ± 60 and 1360 ± 240 µm) were prepared by spraying spherical water droplets into liquid N 2 or crushing and/or sieving the ice in a freezer (Figure 1a-f).H 2 O-ice particles of ~300 µm were produced with the two methods in order to compare spectra of spherical and irregular shaped particles.Spherical ice particles 70 ± 30 µm were prepared with the setup for the production of icy planetary analogs-B (SPIPA-B), made of a peristatic pump and an ultrasonic nebulizer [14].This setup allows for producing in a controlled and reproducible way spherical particles of water ice of mean diameter of 67 ± 31 µm, as measured by cryo-scanning electron microscopy.Larger spherical particles were produced with a garden sprayer, and the particles crystallized in liquid nitrogen were sieved using 200 and 400 µm sieves.Except for the SPIPA-B particles, the size of all other particles was selected and controlled by sieving in the freezer using sieves for particles of 200, 400, 560, 800, 1000, 1120, and 1600 µm in size, followed by a washing of the samples with liquid nitrogen to remove the remaining smaller ice particles (which otherwise stick to the larger ice grains).In addition, three mixtures of ~70 µm and ~1060 µm particles were prepared, with mass mixing ratios of 27%/73%, 52%/48%, and 77%/23% measured during the preparation (Figure 1g,h).These mixtures were made by weighting each type of H 2 O-ice particles in separate aluminum bottles cooled with liquid nitrogen and mixing them in a single bottle put on a vortex mixer (two times during 10 s) to homogenize the mixture.To avoid the metamorphization (change of size/shape) of the ice particles after their crystallization and before their measurements, all bottles, sieves, and other tools used to prepare the samples were filled and thus regularly cooled with liquid nitrogen and manipulated at the bottom of the freezer (dimensions: H/W/D 90.5/137.9/78.7 cm).The freezer temperature was set at a temperature of −30 • C, which was continuously monitored.As the liquid nitrogen evaporates with time, the bottles are manually refilled with liquid nitrogen, approximately every 5 min.
After production, the H 2 O-ice particles were poured in a cylindrical aluminum sample holder of 48 mm in diameter and 10 mm deep (for samples containing 70 ± 30, 300 ± 100, or 680 ± 120 µm particles) or 20 mm deep (for samples containing 1060 ± 60 or 1360 ± 240 µm particles) in order to provide sufficient optical thickness for reflectance measurements.Moreover, the bottom of the sample holder was covered with a dark opaque layer (black aluminum tape) to limit potential secondary reflections.The samples were finally transferred (under dry and cold conditions, using a box containing liquid nitrogen) to the CarboN-IR environmental chamber.Liquid nitrogen keeps the ice at low temperature, while the outgassing of N 2 keeps the atmosphere dry, avoiding the formation of frost on the ice particles during their transportation.

Experiment Setup and Measurement Procedure
Measurements were performed with the SpectropHotometer with variable INcidence and Emergence (SHINE) coupled to the CarboN-IR environmental chamber at the Cold Surfaces Spectroscopy (CSS) facility (https://cold-spectro.sshade.eu/accessed on 22 November 2021) of the Institut de Planétologie et d'Astrophysique de Grenoble (IPAG) (Grenoble, France).CarboN-IR is a cylindrical stainless-steel chamber containing an inner cell made of copper where the sample holder is placed [15].The main cylinder of the chamber is pumped down to high vacuum (HV).The temperature and pressure inside the cell are controlled with a helium cryostat, a heater and a proportional-integral-derivative (PID) temperature controller, and a turbopump system.There is no liquid N 2 in this chamber.After the introduction of the sample and closing of the inner cell and main chamber, the inner cell is first pumped down to primary vacuum, then 5-10 mbar of pure N 2 from a gas bottle (>99.99(9)%N 2 ; the stated impurity of the N 2 bottle is <5.5 ppm with <3 ppm of H 2 O) is introduced in this inner cell to allow a good thermalization of the sample and avoid the sublimation of the water ice (and consecutive change of particles size/shape).Therefore, there is no issue related to H 2 O as a contaminant in the inner-cell where the sample is kept.Residual water potentially present in the main cylinder of the chamber pumped to HV condenses on the cold point of the cryostat at the bottom of the chamber; therefore, it does not deposit on the windows and does not interfere with the reflectance measurements.
All samples were measured at least at 8 different temperatures from 70 to 150 K and under about 1 × 10 −6 mbar of dinitrogen to ensure good thermalization.The irregularly shaped ~300 µm particles were measured twice in order to detect possible irreversible spectral changes due to cooling and heating of the sample.In contrast to the other measurement sequences, the sample was first cooled from 150 to 70 K and then again heated to 220 K.The top of the CarboN-IR chamber and the inner cell are closed with lids equipped with sapphire windows, enabling the illumination of the sample surface and the measurement of its reflectance over the whole visible and near-infrared ranges by the SHINE spectro-gonio radiometer.
SHINE consists of a goniometer arm illuminating the sample with a 200 mm wide collimated monochromatic beam, provided by a monochromator, and a second goniometer arm holds two detectors covering the spectral range from 0.35 to 5.0 µm [16].Spectra were acquired from 0.4 to 4.2 µm using an incidence angle of 0 • and an emission angle of 30 • .The spectral sampling was of 0.050 µm from 0.4 to 1.0 µm, 0.010 µm from 1.0 to 3.6 µm, and 0.050 µm from 3.6 to 4.2 µm.The spectral resolution was 0.003 µm in the wavelength range of 0.4-0.7 µm, 0.006 µm in the wavelength range of 0.7-1.4µm, 0.013 µm in the range of 1.4-3.0µm, and 0.026 µm in the range of 3.0-4.2µm.At each wavelength, the intensity was recorded at 1 s, repeated 10 times, and averaged.Spectralon and Infragold (Labsphere Inc., North Sutton, NH, USA) were used as references, and the nonperfect Lambertian behavior of the Spectralon was considered in the calibration [17].Finally, the collected spectra were corrected for effects of transmission and multiple reflections due to the sapphire windows of the CarboN-IR chamber, as described in [18].The final absolute radiometric accuracy of the laboratory measurement was better than 3%.

Spectral Parameter Analysis
The spectral signature of H 2 O ice and its variation highly depend on the refractive index and absorption coefficient of H 2 O as a function of wavelength.For detailed information, the reader is referred to excellent reviews by [19][20][21].Particularly in the near-(0.78-1.4µm) and shortwave-infrared (1.4-3 µm) wavelength region, H 2 O-ice exhibits diagnostic absorptions caused by overlapping overtones and combinations of the fundamental vibration modes of the H 2 O molecule [19].However, from the UV through the IR, the H 2 O-ice spectral signature is very sensitive to various physical parameters such as particle size, temperature (including thermal history), and crystallinity [21], which can significantly affect the depth or intensity, wavelength position (or band center), width, and shape of the specific absorption [21].Although these effects make the interpretation of the H 2 O-ice spectra more challenging, they also offer valuable indicators for the chemical and physical conditions of H 2 O ice on icy satellites' surfaces.Therefore, the best possible knowledge how these spectral parameters vary between the collected spectra in relation with particle size and temperature is essential for future investigation, and thus is the focus of this study.
Figure 1 shows the spectra of all samples measured for this study at 70 K.All spectra show the typical spectral signature of H 2 O ice with a relatively high reflectance in the visible spectral region, which decreases with increasing wavelength.The major H 2 O-ice absorptions centered at ~1.04, ~1.25, ~1.5, and ~2 µm appear in every spectrum.The band depth of these absorptions has been often used before as an indicator for abundance and particle size [4], whereas the wavelength position of the band center could give indications of varying crystallinity [6].In addition, absorptions of H 2 O ice at 1.31, 1.57, and 1.65 µm can be recognized in some spectra, which correspond to additional H 2 O-ice absorptions previously described by [6] and interpreted to be related to variations in the temperature and crystallinity of H 2 O ice.At wavelengths larger than 3 µm, the reflectance is usually very low as the broad and very strong fundamental band of H 2 O ice saturates in reflectance.Nevertheless, the characteristic Fresnel reflection peak at 3.1 µm sometimes accompanied by a second peak at 3.2 µm, which was previously used as an indicator for the characterization of H 2 O-ice crystallinity [2], is present in all spectra.In this wavelength region, the absorption of light is so large that the particles reflect the light like polished metal spheres [19].The shallower reflection peak at 3.6 µm, which has been used as an indicator of the ice temperature [22,23], is hardly visible due to the very low reflectance in this wavelength region in the particle size range studied in this work.
Figure 2 presents the spectra measured for the five H 2 O-ice particle sizes (Figure 2a) and the second sample for the 300 µm-sized particles (Figure 2b) separately for the major spectral signatures of H 2 O ice.At this closer inspection, distinct variations in the band center, depth or intensity, and shape of the individual H 2 O-ice absorptions and of the Fresnel reflection peak, respectively, become apparent.In order to analyze the observed variations in detail, we derived the band center (BC), band depth (BD), or intensity (I), the full-width at half maximum (FWHM) and the asymmetry (A) of the major H 2 O-ice absorptions at 1.04, 1.25, 1.5, and 2 µm, the additional H 2 O-ice absorptions at 1.31, 1.57, and 1.65 µm, and the two peaks related to the Fresnel reflection centered at 3.1 and 3.2 µm.The band center (BC) of the H 2 O-ice absorptions is defined as the position of the band minimum after spectral continuum removal, whereas the band depth (BD) is defined as 1-R b /R c [24,25], with R b and R c representing, respectively, the reflectance of the absorption band and the interpolated spectral continuum at the BC wavelength.As a descriptor of the shape of the individual H 2 O-ice absorptions, we computed the full-width at half maximum (FWHM), defined as the distance between the two positions on the band wings at halfway reflectance between the band minimum and the continuum [26].Further, the ratio between the HWHM (half-width at half maximum) derived separately for the two wings works as a potential indicator for any arising asymmetry (A) in the shape of the absorption band.In order to also characterize variations related to the Fresnel reflection peaks at 3.1 and 3.2 µm, the intensity (I) of its maximum reflectance and the corresponding band center (BC) were measured.Despite possible small variations in the BCs of the investigated absorptions, in the following section, we refer to the individual absorptions using their generalized wavelength position, such as 1.04, 1.25, 1.31, 1.5, 1.57, 1.65, and 2.0 µm.The same applies to the Fresnel reflection peaks at 3.1 and 3.2 µm.

Spectral Changes with Particle Size
Figure 3 presents the derived BD and I variations depending on the PS of the H2Oice samples measured at the lowest (70 K) and highest temperature (150 K).Please note that the radiometric accuracy of <3% usually does not exceed the used symbol size and therefore is not indicated in the plots.Variations in the BDs of the H2O-ice absorptions at 1.04, 1.25, 1.5, and 2 µm show a strong relationship to the different particle sizes of the samples and follow the general trends highlighted in previous works [27 ,28] with larger particles more absorptive than smaller ones due to the longer distance a photon travels through the ice between the scattering opportunities than in smaller particles [19].At 70 K, BD1.04 increases from ~9% to ~17% monotonically with particle size (Figure 3a).Similarly, at 70 K, BD1.25 increases from ~21% to ~40% between the particle sizes of ~70 and ~1060 µm and reaches a plateau for the two largest particle sizes.A similar trend can be observed in the case of the spectra derived at higher temperatures (Figure 3b).At 150 K, BD1.04 and BD1.25 reach their maxima at a particle size of ~1060 µm and then slightly decrease by about 3%.It has to be noted, however, that these variations lie within the expected range of accuracy of the radiometric calibration (see above).

Spectral Changes with Particle Size
Figure 3 presents the derived BD and I variations depending on the PS of the H 2 O-ice samples measured at the lowest (70 K) and highest temperature (150 K).Please note that the radiometric accuracy of <3% usually does not exceed the used symbol size and therefore is not indicated in the plots.Variations in the BDs of the H 2 O-ice absorptions at 1.04, 1.25, 1.5, and 2 µm show a strong relationship to the different particle sizes of the samples and follow the general trends highlighted in previous works [27 ,28] with larger particles more absorptive than smaller ones due to the longer distance a photon travels through the ice between the scattering opportunities than in smaller particles [19].At 70 K, BD1.04 increases from ~9% to ~17% monotonically with particle size (Figure 3a).Similarly, at 70 K, BD1.25 increases from ~21% to ~40% between the particle sizes of ~70 and ~1060 µm and reaches a plateau for the two largest particle sizes.A similar trend can be observed in the case of the spectra derived at higher temperatures (Figure 3b).At 150 K, BD1.04 and BD1.25 reach their maxima at a particle size of ~1060 µm and then slightly decrease by about 3%.It has to be noted, however, that these variations lie within the expected range of accuracy of the radiometric calibration (see above).In the case of the stronger H2O-ice absorptions at 1.5 and 2 µm, the BDs reach their maxima at even smaller particle sizes (Figure 3a).At 70 K, BD1.5 increases from ~85% to 89% for particle size from ~70 up to ~680 µm but then decreases to ~81% for larger particles.The BD2.0 increases between 70 and 300 µm-sized particles (~89%) and then decreases with increasing particle size even more strongly down to 68%.At 150 K, similar trends can be observed.BD1.5 increases slightly from ~85% to ~88%, until a plateau is reached for particles with the size of ~680 and ~1060 µm.At the largest particles, the BD1.5 is slightly lower (~84%).BD2.0 exhibits the highest values of ~88% for 70 µm and 300 µm particles, which decrease to ~1060 µm-particles and remain at ~78% also for ~1360 µmsized particles (Figure 3b).The decrease in the BDs of H2O-ice absorptions at larger particles is a well-known characteristic of the H2O-ice spectral properties and is believed to be caused by band saturation.The reflectance in the band center at 1.5 and 2 µm reaches a minimum beginning at particle sizes of ~500 and ~200 µm, respectively.A similar trend has been already described by [4], based on model H2O-ice spectra.In the case of saturation, a further increase in particle size cannot further decrease the reflectance in the band In the case of the stronger H 2 O-ice absorptions at 1.5 and 2 µm, the BDs reach their maxima at even smaller particle sizes (Figure 3a).At 70 K, BD1.5 increases from ~85% to 89% for particle size from ~70 up to ~680 µm but then decreases to ~81% for larger particles.The BD2.0 increases between 70 and 300 µm-sized particles (~89%) and then decreases with increasing particle size even more strongly down to 68%.At 150 K, similar trends can be observed.BD1.5 increases slightly from ~85% to ~88%, until a plateau is reached for particles with the size of ~680 and ~1060 µm.At the largest particles, the BD1.5 is slightly lower (~84%).BD2.0 exhibits the highest values of ~88% for 70 µm and 300 µm particles, which decrease to ~1060 µm-particles and remain at ~78% also for ~1360 µm-sized particles (Figure 3b).The decrease in the BDs of H 2 O-ice absorptions at larger particles is a wellknown characteristic of the H 2 O-ice spectral properties and is believed to be caused by band saturation.The reflectance in the band center at 1.5 and 2 µm reaches a minimum beginning at particle sizes of ~500 and ~200 µm, respectively.A similar trend has been already described by [4], based on model H 2 O-ice spectra.In the case of saturation, a further increase in particle size cannot further decrease the reflectance in the band center below these residuals.Due to the progressive reduction of the reflectance of the continuum with increasing particle size, the band depth decreases again.
At least partly, a similar relationship to the H 2 O-ice particle size can be also observed for the BDs of the weaker H 2 O-ice absorptions at 1.31, 1.57, and 1.65 µm at a given temperature (Figure 3c,d).At 70 K, BD1.31 increases from ~7% to ~19% for particle sizes between ~70 and ~1060 µm and drops to 16% for 1360 µm-sized particles.On the contrary, BD1.57 and BD1.65 increase from ~29% to ~38% and ~50% to 61.3%, respectively, between ~70 and ~300 µm-sized particles and then continuously decrease again toward larger particles to ~14% and ~41%, respectively.Although the overall trend of the BD variations is with particle sizes more or less similar at 150 K, the BD values differ between 70 and 150 K up to ~25%, i.e., more strongly than between the particle sizes.Thus, the BD variations of these weaker absorptions cannot be attributed exclusively to particle size, which is discussed in more detail in Section 3.2.
The same applies for the Fresnel reflection peaks at 3.1 and 3.2 µm (Figure 3e,f).Both peaks show a comparable trend, which is also similar for 70 and 150 K.At least for particle sizes from ~300 to ~1060 µm, an increasing in the intensity (I) of both peaks with increasing particle size can be observed.In addition, a distinct drop in I between the samples of ~70 and ~300 µm-particles can be seen.If particles are larger than a few microns, only a single surface particle reflection is expected at the wavelengths of the Fresnel peak, which is then poorly dependent of particle size [2].According to Figure 1, all particle sizes exhibit a similar very low (nearly 0) reflectance around 2.9-3.5 µm, where the 3.1 µm-peak occurs.In this wavelength range, absorption dominates over scattering (except first-reflection scattering causing the peak).Because, there is no increase in spectral continuum for the small particles (~300 and ~70 µm), and because the Fresnel reflection is a surface effect, the decrease in the peak for 300 and 70 µm particles could probably be due to the increased roughness at small scale for these samples compared to the samples made of larger particles.A higher roughness would produce some shadowing effects, mimic a larger porosity of the first layer, and/or induce multiple different directions of first reflections that could decrease the intensity of the Fresnel peak in the averaged-area of the beam light compared to flatter surfaces at the same scale.The BD3.1 and BD3.2 values slightly decrease again between particles with a size of 1060 and 1360 µm.However, as stated above, the overall variations in the intensity I of both peaks over the range of the measured particle sizes do not exceed the I variations between 70 and 150 K, which are discussed in more detail below.
Generally, the BCs of the individual H 2 O-ice absorptions measured at a given temperature remain more or less stable for the studied range of particle sizes (Figure 4).Only the strong H 2 O-ice absorptions at 1.5 and 2 µm show a slight dispersion of BC values over the entire particle size range, which could possibly imply a relation to the saturation of these absorptions.Slight dispersions in the BC also occur with respect to the Fresnel reflection peaks at 3.1 and 3.2 µm (Figure 4e,f).Although, at 70 K, a slight shift of the BC3.2 to shorter wavelengths (from 3.22 to 3.20) with increasing particle size could be possible, such a trend is difficult to see for the BC3.1 and does not exist for the spectra measured at 150 K. Variations with H2O-ice particle size could be only partly recognized for the FWHM derived for the H2O-ice absorptions (Figure 5).At least, the FWHM of the absorptions at 1.5 and 2 µm seem to increase with particle size followed by a slight drop toward ~1360 µm-sized particles for FWHM2.0.Since this trend occurs without any significant shift in the BC (Figure 4) and no apparent relationship between A and particle size of the samples and/or FWHM, the trend between FWHM and particle size could be explained by a broadening of these absorptions due to their progressive saturation.The weaker absorptions do not show significant FWHM variations with particle size.This trend is even more pronounced at 150 K, which is discussed in the next section.Variations with H 2 O-ice particle size could be only partly recognized for the FWHM derived for the H 2 O-ice absorptions (Figure 5).At least, the FWHM of the absorptions at 1.5 and 2 µm seem to increase with particle size followed by a slight drop toward ~1360 µmsized particles for FWHM2.0.Since this trend occurs without any significant shift in the BC (Figure 4) and no apparent relationship between A and particle size of the samples and/or FWHM, the trend between FWHM and particle size could be explained by a broadening of these absorptions due to their progressive saturation.The weaker absorptions do not show significant FWHM variations with particle size.This trend is even more pronounced at 150 K, which is discussed in the next section.

Spectral Changes with Temperature
Figures 6 and 7 show, how the spectral parameters vary depending on the tempera ture exemplarily for the smallest (~70 µm) and largest (~1360 µm) H2O-ice particles of thi study.Generally, the BDs of the H2O-ice absorptions at 1.04, 1.25, 1.5, and 2 µm only vary within ±2% (which is even lower than the derived measurement accuracy) over the stud ied range of temperatures and thus remain fairly constant for a specific particle size (Fig ure 6a,b) and therefore support that the BDs of these absorptions are useful indicators o particle size variations [3,4].As stated above, the deviation (± 10%) of BD2.0 values in th case of large particles (1360 µm) are expected to be caused by band saturation and canno be related to the temperature.
Variations in the spectral parameters depending on temperature are mainly related to the appearance of the additional H2O-ice absorptions such as a distinct kink at ~1.3 µm, close to the center of the absorption at 1.25 µm, and two additional absorption fea tures at ~1.57 and ~1.65 µm, within the high wavelength wing of the absorption at 1.5 µm (Figure 2), which corresponds to the absorptions described for crystalline H2O ice by [6] The H2O-ice absorption at 1.04 µm may also host an additional absorption, but it is too weak to be analyzed in detail.The variations in the BDs of these absorptions with temper ature show a trend (Figure 6c,d), which mostly appear to dominate over the variation due to particle size (Figure 3c,d).BDs of these additional absorptions are strongest at a temperature of 85 K and then weaken progressively with increasing temperature (Figur 6c,d).It has to be noted, however, that the trend weakens at large particles and thus im plies no exclusive trend with temperature independent of the particle size of the H2O-ic sample.By contrast, in the case of the Fresnel reflection peaks (Figure 6e,f), the peaks re main constant within the expected measurement accuracy of 3% between 70 and 150 K but show a trend of weakening intensity I of the peaks with increasing temperature fo large ~1360 µm-sized particles.The difference in the strength and shape of the peak be tween large and small particle sizes diminishes with increasing temperature.

Spectral Changes with Temperature
Figures 6 and 7 show, how the spectral parameters vary depending on the temperature exemplarily for the smallest (~70 µm) and largest (~1360 µm) H 2 O-ice particles of this study.Generally, the BDs of the H 2 O-ice absorptions at 1.04, 1.25, 1.5, and 2 µm only vary within ±2% (which is even lower than the derived measurement accuracy) over the studied range of temperatures and thus remain fairly constant for a specific particle size (Figure 6a,b) and therefore support that the BDs of these absorptions are useful indicators of particle size variations [3,4].As stated above, the deviation (±10%) of BD2.0 values in the case of large particles (1360 µm) are expected to be caused by band saturation and cannot be related to the temperature.Despite that the spectra shown in Figure 2 imply a possible slight shift of the majo H2O-ice absorptions to shorter wavelengths with increasing temperature, their BCs versu temperature do not (Figure 7a,b).However, in the case of the BCs of the additional ab sorptions at 1.31, 1.57, and 1.65 µm (Figure 7c,d) as well as the Fresnel reflection peak (Figure 7e,f), no significant variations appear clearly be related to changes in the temper ature of the H2O-ice samples.Although, however, the BD variations of the additional ab sorptions do not affect their own BC, they apparently affect the BC of the 'parent' absorp tion (Figure 8).Minor BC shifts may already affect the H2O-ice absorption at 1.04 µm a shown in Figure 2, but variations in BC1.04 are difficult to assess, given the spectral reso lution and SNR of the measured spectra.By contrast, Figure 8a illustrates that BC1.25 gen erally lies at shorter wavelengths at 150 K than measured at 70 K.As BD1.31 µm decrease with increasing temperature, the BC of the H2O-ice absorption at 1.25 µm seems to slightly shift to shorter wavelength, i.e., from ~1.27 to ~1.25 µm (BC1.25).A similar trend can be observed with respect to the H2O-ice absorption at ~1.5 µm affected by the absorptions a 1.57 and 1.65 µm (Figure 8b,c).As mentioned above, BC1.5 (and BC2.0) are sometime difficult to measure in the case of larger particle sizes, i.e., ~1060 and ~1360 µm due to the is relatively strong, BC1.5 cannot be measured independently of BC1.57, which could explain the BC1.5 sometimes reaching ~1.57µm as seen for ~1060 µm-sized particles (Figure 8b,c).Although the relationship between BD1.57 and BC1.5 is more difficult to detect, in the case of BD1.65 (Figure 8c), BC1.5 is located at shorter wavelength in the spectra measured at 70 K than at 150 K and at least varies between ~1.52 and ~1.49µm, respectively, while BD1.65 decreases from ~0.60 to ~0.15 across the entire particle size range investigated.Despite a minor increase in the FWHM1.5 for the 70 µm-sized particles and a more pronounced increasing for the larger particles (~1360 µm) with temperature (Figure 9), none of the other studied H2O-ice absorptions show a clear FWHM temperature dependent trend.This broadening of the 1.5 µm-absorption could possibly be caused by a change of the slope due to an additional absorption at 1.8 µm (Figure 2) growing in the continuum, which is too weak to be analyzed in this study but previously described in [6].Variations in the spectral parameters depending on temperature are mainly related to the appearance of the additional H 2 O-ice absorptions such as a distinct kink at ~1.31 µm, close to the center of the absorption at 1.25 µm, and two additional absorption features at ~1.57 and ~1.65 µm, within the high wavelength wing of the absorption at 1.5 µm (Figure 2), which corresponds to the absorptions described for crystalline H 2 O ice by [6].The H 2 O-ice absorption at 1.04 µm may also host an additional absorption, but it is too weak to be analyzed in detail.The variations in the BDs of these absorptions with temperature show a trend (Figure 6c,d), which mostly appear to dominate over the variations due to particle size (Figure 3c,d).BDs of these additional absorptions are strongest at a temperature of 85 K and then weaken progressively with increasing temperature (Figure 6c,d).It has to be noted, however, that the trend weakens at large particles and thus implies no exclusive trend with temperature independent of the particle size of the H 2 O-ice sample.By contrast, in the case of the Fresnel reflection peaks (Figure 6e,f), the peaks remain constant within the expected measurement accuracy of 3% between 70 and 150 K, but show a trend of weakening intensity I of the peaks with increasing temperature for large ~1360 µm-sized particles.The difference in the strength and shape of the peak between large and small particle sizes diminishes with increasing temperature.
Despite that the spectra shown in Figure 2 imply a possible slight shift of the major H 2 O-ice absorptions to shorter wavelengths with increasing temperature, their BCs versus temperature do not (Figure 7a,b).However, in the case of the BCs of the additional absorptions at 1.31, 1.57, and 1.65 µm (Figure 7c,d) as well as the Fresnel reflection peaks (Figure 7e,f), no significant variations appear clearly be related to changes in the temperature of the H 2 O-ice samples.Although, however, the BD variations of the additional absorptions do not affect their own BC, they apparently affect the BC of the 'parent' absorption (Figure 8).Minor BC shifts may already affect the H 2 O-ice absorption at 1.04 µm as shown in Figure 2, but variations in BC1.04 are difficult to assess, given the spectral resolution and SNR of the measured spectra.By contrast, Figure 8a illustrates that BC1.25 generally lies at shorter wavelengths at 150 K than measured at 70 K.As BD1.31 µm decreases with increasing temperature, the BC of the H 2 O-ice absorption at 1.25 µm seems to slightly shift to shorter wavelength, i.e., from ~1.27 to ~1.25 µm (BC1.25).A similar trend can be observed with respect to the H 2 O-ice absorption at ~1.5 µm affected by the absorptions at 1.57 and 1.65 µm (Figure 8b,c).As mentioned above, BC1.5 (and BC2.0) are sometimes difficult to measure in the case of larger particle sizes, i.e., ~1060 and ~1360 µm due to the starting influence of saturation of the absorption.If, in addition, the absorption at 1.57 µm is relatively strong, BC1.5 cannot be measured independently of BC1.57, which could explain the BC1.5 sometimes reaching ~1.57µm as seen for ~1060 µm-sized particles (Figure 8b,c).Although the relationship between BD1.57 and BC1.5 is more difficult to detect, in the case of BD1.65 (Figure 8c), BC1.5 is located at shorter wavelength in the spectra measured at 70 K than at 150 K and at least varies between ~1.52 and ~1.49µm, respectively, while BD1.65 decreases from ~0.60 to ~0.15 across the entire particle size range investigated.Despite a minor increase in the FWHM1.5 for the 70 µm-sized particles and a more pronounced increasing for the larger particles (~1360 µm) with temperature (Figure 9), none of the other studied H 2 O-ice absorptions show a clear FWHM temperature dependent trend.This broadening of the 1.5 µm-absorption could possibly be caused by a change of the slope due to an additional absorption at 1.8 µm (Figure 2) growing in the continuum, which is too weak to be analyzed in this study but previously described in [6].

Spectral Changes between Samples with Different Particle Shape
Figure 10 presents the comparison of the spectral parameters derived for the two samples with ~300 µm-sized particles, one with spherical and one with irregular particles.Whereas the sample with the spherical particles were measured during warming from 70 to 150 K, the sample with irregular particles were measured twice, first during cooling from 150 to 70 K followed by heating again, which ended at 220 K.No major changes in the trends of the spectral parameters of the H2O-ice absorptions (BD, FWHM, BC) other than described above could be observed between the two H2O-ice samples as well as at the same temperature during ascending (sample heating) and descending (sample cooling) temperature ramps of the irregular particles.The only difference in the spectral variations of the spherical or irregular shape of the ~300 µm-sized H2O-ice particles occurs with respect to the BDs of the major H2O-ice absorptions.The BDs of the H2O-ice absorptions at 1.04, 1.25, 1.5, and 2 µm measured for the spherical 300 µm-particles are slightly larger for some temperatures (up to ~10% for BD1.04 and BD1.25 and up to ~5% for BD1.5

Spectral Changes between Samples with Different Particle Shape
Figure 10 presents the comparison of the spectral parameters derived for the two samples with ~300 µm-sized particles, one with spherical and one with irregular particles.Whereas the sample with the spherical particles were measured during warming from 70 to 150 K, the sample with irregular particles were measured twice, first during cooling from 150 to 70 K followed by heating again, which ended at 220 K.No major changes in the trends of the spectral parameters of the H 2 O-ice absorptions (BD, FWHM, BC) other than described above could be observed between the two H 2 O-ice samples as well as at the same temperature during ascending (sample heating) and descending (sample cooling) temperature ramps of the irregular particles.The only difference in the spectral variations of the spherical or irregular shape of the ~300 µm-sized H 2 O-ice particles occurs with respect to the BDs of the major H 2 O-ice absorptions.The BDs of the H 2 O-ice absorptions at 1.04, 1.25, 1.5, and 2 µm measured for the spherical 300 µm-particles are slightly larger for some temperatures (up to ~10% for BD1.04 and BD1.25 and up to ~5% for BD1.5 and BD2.0) at temperatures between 70 and 150 K.The same trend is often visible for the BDs of the absorptions at 1.57 and 1.65 µm but difficult to recognize for BD1.31.Particularly, in the case of BD1.65, the difference reaches a maximum of ~21%.This difference, however, is believed to be rather caused by the different way of the sample preparation than the shape of the particles.We may expect that irregular H 2 O-ice particles have a slightly smaller effective size than the spherical particle sample as a direct result of the sieving process.During the sieving process, elongated particles have the tendency to stay flat on the sieves and have difficulties to go through it.On the other hand, particles with their largest size below the sieve size go, but in average, they have their smallest size smaller than a spherical particle.Previous works already describe that optically irregular particles behave more like particles with a size close to the smallest dimension of their particles than to their largest.Reference [12] showed that the effective optical particle size (D eff ) decreases for oblate and prolate particles compared to the average particle size D (defined from the average of particle cross section σ: = 2 √ σ/π) by a factor of 3 relative to equant particles for aspect ratio around 0.2.In addition, few small particles could have sintered on the largest one with them not well washed during the sample preparation and/or the rougher particle surface of the grinded grains may play a role, with surface asperities acting as small independent scatterers.Finally, the possibly rougher surface of the grinded particles also might cause slight differences in the spectral parameters.By comparison, fewer surface irregularities in the spherical (sprayed) particles could enable the light to penetrate more efficiently in the particle, favoring forward scattering in the medium and requiring a larger number of interactions before emerging, which would result in stronger absorption.Thus, the spectral parameters of the H 2 O-ice absorptions are apparently mainly independent of the particle shape and whether the corresponding temperature of the H 2 O-ice sample is reached either by heating or cooling.
sieving process.During the sieving process, elongated particles have the tendency to stay flat on the sieves and have difficulties to go through it.On the other hand, particles with their largest size below the sieve size go, but in average, they have their smallest size smaller than a spherical particle.Previous works already describe that optically irregular particles behave more like particles with a size close to the smallest dimension of their particles than to their largest.Reference [12] showed that the effective optical particle size (Deff) decreases for oblate and prolate particles compared to the average particle size D (defined from the average of particle cross section σ: 2 /) by a factor of 3 relative to equant particles for aspect ratio around 0.2.In addition, few small particles could have sintered on the largest one with them not well washed during the sample preparation and/or the rougher particle surface of the grinded grains may play a role, with surface asperities acting as small independent scatterers.Finally, the possibly rougher surface of the grinded particles also might cause slight differences in the spectral parameters.By comparison, fewer surface irregularities in the spherical (sprayed) particles could enable the light to penetrate more efficiently in the particle, favoring forward scattering in the medium and requiring a larger number of interactions before emerging, which would result in stronger absorption.Thus, the spectral parameters of the H2O-ice absorptions are apparently mainly independent of the particle shape and whether the corresponding temperature of the H2O-ice sample is reached either by heating or cooling.No significant difference in strength, shape, and position of the Fresnel reflection peaks can be recognized for spherical and irregular 300-µm sized particles, when both samples were heated from 70 to 150 K (Figure 10c,f).The diminishing intensity I of the No significant difference in strength, shape, and position of the Fresnel reflection peaks can be recognized for spherical and irregular 300-µm sized particles, when both samples were heated from 70 to 150 K (Figure 10c,f).The diminishing intensity I of the Fresnel reflection peaks at 3.1 and 3.2 µm together with their slight shift to shorter wavelengths with increasing temperature described above being recognized more clearly due to the wider temperature range of the irregular particles.The Fresnel reflection peak at 3.1 µm and, to a lesser extent, the 3.2 µm peak weaken with increasing temperature (Figure 10c).In addition, a slight shift of the 3.1 µm peak to shorter wavelengths, i.e., from 3.12 to 3.08 µm is observed with increasing temperature from 70 to 220 K (Figure 10f).The same trend, although weak, can also be partly seen for the 3.2 µm-peak, which shifts from 3.24 to 3.19 µm, but with a BC at the shortest wavelength already reached near ~150 K. Only, in the case of the irregular 300 µm-particles the peak at 3.1 µm appears to be stronger at 150 K when measured at the starting point of the cooling sequence from 150 to 70 K than after heating from 70 to 150 K again (Figure 10c).This may be due to some cracking of the crystal due to thermal stress, which reduces the effective grain size and consequently reduces the peak.

Fresnel Reflection-Particle Size, Crystallinity, or Temperature Related?
In previous works, the Fresnel reflection peak at 3.1 µm together with other spectral signatures, i.e., the BCs of the H 2 O-ice absorptions at 1.5 and 2 µm and the existence and strength of the additional absorptions at 1.31, 1.57, and 1.65 µm, has often been used as indicator for the crystallinity of H 2 O ice, which can occur in amorphous and crystalline (cubic and hexagonal) phase depending on the condensation temperature, condensation rate, and temperature history [2,[6][7][8]29].We therefore shortly discuss the results of our measurements here in this context.Since, however, the spectral properties of cubic and hexagonal H 2 O ice do not differ at wavelengths < 70 µm [30], we only distinguish here between crystalline and amorphous H 2 O ice.
Any H 2 O ice formed from the liquid phase is expected to be crystalline [6] except in special conditions of ultrafast cooling of submicron particles, which is not our case.The spherical H 2 O-ice particles were produced by spraying liquid water droplets on liquid nitrogen.The freezing of the droplet is controlled by the Leidenfrost effect: the relatively warm droplet induces the evaporation of the N 2 , allowing the levitation of the droplet above the liquid nitrogen for a few seconds.The duration of the levitation allows the progressive crystallization of the liquid water from the surface to the interior of the spherical droplet and the resulting ice to be fully crystalline [31].Irregular particles were retrieved from a perfectly crystalline ice slab derived from a relatively large mass of water put in the fridge to slowly freeze overnight and then crushed and sieved in the fridge.The temperature of the freezer lies at −30 • C (~243 K), which is much higher than the temperature required to build amorphous H 2 O ice and rather offer perfect conditions to form hexagonal H 2 O ice [8].Therefore, according to the way and the conditions under which the samples were prepared for this study, we expect our samples to consist of crystalline probably hexagonal H 2 O ice.
The typical spectral characteristics of crystalline H 2 O ice as described by [13] can be seen in most of the spectra measured at 70 K, which show a pronounced Fresnel reflection peak at 3.1 µm with a complex shape including an additional shoulder at 3.2 µm, strong major H 2 O-ice absorptions centered at ~1.5 and ~2.02 µm, and distinct and relatively strong additional H 2 O-ice absorptions such as the one centered at 1.65 µm [6-8].Nevertheless, spectral variations depending on particle size and/or temperature as described above, were also used as indicators for changing in the H 2 O-ice crystallinity on icy satellites before.This, a simple-shaped and relatively weak Fresnel reflection and non-existent or weak additional absorptions at 1.31, 1.57, and 1.65 µm observed in spectra of relatively small ~70 and 300 µm-sized H 2 O-ice particles or measured a relatively high temperature (>100 K), can be also observed in the case of amorphous H 2 O ice [2,6,8].What is more, the observed BC variations of the major H 2 O-ice absorptions at 1.5 and 2 µm in the studied particle size and temperature range are also problematic.Amorphous and crystalline H 2 O ice differ with respect to BC1.5 and BC2.0, which shift from 1.492 to 1.507 µm and 1.996 to 2.024 µm [7,8,32], respectively.Probably, as mentioned above, due to a mixture of variations in temperature as well as saturation effects on the BCs of the H 2 O-ice absorptions at ~1.5 and ~2 µm in the case of large particles, the spectra derived in this study exhibit even significant larger BC variations, i.e., from ~1.49 to ~1.57µm (BC1.5) and ~1.95 to ~2.06 µm (BC2.0) or saturation effects prevent an accurate determination of the band center (BC) at all.Only the BC of the 3.1-µm Fresnel reflection peak measured in this work can be clearly associated to crystalline H 2 O ice.Although the peak at 3.1 µm shifts to shorter wavelengths, i.e., from ~3.11 to ~3.09 µm between 120 and 80 K, in the case of amorphous H 2 O ice the peak should be located at even shorter wavelengths (between ~3.07 and ~3.08 µm) as presented by [8] for the same temperature range and thus supports that the H 2 O-ice samples investigated in this study represent crystalline H 2 O ice.

Spectral Changes in Particle Size Mixtures
For the investigation of icy surfaces by telescopes and spacecrafts, it is interesting to characterize how the spectral properties of intimate particle size mixtures would differ depending on the particle size endmembers.Figure 11a presents the spectra measured for the three intimate mixtures prepared with the ~70 and ~1060 µm-sized H 2 O-ice particles separately for the major spectral signatures of H 2 O ice in order to study the spectral effects of either small or large particles onto the spectra of the mixture.All measured spectra nicely exhibit the spectral properties intermediate between the two particle size endmembers.In addition, the variations in the spectral parameters with particle size and temperature exhibit similar trends as described above (Sections 3.1 and 3.2).The H 2 O-ice absorptions at 1.04 and 1.25 µm deepen with increasing portion of the 1060 µm-sized particles.A similar trend is still visible in the case of the H 2 O-ice absorption at 1.5 µm, whereas the BD of the 2 µm-absorption decreases by a larger portion of the 1060 µm-sized particles due to saturation effects.Spectral changes with temperature again rather affect the absorptions at 1.31, 1.57, and 1.65 µm, which weaken with increasing temperature.Although the Fresnel reflection peaks at 3.1 and 3.2 µm also become stronger with growing portion of the larger particles, their intensity (I) also weakens with increasing temperature.
Since the investigation of icy surfaces by telescopes and spacecrafts still largely depends on comparisons with spectra derived from radiative transfer models, it is also interesting how the spectral properties of intimate particle size mixtures would differ from simulated predictions representing linear (areal or macroscopic) mixtures, i.e., a combined spectrum representing the sum of the fractional area times the spectrum of each component or particle size endmember [33].Although such an approach is not theoretically rigorous, when applied to a mixture of small H 2 O-ice particles sticking onto the surface of larger ones, it would have the advantage of simplifying the estimation of important parameters of icy satellite surfaces with respect to the application of more complex and computer-intensive theoretical modelling of intimate mixtures [34,35], which also have other inherent assumptions and simplifications (see Section 1).Therefore, three areal particle size mixtures have been simulated according to the three measured particle size mixtures (Figure 11a) using the measured spectra of the ~70 and the ~1060 µm-sized particles and with varying mass mixing ratios of these two endmembers measured during the preparation of the intimate particle mixtures (27%/73%, 52%/48%, and 77%/23%).A comparison of the simulated and the measured spectra separately for the major H 2 O-ice features and exemplary for the temperatures 70 and 150 K are shown in Figure 11b,c, respectively.The changes in the BD of the H 2 O-ice absorptions and I of the Fresnel reflection peaks depending on particle size and/or temperature derived for the simulated mixtures follow the ones described for the measured particle size mixtures.Only small differences in the spectral signature of H 2 O ice can be recognized between the spectra of the measured intimate mixtures and the simulated mixtures.The only major difference apparently occurs with respect to the strength of the Fresnel reflection peak at 3.1 µm.Intriguingly, these differences appear to be less for the measured and simulated spectra with the largest amount of ~70 µm-sized particles at low temperatures (70 K) or the largest amount of ~1060 µm-sized particles at high temperature (150 K).
This can also be seen in Figures 12-14, which directly compare the spectral parameters derived from the measured (orange) and simulated spectra (black) in comparison to the spectral parameters derived for the spectral endmembers (gray) depending on average particle size or temperature, respectively.The plots support that differences appear to be relatively small for the major H 2 O-ice absorptions at 1.04, 1.25, 1.5, and 2 µm but reveal that differences can be sometimes significant for all three mixtures in the case of the temperature-sensitive additional absorptions at 1.31, 1.57, and 1.65 as well as the Fresnel reflection peak at 3.1 µm.Figures 13 and 14 illustrate the spectral differences between the mixtures dominated either by ~70 or ~1060 µm-sized particles.In the case of the smaller particles exhibiting ~77% of the mixture (Figure 13), the H 2 O-ice absorptions at 1.31, 1.57, and 1.65 appear to differ most at low temperatures and become more similar at higher temperatures.By contrast, the Fresnel reflection peaks are more similar at low temperatures.When the ~1060 µm-sized particles make up ~73% of the mixture (Figure 14), the additional  This can also be seen in Figures 12-14, which directly compare the spectral parameters derived from the measured (orange) and simulated spectra (black) in comparison to the spectral parameters derived for the spectral endmembers (gray) depending on average particle size or temperature, respectively.The plots support that differences appear to particles exhibiting ~77% of the mixture (Figure 13), the H2O-ice absorptions at 1.31, 1.57, and 1.65 appear to differ most at low temperatures and become more similar at higher temperatures.By contrast, the Fresnel reflection peaks are more similar at low temperatures.When the ~1060 µm-sized particles make up ~73% of the mixture (Figure 14), the additional H2O-ice absorptions at 1.31, 1.57, and 1.65 µm are relatively similar in strength across the entire temperature range, but the Fresnel reflection peaks still differ at low temperatures.
Figure 12.Comparison of the (a-g) BDs of the H2O-ice absorption bands and (h,i) the strength of the Fresnel reflection peak at 3.1 and 3.2 µm between the measured (orange) and simulated (black) H2O-ice spectra of three mixtures (27%/73%, 52%/48%, and 77%/23% depending on particle size (PS), i.e., the ~70 and ~1060 µm-sized particles and the average particle size of the mixtures being about ~800, ~550, and ~300 µm, respectively.The measured mixtures are intimate (wt%) while the simulated ones are areal (fractional area in %).The two measured endmember spectra with 70 and 1060 µm-sized  and (h,i) the strength of the Fresnel reflection peak at 3.1 and 3.2 µm between the measured (orange) and simulated (black) H 2 O-ice spectra of three mixtures (27%/73%, 52%/48%, and 77%/23% depending on particle size (PS), i.e., the ~70 and ~1060 µm-sized particles and the average particle size of the mixtures being about ~800, ~550, and ~300 µm, respectively.The measured mixtures are intimate (wt%) while the simulated ones are areal (fractional area in %).The two measured endmember spectra with 70 and 1060 µm-sized particles are also shown as references (gray).Please note that the spectral parameters are exemplarily shown for the lowest (70 K) and highest (150 K) temperature.
particles are also shown as references (gray).Please note that the spectral parameters are exemplarily shown for the lowest (70 K) and highest (150 K) temperature.In order to see how these differences can affect the possibility to distinguish H 2 O-ice particle sizes based on their spectral properties, we compared the differences between measured and simulated spectra to the spectral differences between the average particle sizes of the three mixtures (BD DIFF (PS)).The difference in BD between measured and simulated spectra (BD(mix)-BD(sim)) as well as the percentage of this difference (%BD(mix)-BD(sim)) on the entire band depth of the corresponding spectral features were calculated.Figure 15 shows the results with respect to the major H 2 O-ice absorptions at 1.04, 1.25, 1.5, and 2 µm, and Figure 16 illustrates the effects onto the additional H 2 O-ice absorptions at 1.31, 1.57, and 1.65 µm.Generally, the differences in the BDs of a measured and simulated spectrum do not exceed ~20% for the weaker H 2 O-ice absorptions at 1.04 and 1.25 µm and are less than 5% for the stronger absorptions at 1.5 and 2 µm.Often, however, these differences are comparable with the BD differences between the different measured particle size mixtures (Figure 15i-l) and thus only allow at best to distinguish the varying particle size of the mixture by their BD variations between the mixtures with the smallest (~300 µm) and the largest average particle size (800 µm).The largest differences of up to 40% between the measured and the simulated mixtures occur with respect to the BDs of the absorptions at 1.31, 1.57, and 1.65 µm (Figure 16).However, it is not clear if these differences are related to variations in particle size and/or temperature.

Figure 14.
Comparison of (a-g) the BDs of the H2O-ice absorption bands and (h,i) the strength of the Fresnel reflection peak at 3.1 and 3.2 µm between the measured (orange) and simulated (black) H2O-ice spectra of the mixture with (27%/73% of the ~70 and ~1060 µm-sized particles and an average particle size being about ~800 µm depending on temperature.The measured mixtures are intimate (wt%), while the simulated ones are areal (fractional area in %).
In order to see how these differences can affect the possibility to distinguish H2O-ice particle sizes based on their spectral properties, we compared the differences between measured and simulated spectra to the spectral differences between the average particle sizes of the three mixtures (BD DIFF (PS)).The difference in BD between measured and simulated spectra (BD(mix)-BD(sim)) as well as the percentage of this difference (%BD(mix)-BD(sim)) on the entire band depth of the corresponding spectral features were calculated.Figure 15 shows the results with respect to the major H2O-ice absorptions at 1.04, 1.25, 1.5, and 2 µm, and Figure 16 illustrates the effects onto the additional H2Oice absorptions at 1.31, 1.57, and 1.65 µm.Generally, the differences in the BDs of a measured and simulated spectrum do not exceed ~20% for the weaker H2O-ice absorptions at 1.04 and 1.25 µm and are less than 5% for the stronger absorptions at 1.5 and 2 µm.Often, however, these differences are comparable with the BD differences between the different measured particle size mixtures (Figure 15i-l) and thus only allow at best to distinguish the varying particle size of the mixture by their BD variations between the mixtures with the smallest (~300 µm) and the largest average particle size (800 µm).The largest differences of up to 40% between the measured and the simulated mixtures occur with respect to the BDs of the absorptions at 1.31, 1.57, and 1.65 µm (Figure 16).However, it is not clear if these differences are related to variations in particle size and/or temperature.(b,d,f,h) the % of the BD diff of BD (mix) of three particles size mixtures (27%/73%, 52%/48%, and 77%/23% of the ~70 and ~1060 µm-sized particles) with the average particle size (PS) being about ~800, ~550, and ~300 µm, respectively, in comparison with (i-l) the difference in the BDs (BD DIFF) between the measured particle size mixtures.The measured mixtures are intimate (wt%), while the simulated ones are areal (fractional area in %).(b,d,f,h) the % of the BD diff of BD (mix) of three particles size mixtures (27%/73%, 52%/48%, and 77%/23% of the ~70 and ~1060 µm-sized particles) with the average particle size (PS) being about ~800, ~550, and ~300 µm, respectively, in comparison with (i-l) the difference in the BDs (BD DIFF) between the measured particle size mixtures.The measured mixtures are intimate (wt%), while the simulated ones are areal (fractional area in %).

Discussion
The major results from the presented reflectance measurements on H2O-ice particulate samples for different particle sizes and temperatures can be summarized as follows: • the band depths (BDs) of the H2O-ice absorptions at 1.04, 1.25, 1.5, and 2 µm vary depending on the size of the H2O-ice particles but are fairly stable within the studied range of temperatures, thus supporting that the BDs of these absorptions are a good indicator for particle size variations in the case of pure H2O ice; • the BDs of the H2O-ice absorptions at 1.31, 1.57, and 1.65 µm show a stronger relationship to the temperature but are still to some degree affected by the particle size of the H2O-ice samples; • as the BDs of the H2O-ice absorptions at 1.31, 1.57, and 1.65 µm increase with decreasing temperature, they can affect the band centers (BCs) of the H2O-ice absorptions at 1.25 µm and possibly also at 1.5 and 2 µm; • similarly, also the Fresnel reflection peak at 3.1 µm with its associated peak at 3.2 µm also weakens with increasing temperature and the BC slightly shifts to shorter wavelengths;

Discussion
The major results from the presented reflectance measurements on H 2 O-ice particulate samples for different particle sizes and temperatures can be summarized as follows: similarly, also the Fresnel reflection peak at 3.1 µm with its associated peak at 3.2 µm also weakens with increasing temperature and the BC slightly shifts to shorter wavelengths; • the intensity of the Fresnel reflection peak also varies with H 2 O-ice particle size at low temperatures; • the trends in the spectral parameters (BD, BC, FWHM) with temperature are the same for cooling or heating the samples; • the trends in the spectral parameters (BD, BC, FWHM) with temperature are the same for particles with a spherical and an irregular shape; • the shape and the BC of the H 2 O-ice absorption bands at 1.5 and 2 µm can be affected by the saturation of the bands in the case of large H 2 O-ice particles (>~680 µm).
The results of our measurements support that mapping particle size variations on icy satellites, and in particular on Ganymede's and Callisto's surface, can be performed independently from temperature effects by using the BDs of the major H 2 O-ice absorptions at least in the case of H 2 O ice expected to be more or less pure [4].This is important when studying global variations and local surface features such as impact craters at different geographic locations and/or on different icy bodies.Special care, however, should be taken in the case of large H 2 O-ice particles and major H 2 O-ice absorptions affected by saturation.Moreover, the measured spectra provide an online spectral library [36] for future studies of icy satellites surface and the quantification of H 2 O-ice particle size on these bodies.Since the BDs of the absorptions at 1.04, 1.25, 1.5, and 2 µm inferred from linear mixtures of spectra of individual particles sizes provide a relatively close match to the BDs inferred from measured spectra of particle size binary mixtures, when relatively large variations in the size of the H 2 O-ice particles are expected, linear combinations of measured spectra can prove to be useful for providing a first characterization of the average particle size on these bodies.Although, it has to be noted that this assumption only covers the particle size range discussed here and might not be valid in the case of micrometer and submicrometer-sized H 2 O-ice particles, which are supposed to cover some of the Saturnian icy satellites [37].Areal mixtures have been claimed to be sufficient to describe the mixtures of H 2 O-ice and dark non-ice material(s) on Ganymede and Callisto due to thermal segregation of both materials over time [38].However, [3,4] advocated that the ratios of BDs of the H 2 O-ice absorptions are usable as indicators for particle size variations on icy bodies, even when H 2 O ice is mixed with most non-ice material(s) supposed to exist on the bodies.This has yet to be verified in the laboratory for the particle size ranges and temperatures on these bodies and will be the focus of future works.
Contrary to the H 2 O-ice particle size measurements, the mapping of temperature effects and/or crystallinity of H 2 O ice on icy satellites' surfaces could be more challenging, particularly when the H 2 O ice is mixed with non-icy strongly absorbing compounds.Spectral indicators used in previous studies such as [23,37] to map differences in the surface temperature on icy Saturnian satellites and Europa could be difficult to apply on Ganymede and Callisto, because they make use of the 3.6 µm-peak of H 2 O ice, which is too weak in the particle size range expected on these bodies and studied in this work (Figure 1).Nevertheless, Ref. [39] exploited the relative intensities of the 1.65 and 1.5 µm bands to infer surface temperature but using numerical modeling.They noted that the temperature can be well constrained when the surface is dominated by crystalline H 2 O ice, independently of the presence of other absorbers, ice-free regions, and unusual scattering effects but "errors in derived temperatures can result if the ice is contaminated with other H 2 O ice phases, with bound water in hydrated minerals, by radiation damage, or with other unanticipated absorbers with spectral features coincident with the H 2 O absorptions", such as for Europa and Enceladus.
Although, the samples studied in this work have be shown to be dominated by crystalline (hexagonal) H 2 O ice, spectral variations have been identified that were previously used as an indicator for variations in the crystallinity of H 2 O ice [32].The appearance and strength of the absorptions at 1.31, 1.57, and 1.65 µm and the Fresnel reflection peak of H 2 O ice at 3.1 µm, which have been used in previous works to distinguish between amorphous and crystalline H 2 O ice [2,6], have been found to vary with temperature (and particle size) as well.Further, the H 2 O-ice absorptions 1.5 and 2 µm can exhibit BCs sim-ilar to amorphous H 2 O ice but are rather thought to be the result of saturation affecting these absorptions in the case of large particle sizes.Particularly on icy satellites such as Ganymede, although crystalline H 2 O ice is dominant on its surface, amorphous H 2 O ice is also assumed to be present [2,40] due to radiation-or impact-induced amorphization.In this case, a combined analysis of the spectral parameters of the spectral features of H 2 O-ice including the Fresnel reflection peaks in combination with their geographical and geological context might be necessary to disentangle the influences of the local surface temperature and the radiation environment and to evaluate separately the particles size, crystallinity, and surface temperature of H 2 O ice.

Figure 2 .
Figure 2. Close view onto the individual spectral signatures of the H2O-ice spectra varying with temperature (a) spectra of the samples with the five different particle sizes measured between 70 and 150 K, (b) the irregular 300 µm particles measured during cooling from 150 to 70 K and heating from 70 to 150 K and then 220 K. From left to right: the band depth (BD) of the H2O-ice absorptions at 1.04 and 1.25 µm with the 1.31-µm absorption, at 1.5 µm with the absorption at 1.57 and 1.65 µm, and at 2 µm, as well as the Fresnel reflection peaks at 3.1 and 3.2 µm.All spectra are normalized to the spectral continuum (except for the 3.1 and 3.2 µm peak) and offset for clarity.

Figure 2 .
Figure 2. Close view onto the individual spectral signatures of the H 2 O-ice spectra varying with temperature (a) spectra of the samples with the five different particle sizes measured between 70 and 150 K, (b) the irregular 300 µm particles measured during cooling from 150 to 70 K and heating from 70 to 150 K and then 220 K. From left to right: the band depth (BD) of the H 2 O-ice absorptions at 1.04 and 1.25 µm with the 1.31-µm absorption, at 1.5 µm with the absorption at 1.57 and 1.65 µm, and at 2 µm, as well as the Fresnel reflection peaks at 3.1 and 3.2 µm.All spectra are normalized to the spectral continuum (except for the 3.1 and 3.2 µm peak) and offset for clarity.

Figure 3 .
Figure 3.The derived band depth (BD) of (a,b) the H2O-ice absorptions at 1.04, 1.25, 1.5, and 2 µm and (c,d) at 1.31, 1.57, and 1.65 µm, as well as (e,f) the intensity (I) of the Fresnel reflection peak separately measured at 3.1 and 3.2 µm versus the average particle size (PS) of the samples as derived from spectra measured at a temperature of 70 (left column) and 150 K (right column).The studied H2O-ice features are indicated by different colors.Error bars indicate the particle size range given by the sieves to prepare the samples.

Figure 3 .
Figure 3.The derived band depth (BD) of (a,b) the H 2 O-ice absorptions at 1.04, 1.25, 1.5, and 2 µm and (c,d) at 1.31, 1.57, and 1.65 µm, as well as (e,f) the intensity (I) of the Fresnel reflection peak separately measured at 3.1 and 3.2 µm versus the average particle size (PS) of the samples as derived from spectra measured at a temperature of 70 (left column) and 150 K (right column).The studied H 2 O-ice features are indicated by different colors.Error bars indicate the particle size range given by the sieves to prepare the samples.

Figure 4 .
Figure 4.The derived band center (BC) of (a,b) the H2O-ice absorptions at 1.04, 1.25, 1.5, and 2 µm and (c,d) at 1.31, 1.57, and 1.65 µm, as well as (e,f) the Fresnel reflection peak separately measured at 3.1 and 3.2 µm versus the average particle size (PS) of the samples as derived from spectra measured at a temperature of 70 (left column) and 150 K (right column).The studied H2O-ice features are indicated by different colors.Error bars indicate the particle size range given by the sieves to prepare the samples.

Figure 4 .
Figure 4.The derived band center (BC) of (a,b) the H 2 O-ice absorptions at 1.04, 1.25, 1.5, and 2 µm and (c,d) at 1.31, 1.57, and 1.65 µm, as well as (e,f) the Fresnel reflection peak separately measured at 3.1 and 3.2 µm versus the average particle size (PS) of the samples as derived from spectra measured at a temperature of 70 (left column) and 150 K (right column).The studied H 2 O-ice features are indicated by different colors.Error bars indicate the particle size range given by the sieves to prepare the samples.

Figure 5 .
Figure 5.The derived full width at half maximum (FWHM) of (a,b) the H2O-ice absorptions at 1.04 1.25, 1.5, and 2 µm and (c,d) at 1.31, 1.57, and 1.65 µm.The studied H2O-ice features are indicated by different colors.Error bars indicate the particle size range given by the sieves to prepare th samples.

Figure 5 .
Figure 5.The derived full width at half maximum (FWHM) of (a,b) the H 2 O-ice absorptions at 1.04, 1.25, 1.5, and 2 µm and (c,d) at 1.31, 1.57, and 1.65 µm.The studied H 2 O-ice features are indicated by different colors.Error bars indicate the particle size range given by the sieves to prepare the samples.

Figure 6 .
Figure 6.The derived band depth (BD) of (a,b) the H2O-ice absorptions at 1.04, 1.25, 1.5, and 2 µm and (c,d) at 1.31, 1.57, and 1.65 µm, as well as (e,f) the intensity (I) of the Fresnel reflection peak separately measured at 3.1 and 3.2 µm versus the temperature (T) of the samples as derived from spectra of the smallest and largest studied particle size (PS) of ~70 µm (left column) and ~1360 µm (right column).The studied H2O-ice features are indicated by different colors.

Figure 6 .
Figure 6.The derived band depth (BD) of (a,b) the H 2 O-ice absorptions at 1.04, 1.25, 1.5, and 2 µm and (c,d) at 1.31, 1.57, and 1.65 µm, as well as (e,f) the intensity (I) of the Fresnel reflection peak separately measured at 3.1 and 3.2 µm versus the temperature (T) of the samples as derived from spectra of the smallest and largest studied particle size (PS) of ~70 µm (left column) and ~1360 µm (right column).The studied H 2 O-ice features are indicated by different colors.

Figure 7 .
Figure 7.The derived band center (BC) of (a,b) the H2O-ice absorptions at 1.04, 1.25, 1.5, and 2 µm and (c,d) at 1.31, 1.57, and 1.65 µm, as well as (e,f) the Fresnel reflection peak separately measured at 3.1 and 3.2 µm versus the temperature (T) of the samples as derived from spectra of the smallest and largest studied particle size (PS) of ~70 (left column) and ~1360 µm (right column).The studied H2O-ice features are indicated by different colors.

Figure 7 .
Figure 7.The derived band center (BC) of (a,b) the H 2 O-ice absorptions at 1.04, 1.25, 1.5, and 2 µm and (c,d) at 1.31, 1.57, and 1.65 µm, as well as (e,f) the Fresnel reflection peak separately measured at 3.1 and 3.2 µm versus the temperature (T) of the samples as derived from spectra of the smallest and largest studied particle size (PS) of ~70 (left column) and ~1360 µm (right column).The studied H 2 O-ice features are indicated by different colors.

Minerals 2021, 11 , 1328 13 of 26 Figure 8 .
Figure 8. Band depth (BD) variations of the H2O-ice absorptions at (a) 1.31, (b) 1.57, and (c) 1.65 µm derived for the different particles sizes (PS) at temperatures of 70 and 150 K versus the band center (BC) of the associated H2O-ice absorptions at 1.25 and 1.5 µm, respectively.The temperature is indicated by different colors.The different symbols indicate the particle sizes (PS) of the measured samples.

Figure 8 .
Figure 8. Band depth (BD) variations of the H 2 O-ice absorptions at (a) 1.31, (b) 1.57, and (c) 1.65 µm derived for the different particles sizes (PS) at temperatures of 70 and 150 K versus the band center (BC) of the associated H 2 O-ice absorptions at 1.25 and 1.5 µm, respectively.The temperature is indicated by different colors.The different symbols indicate the particle sizes (PS) of the measured samples.

Figure 8 .
Figure 8. Band depth (BD) variations of the H2O-ice absorptions at (a) 1.31, (b) 1.57, and (c) 1.65 µm derived for the different particles sizes (PS) at temperatures of 70 and 150 K versus the band center (BC) of the associated H2O-ice absorptions at 1.25 and 1.5 µm, respectively.The temperature is indicated by different colors.The different symbols indicate the particle sizes (PS) of the measured samples.

Figure 9 .
Figure 9.The derived full width at half maximum (FWHM) of (a,b) the H2O-ice absorptions at 1.04, 1.25, 1.5, and 2 µm and (c,d) at 1.31, 1.57, and 1.65 µm versus the temperature (T) of the samples as derived from spectra of the smallest and largest studied particle size (PS) of ~70 (left column) and ~1360 µm (right column).The studied H2O-ice features are indicated by different colors.

Figure 9 .
Figure 9.The derived full width at half maximum (FWHM) of (a,b) the H 2 O-ice absorptions at 1.04, 1.25, 1.5, and 2 µm and (c,d) at 1.31, 1.57, and 1.65 µm versus the temperature (T) of the samples as derived from spectra of the smallest and largest studied particle size (PS) of ~70 (left column) and ~1360 µm (right column).The studied H 2 O-ice features are indicated by different colors.

Figure 10 .
Figure 10.Spectral parameters of both spherical and irregular ~300 µm-sized H2O-ice particles during cooling from 150 to 70 K and heating (only for irregular particles) from 70 to 220 K: band depths (BDs) of the H2O-ice absorptions at (a) 1.04, 1.25, 1.5, and 2 µm, (b) at 1.31, 1.57, and 1.65 µm, (c) the intensity (I) of the Fresnel reflection peak at 3.1 and 3.2 µm, as well as (d) the FWHM and (e) the band center (BC) of all H2O-ice absorptions, and (f) the position (BC) of the Fresnel reflection peaks at 3.1 and 3.2 µm versus temperature of the measured samples.Please note that the irregular 300 µmsized particles were the only ones that could be measured up to 220 K.

Figure 10 .
Figure 10.Spectral parameters of both spherical and irregular ~300 µm-sized H 2 O-ice particles during cooling from 150 to 70 K and heating (only for irregular particles) from 70 to 220 K: band depths (BDs) of the H 2 O-ice absorptions at (a) 1.04, 1.25, 1.5, and 2 µm, (b) at 1.31, 1.57, and 1.65 µm, (c) the intensity (I) of the Fresnel reflection peak at 3.1 and 3.2 µm, as well as (d) the FWHM and (e) the band center (BC) of all H 2 O-ice absorptions, and (f) the position (BC) of the Fresnel reflection peaks at 3.1 and 3.2 µm versus temperature of the measured samples.Please note that the irregular 300 µm-sized particles were the only ones that could be measured up to 220 K.

H 2 O
-ice absorptions at 1.31, 1.57, and 1.65 µm are relatively similar in strength across the entire temperature range, but the Fresnel reflection peaks still differ at low temperatures.Minerals 2021, 11, 1328 17 of 26 µm-sized particles at low temperatures (70 K) or the largest amount of ~1060 µm-sized particles at high temperature (150 K).

Figure 11 .
Figure 11.Close view of the individual spectral signatures of the H2O-ice spectra of (a) the three intimate particle size mixtures measured at temperatures between 70 (colored) and 150 K (gray) compared to the spectra of the ~70 and ~1060 µm endmember samples, and (b,c) the three measured intimate mixtures (colored) compared to the corresponding simulated spectra of areal mixtures (gray) of the end-members at (b) 70 K and (c) 150 K. From left to right: the absorptions of the H2O ice at 1.04 and 1.25, 1.5, and 2 µm including the absorptions at 1.31, 1.57, and 1.65 µm, respectively, as well as the Fresnel reflection peak at 3.1 µm.All spectra are normalized (continuum removed, except for the 3.1-µm peak) and offset for clarity.

Figure 11 .
Figure 11.Close view of the individual spectral signatures of the H 2 O-ice spectra of (a) the three intimate particle size mixtures measured at temperatures between 70 (colored) and 150 K (gray) compared to the spectra of the ~70 and ~1060 µm endmember samples, and (b,c) the three measured intimate mixtures (colored) compared to the corresponding simulated spectra of areal mixtures (gray) of the end-members at (b) 70 K and (c) 150 K. From left to right: the absorptions of the H 2 O ice at 1.04 and 1.25, 1.5, and 2 µm including the absorptions at 1.31, 1.57, and 1.65 µm, respectively, as well as the Fresnel reflection peak at 3.1 µm.All spectra are normalized (continuum removed, except for the 3.1-µm peak) and offset for clarity.

Figure 12 .
Figure 12.Comparison of the (a-g) BDs of the H 2 O-ice absorption bands and(h,i) the strength of the Fresnel reflection peak at 3.1 and 3.2 µm between the measured (orange) and simulated (black) H 2 O-ice spectra of three mixtures (27%/73%, 52%/48%, and 77%/23% depending on particle size (PS), i.e., the ~70 and ~1060 µm-sized particles and the average particle size of the mixtures being about ~800, ~550, and ~300 µm, respectively.The measured mixtures are intimate (wt%) while the simulated ones are areal (fractional area in %).The two measured endmember spectra with 70 and 1060 µm-sized particles are also shown as references (gray).Please note that the spectral parameters are exemplarily shown for the lowest (70 K) and highest (150 K) temperature.

Figure 13 .
Figure13.Comparison of (a-g) the BDs of the H2O-ice absorption bands and (h,i) the strength of the Fresnel reflection peak at 3.1 and 3.2 µm between the measured (orange) and simulated (black) H2O-ice spectra of mixture with 77%/23% of the ~70 and ~1060 µm-sized particles and an average particle size being about ~300 µm depending on temperature.The measured mixtures are intimate (wt%), while the simulated ones are areal (fractional area in %).

Figure 13 .
Figure 13.Comparison of (a-g) the BDs of the H 2 O-ice absorption bands and (h,i) the strength of the Fresnel reflection peak at 3.1 and 3.2 µm between the measured (orange) and simulated (black) H 2 O-ice spectra of mixture with 77%/23% of the ~70 and ~1060 µm-sized particles and an average particle size being about ~300 µm depending on temperature.The measured mixtures are intimate (wt%), while the simulated ones are areal (fractional area in %).

Figure 14 .
Figure 14.Comparison of (a-g) the BDs of the H 2 O-ice absorption bands and (h,i) the strength of the Fresnel reflection peak at 3.1 and 3.2 µm between the measured (orange) and simulated (black) H 2 O-ice spectra of the mixture with (27%/73% of the ~70 and ~1060 µm-sized particles and an average particle size being about ~800 µm depending on temperature.The measured mixtures are intimate (wt%), while the simulated ones are areal (fractional area in %).

Figure 16 .
Figure 16.Comparison of (a,d,g) the difference in the BDs (BD DIFF) of the H 2 O-ice absorption bands at 1.31, 1.57, and 1.65 µm between the measured (BD (mix)) and simulated (BD (sim)) H 2 O-ice spectra, (b,e,h) the % of the BD diff of BD (mix) of three particles size mixtures (27%/73%, 52%/48%, and 77%/23% of the ~70 and ~1060 µm-sized particles) with the average particle size (PS) being about ~800, ~550, and ~300 µm, respectively, in comparison with (c,f,i) the difference in the BDs (BD DIFF) between the measured particle size mixtures.The measured mixtures are intimate (wt%), while the simulated ones are areal (fractional area in %).

•
the band depths (BDs) of the H 2 O-ice absorptions at 1.04, 1.25, 1.5, and 2 µm vary depending on the size of the H 2 O-ice particles but are fairly stable within the studied range of temperatures, thus supporting that the BDs of these absorptions are a good indicator for particle size variations in the case of pure H 2 O ice; • the BDs of the H 2 O-ice absorptions at 1.31, 1.57, and 1.65 µm show a stronger relationship to the temperature but are still to some degree affected by the particle size of the H 2 O-ice samples; • as the BDs of the H 2 O-ice absorptions at 1.31, 1.57, and 1.65 µm increase with decreasing temperature, they can affect the band centers (BCs) of the H 2 O-ice absorptions at 1.25 µm and possibly also at 1.5 and 2 µm; •