The evolution of a Pluto-like system during the migration of the ice giants

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Atmospheric radiative transfer in multiple scattering conditions. Application to NWP models.

High spectral resolution radiative transfer (RT) codes are essential tools in the study of the radiative energy transfer in the Earth atmosphere and a support for the development of parameterizations for fast RT codes used in climate and weather prediction models. Cirrus clouds cover permanently 30% of the Earth's surface, representing an important contribution to the Earth-atmosphere radiation balance. The work has been focussed on the development of the RT model LBLMS. The model, widely tested in the infra-red spectral range, has been extended to the short wave spectrum and it has been used in comparison with airborne and satellite measurements to study the optical properties of cirrus clouds. A new database of single scattering properties has been developed for mid latitude cirrus clouds. Ice clouds are treated as a mixture of ice crystals with various habits. The optical properties of the mixture are tested in comparison to radiometric measurements in selected case studies. Finally, a parameterization of the mixture for application to weather prediction and global circulation models has been developed. The bulk optical properties of ice crystals are parameterized as functions of the effective dimension of measured particle size distributions that are representative of mid latitude cirrus clouds. Tests with the Limited Area Weather Prediction model COSMO have shown the impact of the new parameterization with respect to cirrus cloud optical properties based on ice spheres.

Influence of ice crystal habit and cirrus spatial inhomogeneities on the retrieval of cirrus optical thickness and effective radius

Diese Dissertation untersucht den Einfluss von Eiskristallform und räumlicher Inhomogenität von Zirren auf das Retrieval von optischer Wolkendicke und effektivem Eispartikelradius. Zu diesem Zweck werden flugzeuggetragene spektrale Messungen solarer Strahlung sowie solare und langwellige Strahlungstransfersimulationen durchgeführt. Flugzeuggetragene spektrale aufwärtsgerichtete Radianzen (Strahldichten) sind mit dem SMART-Albedometer (Spectral Modular Airborne Radiation measurement sysTem) während des CIRCLE-2 (CIRrus CLoud Experiment-2) Feldexperiments im Mai 2007 gemessen worden. Basierend auf diesen Radianzdaten werden mittels eines Wolkenretrievalalgorithmus optische Wolkendicken und effektive Eispartikelradien anhand von eindimensionalen Strahlungstransferrechnungen bestimmt. Die Auswirkung der Annahme unterschiedlicher Eiskristallformen auf die retrievten Parameter wird durch Variation der Einfachstreueigenschaften der Eispartikel untersucht. Darüber hinaus wird mittels Strahlungstransferrechnungen auch der Einfluss der Eiskristallform auf den Strahlungsantrieb von Eiswolken ermittelt. Die Frage nach dem relativen Einfluss von räumlicher Wolkeninhomogenität und Eiskristallform wird anhand von dreidimensionalen und independent pixel approximation (IPA) Strahlungssimulationen untersucht. Die Analyse basiert auf einer Modelleiswolke, die aus Daten des NASA (National Aeronautics and Space Administration) TC4 (Tropical Composition, Cloud, and Climate Coupling) Feldexperiments im Sommer 2007 in Costa Rica erzeugt wurde. Lokal gesehen können beide Effekte - Eiskristallform und räumliche Eiswolkeninhomogenität - die gleiche Grössenordnung haben und zu einer Unter- bzw. Überschätzung der retrievten Parameter um 40 – 60% führen. Gemittelt über die ganze Wolke ist jedoch der Einfluss der Eiskristallform viel bedeutender als der von räumlichen Inhomogenitäten.

Airborne in situ measurements of ice particles in the tropical tropopause layer

Ice clouds have a strong effect on the Earth-atmosphere radiative energy balance, on the distribution of condensable gases in the atmosphere, as well as on the chemical composition of the air. The ice particles in these clouds can take on a variety of shapes which makes the description of the cloud microphysical properties more difficult. In the tropical upper troposphere/lower stratosphere (UTLS), a region where ice cloud abundance is relatively high, different types of ice clouds can be observed. However, in situ measurements are rare due to the high altitude of these clouds and the few available research aircraft, only three worldwide, that can fly at such altitudes.rnThis work focuses on in situ measurements of the tropical UTLS clouds performedrnwith a Cloud Imaging Probe (CIP) and a Forward Scattering Spectrometer Probern(FSSP-100), whereof the CIP is the key instrument of this thesis. The CIP is anrnairborne in situ instrument that obtains two-dimensional shadow images of cloud particles. Several cloud microphysical parameters can be derived from these measurements, e.g. number concentrations and size distributions. In order to obtain a high quality data set, a careful image analysis and several corrections need to be applied to the CIP observations. These methods are described in detail.rnMeasurements within the tropical UTLS have been performed during two campaigns:rnSCOUT-O3, 2005 in Northern Australia and SCOUT-AMMA, 2006 inWest Africa. Thernobtained data set includes first observations of subvisible cirrus clouds over a continental area and observations of the anvils of deep convective clouds. The latter can be further divided into clouds in mesoscale convective system outflows of different ages and clouds in overshooting cloud turrets that even penetrated the stratosphere. The microphysical properties of these three cloud types are discussed in detail. Furthermore, the vertical structure of the ice clouds in the UTLS is investigated. The values of the microphysical parameters were found to decrease with increasing altitude in the upper troposphere. Particle numbers and maximum sizes were also decreasing with increasing age of the outflow clouds. Further differences between the deep convective clouds and subvisible cirrus were found in the particle morphology as well as in the ratio of the observed aerosol particles to cloud particles which indicates that the different freezing processes (deposition, contact, immersion freezing) play different roles in the formation of the respective clouds. For the achievementrnof a better microphysical characterisation and description numerical fits have been adjusted onto the cloud particle size distributions of the subvisible cirrus as well as on the size distributions of the clouds at different altitudes in the UTLS.

Fast and accurate numerical solutions in some problems of particle and radiation transport: synthetic acceleration for the method of short characteristics, Doppler-broadened scattering kernel, remote sensing of the cryosphere

The aim of this work is to present various aspects of numerical simulation of particle and radiation transport for industrial and environmental protection applications, to enable the analysis of complex physical processes in a fast, reliable, and efficient way. In the first part we deal with speed-up of numerical simulation of neutron transport for nuclear reactor core analysis. The convergence properties of the source iteration scheme of the Method of Characteristics applied to be heterogeneous structured geometries has been enhanced by means of Boundary Projection Acceleration, enabling the study of 2D and 3D geometries with transport theory without spatial homogenization. The computational performances have been verified with the C5G7 2D and 3D benchmarks, showing a sensible reduction of iterations and CPU time. The second part is devoted to the study of temperature-dependent elastic scattering of neutrons for heavy isotopes near to the thermal zone. A numerical computation of the Doppler convolution of the elastic scattering kernel based on the gas model is presented, for a general energy dependent cross section and scattering law in the center of mass system. The range of integration has been optimized employing a numerical cutoff, allowing a faster numerical evaluation of the convolution integral. Legendre moments of the transfer kernel are subsequently obtained by direct quadrature and a numerical analysis of the convergence is presented. In the third part we focus our attention to remote sensing applications of radiative transfer employed to investigate the Earth's cryosphere. The photon transport equation is applied to simulate reflectivity of glaciers varying the age of the layer of snow or ice, its thickness, the presence or not other underlying layers, the degree of dust included in the snow, creating a framework able to decipher spectral signals collected by orbiting detectors.

Observations and modeling of the Marginal Ice Zone

Global climate change in recent decades has strongly influenced the Arctic generating pronounced warming accompanied by significant reduction of sea ice in seasonally ice-covered seas and a dramatic increase of open water regions exposed to wind [Stephenson et al., 2011]. By strongly scattering the wave energy, thick multiyear ice prevents swell from penetrating deeply into the Arctic pack ice. However, with the recent changes affecting Arctic sea ice, waves gain more energy from the extended fetch and can therefore penetrate further into the pack ice. Arctic sea ice also appears weaker during melt season, extending the transition zone between thick multi-year ice and the open ocean. This region is called the Marginal Ice Zone (MIZ). In the Arctic, the MIZ is mainly encountered in the marginal seas, such as the Nordic Seas, the Barents Sea, the Beaufort Sea and the Labrador Sea. Formed by numerous blocks of sea ice of various diameters (floes) the MIZ, under certain conditions, allows maritime transportation stimulating dreams of industrial and touristic exploitation of these regions and possibly allowing, in the next future, a maritime connection between the Atlantic and the Pacific. With the increasing human presence in the Arctic, waves pose security and safety issues. As marginal seas are targeted for oil and gas exploitation, understanding and predicting ocean waves and their effects on sea ice become crucial for structure design and for real time safety of operations. The juxtaposition of waves and sea ice represents a risk for personnel and equipment deployed on ice, and may complicate critical operations such as platform evacuations. The risk is difficult to evaluate because there are no long-term observations of waves in ice, swell events are difficult to predict from local conditions, ice breakup can occur on very short time-scales and wave-ice interactions are beyond the scope of current forecasting models [Liu and Mollo-Christensen, 1988,Marko, 2003]. In this thesis, a newly developed Waves in Ice Model (WIM) [Williams et al., 2013a,Williams et al., 2013b] and its related Ocean and Sea Ice model (OSIM) will be used to study the MIZ and the improvements of wave modeling in ice infested waters. The following work has been conducted in collaboration with the Nansen Environmental and Remote Sensing Center and within the SWARP project which aims to extend operational services supporting human activity in the Arctic by including forecast of waves in ice-covered seas, forecast of sea-ice in the presence of waves and remote sensing of both waves and sea ice conditions. The WIM will be included in the downstream forecasting services provided by Copernicus marine environment monitoring service.

Micrometer-sized ice particles for planetary-science experiments – II. Bidirectional reflectance

We have measured the bidirectional reflectance of spherical micrometer-sized water-ice particles in the visible spectral range over a wide range of incidence and emission angles. The small ice spheres were produced by spraying fine water droplets directly into liquid nitrogen. The resulting mean particle radii are 1.47 + 0.96 - 0.58 μm. Such a material shares many properties with ice in comets and at the surface of icy satellites. Measurements show that the fresh sample material is highly backscattering, contrasting with natural terrestrial snow and frost. The formation of agglomerates of particles during the sample production results in a noticeable variability of the photometric properties of the samples in their initial state. We have also observed significant temporal evolutions of the scattering behavior of the samples, shifting towards more forward scattering within some tens of hours, resulting most likely from sintering processes. All reflectance data are fitted by the Hapke photometric model (1993 and 2002 formulation) with a one/two/three-parameter Henyey-Greenstein phase function and the resulting Hapke parameters are provided. These parameters can be used to compare laboratory results with the observed photometric behaviors of astronomical objects. We show, in particular, that the optical properties of the fresh micrometer-sized ice samples can be used to reproduce the predominant backscattering in the phase curves of Enceladus and Europa.

Solar radiation over and under sea ice and photographic quantification of Ice algal aggregates during the POLARSTERN cruise ARK-XXVII/3 (IceArc) in summer 2012

The ice cover of the Arctic Ocean has been changing dramatically in the last decades and the consequences for the sea-ice associated ecosystem remain difficult to assess. Algal aggregates underneath sea ice have been described sporadically but the frequency and distribution of their occurrence is not well quantified. We used upward looking images obtained by a remotely operated vehicle (ROV) to derive estimates of ice algal aggregate biomass and to investigate their spatial distribution. During the IceArc expedition (ARK-XXVII/3) of RV Polarstern in late summer 2012, different types of algal aggregates were observed floating underneath various ice types in the Central Arctic basins. Our results show that the floe scale distribution of algal aggregates in late summer is very patchy and determined by the topography of the ice underside, with aggregates collecting in dome shaped structures and at the edges of pressure ridges. The buoyancy of the aggregates was also evident from analysis of the aggregate size distribution. Different approaches used to estimate aggregate biomass yield a wide range of results. This highlights that special care must be taken when upscaling observations and comparing results from surveys conducted using different methods or on different spatial scales.

(Table 1) Attenuation coefficient of lateral propagating light in sea ice, southern Beaufort Sea

The attenuation property of a lateral propagating light (LPL) in sea ice was measured using an artificial lamp in the Canadian Arctic during the 2007/2008 winter. A measurement method is proposed and applied whereby a recording instrument is buried in the sea ice and an artificial lamp is moved across the instrument. The apparent attenuation coefficient µ(lamda) for the lateral propagating light is obtained from the measured logarithmic relative variation rate. With the exception of blue and red lights, the attenuation coefficient changed little with wavelength, but changed considerably with depth. The vertical decrease of the attenuation coefficient was found to be correlated with salinity: the greater the salinity, the greater the attenuation coefficient. A clear linear relation of salinity and the lateral attenuation coefficient with R2 = 0.939 exists to address the close correlation of the attenuation of LPL with the scattering from the brine. The observed attenuation coefficient of LPL is much larger than that of the vertical propagation light, which we speculate to be caused by scattering. Part of this scattered component is transmitted out of the sea ice from the upper and lower surfaces.

(Table 1) White ice and melt pond characteristics in Darnley and Franklin Bay, Canada

Melt pond covered sea ice is a ubiquitous feature of the summertime Arctic Ocean when meltwater collects in lower-lying areas of ice surfaces. Horizontal transects were conducted during June 2008 above and below landfast sea ice with melt ponds to characterize surface and bottom topography together with variations in transmitted spectral irradiance. We captured a rapid progression from a highly flooded sea ice surface with lateral drainage toward flaws and seal breathing holes to the formation of distinct melt ponds with steep edges. As the mass of the ice cover decreased due to meltwater drainage and rose upward with respect to the seawater level, the high-scattering properties of ice above the water level (i.e., white ice) were continuously regenerated, while pond waters remained transparent compared to underlying ice. The relatively stable albedos observed throughout the study, even as ice thickness decreased, were directly related to these surface processes. Transmission through the ice cover of incident irradiance in the 400-700 nm wave band ranged from 38% to 67% and from 5% to 16% beneath ponded and white ice, respectively. Our results show that this transmission varied not only as a function of surface type (melt ponds or white ice) areal coverage but also in relation to ice thickness and proximity to other surface types through the influence of horizontal spreading of light. Thus, in contrast to albedo, this implies that regional transmittance estimates need to consider melt pond size and shape distributions and variations in optical properties and thickness of the ice cover.

(Table 3) Air temperatures, and snow and ice characteristics of fresh and brackish water lakes in the Northwest Territories, Canada

The algorithms designed to estimate snow water equivalent (SWE) using passive microwave measurements falter in lake-rich high-latitude environments due to the emission properties of ice covered lakes on low frequency measurements. Microwave emission models have been used to simulate brightness temperatures (Tbs) for snowpack characteristics in terrestrial environments but cannot be applied to snow on lakes because of the differing subsurface emissivities and scattering matrices present in ice. This paper examines the performance of a modified version of the Helsinki University of Technology (HUT) snow emission model that incorporates microwave emission from lake ice and sub-ice water. Inputs to the HUT model include measurements collected over brackish and freshwater lakes north of Inuvik, Northwest Territories, Canada in April 2008, consisting of snowpack (depth, density, and snow water equivalent) and lake ice (thickness and ice type). Coincident airborne radiometer measurements at a resolution of 80x100 m were used as ground-truth to evaluate the simulations. The results indicate that subsurface media are simulated best when utilizing a modeled effective grain size and a 1 mm RMS surface roughness at the ice/water interface compared to using measured grain size and a flat Fresnel reflective surface as input. Simulations at 37 GHz (vertical polarization) produce the best results compared to airborne Tbs, with a Root Mean Square Error (RMSE) of 6.2 K and 7.9 K, as well as Mean Bias Errors (MBEs) of -8.4 K and -8.8 K for brackish and freshwater sites respectively. Freshwater simulations at 6.9 and 19 GHz H exhibited low RMSE (10.53 and 6.15 K respectively) and MBE (-5.37 and 8.36 K respectively) but did not accurately simulate Tb variability (R= -0.15 and 0.01 respectively). Over brackish water, 6.9 GHz simulations had poor agreement with airborne Tbs, while 19 GHz V exhibited a low RMSE (6.15 K), MBE (-4.52 K) and improved relative agreement to airborne measurements (R = 0.47). Salinity considerations reduced 6.9 GHz errors substantially, with a drop in RMSE from 51.48 K and 57.18 K for H and V polarizations respectively, to 26.2 K and 31.6 K, although Tb variability was not well simulated. With best results at 37 GHz, HUT simulations exhibit the potential to track Tb evolution, and therefore SWE through the winter season.

Radio-echo sounding and ice volume estimates of western Nordenskiöld Land glaciers, Svalbard

As part of ongoing work to obtain a reliable estimate of the total ice volume of Svalbard glaciers and their potential contribution to sea-level rise, we present here volume calculations, with detailed error estimates, for ten glaciers on western Nordenskiöld Land, central Spitsbergen, Svalbard. The volume estimates are based upon a dense net of GPR-retrieved ice thickness data collected over several field campaigns spanning the period 1999-2012. On the basis of the pattern of scattering in theradargrams, we also analyse the hydrothermal structure of these glaciers.

Ice volume estimates from ground-penetrating radar surveys, western Nordenskiöld Land glaciers, Svalbard

As part of ongoing work within the SvalGlac project aimed to obtain a reliable estimate of the total ice volume of Svalbard glaciers and their potential contribution to sea level rise, in this contribution we present volume calculations, with detailed error estimates, for ten glaciers on western Nordenskiöld Land, central Spitsbergen, Svalbard. The volume estimates are based upon a dense net of GPR-retrieved ice thickness data collected over several field campaigns spanning the period 1999-2012, all of them except one within 2010-2012. The total area and volume of the ensemble are 113.38±0.09 km2 and 10.439±0.185 km3, respectively, while the individual areas, volumes and average ice thickness lie within 2.5-49.1 km2, 0.08-5.48 km3 and 29-108 m, respectively. The maximum recorded ice thickness, 265±15 m, corresponds to Fridtjovbreen, which has also the largest average thickness (108±1m). Available empirical formulae for Svalbard glaciers overestimate the total volume of these glaciers by 24% with respect to our calculation. On the basis of the pattern of scattering in the radargrams, we also analyse the hydrothermal structure of these glaciers. Nine out of ten are polythermal, while only one is entirely cold.