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.

Mixed pixel analysis for flood mapping using extended support vector machine

This paper addresses the challenges of flood mapping using multispectral images. Quantitative flood mapping is critical for flood damage assessment and management. Remote sensing images obtained from various satellite or airborne sensors provide valuable data for this application, from which the information on the extent of flood can be extracted. However the great challenge involved in the data interpretation is to achieve more reliable flood extent mapping including both the fully inundated areas and the 'wet' areas where trees and houses are partly covered by water. This is a typical combined pure pixel and mixed pixel problem. In this paper, an extended Support Vector Machines method for spectral unmixing developed recently has been applied to generate an integrated map showing both pure pixels (fully inundated areas) and mixed pixels (trees and houses partly covered by water). The outputs were compared with the conventional mean based linear spectral mixture model, and better performance was demonstrated with a subset of Landsat ETM+ data recorded at the Daly River Basin, NT, Australia, on 3rd March, 2008, after a flood event.

Spectral unmixing with estimated adaptive endmember index using Extended Support Vector Machine

The most difficult operation in the flood inundation mapping using optical flood images is to separate fully inundated areas from the ‘wet’ areas where trees and houses are partly covered by water. This can be referred as a typical problem the presence of mixed pixels in the images. A number of automatic information extraction image classification algorithms have been developed over the years for flood mapping using optical remote sensing images. Most classification algorithms generally, help in selecting a pixel in a particular class label with the greatest likelihood. However, these hard classification methods often fail to generate a reliable flood inundation mapping because the presence of mixed pixels in the images. To solve the mixed pixel problem advanced image processing techniques are adopted and Linear Spectral unmixing method is one of the most popular soft classification technique used for mixed pixel analysis. The good performance of linear spectral unmixing depends on two important issues, those are, the method of selecting endmembers and the method to model the endmembers for unmixing. This paper presents an improvement in the adaptive selection of endmember subset for each pixel in spectral unmixing method for reliable flood mapping. Using a fixed set of endmembers for spectral unmixing all pixels in an entire image might cause over estimation of the endmember spectra residing in a mixed pixel and hence cause reducing the performance level of spectral unmixing. Compared to this, application of estimated adaptive subset of endmembers for each pixel can decrease the residual error in unmixing results and provide a reliable output. In this current paper, it has also been proved that this proposed method can improve the accuracy of conventional linear unmixing methods and also easy to apply. Three different linear spectral unmixing methods were applied to test the improvement in unmixing results. Experiments were conducted in three different sets of Landsat-5 TM images of three different flood events in Australia to examine the method on different flooding conditions and achieved satisfactory outcomes in flood mapping.

Mixed pixel analysis and assessment for flood mapping

The most difficult operation in flood inundation mapping using optical flood images is to map the ‘wet’ areas where trees and houses are partly covered by water. This can be referred to as a typical problem of the presence of mixed pixels in the images. A number of automatic information extracting image classification algorithms have been developed over the years for flood mapping using optical remote sensing images, with most labelling a pixel as a particular class. However, they often fail to generate reliable flood inundation mapping because of the presence of mixed pixels in the images. To solve this problem, spectral unmixing methods have been developed. In this thesis, methods for selecting endmembers and the method to model the primary classes for unmixing, the two most important issues in spectral unmixing, are investigated. We conduct comparative studies of three typical spectral unmixing algorithms, Partial Constrained Linear Spectral unmixing, Multiple Endmember Selection Mixture Analysis and spectral unmixing using the Extended Support Vector Machine method. They are analysed and assessed by error analysis in flood mapping using MODIS, Landsat and World View-2 images. The Conventional Root Mean Square Error Assessment is applied to obtain errors for estimated fractions of each primary class. Moreover, a newly developed Fuzzy Error Matrix is used to obtain a clear picture of error distributions at the pixel level. This thesis shows that the Extended Support Vector Machine method is able to provide a more reliable estimation of fractional abundances and allows the use of a complete set of training samples to model a defined pure class. Furthermore, it can be applied to analysis of both pure and mixed pixels to provide integrated hard-soft classification results. Our research also identifies and explores a serious drawback in relation to endmember selections in current spectral unmixing methods which apply fixed sets of endmember classes or pure classes for mixture analysis of every pixel in an entire image. However, as it is not accurate to assume that every pixel in an image must contain all endmember classes, these methods usually cause an over-estimation of the fractional abundances in a particular pixel. In this thesis, a subset of adaptive endmembers in every pixel is derived using the proposed methods to form an endmember index matrix. The experimental results show that using the pixel-dependent endmembers in unmixing significantly improves performance.

Hall thruster sequencer using Intel-8051 microcontroller for ion thrust measurement

Hall thrusters, such as Stationary Plasma Thruster (SPT), have been widely used on board modern satellites placed in geo-synchronows orbits for reasons such as orbit maintenance, repositioning and attitude control. In order to study the performance of the stationary plasma thruster, the thrust produced by it has been measured, using a thrust balance with strain gauge sensors under vacuum conditions, by activating the thruster. This activation of thruster has been carried out by switching ON and switching OFF of the necessary power supplies and control of other feed system such as the propellant flow in a particular sequence. Hitherto, these operations were done manually in the required sequence. This paper reports the attempt made to automate the sequential operation of the power supplies and the necessary control valves of the feed system using Intel 8051 microcontroller. This automation has made thrust measurements easier and more sophisticated.

Variations in the satellited chromosomes ofCicer arietinum Linn.

The diversity in the morphology of the SAT-chromosomes of C. arietinum in different roots are analysed in the context of similar variations observed in the different cells of one of the roots. Intergrades between the two pairs of chromosomes with normal satellites and those with tandem ones were observed in anaphases. The differences in the number as well as the size of the SAT-grains of replicated chromosomes indicate that non-reciprocal translocations between sister chromosomes may be of parts of SAT-grains and not of entire ones.

Orbital evolution and orbital phase resolved spectroscopy of the HMXB pulsar 4U 1538-52 with RXTE-PCA and BeppoSAX

We report here results from detailed timing and spectral studies of the high mass X-ray binary pulsar 4U 1538-52 over several binary periods using observations made with the Rossi X-ray Timing Explorer (RXTE) and BeppoSAX satellites. Pulse timing analysis with the 2003 RXTE data over two binary orbits confirms an eccentric orbit of the system. Combining the orbitial parameters determined from this observation with the earlier measurements we did not find any evidence of orbital decay in this X-ray binary. We have carried out orbital phase resolved spectroscopy to measure changes in the spectral parameters with orbital phase, particularly the absorption column density and the iron line flux. The RXTE-PCA spectra in the 3-20 keV energy range were fitted with a power law and a high energy cut-off along with a Gaussian line at similar to 6.4 keV, whereas the BeppoSAX spectra needed only a power law and Gaussian emission line at similar to 6.4keV in the restricted energy range of 0.3-10.0 keV. An absorption along the line of sight was included for both the RXTE and BeppoSAX data. The variation of the free spectral parameters over the binary orbit was investigated and we found that the variation of the column density of absorbing material in the line of sight with orbital phase is in reasonable agreement with a simple model of a spherically symmetric stellar wind from the companion star.

Analysis of solar radiation pressure induced coupled liberations of a gravity stabilized axisymmetric satellite

This paper presents an analysis of solar radiation pressure induced coupled librations of gravity stabilized cylindrical spacecraft with a special reference to geostationary communication satellites. The Lagrangian approach is used to obtain the corresponding equations of motion. The solar induced torques are assumed to be free of librational angles and are represented by their Fourier expansion. The response and periodic solutions are obtained through linear and nonlinear analyses, using the method of harmonic balance in the latter case. The stability conditions are obtained using Routh-Hurwitz criteria. To establish the ranges of validity the analytic response is compared with the numerical solution. Finally, values of the system parameters are suggested to make the satellite behave as desired. Among these is a possible approach to subdue the solar induced roll resonance. It is felt that the approximate analysis presented here should significantly reduce the computational efforts involved in the design and stability analysis of the systems.

Enantiodiscrimination and extraction of short and long range homo- and hetero-nuclear residual dipolar couplings by a spin selective correlation experiment

A two dimensional correlation experiment for the measurement of short and long range homo- and hetero- nuclear residual dipolar couplings (RDCs) from the broad and featureless proton NMR spectra including C-13 satellites is proposed. The method employs a single natural abundant C-13 spin as a spy nucleus to probe all the coupled protons and permits the determination of RDCs of negligible strengths. The technique has been demonstrated for the study of organic chiral molecules aligned in chiral liquid crystal, where additional challenge is to unravel the overlapped spectrum of enantiomers. The significant advantage of the method is demonstrated in better chiral discrimination using homonuclear RDCs as additional parameters. (C) 2010 Elsevier B.V. All rights reserved.

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.

Light scattering in random media with wavelength-scale structures : astronomical and industrial applications

Light scattering, or scattering and absorption of electromagnetic waves, is an important tool in all remote-sensing observations. In astronomy, the light scattered or absorbed by a distant object can be the only source of information. In Solar-system studies, the light-scattering methods are employed when interpreting observations of atmosphereless bodies such as asteroids, atmospheres of planets, and cometary or interplanetary dust. Our Earth is constantly monitored from artificial satellites at different wavelengths. With remote sensing of Earth the light-scattering methods are not the only source of information: there is always the possibility to make in situ measurements. The satellite-based remote sensing is, however, superior in the sense of speed and coverage if only the scattered signal can be reliably interpreted. The optical properties of many industrial products play a key role in their quality. Especially for products such as paint and paper, the ability to obscure the background and to reflect light is of utmost importance. High-grade papers are evaluated based on their brightness, opacity, color, and gloss. In product development, there is a need for computer-based simulation methods that could predict the optical properties and, therefore, could be used in optimizing the quality while reducing the material costs. With paper, for instance, pilot experiments with an actual paper machine can be very time- and resource-consuming. The light-scattering methods presented in this thesis solve rigorously the interaction of light and material with wavelength-scale structures. These methods are computationally demanding, thus the speed and accuracy of the methods play a key role. Different implementations of the discrete-dipole approximation are compared in the thesis and the results provide practical guidelines in choosing a suitable code. In addition, a novel method is presented for the numerical computations of orientation-averaged light-scattering properties of a particle, and the method is compared against existing techniques. Simulation of light scattering for various targets and the possible problems arising from the finite size of the model target are discussed in the thesis. Scattering by single particles and small clusters is considered, as well as scattering in particulate media, and scattering in continuous media with porosity or surface roughness. Various techniques for modeling the scattering media are presented and the results are applied to optimizing the structure of paper. However, the same methods can be applied in light-scattering studies of Solar-system regoliths or cometary dust, or in any remote-sensing problem involving light scattering in random media with wavelength-scale structures.

Discriminating between natural and anthropogenic aerosol contributions to optical depth over Afro-Asian region

To accurately assess the impact of anthropogenic aerosols on climate, spatial and temporal distribution of its radiative properties is essential. The first step towards separating the radiative impact of natural aerosol from its anthropogenic counterparts is to gather information on natural aerosols. In this paper, we have used data from multiple satellites to derive the anthropogenic aerosol fraction (AAF) over the Afro-Asian region. The AAF was largest during the pre-monsoon season (May-June) and lowest during winter. We have shown that over desert locations the AAF was unexpectedly large (>0.4) and the regionally (and annually) averaged anthropogenic fraction over the Afro-Asian region was 0.54 +/- 0.12. Copyright (C) 2010 Royal Meteorological Society

Demixing of severely overlapped H-1 NMR resonances and interpretation of anomalous intensity pattern of dipolar coupled A(3) spins in a weakly aligning medium

We report a single C-13 spin edited selective proton-proton correlation experiment to decipher overcrowded 13C coupled proton NMR spectra of weakly dipolar coupled spin systems. The experiment unravels the masked C-13 satellites in proton spectrum and permits the measurement of one bond carbon-proton residual dipolar couplings in I3S and for each diastereotopic proton in I2S groups. It also provides all the possible homonuclear proton-proton residual couplings which are otherwise difficult to extract from the broad and featureless one dimensional H-1 spectrum, in addition to enantiodifferentiation in a chiral molecule. Employment of heteronuclear (C-13) decoupling in the evolution period results in complete demixing of overlapped signals from enantiomers. The observed anomalous intensity pattern in strongly dipolar coupled methyl protons in methyl selective correlation experiment has been interpreted using polarization operator formalism. (C) 2010 Elsevier Inc. All rights reserved.

The Radiation, snow characteristics and albedo at summit (RASCALS) expedition report

The RASCALS expedition spent over three weeks at the Summit camp research station near the top of the Greenland Ice Sheet during polar summer 2010. During this time, detailed measurements of the physical and optical properties of Arctic perennial snow were carried out concurrently with snow albedo and reflectance measurements. Favorable weather conditions during the campaign enabled the collection of a large dataset on Arctic snow albedo and associated quantities for use in developing and validating remote sensing algorithms for snow albedo using satellites. This report provides a description of the measurements and conditions during the campaign. RASCALS-retkikunnan tehtävä oli tutkia Grönlannin mannerjäätikön lumen fysikaalisia ja optisia ominaisuuksia sekä Auringon valon vuorovaikutusta lumen kanssa. Retikunta vietti hieman yli kolme viikkoa mannerjäätikön keskellä sijaitsevalla Summit Camp-tutkimusasemalla tehden mittauksia. Sääolot suosivat kampanjaa, jonka seurauksena onnistuttiin keräämään laaja ja monipuolinen tietoaineisto mannerjäätikön lumen pintakerroksesta ja eritoten lumen heijastavuuden (albedon)käyttäytymisestä. Aineisto on hyödyllinen kehitettäessä ja varmennettaessa lumen albedon kaukokartoitusmenetelmiä satelliiteilla.