Enhanced processing of SPOT multispectral satellite imagery for environmental monitoring and modelling

The Taita Hills in southeastern Kenya form the northernmost part of Africa’s Eastern Arc Mountains, which have been identified by Conservation International as one of the top ten biodiversity hotspots on Earth. As with many areas of the developing world, over recent decades the Taita Hills have experienced significant population growth leading to associated major changes in land use and land cover (LULC), as well as escalating land degradation, particularly soil erosion. Multi-temporal medium resolution multispectral optical satellite data, such as imagery from the SPOT HRV, HRVIR, and HRG sensors, provides a valuable source of information for environmental monitoring and modelling at a landscape level at local and regional scales. However, utilization of multi-temporal SPOT data in quantitative remote sensing studies requires the removal of atmospheric effects and the derivation of surface reflectance factor. Furthermore, for areas of rugged terrain, such as the Taita Hills, topographic correction is necessary to derive comparable reflectance throughout a SPOT scene. Reliable monitoring of LULC change over time and modelling of land degradation and human population distribution and abundance are of crucial importance to sustainable development, natural resource management, biodiversity conservation, and understanding and mitigating climate change and its impacts. The main purpose of this thesis was to develop and validate enhanced processing of SPOT satellite imagery for use in environmental monitoring and modelling at a landscape level, in regions of the developing world with limited ancillary data availability. The Taita Hills formed the application study site, whilst the Helsinki metropolitan region was used as a control site for validation and assessment of the applied atmospheric correction techniques, where multiangular reflectance field measurements were taken and where horizontal visibility meteorological data concurrent with image acquisition were available. The proposed historical empirical line method (HELM) for absolute atmospheric correction was found to be the only applied technique that could derive surface reflectance factor within an RMSE of < 0.02 ps in the SPOT visible and near-infrared bands; an accuracy level identified as a benchmark for successful atmospheric correction. A multi-scale segmentation/object relationship modelling (MSS/ORM) approach was applied to map LULC in the Taita Hills from the multi-temporal SPOT imagery. This object-based procedure was shown to derive significant improvements over a uni-scale maximum-likelihood technique. The derived LULC data was used in combination with low cost GIS geospatial layers describing elevation, rainfall and soil type, to model degradation in the Taita Hills in the form of potential soil loss, utilizing the simple universal soil loss equation (USLE). Furthermore, human population distribution and abundance were modelled with satisfactory results using only SPOT and GIS derived data and non-Gaussian predictive modelling techniques. The SPOT derived LULC data was found to be unnecessary as a predictor because the first and second order image texture measurements had greater power to explain variation in dwelling unit occurrence and abundance. The ability of the procedures to be implemented locally in the developing world using low-cost or freely available data and software was considered. The techniques discussed in this thesis are considered equally applicable to other medium- and high-resolution optical satellite imagery, as well the utilized SPOT data.

Modeling of Surface UV Radiation Using Satellite Data

Solar UV radiation is harmful for life on planet Earth, but fortunately the atmospheric oxygen and ozone absorb almost entirely the most energetic UVC radiation photons. However, part of the UVB radiation and much of the UVA radiation reaches the surface of the Earth, and affect human health, environment, materials and drive atmospheric and aquatic photochemical processes. In order to quantify these effects and processes there is a need for ground-based UV measurements and radiative transfer modeling to estimate the amounts of UV radiation reaching the biosphere. Satellite measurements with their near-global spatial coverage and long-term data conti-nuity offer an attractive option for estimation of the surface UV radiation. This work focuses on radiative transfer theory based methods used for estimation of the UV radiation reaching the surface of the Earth. The objectives of the thesis were to implement the surface UV algorithm originally developed at NASA Goddard Space Flight Center for estimation of the surface UV irradiance from the meas-urements of the Dutch-Finnish built Ozone Monitoring Instrument (OMI), to improve the original surface UV algorithm especially in relation with snow cover, to validate the OMI-derived daily surface UV doses against ground-based measurements, and to demonstrate how the satellite-derived surface UV data can be used to study the effects of the UV radiation. The thesis consists of seven original papers and a summary. The summary includes an introduction of the OMI instrument, a review of the methods used for modeling of the surface UV using satellite data as well as the con-clusions of the main results of the original papers. The first two papers describe the algorithm used for estimation of the surface UV amounts from the OMI measurements as well as the unique Very Fast Delivery processing system developed for processing of the OMI data received at the Sodankylä satellite data centre. The third and the fourth papers present algorithm improvements related to the surface UV albedo of the snow-covered land. Fifth paper presents the results of the comparison of the OMI-derived daily erythemal doses with those calculated from the ground-based measurement data. It gives an estimate of the expected accuracy of the OMI-derived sur-face UV doses for various atmospheric and other conditions, and discusses the causes of the differences between the satellite-derived and ground-based data. The last two papers demonstrate the use of the satellite-derived sur-face UV data. Sixth paper presents an assessment of the photochemical decomposition rates in aquatic environment. Seventh paper presents use of satellite-derived daily surface UV doses for planning of the outdoor material weathering tests.

Categorical Meteorological Products: Evaluation and Analysis

In meteorology, observations and forecasts of a wide range of phenomena for example, snow, clouds, hail, fog, and tornados can be categorical, that is, they can only have discrete values (e.g., "snow" and "no snow"). Concentrating on satellite-based snow and cloud analyses, this thesis explores methods that have been developed for evaluation of categorical products and analyses. Different algorithms for satellite products generate different results; sometimes the differences are subtle, sometimes all too visible. In addition to differences between algorithms, the satellite products are influenced by physical processes and conditions, such as diurnal and seasonal variation in solar radiation, topography, and land use. The analysis of satellite-based snow cover analyses from NOAA, NASA, and EUMETSAT, and snow analyses for numerical weather prediction models from FMI and ECMWF was complicated by the fact that we did not have the true knowledge of snow extent, and we were forced simply to measure the agreement between different products. The Sammon mapping, a multidimensional scaling method, was then used to visualize the differences between different products. The trustworthiness of the results for cloud analyses [EUMETSAT Meteorological Products Extraction Facility cloud mask (MPEF), together with the Nowcasting Satellite Application Facility (SAFNWC) cloud masks provided by Météo-France (SAFNWC/MSG) and the Swedish Meteorological and Hydrological Institute (SAFNWC/PPS)] compared with ceilometers of the Helsinki Testbed was estimated by constructing confidence intervals (CIs). Bootstrapping, a statistical resampling method, was used to construct CIs, especially in the presence of spatial and temporal correlation. The reference data for validation are constantly in short supply. In general, the needs of a particular project drive the requirements for evaluation, for example, for the accuracy and the timeliness of the particular data and methods. In this vein, we discuss tentatively how data provided by general public, e.g., photos shared on the Internet photo-sharing service Flickr, can be used as a new source for validation. Results show that they are of reasonable quality and their use for case studies can be warmly recommended. Last, the use of cluster analysis on meteorological in-situ measurements was explored. The Autoclass algorithm was used to construct compact representations of synoptic conditions of fog at Finnish airports.

Burnt area mapping in insular Southeast Asia using medium resolution satellite imagery

Burnt area mapping in humid tropical insular Southeast Asia using medium resolution (250-500m) satellite imagery is characterized by persisting cloud cover, wide range of land cover types, vast amount of wetland areas and highly varying fire regimes. The objective of this study was to deepen understanding of three major aspects affecting the implementation and limits of medium resolution burnt area mapping in insular Southeast Asia: 1) fire-induced spectral changes, 2) most suitable multitemporal compositing methods and 3) burn scars patterns and size distribution. The results revealed a high variation in fire-induced spectral changes depending on the pre-fire greenness of burnt area. It was concluded that this variation needs to be taken into account in change detection based burnt area mapping algorithms in order to maximize the potential of medium resolution satellite data. Minimum near infrared (MODIS band 2, 0.86μm) compositing method was found to be the most suitable for burnt area mapping purposes using Moderate Resolution Imaging Spectroradiometer (MODIS) data. In general, medium resolution burnt area mapping was found to be usable in the wetlands of insular Southeast Asia, whereas in other areas the usability was seriously jeopardized by the small size of burn scars. The suitability of medium resolution data for burnt area mapping in wetlands is important since recently Southeast Asian wetlands have become a major point of interest in many fields of science due to yearly occurring wild fires that not only degrade these unique ecosystems but also create regional haze problem and release globally significant amounts of carbon into the atmosphere due to burning peat. Finally, super-resolution MODIS images were tested but the test failed to improve the detection of small scars. Therefore, super-resolution technique was not considered to be applicable to regional level burnt area mapping in insular Southeast Asia.

Observations of the solar wind-magnetosphere-ionosphere coupling

In this thesis, the solar wind-magnetosphere-ionosphere coupling is studied observationally, with the main focus on the ionospheric currents in the auroral region. The thesis consists of five research articles and an introductory part that summarises the most important results reached in the articles and places them in a wider context within the field of space physics. Ionospheric measurements are provided by the International Monitor for Auroral Geomagnetic Effects (IMAGE) magnetometer network, by the low-orbit CHAllenging Minisatellite Payload (CHAMP) satellite, by the European Incoherent SCATter (EISCAT) radar, and by the Imager for Magnetopause-to-Aurora Global Exploration (IMAGE) satellite. Magnetospheric observations, on the other hand, are acquired from the four spacecraft of the Cluster mission, and solar wind observations from the Advanced Composition Explorer (ACE) and Wind spacecraft. Within the framework of this study, a new method for determining the ionospheric currents from low-orbit satellite-based magnetic field data is developed. In contrast to previous techniques, all three current density components can be determined on a matching spatial scale, and the validity of the necessary one-dimensionality approximation, and thus, the quality of the results, can be estimated directly from the data. The new method is applied to derive an empirical model for estimating the Hall-to-Pedersen conductance ratio from ground-based magnetic field data, and to investigate the statistical dependence of the large-scale ionospheric currents on solar wind and geomagnetic parameters. Equations describing the amount of field-aligned current in the auroral region, as well as the location of the auroral electrojets, as a function of these parameters are derived. Moreover, the mesoscale (10-1000 km) ionospheric equivalent currents related to two magnetotail plasma sheet phenomena, bursty bulk flows and flux ropes, are studied. Based on the analysis of 22 events, the typical equivalent current pattern related to bursty bulk flows is established. For the flux ropes, on the other hand, only two conjugate events are found. As the equivalent current patterns during these two events are not similar, it is suggested that the ionospheric signatures of a flux rope depend on the orientation and the length of the structure, but analysis of additional events is required to determine the possible ionospheric connection of flux ropes.

Ionosphere-atmosphere interaction during solar proton events

Among the most striking natural phenomena affecting ozone are solar proton events (SPE), during which high-energy protons precipitate into the middle atmosphere in the polar regions. Ionisation caused by the protons results in changes in the lower ionosphere, and in production of neutral odd nitrogen and odd hydrogen species which then destroy ozone in well-known catalytic chemical reaction chains. Large SPEs are able to decrease the ozone concentration of upper stratosphere and mesosphere, but are not expected to significantly affect the ozone layer at 15--30~km altitude. In this work we have used the Sodankylä Ion and Neutral Chemistry Model (SIC) in studies of the short-term effects caused by SPEs. The model results were found to be in a good agreement with ionospheric observations from incoherent scatter radars, riometers, and VLF radio receivers as well as with measurements from the GOMOS/Envisat satellite instrument. For the first time, GOMOS was able to observe the SPE effects on odd nitrogen and ozone in the winter polar region. Ozone observations from GOMOS were validated against those from MIPAS/Envisat instrument, and a good agreement was found throughout the middle atmosphere. For the case of the SPE of October/November 2003, long-term ozone depletion was observed in the upper stratosphere. The depletion was further enhanced by the descent of odd nitrogen from the mesosphere inside the polar vortex, until the recovery occurred in late December. During the event, substantial diurnal variation of ozone depletion was seen in the mesosphere, caused mainly by the the strong diurnal cycle of the odd hydrogen species. In the lower ionosphere, SPEs increase the electron density which is very low in normal conditions. Therefore, SPEs make radar observations easier. In the case of the SPE of October, 1989, we studied the sunset transition of negative charge from electrons to ions, a long-standing problem. The observed phenomenon, which is controlled by the amount of solar radiation, was successfully explained by considering twilight changes in both the rate of photodetachment of negative ions and concentrations of minor neutral species. Changes in the magnetic field of the Earth control the extent of SPE-affected area. For the SPE of November 2001, the results indicated that for low and middle levels of geomagnetic disturbance the estimated cosmic radio noise absorption levels based on a magnetic field model are in a good agreement with ionospheric observations. For high levels of disturbance, the model overestimates the stretching of the geomagnetic field and the geographical extent of SPE-affected area. This work shows the importance of ionosphere-atmosphere interaction for SPE studies. By using both ionospheric and atmospheric observations, we have been able to cover for the most part the whole chain of SPE-triggered processes, from proton-induced ionisation to depletion of ozone.

Aerosol Particle Formation: Meteorological and Synoptic Processes behind the Event

Aerosol particles in the atmosphere are known to significantly influence ecosystems, to change air quality and to exert negative health effects. Atmospheric aerosols influence climate through cooling of the atmosphere and the underlying surface by scattering of sunlight, through warming of the atmosphere by absorbing sun light and thermal radiation emitted by the Earth surface and through their acting as cloud condensation nuclei. Aerosols are emitted from both natural and anthropogenic sources. Depending on their size, they can be transported over significant distances, while undergoing considerable changes in their composition and physical properties. Their lifetime in the atmosphere varies from a few hours to a week. New particle formation is a result of gas-to-particle conversion. Once formed, atmospheric aerosol particles may grow due to condensation or coagulation, or be removed by deposition processes. In this thesis we describe analyses of air masses, meteorological parameters and synoptic situations to reveal conditions favourable for new particle formation in the atmosphere. We studied the concentration of ultrafine particles in different types of air masses, and the role of atmospheric fronts and cloudiness in the formation of atmospheric aerosol particles. The dominant role of Arctic and Polar air masses causing new particle formation was clearly observed at Hyytiälä, Southern Finland, during all seasons, as well as at other measurement stations in Scandinavia. In all seasons and on multi-year average, Arctic and North Atlantic areas were the sources of nucleation mode particles. In contrast, concentrations of accumulation mode particles and condensation sink values in Hyytiälä were highest in continental air masses, arriving at Hyytiälä from Eastern Europe and Central Russia. The most favourable situation for new particle formation during all seasons was cold air advection after cold-front passages. Such a period could last a few days until the next front reached Hyytiälä. The frequency of aerosol particle formation relates to the frequency of low-cloud-amount days in Hyytiälä. Cloudiness of less than 5 octas is one of the factors favouring new particle formation. Cloudiness above 4 octas appears to be an important factor that prevents particle growth, due to the decrease of solar radiation, which is one of the important meteorological parameters in atmospheric particle formation and growth. Keywords: Atmospheric aerosols, particle formation, air mass, atmospheric front, cloudiness

Development of biogenic VOC emission inventories for the boreal forest

Luonnosta haihtuvat orgaaniset yhdisteet, joita pääsee ilmaan etenkin metsistä, voivat vaikuttaa paikalliseen ja alueelliseen ilmanlaatuun, koska ne reagoivat ilmakehässä. Niiden reaktiotuotteet voivat myös osallistua uusien hiukkasten muodostumiseen ja kasvuun, millä voi olla vaikutusta ilmakehän säteilytaseeseen ja tätä kautta myös ilmastoon. Hiukkaset absorboivat ja sirottavat auringon säteilyä ja maapallon lämpösäteilyä minkä lisäksi ne vaikuttavat pilvien säteilyominaisuuksiin, määrään ja elinikään. Koko maapallon mittakaavassa luonnosta tulevat hiilivetypäästöt ylittävät ihmistoiminnan aiheuttamat päästöt moninkertaisesti. Tämän vuoksi luonnon päästöjen arviointi on tärkeää kun halutaan kehittää tehokkaita ilmanlaatu- ja ilmastostrategioita. Tämä tutkimus käsittelee boreaalisen metsän hiilivetypäästöjä. Boreaalinen metsä eli pohjoinen havumetsä on suurin maanpäällinen ekosysteemi, ja se ulottuu lähes yhtenäisenä nauhana koko pohjoisen pallonpuoliskon ympäri. Sille on tyypillistä puulajien suhteellisen pieni kirjo sekä olosuhteiden ja kasvun voimakkaat vuodenaikaisvaihtelut. Työssä on tutkittu Suomen yleisimmän boreaalisen puun eli männyn hiilivetypäästöjen vuodenaikaisvaihtelua sekä päästöjen riippuvuutta lämpötilasta ja valosta. Saatuja tuloksia on käytetty yhdessä muiden boreaalisilla puilla tehtyjen päästömittaustulosten kanssa Suomen metsiä varten kehitetyssä päästömallissa. Malli perustuu lisäksi maankäyttötietoihin, suomen metsille kehitettyyn luokitukseen ja meteorologisiin tietoihin, joiden avulla se laskee metsien hiilivetypäästöt kasvukauden aikana. Suomen metsien päästöt koostuvat koko kasvukauden ajan suurelta osin alfa- ja beta-pineenistä sekä delta-kareenista. Kesällä ja syksyllä päästöissä on myös paljon sabineenia, jota tulee etenkin lehtipuista. Päästöt seuraavat lämpötilan keskimääräistä vaihtelua, ovat suurimmillaan maan eteläosissa ja laskevat tasaisesti pohjoiseen siirryttäessä. Metsän isopreenipäästö on suhteellisen pieni – Suomessa tärkein isopreeniä päästävä puu on vähäpäästöinen kuusi, koska runsaspäästöisten pajun ja haavan osuus metsän lehtimassasta on hyvin pieni. Tässä työssä on myös laskettu ensimmäinen arvio metsän seskviterpeenipäästöistä. Seskviterpeenipäästöt alkavat Juhannuksen jälkeen ja ovat kasvukauden aikana samaa suuruusluokkaa kuin isopreenipäästöt. Vuositasolla Suomen metsien hiilivetypäästöt ovat noin kaksinkertaiset ihmistoiminnasta aiheutuviin päästöihin verrattuna.

Observations of production and transport of NOx formed by energetic particle precipitation in the polar night atmosphere

This work is focused on the effects of energetic particle precipitation of solar or magnetospheric origin on the polar middle atmosphere. The energetic charged particles have access to the atmosphere in the polar areas, where they are guided by the Earth's magnetic field. The particles penetrate down to 20-100 km altitudes (stratosphere and mesosphere) ionising the ambient air. This ionisation leads to production of odd nitrogen (NOx) and odd hydrogen species, which take part in catalytic ozone destruction. NOx has a very long chemical lifetime during polar night conditions. Therefore NOx produced at high altitudes during polar night can be transported to lower stratospheric altitudes. Particular emphasis in this work is in the use of both space and ground based observations: ozone and NO2 measurements from the GOMOS instrument on board the European Space Agency's Envisat-satellite are used together with subionospheric VLF radio wave observations from ground stations. Combining the two observation techniques enabled detection of NOx enhancements throughout the middle atmosphere, including tracking the descent of NOx enhancements of high altitude origin down to the stratosphere. GOMOS observations of the large Solar Proton Events of October-November 2003 showed the progression of the SPE initiated NOx enhancements through the polar winter. In the upper stratosphere, nighttime NO2 increased by an order of magnitude, and the effect was observed to last for several weeks after the SPEs. Ozone decreases up to 60 % from the pre-SPE values were observed in the upper stratosphere nearly a month after the events. Over several weeks the GOMOS observations showed the gradual descent of the NOx enhancements to lower altitudes. Measurements from years 2002-2006 were used to study polar winter NOx increases and their connection to energetic particle precipitation. NOx enhancements were found to occur in a good correlation with both increased high-energy particle precipitation and increased geomagnetic activity. The average wintertime polar NOx was found to have a nearly linear relationship with the average wintertime geomagnetic activity. The results from this thesis work show how important energetic particle precipitation from outside the atmosphere is as a source of NOx in the middle atmosphere, and thus its importance to the chemical balance of the atmosphere.

Hybrid Modelling of Titan's Interaction with the Magnetosphere of Saturn

In this dissertation we study the interaction between Saturn's moon Titan and the magnetospheric plasma and magnetic field. The method of research is a three-dimensional computer simulation model, that is used to simulate this interaction. The simulation model used is a hybrid model. Hybrid models enable individual tracking or tracing of ions and also take into account the particle motion in the propagation of the electromagnetic fields. The hybrid model has been developed at the Finnish Meteorological Institute. This thesis gives a general description of the effects that the solar wind has on Earth and other planets of our solar system. Planetary satellites can also have similar interactions with the solar wind but also with the plasma flows of planetary magnetospheres. Titan is clearly the largest among the satellites of Saturn and also the only known satellite with a dense atmosphere. It is the atmosphere that makes Titan's plasma interaction with the magnetosphere of Saturn so unique. Nevertheless, comparisons with the plasma interactions of other solar system bodies are valuable. Detecting charged plasma particles requires in situ measurements obtainable through scientific spacecraft. The Cassini mission has been one of the most remarkable international efforts in space science. Since 2004 the measurements and images obtained from instruments onboard the Cassini spacecraft have increased the scientific knowledge of Saturn as well as its satellites and magnetosphere in a way no one was probably able to predict. The current level of science on Titan is practically unthinkable without the Cassini mission. Many of the observations by Cassini instrument teams have influenced this research both the direct measurements of Titan as well as observations of its plasma environment. The theoretical principles of the hybrid modelling approach are presented in connection to the broader context of plasma simulations. The developed hybrid model is described in detail: e.g. the way the equations of the hybrid model are solved is shown explicitly. Several simulation techniques, such as the grid structure and various boundary conditions, are discussed in detail as well. The testing and monitoring of simulation runs is presented as an essential routine when running sophisticated and complex models. Several significant improvements of the model, that are in preparation, are also discussed. A main part of this dissertation are four scientific articles based on the results of the Titan model. The Titan model developed during the course of the Ph.D. research has been shown to be an important tool to understand Titan's plasma interaction. One reason for this is that the structures of the magnetic field around Titan are very much three-dimensional. The simulation results give a general picture of the magnetic fields in the vicinity of Titan. The magnetic fine structure of Titan's wake as seen in the simulations seems connected to Alfvén waves an important wave mode in space plasmas. The particle escape from Titan is also a major part of these studies. Our simulations show a bending or turning of Titan's ionotail that we have shown to be a direct result of the basic principles in plasma physics. Furthermore, the ion flux from the magnetosphere of Saturn into Titan's upper atmosphere has been studied. The modelled ion flux has asymmetries that would likely have a large impact in the heating in different parts of Titan's upper atmosphere.

Assessment of insect occurrence in boreal forests based on satellite imagery and field measurements.

The presence/absence data of twenty-seven forest insect taxa (e.g. Retinia resinella, Formica spp., Pissodes spp., several scolytids) and recorded environmental variation were used to investigate the applicability of modelling insect occurrence based on satellite imagery. The sampling was based on 1800 sample plots (25 m by 25 m) placed along the sides of 30 equilateral triangles (side 1 km) in a fragmented forest area (approximately 100 km2) in Evo, S Finland. The triangles were overlaid on land use maps interpreted from satellite images (Landsat TM 30 m multispectral scanner imagery 1991) and digitized geological maps. Insect occurrence was explained using either environmental variables measured in the field or those interpreted from the land use and geological maps. The fit of logistic regression models varied between species, possibly because some species may be associated with the characteristics of single trees while other species with stand characteristics. The occurrence of certain insect species at least, especially those associated with Scots pine, could be relatively accurately assessed indirectly on the basis of satellite imagery and geological maps. Models based on both remotely sensed and geological data better predicted the distribution of forest insects except in the case of Xylechinus pilosus, Dryocoetes sp. and Trypodendron lineatum, where the differences were relatively small in favour of the models based on field measurements. The number of species was related to habitat compartment size and distance from the habitat edge calculated from the land use maps, but logistic regressions suggested that other environmental variables in general masked the effect of these variables in species occurrence at the present scale.