Aumento della risoluzione spaziale per il sondaggio di temperatura e umidità da satellite geostazionario mediante radiometria ad onde millimetriche e submillimetriche

The purpose of this study is to develop and evaluate techniques that improve the spatial resolution of the channels already selected in the preliminary studies for Geostationary Observatory for Microwave Atmospheric Soundings (GOMAS). Reference high resolution multifrequency brightness temperatures scenarios have been derived by applying radiative transfer calculation to the spatially and microphysically detailed output of meteorological events simulated by the University of Wisconsin - Non-hydrostatic Model System (UW-NMS). Three approaches, Wiener filter, Super-Resolution and Image Fusion have been applied to some representative GOMAS frequency channels to enhance the resolution of antenna temperatures. The Wiener filter improved resolution of the largely oversampled images by a factor 1.5- 2.0 without introducing any penalty in the radiometric accuracy. Super-resolution, suitable for not largely oversampled images, improved resolution by a factor ~1.5 but introducing an increased radiometric noise by a factor 1.4-2.5. The image fusion allows finally to further increase the spatial frequency of the images obtained by the Wiener filter increasing the total resolution up to a factor 5.0 with an increased radiometric noise closely linked to the radiometric frequency and to the examined case study.

Entwicklung und Evaluierung einer Modellhierarchie zur Simulation der aktinischen Strahlung in aerosolbelasteter und bewölkter Atmosphäre

Das Verstaendnis der in der Atmosphaere stattfindendenOzonproduktions- und Ozonabbaumechanismen ist eines derwichtigsten Ziele der aktuellen Umweltforschung. Da Ozon inder Atmosphaere ausschliesslich durch photochemischeProzesse gebildet wird, ist die Kenntnis der aktinischenStrahlung eine grundlegende Voraussetzung fuer die Prognosevon Ozonkonzentrationen. Die vorliegende Arbeit befasst sichdaher mit der Modellierung der aktinischen Strahlung alsBasis zur Untersuchung photochemischer Prozesse. Dazu wurde eine Modellhierarchie entwickelt, welcheausgehend von der Berechnung der optischen Eigenschaften deratmosphaerischen Extingenten den Strahlungstransport in derAtmosphaere bestimmt. Dazu wurden verschiedeneStrahlungstransportmodelle verwendet, welche sich durch dieLoesungsmethode der Strahlungstransportgleichung und dieAnzahl der raeumlichen Dimensionen unterscheiden. Zur Ueberpruefung der entwickelten Modellhierarchie und zurUntersuchung des Einflusses der verschiedenen Extingentenauf das atmosphaerische Strahlungsfeld wurden zahlreicheSensitivitaetsstudien durchgefuehrt. Dabei zeigte sich, dasssich der Einfluss von Ozon fast ausschliesslich auf denAnteil der Strahlung mit Wellenlaengen kleiner 340 nmauswirkt. Bei der Untersuchung des Einflusses vonAerosolpartikeln und Wolkentropfen auf die aktinischeStrahlung wurde festgestellt, dass die chemischeZusammensetzung, das Vertikalprofil, die Groesse und dieGroessenverteilung der Partikel einen entscheidendenEinfluss nehmen. Eine Evaluierung des 1-dimensionalenStrahlungstransportmodells an Hand von Vergleichen mit denErgebnissen anderen Strahlungstransportmodelle undverschiedenen Messungen ergab, dass die entwickelteModellhierarchie sehr zuverlaessige Resultate liefert. Die Effekte durchbrochener Bewoelkung konnten mit einer2-dimensionalen Simulation untersucht werden. DieUntersuchungen zeigen, dass die haeufig durchgefuehrteapproximative Betrachtung von partieller Bewoelkung mit1-dimensionalen Strahlungstransportmodellen im Vergleich zuden Ergebnissen der 2-dimensionalen Simulation grosse Fehlerverursacht, da die Effekte an den seitlichen Wolkenraendernvon den 1-dimensionalen Strahlungstransportmodellen nichtberuecksichtigt werden koennen.

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.

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.

Antarctic cloud spectral emission from ground-based measurements, a focus on far infrared signatures

The present work belongs to the PRANA project, the first extensive field campaign of observation of atmospheric emission spectra covering the Far InfraRed spectral region, for more than two years. The principal deployed instrument is REFIR-PAD, a Fourier transform spectrometer used by us to study Antarctic cloud properties. A dataset covering the whole 2013 has been analyzed and, firstly, a selection of good quality spectra is performed, using, as thresholds, radiance values in few chosen spectral regions. These spectra are described in a synthetic way averaging radiances in selected intervals, converting them into BTs and finally considering the differences between each pair of them. A supervised feature selection algorithm is implemented with the purpose to select the features really informative about the presence, the phase and the type of cloud. Hence, training and test sets are collected, by means of Lidar quick-looks. The supervised classification step of the overall monthly datasets is performed using a SVM. On the base of this classification and with the help of Lidar observations, 29 non-precipitating ice cloud case studies are selected. A single spectrum, or at most an average over two or three spectra, is processed by means of the retrieval algorithm RT-RET, exploiting some main IR window channels, in order to extract cloud properties. Retrieved effective radii and optical depths are analyzed, to compare them with literature studies and to evaluate possible seasonal trends. Finally, retrieval output atmospheric profiles are used as inputs for simulations, assuming two different crystal habits, with the aim to examine our ability to reproduce radiances in the FIR. Substantial mis-estimations are found for FIR micro-windows: a high variability is observed in the spectral pattern of simulation deviations from measured spectra and an effort to link these deviations to cloud parameters has been performed.

Stratospheric aerosol extinction retrieved from SCIAMACHY measurements in limb geometry

Stratosphärische Partikel sind typischerweise mit dem bloßen Auge nicht wahrnehmbar. Dennoch haben sie einen signifikanten Einfluss auf die Strahlungsbilanz der Erde und die heteorogene Chemie in der Stratosphäre. Kontinuierliche, vertikal aufgelöste, globale Datensätze sind daher essenziell für das Verständnis physikalischer und chemischer Prozesse in diesem Teil der Atmosphäre. Beginnend mit den Messungen des zweiten Stratospheric Aerosol Measurement (SAM II) Instruments im Jahre 1978 existiert eine kontinuierliche Zeitreihe für stratosphärische Aerosol-Extinktionsprofile, welche von Messinstrumenten wie dem zweiten Stratospheric Aerosol and Gas Experiment (SAGE II), dem SCIAMACHY, dem OSIRIS und dem OMPS bis heute fortgeführt wird. rnrnIn dieser Arbeit wird ein neu entwickelter Algorithmus vorgestellt, der das sogenannte ,,Zwiebel-Schäl Prinzip'' verwendet, um Extinktionsprofile zwischen 12 und 33 km zu berechnen. Dafür wird der Algorithmus auf Radianzprofile einzelner Wellenlängen angewandt, die von SCIAMACHY in der Limb-Geometrie gemessen wurden. SCIAMACHY's einzigartige Methode abwechselnder Limb- und Nadir-Messungen bietet den Vorteil, hochaufgelöste vertikale und horizontale Messungen mit zeitlicher und räumlicher Koinzidenz durchführen zu können. Die dadurch erlangten Zusatzinformationen können verwendet werden, um die Effekte von horizontalen Gradienten entlang der Sichtlinie des Messinstruments zu korrigieren, welche vor allem kurz nach Vulkanausbrüchen und für polare Stratosphärenwolken beobachtet werden. Wenn diese Gradienten für die Berechnung von Extinktionsprofilen nicht beachtet werden, so kann dies dazu führen, dass sowohl die optischen Dicke als auch die Höhe von Vulkanfahnen oder polarer Stratosphärenwolken unterschätzt werden. In dieser Arbeit wird ein Verfahren vorgestellt, welches mit Hilfe von dreidimensionalen Strahlungstransportsimulationen und horizontal aufgelösten Datensätzen die berechneten Extinktionsprofile korrigiert.rnrnVergleichsstudien mit den Ergebnissen von Satelliten- (SAGE II) und Ballonmessungen zeigen, dass Extinktionsprofile von stratosphärischen Partikeln mit Hilfe des neu entwickelten Algorithmus berechnet werden können und gut mit bestehenden Datensätzen übereinstimmen. Untersuchungen des Nabro Vulkanausbruchs 2011 und des Auftretens von polaren Stratosphärenwolken in der südlichen Hemisphäre zeigen, dass das Korrekturverfahren für horizontale Gradienten die berechneten Extinktionsprofile deutlich verbessert.

In-situ measurement and characterization of cloud particles using digital in-line holography

Satellite measurement validations, climate models, atmospheric radiative transfer models and cloud models, all depend on accurate measurements of cloud particle size distributions, number densities, spatial distributions, and other parameters relevant to cloud microphysical processes. And many airborne instruments designed to measure size distributions and concentrations of cloud particles have large uncertainties in measuring number densities and size distributions of small ice crystals. HOLODEC (Holographic Detector for Clouds) is a new instrument that does not have many of these uncertainties and makes possible measurements that other probes have never made. The advantages of HOLODEC are inherent to the holographic method. In this dissertation, I describe HOLODEC, its in-situ measurements of cloud particles, and the results of its test flights. I present a hologram reconstruction algorithm that has a sample spacing that does not vary with reconstruction distance. This reconstruction algorithm accurately reconstructs the field to all distances inside a typical holographic measurement volume as proven by comparison with analytical solutions to the Huygens-Fresnel diffraction integral. It is fast to compute, and has diffraction limited resolution. Further, described herein is an algorithm that can find the position along the optical axis of small particles as well as large complex-shaped particles. I explain an implementation of these algorithms that is an efficient, robust, automated program that allows us to process holograms on a computer cluster in a reasonable time. I show size distributions and number densities of cloud particles, and show that they are within the uncertainty of independent measurements made with another measurement method. The feasibility of another cloud particle instrument that has advantages over new standard instruments is proven. These advantages include a unique ability to detect shattered particles using three-dimensional positions, and a sample volume size that does not vary with particle size or airspeed. It also is able to yield two-dimensional particle profiles using the same measurements.

Near-infrared digital photography to estimate snow correlation length for microwave emission modeling

The study is based on experimental work conducted in alpine snow. We made microwave radiometric and near-infrared reflectance measurements of snow slabs under different experimental conditions. We used an empirical relation to link near-infrared reflectance of snow to the specific surface area (SSA), and converted the SSA into the correlation length. From the measurements of snow radiances at 21 and 35 GHz , we derived the microwave scattering coefficient by inverting two coupled radiative transfer models (the sandwich and six-flux model). The correlation lengths found are in the same range as those determined in the literature using cold laboratory work. The technique shows great potential in the determination of the snow correlation length under field conditions.

Modeling the heterogeneous ice and gas coma of Comet 103P/Hartley 2

The spectacular images of Comet 103P/Hartley 2 recorded by the Medium Resolution Instrument (MRI) and High Resolution Instrument (HRI) on board of the Extrasolar Planet Observation and Deep Impact Extended Investigation (EPOXI) spacecraft, as the Deep Impact extended mission, revealed that its bi-lobed very active nucleus outgasses volatiles heterogeneously. Indeed, CO2 is the primary driver of activity by dragging out chunks of pure ice out of the nucleus from the sub-solar lobe that appear to be the main source of water in Hartley 2's coma by sublimating slowly as they go away from the nucleus. However, water vapor is released by direct sublimation of the nucleus at the waist without any significant amount of either CO2 or icy grains. The coma structure for a comet with such areas of diverse chemistry differs from the usual models where gases are produced in a homogeneous way from the surface. We use the fully kinetic Direct Simulation Monte Carlo model of Tenishev et al. (Tenishev, V.M., Combi, M.R., Davidsson, B. [2008]. Astrophys. J. 685, 659-677; Tenishev, V.M., Combi, M.R., Rubin, M. [2011]. Astrophys. J. 732, 104-120) applied to Comet 103P/Hartley 2 including sublimating icy grains to reproduce the observations made by EPOXI and ground-based measurements. A realistic bi-lobed nucleus with a succession of active areas with different chemistry was included in the model enabling us to study in details the coma of Hartley 2. The different gas production rates from each area were found by fitting the spectra computed using a line-by-line non-LTE radiative transfer model to the HRI observations. The presence of icy grains with long lifetimes, which are pushed anti-sunward by radiation pressure, explains the observed OH asymmetry with enhancement on the night side of the coma.

The Value of Adding Optics to Ecosystem Models: A Case Study

Many ecosystem models have been developed to study the ocean's biogeochemical properties, but most of these models use simple formulations to describe light penetration and spectral quality. Here, an optical model is coupled with a previously published ecosystem model that explicitly represents two phytoplankton (picoplankton and diatoms) and two zooplankton functional groups, as well as multiple nutrients and detritus. Surface ocean color fields and subsurface light fields are calculated by coupling the ecosystem model with an optical model that relates biogeochemical standing stocks with inherent optical properties (absorption, scattering); this provides input to a commercially available radiative transfer model (Ecolight). We apply this bio-optical model to the equatorial Pacific upwelling region, and find the model to be capable of reproducing many measured optical properties and key biogeochemical processes in this region. Our model results suggest that non-algal particles largely contribute to the total scattering or attenuation (> 50% at 660 nm) but have a much smaller contribution to particulate absorption (< 20% at 440 nm), while picoplankton dominate the total phytoplankton absorption (> 95% at 440 nm). These results are consistent with the field observations. In order to achieve such good agreement between data and model results, however, key model parameters, for which no field data are available, have to be constrained. Sensitivity analysis of the model results to optical parameters reveals a significant role played by colored dissolved organic matter through its influence on the quantity and quality of the ambient light. Coupling explicit optics to an ecosystem model provides advantages in generating: (1) a more accurate subsurface light-field, which is important for light sensitive biogeochemical processes such as photosynthesis and photo-oxidation, (2) additional constraints on model parameters that help to reduce uncertainties in ecosystem model simulations, and (3) model output which is comparable to basic remotely-sensed properties. In addition, the coupling of biogeochemical models and optics paves the road for future assimilation of ocean color and in-situ measured optical properties into the models.

THE EFFECTS OF IRRADIATION ON HOT JOVIAN ATMOSPHERES: HEAT REDISTRIBUTION AND ENERGY DISSIPATION

Hot Jupiters, due to the proximity to their parent stars, are subjected to a strong irradiating flux that governs their radiative and dynamical properties. We compute a suite of three-dimensional circulation models with dual-band radiative transfer, exploring a relevant range of irradiation temperatures, both with and without temperature inversions. We find that, for irradiation temperatures T irr lsim 2000 K, heat redistribution is very efficient, producing comparable dayside and nightside fluxes. For T irr ≈ 2200-2400 K, the redistribution starts to break down, resulting in a high day-night flux contrast. Our simulations indicate that the efficiency of redistribution is primarily governed by the ratio of advective to radiative timescales. Models with temperature inversions display a higher day-night contrast due to the deposition of starlight at higher altitudes, but we find this opacity-driven effect to be secondary compared to the effects of irradiation. The hotspot offset from the substellar point is large when insolation is weak and redistribution is efficient, and decreases as redistribution breaks down. The atmospheric flow can be potentially subjected to the Kelvin-Helmholtz instability (as indicated by the Richardson number) only in the uppermost layers, with a depth that penetrates down to pressures of a few millibars at most. Shocks penetrate deeper, down to several bars in the hottest model. Ohmic dissipation generally occurs down to deeper levels than shock dissipation (to tens of bars), but the penetration depth varies with the atmospheric opacity. The total dissipated Ohmic power increases steeply with the strength of the irradiating flux and the dissipation depth recedes into the atmosphere, favoring radius inflation in the most irradiated objects. A survey of the existing data, as well as the inferences made from them, reveals that our results are broadly consistent with the observational trends.

Lake surface water temperatures of European Alpine lakes (1989–2013) based on the Advanced Very High Resolution Radiometer (AVHRR) 1 km data set

Lake water temperature (LWT) is an important driver of lake ecosystems and it has been identified as an indicator of climate change. Consequently, the Global Climate Observing System (GCOS) lists LWT as an essential climate variable. Although for some European lakes long in situ time series of LWT do exist, many lakes are not observed or only on a non-regular basis making these observations insufficient for climate monitoring. Satellite data can provide the information needed. However, only few satellite sensors offer the possibility to analyse time series which cover 25 years or more. The Advanced Very High Resolution Radiometer (AVHRR) is among these and has been flown as a heritage instrument for almost 35 years. It will be carried on for at least ten more years, offering a unique opportunity for satellite-based climate studies. Herein we present a satellite-based lake surface water temperature (LSWT) data set for European water bodies in or near the Alps based on the extensive AVHRR 1 km data record (1989–2013) of the Remote Sensing Research Group at the University of Bern. It has been compiled out of AVHRR/2 (NOAA-07, -09, -11, -14) and AVHRR/3 (NOAA-16, -17, -18, -19 and MetOp-A) data. The high accuracy needed for climate related studies requires careful pre-processing and consideration of the atmospheric state. The LSWT retrieval is based on a simulation-based scheme making use of the Radiative Transfer for TOVS (RTTOV) Version 10 together with ERA-interim reanalysis data from the European Centre for Medium-range Weather Forecasts. The resulting LSWTs were extensively compared with in situ measurements from lakes with various sizes between 14 and 580 km2 and the resulting biases and RMSEs were found to be within the range of −0.5 to 0.6 K and 1.0 to 1.6 K, respectively. The upper limits of the reported errors could be rather attributed to uncertainties in the data comparison between in situ and satellite observations than inaccuracies of the satellite retrieval. An inter-comparison with the standard Moderate-resolution Imaging Spectroradiometer (MODIS) Land Surface Temperature product exhibits RMSEs and biases in the range of 0.6 to 0.9 and −0.5 to 0.2 K, respectively. The cross-platform consistency of the retrieval was found to be within ~ 0.3 K. For one lake, the satellite-derived trend was compared with the trend of in situ measurements and both were found to be similar. Thus, orbital drift is not causing artificial temperature trends in the data set. A comparison with LSWT derived through global sea surface temperature (SST) algorithms shows lower RMSEs and biases for the simulation-based approach. A running project will apply the developed method to retrieve LSWT for all of Europe to derive the climate signal of the last 30 years. The data are available at doi:10.1594/PANGAEA.831007.

Zeeman effect in atmospheric O₂ measured by ground-based microwave radiometry

In this work we study the Zeeman effect on stratospheric O₂ using ground-based microwave radiometer measurements. The interaction of the Earth magnetic field with the oxygen dipole leads to a splitting of O₂ energy states, which polarizes the emission spectra. A special campaign was carried out in order to measure this effect in the oxygen emission line centered at 53.07 GHz. Both a fixed and a rotating mirror were incorporated into the TEMPERA (TEMPERature RAdiometer) in order to be able to measure under different observational angles. This new configuration allowed us to change the angle between the observational path and the Earth magnetic field direction. Moreover, a high-resolution spectrometer (1 kHz) was used in order to measure for the first time the polarization state of the radiation due to the Zeeman effect in the main isotopologue of oxygen from ground-based microwave measurements. The measured spectra showed a clear polarized signature when the observational angles were changed, evidencing the Zeeman effect in the oxygen molecule. In addition, simulations carried out with the Atmospheric Radiative Transfer Simulator (ARTS) allowed us to verify the microwave measurements showing a very good agreement between model and measurements. The results suggest some interesting new aspects for research of the upper atmosphere.

MEMLS3&a: Microwave Emission Model of Layered Snowpacks adapted to include backscattering

The Microwave Emission Model of Layered Snowpacks (MEMLS) was originally developed for microwave emissions of snowpacks in the frequency range 5–100 GHz. It is based on six-flux theory to describe radiative transfer in snow including absorption, multiple volume scattering, radiation trapping due to internal reflection and a combination of coherent and incoherent superposition of reflections between horizontal layer interfaces. Here we introduce MEMLS3&a, an extension of MEMLS, which includes a backscatter model for active microwave remote sensing of snow. The reflectivity is decomposed into diffuse and specular components. Slight undulations of the snow surface are taken into account. The treatment of like- and cross-polarization is accomplished by an empirical splitting parameter q. MEMLS3&a (as well as MEMLS) is set up in a way that snow input parameters can be derived by objective measurement methods which avoid fitting procedures of the scattering efficiency of snow, required by several other models. For the validation of the model we have used a combination of active and passive measurements from the NoSREx (Nordic Snow Radar Experiment) campaign in Sodankylä, Finland. We find a reasonable agreement between the measurements and simulations, subject to uncertainties in hitherto unmeasured input parameters of the backscatter model. The model is written in Matlab and the code is publicly available for download through the following website: http://www.iapmw.unibe.ch/research/projects/snowtools/memls.html.