Topographical, structural and geophysical characterization of fracture zones: implications for groundwater flow and vulnerability

The main objective of this study is to evaluate selected geophysical, structural and topographic methods on regional, local, and tunnel and borehole scales, as indicators of the properties of fracture zones or fractures relevant to groundwater flow. Such information serves, for example, groundwater exploration and prediction of the risk of groundwater inflow in underground construction. This study aims to address how the features detected by these methods link to groundwater flow in qualitative and semi-quantitative terms and how well the methods reveal properties of fracturing affecting groundwater flow in the studied sites. The investigated areas are: (1) the Päijänne Tunnel for water-conveyance whose study serves as a verification of structures identified on regional and local scales; (2) the Oitti fuel spill site, to telescope across scales and compare geometries of structural assessment; and (3) Leppävirta, where fracturing and hydrogeological environment have been studied on the scale of a drilled well. The methods applied in this study include: the interpretation of lineaments from topographic data and their comparison with aeromagnetic data; the analysis of geological structures mapped in the Päijänne Tunnel; borehole video surveying; groundwater inflow measurements; groundwater level observations; and information on the tunnel s deterioration as demonstrated by block falls. The study combined geological and geotechnical information on relevant factors governing groundwater inflow into a tunnel and indicators of fracturing, as well as environmental datasets as overlays for spatial analysis using GIS. Geophysical borehole logging and fluid logging were used in Leppävirta to compare the responses of different methods to fracturing and other geological features on the scale of a drilled well. Results from some of the geophysical measurements of boreholes were affected by the large diameter (gamma radiation) or uneven surface (caliper) of these structures. However, different anomalies indicating more fractured upper part of the bedrock traversed by well HN4 in Leppävirta suggest that several methods can be used for detecting fracturing. Fracture trends appear to align similarly on different scales in the zone of the Päijänne Tunnel. For example, similarities of patterns were found between the regional magnetic trends, correlating with orientations of topographic lineaments interpreted as expressions of fracture zones. The same structural orientations as those of the larger structures on local or regional scales were observed in the tunnel, even though a match could not be made in every case. The size and orientation of the observation space (patch of terrain at the surface, tunnel section, or borehole), the characterization method, with its typical sensitivity, and the characteristics of the location, influence the identification of the fracture pattern. Through due consideration of the influence of the sampling geometry and by utilizing complementary fracture characterization methods in tandem, some of the complexities of the relationship between fracturing and groundwater flow can be addressed. The flow connections demonstrated by the response of the groundwater level in monitoring wells to pressure decrease in the tunnel and the transport of MTBE through fractures in bedrock in Oitti, highlight the importance of protecting the tunnel water from a risk of contamination. In general, the largest values of drawdown occurred in monitoring wells closest to the tunnel and/or close to the topographically interpreted fracture zones. It seems that, to some degree, the rate of inflow shows a positive correlation with the level of reinforcement, as both are connected with the fracturing in the bedrock. The following geological features increased the vulnerability of tunnel sections to pollution, especially when several factors affected the same locations: (1) fractured bedrock, particularly with associated groundwater inflow; (2) thin or permeable overburden above fractured rock; (3) a hydraulically conductive layer underneath the surface soil; and (4) a relatively thin bedrock roof above the tunnel. The observed anisotropy of the geological media should ideally be taken into account in the assessment of vulnerability of tunnel sections and eventually for directing protective measures.

Intensity Oscillation in the Corona as Observed during the Total Solar Eclipse of 29 March 2006

We obtained the images of the eastern part of the solar corona in the Fe xiv 530.3 nm (green) and Fe x 637.4 nm (red) coronal emission lines during the total solar eclipse of 29 March 2006 at Manavgat, Antalya, Turkey. The images were obtained using a 35 cm Meade telescope equipped with a Peltier-cooled 2k x 2k CCD and 0.3 nm pass-band interference filters at the rates of 2.95 s (exposure times of 100 ms) and 2.0 s (exposure times of 300 ms) in the Fe xiv and Fe x emission lines,respectively. The analysis of the data indicates intensity variations at some locations with period of strongest power around 27 s for the green line and 20 s for the red line. These results confirm earlier findings of variations in the continuum intensity with periods in the range of 5 to 56 s by Singh et al. (Solar Phys. 170, 235, 1997). The wavelet analysis has been used to identify significant intensity oscillations at all pixels within our field of view. Significant oscillations with high probability estimates were detected for some locations only. These locations seem to follow the boundary of an active region and in the neighborhood, rather than within the loops themselves. These intensity oscillations may be caused by fast magneto-sonic waves in the solar corona and partly account for heating of the plasma in the corona.

Contribution from unresolved discrete sources to the extragalactic gamma-ray background (EGRB)

The origin of the extragalactic gamma-ray background (EGRB) is still an open question, even nearly forty years after its discovery. The emission could originate either from truly diffuse processes or from unresolved point sources. Although the majority of the 271 point sources detected by EGRET (Energetic Gamma Ray Experiment Telescope) are unidentified, of the identified sources, blazars are the dominant candidates. Therefore, unresolved blazars may be considered the main contributor to the EGRB, and many studies have been carried out to understand their distribution, evolution and contribution to the EGRB. Considering that gamma-ray emission comes mostly from jets of blazars and that the jet emission decreases rapidly with increasing jet to line-of-sight angle, it is not surprising that EGRET was not able to detect many large inclination angle active galactic nuclei (AGNs). Though Fermi could only detect a few large inclination angle AGNs during the first three months of its survey, it is expected to detect many such sources in the near future. Since non-blazar AGNs are expected to have higher density as compared to blazars, these could also contribute significantly to the EGRB. In this paper, we discuss contributions from unresolved discrete sources including normal galaxies, starburst galaxies, blazars and off-axis AGNs to the EGRB.

Geant4 Hadronic Cascade Models and CMS Data Analysis : Computational Challenges in the LHC era

This work belongs to the field of computational high-energy physics (HEP). The key methods used in this thesis work to meet the challenges raised by the Large Hadron Collider (LHC) era experiments are object-orientation with software engineering, Monte Carlo simulation, the computer technology of clusters, and artificial neural networks. The first aspect discussed is the development of hadronic cascade models, used for the accurate simulation of medium-energy hadron-nucleus reactions, up to 10 GeV. These models are typically needed in hadronic calorimeter studies and in the estimation of radiation backgrounds. Various applications outside HEP include the medical field (such as hadron treatment simulations), space science (satellite shielding), and nuclear physics (spallation studies). Validation results are presented for several significant improvements released in Geant4 simulation tool, and the significance of the new models for computing in the Large Hadron Collider era is estimated. In particular, we estimate the ability of the Bertini cascade to simulate Compact Muon Solenoid (CMS) hadron calorimeter HCAL. LHC test beam activity has a tightly coupled cycle of simulation-to-data analysis. Typically, a Geant4 computer experiment is used to understand test beam measurements. Thus an another aspect of this thesis is a description of studies related to developing new CMS H2 test beam data analysis tools and performing data analysis on the basis of CMS Monte Carlo events. These events have been simulated in detail using Geant4 physics models, full CMS detector description, and event reconstruction. Using the ROOT data analysis framework we have developed an offline ANN-based approach to tag b-jets associated with heavy neutral Higgs particles, and we show that this kind of NN methodology can be successfully used to separate the Higgs signal from the background in the CMS experiment.

Multiwavelength studies of regolith effects in planetary remote sensing

A large proportion of our knowledge about the surfaces of atmosphereless solar-system bodies is obtained through remote-sensing measurements. The measurements can be carried out either as ground-based telescopic observations or space-based observations from orbiting spacecraft. In both cases, the measurement geometry normally varies during the observations due to the orbital motion of the target body, the spacecraft, etc.. As a result, the data are acquired over a variety of viewing and illumination angles. Surfaces of planetary bodies are usually covered with a layer of loose, broken-up rock material called the regolith whose physical properties affect the directional dependence of remote-sensed measurements. It is of utmost importance for correct interpretation of the remote-sensed data to understand the processes behind this alteration. In the thesis, the multi-angular effects that the physical properties of the regolith have on remote-sensing measurements are studied in two regimes of electromagnetic radiation, visible to near infrared and soft X-rays. These effects are here termed generally the regolith effects in remote sensing. Although the physical mechanisms that are important in these regions are largely different, notable similarities arise in the methodology that is used in the study of the regolith effects, including the characterization of the regolith both in experimental studies and in numerical simulations. Several novel experimental setups have been constructed for the thesis. Alongside the experimental work, theoretical modelling has been carried out, and results from both approaches are presented. Modelling of the directional behaviour of light scattered from a regolith is utilized to obtain shape and spin-state information of several asteroids from telescopic observations and to assess the surface roughness and single-scattering properties of lunar maria from spacecraft observations. One of the main conclusions is that the azimuthal direction is an important factor in detailed studies of planetary surfaces. In addition, even a single parameter, such as porosity, can alter the light scattering properties of a regolith significantly. Surface roughness of the regolith is found to alter the elemental fluorescence line ratios of a surface obtained through planetary soft X-ray spectrometry. The results presented in the thesis are among the first to report this phenomenon. Regolith effects need to be taken into account in the analysis of remote-sensed data, providing opportunities for retrieving physical parameters of the surface through inverse methods.

GMRT observations of the group Holmberg 124: Evolution by tidal forces and ram pressure?

We report new radio continuum and 21 cm HI observations using the Giant Metrewave Radio Telescope (GMRT) of the group Holmberg 124 ( Ho 124) comprising four late-type galaxies, namely NGC 2820, Mrk 108, NGC 2814 and NGC 2805. The three galaxies, NGC 2820, Mrk 108 and NGC 2814 which are closely located in the sky plane have clearly undergone tidal interactions as seen from the various morphological tidal signatures and debris. Moreover we note various features in the group members which we believe might be due to ram pressure. In this paper, we describe four interesting results emerging from our observations: a) detection of the tidal radio continuum bridge at 330 MHz connecting the galaxies NGC 2820+ Mrk 108 with NGC 2814. The radio bridge was discovered at 1465 MHz by van der Hulst & Hummel ( 1985, A& A, 150, 17). We find that the bridge has a fairly steep spectrum with a spectral index alpha(S proportional to nu(alpha)) of - 1.8(-0.2)(+0.3) which is much steeper than the - 0.8 quoted by van der Hulst & Hummel ( 1985); b) detection of other tidal features like the tilted HI and radio continuum disk of NGC 2814, a HI streamer and a radio continuum tail arising from the south of NGC 2814. We also report the detection of a possible tidal dwarf galaxy in HI; c) sharp truncation in the HI distribution in the south of NGC 2820 and in the HI and radio continuum distribution in the north of NGC 2814. The optical disks in both the cases look undisturbed. As pointed out by Davis et al. ( 1997, AJ, 114, 613), ram pressure affects different components of the interstellar medium to varying degrees. Simple estimates of pressure in different components of the interstellar medium ( radio continuum, Ha and HI) in NGC 2820 indicate that ram pressure will significantly influence HI; d) detection of a large one-sided HI loop to the north of NGC 2820. No radio continuum emission or Ha emission is associated with the HI loop. We discuss various scenarios for the origin of this loop including a central starburst, ram pressure stripping and tidal interaction. We do not support the central starburst scenario since the loop is not detected in ionized gas. Using the upper limit on X-ray luminosity of Ho 124 (Mulchaey et al. 2003, ApJS, 145, 39), we estimate an upper limit on the intragroup medium (IGrM) density of 8.8 x 10(-4) cm(-3). For half this electron density, we estimate the ram pressure force of the IGrM to be comparable to the gravitational pull of the disk of NGC 2820. Since tidal interaction has obviously influenced the group, we suggest that the loop could have formed by ram pressure stripping if tidal effects had reduced the surface density of HI in NGC 2820. From the complex observational picture of Ho 124 and the numerical estimates, we suggest that the evolution of the Ho 124 group may be governed by both tidal forces due to the interaction and the ram pressure due to motion of the member galaxies in the IGrM and that the IGrM densities should not be too low (i.e. >= 4 x 10(-4)). However this needs to be verified by further observations.

Starburst in the ultraluminous galaxy Arp 220: Constraints from observations of radio recombination lines and continuum

We present observations of radio recombination lines (RRL) from the starburst galaxy Arp 220 at 8.1 GHz (H92 alpha) and 1.4 GHz (H167 alpha and H165 alpha) and at 84 GHz (H42 alpha), 96 GHz (H40 alpha) and 207 GHz (H31 alpha) using the Very Large Array and the IRAM 30 m telescope, respectively. RRLs were detected at all the frequencies except 1.4 GHz, where a sensitive upper limit was obtained. We also present continuum flux measurements at these frequencies as well as at 327 MHz made with the VLA. The continuum spectrum, which has a spectral index alpha similar to -0.6 (S-nu proportional to nu(alpha)) between 5 and 10 GHz, shows a break near 1.5 GHz, a prominent turnover below 500 MHz, and a flatter spectral index above 50 GHz. We show that a model with three components of ionized gas with different densities and area covering factors can consistently explain both RRL and continuum data. The total mass of ionized gas in the three components is 3.2 x 10(7) M., requiring 3 x 10(5) O5 stars with a total Lyman continuum production rate N-Lyc similar to 1.3 x 10(55) photons s(-1). The ratio of the expected to observed Br alpha and Br gamma fluxes implies a dust extinction A(V) similar to 45 mag. The derived Lyman continuum photon production rate implies a continuous star formation rate (SFR) averaged over the lifetime of OB stars of similar to 240 M yr(-1). The Lyman continuum photon Production rate of similar to 3% associated with the high-density H II regions implies a similar SFR at recent epochs (t < 10(5) yr). An alternative model of high-density gas, which cannot be excluded on the basis of the available data, predicts 10 times higher SFR at recent epochs. If confirmed, this model implies that star formation in Arp 220 consists of multiple starbursts of very high SFR (few times 10(3) M. yr(-1)) and short duration (similar to 10(5) yr). The similarity of IR excess, L-IR/L-Ly alpha similar to 24, in Arp 220 to values observed in starburst galaxies shows that most of the high luminosity of Arp 220 is due to the ongoing starburst rather than to a hidden active galactic nucleus (AGN). A comparison of the IR excesses in Arp 220, the Galaxy, and M33 indicates that the starburst in Arp 220 has an initial mass function that is similar to that in normal galaxies and has a duration longer than 107 yr. If there was no infall of gas during this period, then the star formation efficiency (SFE) in Arp 220 is similar to 50%. The high SFR and SFE in Arp 220 is consistent with their known dependences on mass and density of gas in star-forming regions of normal galaxies.

Position-sensitive silicon strip detector characterization using particle beams

Silicon strip detectors are fast, cost-effective and have an excellent spatial resolution. They are widely used in many high-energy physics experiments. Modern high energy physics experiments impose harsh operation conditions on the detectors, e.g., of LHC experiments. The high radiation doses cause the detectors to eventually fail as a result of excessive radiation damage. This has led to a need to study radiation tolerance using various techniques. At the same time, a need to operate sensors approaching the end their lifetimes has arisen. The goal of this work is to demonstrate that novel detectors can survive the environment that is foreseen for future high-energy physics experiments. To reach this goal, measurement apparatuses are built. The devices are then used to measure the properties of irradiated detectors. The measurement data are analyzed, and conclusions are drawn. Three measurement apparatuses built as a part of this work are described: two telescopes measuring the tracks of the beam of a particle accelerator and one telescope measuring the tracks of cosmic particles. The telescopes comprise layers of reference detectors providing the reference track, slots for the devices under test, the supporting mechanics, electronics, software, and the trigger system. All three devices work. The differences between these devices are discussed. The reconstruction of the reference tracks and analysis of the device under test are presented. Traditionally, silicon detectors have produced a very clear response to the particles being measured. In the case of detectors nearing the end of their lifefimes, this is no longer true. A new method benefitting from the reference tracks to form clusters is presented. The method provides less biased results compared to the traditional analysis, especially when studying the response of heavily irradiated detectors. Means to avoid false results in demonstrating the particle-finding capabilities of a detector are also discussed. The devices and analysis methods are primarily used to study strip detectors made of Magnetic Czochralski silicon. The detectors studied were irradiated to various fluences prior to measurement. The results show that Magnetic Czochralski silicon has a good radiation tolerance and is suitable for future high-energy physics experiments.

Polarization studies in electromagnetic scattering by small Solar-system particles

In remote-sensing studies, particles that are comparable to the wavelength exhibit characteristic features in electromagnetic scattering, especially in the degree of linear polarization. These features vary with the physical properties of the particles, such as shape, size, refractive index, and orientation. In the thesis, the direct problem of computing the unknown scattered quantities using the known properties of the particles and the incident radiation is solved at both optical and radar spectral regions in a unique way. The internal electromagnetic fields of wavelength-scale particles are analyzed by using both novel and established methods to show how the internal fields are related to the scattered fields in the far zone. This is achieved by using the tools and methods that were developed specifically to reveal the internal field structure of particles and to study the mechanisms that relate the structure to the scattering characteristics of those particles. It is shown that, for spherical particles, the internal field is a combination of a forward propagating wave with the apparent wavelength determined by the refractive index of the particle, and a standing wave pattern with the apparent wavelength the same as for the incident wave. Due to the surface curvature and dielectric nature of the particle, the incident wave front undergoes a phase shift, and the resulting internal wave is focused mostly at the forward part of the particle similar to an optical lens. This focusing is also seen for irregular particles. It is concluded that, for both spherical and nonspherical particles, the interference at the far field between the partial waves that originate from these concentrated areas in the particle interior, is responsible for the specific polarization features that are common for wavelength-scale particles, such as negative values and local extrema in the degree of linear polarization, asymmetry of the phase function, and enhancement of intensity near the backscattering direction. The papers presented in this thesis solve the direct problem for particles with both simple and irregular shapes to demonstrate that these interference mechanisms are common for all dielectric wavelength-scale particles. Furthermore, it is shown that these mechanisms can be applied to both regolith particles in the optical wavelengths and hydrometeors at microwave frequencies. An advantage from this kind of study is that it does not matter whether the observation is active (e.g., polarimetric radar) or passive (e.g., optical telescope). In both cases, the internal field is computed for two mutually perpendicular incident polarizations, so that the polarization characteristics can then be analyzed according to the relation between these fields and the scattered far field.

A Multi frequency study of Radio Intensity Variations For Active Galactic Nuclei of Different Optical Classes

A study of radio intensity variations at seven frequencies in the range 0.3 to 90 GHz for compact extragalactic radio sources classified as BL Lacs and high- and low-optical polarization quasars (HPQs and LPQs) is presented. This include the results of flux-density monitoring of 33 compact sources for three years at 327 MHz with the Ooty Synthesis Radio Telescope. The degrees of 'short-term' (tau less than about 1 yr) variability for the three optical types are found to be indistinguishable at low frequencies (less than 1 GHz), pointing to an extrinsic origin for the low-frequency variability. At high frequencies, a distinct dependence on optical type is present, the variability increasing from LPQs, through HPQs to BL Lacs. This trend persists even when only sources with ultra-flat radio spectra (alpha greater than -0.2) are considered. Implications of this for the phenomenon of high-frequency variability and the proposed unification schemes for different optical types of active galactic nuclei are discussed.

Optimum weighting function in bispectral analysis in speckle interferometry - Binary star parity detection

In order to describe the atmospheric turbulence which limits the resolution of long-exposure images obtained using ground-based large telescopes, a simplified model of a speckle pattern, reducing the complexity of calculating field-correlations of very high order, is presented. Focal plane correlations are used instead of correlations in the spatial frequency domain. General tripple correlations for a point source and for a binary are calculated and it is shown that they are not a strong function of the binary separation. For binary separations close to the diffraction limit of the telescope, the genuine triple correlation technique ensures a better SNR than the near-axis Knox-Thompson technique. The simplifications allow a complete analysis of the noise properties at all levels of light.

Gain calibration of CCD systems at VBO

The system gain of two CCD systems in regular use at the Vainu Bappu Observatory, Kavalur, is determined at a few gain settings. The procedure used for the determination of system gain and base-level noise is described in detail. The Photometrics CCD system at the 1-m reflector uses a Thomson-CSF TH 7882 CDA chip coated for increased ultraviolet sensitivity. The gain is programme-selected through the parameter 'cgain' varying between 0 and 4095 in steps of 1. The inverse system gain for this system varies almost linearly from 27.7 electrons DN-1 at cgain = 0 to 1.5 electrons DN-1 at cgain = 500. The readout noise is less than or similar 11 electrons at cgain = 66. The Astromed CCD system at 2.3-m Vainu Bappu Telescope uses a GEC P8603 chip which is also coated for enhanced ultraviolet sensitivity. The amplifier gain is selected in discrete steps using switches in the controller. The inverse system gain is 4.15 electrons DN-1 at the gain setting of 9.2, and the readout noise approximately 8 electrons.

SiO Maser Emission and the Intrinsic Properties of Mira Variables

This paper presents observations of SiO maser emission from 161 Mira variables distributed over a wide range of intrinsic parameters like spectral type, bolometric magnitude and amplitude of pulsation. The observations were made at 86.243 GHz, using the 10.4 m millimeter-wave telescope of the Raman Research Institute at Bangalore, India. These are the first observations made using this telescope. From these observations, we have established that the maser emission is restricted to Miras having mean spectral types between M6 and M10. The infrared period-luminosity relation for Mira variables is used to calculate their distances and hence estimate their maser luminosities from the observed fluxes. The maser luminosity is found to be correlated with the bolometric magnitude of the Mira variable. On an H-R diagram, the masing Mira variables are shown to lie in a region distinct from that for the non-masing ones.

Carbon recombination lines between 34.5 and 770 MHz toward Cassiopeia A

We present observations of low-frequency recombination lines of carbon toward Cas A near 34.5 MHz (n similar to 575) using the Gauribidanur radio telescope and near 560 MHz (n similar to 225) and 770 MHz (n similar to 205) using the NRAO 140 foot (43 m) telescope in Greenbank. We also present high angular resolution (1') observations of the C270 alpha line near 332 MHz using the Very Large Array in B-configuration. A high signal-to-noise ratio spectrum is obtained at 34.5 MHz, which clearly shows a Voigt profile with distinct Lorentzian wings, resulting from significant pressure and radiation broadening at such high quantum numbers. The emission lines detected near 332, 550, and 770 MHz, on the other hand, are narrow and essentially Doppler-broadened. The measured Lorentzian width at 34.5 MHz constrains the allowed combinations of radiation temperature, electron density, and electron temperature in the line-forming region. Radiation broadening at 34.5 MHz places a lower limit of 115 pc on the separation between Cas A and the line-forming clouds. Modeling the variation in the integrated line-to-continuum ratio with frequency indicates that the region is likely to be associated with the cold atomic hydrogen component of the interstellar medium, and the physical properties of this region are likely to be T-e = 75 K, n(e) = 0.02 cm(-3), T-R100 = 3200 K, and n(H) T-e = 10,000 cm(-3) K. Comparison of the distribution of the C270 alpha recombination line emission across Cas A with that of (CO)-C-12 and H I also supports the above conclusion.

Mapping of large scale 158 mu m [CII] line emission: Orion A

We present the first results of an observational programme undertaken to map the fine structure line emission of singly ionized carbon ([ CII] 157 : 7409 mum) over extended regions using a Fabry Perot spectrometer newly installed at the focal plane of a 100 cm balloon- borne far- infrared telescope. This new combination of instruments has a velocity resolution of similar to 200 km s(-1) and an angular resolution of 1.'5. During the first flight, an area of 30' x 15' in Orion A was mapped. These observations extend over a larger area than previous observations, the map is fully sampled and the spectral scanning method used enables reliable estimation of the continuum emission at frequencies adjacent to the [ CII] line. The total [ CII] line luminosity, calculated by considering up to 20% of the maximum line intensity is 0.04% of the luminosity of the far- infrared continuum. We have compared the [ CII] intensity distribution with the velocity- integrated intensity distributions of (CO)-C-13(1- 0), CI(1- 0) and CO( 3- 2) from the literature. Comparison of the [ CII], [ CI] and the radio continuum intensity distributions indicates that the largescale [ CII] emission originates mainly from the neutral gas, except at the position of M 43, where no [ CI] emission corresponding to the [ CII] emission is seen. Substantial part of the [ CII] emission from here originates from the ionized gas. The observed line intensities and ratios have been analyzed using the PDR models by Kaufman et al. ( 1999) to derive the incident UV flux and volume density at a few selected positions. The models reproduce the observations reasonably well at most positions excepting the [ CII] peak ( which coincides with the position of theta(1) Ori C). Possible reason for the failure could be the simplifying assumption of a homogeneous plane parallel slab in place of a more complicated geometry.