Advanced modelling of time domain electromagnetic data with updated hydrogeological interpretations

Several countries have acquired, over the past decades, large amounts of area covering Airborne Electromagnetic data. Contribution of airborne geophysics has dramatically increased for both groundwater resource mapping and management proving how those systems are appropriate for large-scale and efficient groundwater surveying. We start with processing and inversion of two AEM dataset from two different systems collected over the Spiritwood Valley Aquifer area, Manitoba, Canada respectively, the AeroTEM III (commissioned by the Geological Survey of Canada in 2010) and the “Full waveform VTEM” dataset, collected and tested over the same survey area, during the fall 2011. We demonstrate that in the presence of multiple datasets, either AEM and ground data, due processing, inversion, post-processing, data integration and data calibration is the proper approach capable of providing reliable and consistent resistivity models. Our approach can be of interest to many end users, ranging from Geological Surveys, Universities to Private Companies, which are often proprietary of large geophysical databases to be interpreted for geological and\or hydrogeological purposes. In this study we deeply investigate the role of integration of several complimentary types of geophysical data collected over the same survey area. We show that data integration can improve inversions, reduce ambiguity and deliver high resolution results. We further attempt to use the final, most reliable output resistivity models as a solid basis for building a knowledge-driven 3D geological voxel-based model. A voxel approach allows a quantitative understanding of the hydrogeological setting of the area, and it can be further used to estimate the aquifers volumes (i.e. potential amount of groundwater resources) as well as hydrogeological flow model prediction. In addition, we investigated the impact of an AEM dataset towards hydrogeological mapping and 3D hydrogeological modeling, comparing it to having only a ground based TEM dataset and\or to having only boreholes data.

Biogenic volatile organic compounds in tropical, temperate and boreal forest ecosystems

Biogene flüchtige organische Verbindungen (BFOV) werden in großen Mengen aus terrestrischenrnÖkosystemen, insbesondere aus Wäldern und Wiesen, in die untere Troposphäre emittiert. Austausch-rnFlüsse von BFOVs sind in der troposphärischen Chemie wichtig, weil sie eine bedeutende Rolle in derrnOzon- und Aerosolbildung haben. Trotzdem bleiben die zeitliche und räumliche Änderung der BFOVrnEmissionen und ihre Rolle in Bildung und Wachstum von Aerosolen ungewiss.rnDer Fokus dieser Arbeit liegt auf der in-situ Anwendung der Protonen Transfer ReaktionsrnMassenspektrometrie (PTR-MS) und der Messung von biogenen flüchtigen organischen Verbindungen inrnnordländischen, gemäßigten und tropischen Waldökosystemen während drei unterschiedlicherrnFeldmesskampagnen. Der Hauptvorteil der PTR-MS-Technik liegt in der hohen Messungsfrequenz,rnwodurch eine eventuelle Änderung in der Atmosphäre durch Transport, Vermischung und Chemiernonline beobachtet werden kann. Die PTR-MS-Messungen wurden zweimal am Boden aus und einmalrnvon einem Forschungsflugzug durchgeführt.rnIn Kapitel 3 werden die PTR-MS-Daten, gesammelt während der Flugmesskampagne über demrntropischen Regenwald, vorgelegt. Diese Studie zeigt den Belang der Grenzschichtdynamik für diernVerteilung von Spurengasen mittels eines eindimensionalen Säule - Chemie und KlimaModells (SCM).rnDer Tagesablauf von Isopren zeigte zwischen 14:00 und 16:15 Uhr lokaler Zeit einen Mittelwert vonrn5.4 ppbv auf der Höhe der Baumspitzen und von 3.3 ppbv über 300 m Höhe. Dies deutet darauf hin, dassrnsowohl der turbulente Austausch als auch die hohe Reaktionsfähigkeit von Isopren mit den OxidantienrnOH und Ozon eine wichtige Rolle spielen. Nach dem Ausgleich von chemischen Verlusten undrnEntrainment (Ein- und Ausmischung von Luft an der Grenzschicht), wurde ein Fluss vonrn8.4 mg Isopren m-2h-1 unter teilweise bewölkten Bedingungen für den tropischen Regenwald in derrnGuyanregion abgeschätzt. Dies entspricht einem täglichen Emissionsfluss von 28 mg Isopren prornQuadratmeter.rnIm Kapitel 4 werden die Messungen, welche auf einer Hügellage in gemäßigter Breite inrnsüddeutschland stattgefunden haben, diskutiert. Bei diesem Standort ist die Grenzschicht nachts unter diernStandorthöhe abgefallen, was den Einsatzort von Emissionen abgesondert hatte. Während diernGrenzschicht morgens wieder über die Höhe des Einsatzortes anstieg, konnten die eingeschlossenenrnnächtlichen Emissionen innerhalb der bodennahen Schicht beobachtet werden. Außerdem wurde einrndeutlicher Anstieg von flüchtigen organischen Verbindungen gemessen, wenn die Luftmassen überrnMünchen geführt wurden oder wenn verschmutzte Luftmassen aus dem Po-Tal über die Alpen nachrnDeutschland transportiert wurden. Daten von dieser Kampagne wurden genutzt, um die Änderungen inrndem Mischungsverhältnis der flüchtigen organischen Verbindungen, verbunden mit dem Durchfluss vonrnwarmen und kalten Wetterfronten sowie bei Regen zu untersuchen.rnIm Kapitel 5 werden PTR-MS-Messungen aus dem nördlichen Nadelwaldgürtel beschrieben. Starkernnächtliche Inversionen mit einer niedrigen Windgeschwindigkeit fingen die Emissionen vonrnnahegelegenen Kiefernwäldern und andere BFOV-Quellen ab, was zu hohen nächtlichen BFOVMischverhältnissenrnführte. Partikelereignisse wurden für Tag und Nacht detailliert analysiert. Diernnächtlichen Partikelereignisse erfolgten synchron mit starken extremen von Monoterpenen, obwohl dasrnzweite Ereignis Kernbildung einschloss und nicht mit Schwefelsäure korrelierte. Die MonoterpenrnMischungsverhältnisse von über 16 ppbv waren unerwartet hoch für diese Jahreszeit. NiedrigernWindgeschwindigkeiten und die Auswertung von Rückwärtstrajektorien deuten auf eine konzentrierternQuelle in der Nähe von Hyytiälä hin. Die optische Stereoisomerie von Monoterpenen hat bestätigt, dassrndie Quelle unnatürlich ist, da das Verhältnis von [(+)-α-pinen]/[(−)-α-pinen] viel höher ist als dasrnnatürliches Verhältnis der beiden Enantiomeren.

The development and deployment of FOTOS: a fast observation of trace organics system

Flüchtige organische Bestandteile (engl.: VOC) sind in der Atmosphäre in Spuren vorhanden, spielen aber trotzdem eine wichtige Rolle in der Luftchemie: sie beeinflussen das Ozon der Troposphäre, städtischen Smog, Oxidationskapazität und haben direkte und indirekte Auswirkungen auf die globale Klimaveränderung. Eine wichtige Klasse der VOC sind die Nicht-Methan-Kohlenwasserstoffe (engl.: NMHC), die überwiegend von anthropogenen Quellen kommen. Aus diesem Grund ist für Luftchemiker ein Messinstrument nötig, das die VOC, die NMHC eingeschlossen, mit einer höheren Zeitauflösung misst, besonders für Echtzeitmessungen an Bord eines Forschungsflugzeuges. Dafür wurde das System zur schnellen Beobachtung von organischen Spuren (engl.: FOTOS) entworfen, gebaut für den Einsatz in einem neuen Wissenschaftlichen Flugzeug, das in großen Höhen und über weite Strecken fliegt, genannt HALO. In der Folge wurde FOTOS in zwei Messkampagnen am Boden getestet. FOTOS wurde entworfen und gebaut mit einem speziell angefertigten, automatisierten, kryogenen Probensystem mit drei Fallen und einem angepassten, erworbenen schnellen GC-MS. Ziel dieses Aufbaus war es, die Vielseitigkeit zu vergrößern und das Störungspotential zu verringern, deshalb wurden keine chemischen Trocknungsmittel oder adsorbierenden Stoffe verwendet. FOTOS erreichte eine Probenfrequenz von 5.5 Minuten, während es mindestens 13 verschiedene C2- bis C5-NMHC maß. Die Drei-Sigma-Detektionsgrenze für n- und iso-Pentan wurde als 2.6 und 2.0 pptv ermittelt, in dieser Reihenfolge. Labortests bestätigten, dass FOTOS ein vielseitiges, robustes, hochautomatisiertes, präzises, genaues, empfindliches Instrument ist, geeignet für Echtzeitmessungen von VOC in Probenfrequenzen, die angemessen sind für ein Forschungsflugzeug wie HALO. Um die Leistung von FOTOS zu bestätigen, wurde vom 26. Januar bis 4. Februar 2010 ein Zwischenvergleich gemacht mit dem GC-FID-System am Meteorologischen Observatorium Hohenpeißenberg, einer WMO-GAW-globalen Station. Dreizehn verschiedene NMHC wurden innerhalb des Rahmens der GWA Data Quality Objectives (DQO) analysiert und verglichen. Mehr als 80% der Messungen von sechs C3- bis C5-NMHC erfüllten diese DQO. Diese erste Messkampagne im Feld hob die Robustheit und Messgenauigkeit von FOTOS hervor, zusätzlich zu dem Vorteil der höheren Probenfrequenz, sogar in einer Messung am Boden. Um die Möglichkeiten dieses Instrumentes im Feld zu zeigen, maß FOTOS ausgewählte leichte NMHC während einer Messkampagne im Borealen Waldgebiet, HUMPPA-COPEC 2010. Vom 12. Juli bis zum 12. August 2010 beteiligte sich eine internationale Gruppe von Instituten und Instrumenten an Messungen physikalischer und chemischer Größen der Gas- und Partikelphasen der Luft über dem Borealen Wald an der SMEAR II-Station nahe Hyyttiälä, Finnland. Es wurden mehrere Hauptpunkte von Interesse im Mischungsverhältnis der Alkane und im Isomerenverhätnis von Pentan identifiziert, insbesondere sehr unterschiedliche Perioden niedriger und hoher Variabilität, drei Rauchschwaden von Biomassen-Verbrennung von russischen Waldbränden und zwei Tage mit extrem sauberer Luft aus der Polarregion. Vergleiche der NMHC mit anderen anthropogenen Indikatoren zeigten mehrere Quellen anthropogener Einflüsse am Ort auf und erlaubten eine Unterscheidung zwischen lokalen und weiter entfernten Quellen. Auf einen minimalen natürlichen Beitrag zum 24h-Kreislauf von NOx wurde geschlussfolgert aus der Korrelation von NOx mit Alkanen. Altersschätzungen der Luftmassen durch das Isomerenverhältnis von Pentan wurden erschwert durch sich verändernde Verhältnisse der Quellen und durch Besonderheiten der Photochemie während des Sommers im hohen Norden. Diese Messungen zeigten den Wert des Messens leichter NMHC, selbst in abgelegenen Regionen, als einen zusätzlichen spezifischen Marker von anthropogenem Einfluss.

A multi-instrumental approach to the study of equatorial ionosphere over South America

An extensive study of the morphology and the dynamics of the equatorial ionosphere over South America is presented here. A multi parametric approach is used to describe the physical characteristics of the ionosphere in the regions where the combination of the thermospheric electric field and the horizontal geomagnetic field creates the so-called Equatorial Ionization Anomalies. Ground based measurements from GNSS receivers are used to link the Total Electron Content (TEC), its spatial gradients and the phenomenon known as scintillation that can lead to a GNSS signal degradation or even to a GNSS signal ‘loss of lock’. A new algorithm to highlight the features characterizing the TEC distribution is developed in the framework of this thesis and the results obtained are validated and used to improve the performance of a GNSS positioning technique (long baseline RTK). In addition, the correlation between scintillation and dynamics of the ionospheric irregularities is investigated. By means of a software, here implemented, the velocity of the ionospheric irregularities is evaluated using high sampling rate GNSS measurements. The results highlight the parallel behaviour of the amplitude scintillation index (S4) occurrence and the zonal velocity of the ionospheric irregularities at least during severe scintillations conditions (post-sunset hours). This suggests that scintillations are driven by TEC gradients as well as by the dynamics of the ionospheric plasma. Finally, given the importance of such studies for technological applications (e.g. GNSS high-precision applications), a validation of the NeQuick model (i.e. the model used in the new GALILEO satellites for TEC modelling) is performed. The NeQuick performance dramatically improves when data from HF radar sounding (ionograms) are ingested. A custom designed algorithm, based on the image recognition technique, is developed to properly select the ingested data, leading to further improvement of the NeQuick performance.

Retrieval of tropospheric water vapour by using spectra of a 22 GHz radiometer

In this paper, we present an approach to retrieve tropospheric water vapour profiles from pressure broadened emission spectra at 22 GHz, measured by a ground based microwave radiometer installed in the south of Bern at 905 m. Classical microwave instruments concentrating on the troposphere observe several channels in the center and the wings of the water vapour line (20–30 Ghz), whereas our retrieval approach uses spectra with a bandwidth of 1 GHz and a high resolution around the center of the 22 GHz water vapour line. The retrieval is sensitive up to 7 km with a vertical resolution of 3–5 km. Comparisons with profiles from operational balloon soundings, performed at Payerne, 40 km away from the radiometer location, showed a good agreement up to 7 km with a correlation of above 0.8. The retrievals shows a wet bias of 10–20% compared to the sounding.

Gasoline emissions dominate over diesel in formation of secondary organic aerosol mass

Although laboratory experiments have shown that organic compounds in both gasoline fuel and diesel engine exhaust can form secondary organic aerosol (SOA), the fractional contribution from gasoline and diesel exhaust emissions to ambient SOA in urban environments is poorly known. Here we use airborne and ground-based measurements of organic aerosol (OA) in the Los Angeles (LA) Basin, California made during May and June 2010 to assess the amount of SOA formed from diesel emissions. Diesel emissions in the LA Basin vary between weekdays and weekends, with 54% lower diesel emissions on weekends. Despite this difference in source contributions, in air masses with similar degrees of photochemical processing, formation of OA is the same on weekends and weekdays, within the measurement uncertainties. This result indicates that the contribution from diesel emissions to SOA formation is zero within our uncertainties. Therefore, substantial reductions of SOA mass on local to global scales will be achieved by reducing gasoline vehicle emissions.

CH4, CO, and H2O spectroscopy for the Sentinel-5 Precursor mission: an assessment with the Total Carbon Column Observing Network measurements

The TROPOspheric Monitoring Instrument (TROPOMI) will be part of ESA's Sentinel-5 Precursor (S5P) satellite platform scheduled for launch in 2015. TROPOMI will monitor methane and carbon monoxide concentrations in the Earth's atmosphere by measuring spectra of back-scattered sunlight in the short-wave infrared (SWIR). S5P will be the first satellite mission to rely uniquely on the spectral window at 4190–4340 cm−1 (2.3 μm) to retrieve CH4 and CO. In this study, we investigated if the absorption features of the three relevant molecules CH4, CO, and H2O are adequately known. To this end, we retrieved total columns of CH4, CO, and H2O from absorption spectra measured by two ground-based Fourier transform spectrometers that are part of the Total Carbon Column Observing Network (TCCON). The retrieval results from the 4190–4340 cm−1 range at the TROPOMI resolution (0.45 cm−1) were then compared to the CH4 results obtained from the 6000 cm−1 region, and the CO results obtained from the 4190–4340 cm−1 region at the higher TCCON resolution (0.02 cm−1). For TROPOMI-like settings, we were able to reproduce the CH4 columns to an accuracy of 0.3% apart from a constant bias of 1%. The CO retrieval accuracy was, through interference, systematically influenced by the shortcomings of the CH4 and H2O spectroscopy. In contrast to CH4, the CO column error also varied significantly with atmospheric H2O content. Unaddressed, this would introduce seasonal and latitudinal biases to the CO columns retrieved from TROPOMI measurements. We therefore recommend further effort from the spectroscopic community to be directed at the H2O and CH4 spectroscopy in the 4190–4340 cm−1 region.

Improved water vapour spectroscopy in the 4174-4300 cm-1 region and its impact on SCIAMACHY HDO/H2O measurements

The relative abundance of the heavy water isotopologue HDO provides a deeper insight into the atmospheric hydrological cycle. The SCanning Imaging Absorption spectroMeter for Atmospheric CartograpHY (SCIAMACHY) allows for global retrievals of the ratio HDO/H2O in the 2.3 micron wavelength range. However, the spectroscopy of water lines in this region remains a large source of uncertainty for these retrievals. We therefore evaluate and improve the water spectroscopy in the range 4174–4300 cm−1 and test if this reduces systematic uncertainties in the SCIAMACHY retrievals of HDO/H2O. We use a laboratory spectrum of water vapour to fit line intensity, air broadening and wavelength shift parameters. The improved spectroscopy is tested on a series of ground-based high resolution FTS spectra as well as on SCIAMACHY retrievals of H2O and the ratio HDO/H2O. We find that the improved spectroscopy leads to lower residuals in the FTS spectra compared to HITRAN 2008 and Jenouvrier et al. (2007) spectroscopy, and the retrievals become more robust against changes in the retrieval window. For both the FTS and SCIAMACHY measurements, the retrieved total H2O columns decrease by 2–4% and we find a negative shift of the HDO/H2O ratio, which for SCIAMACHY is partly compensated by changes in the retrieval setup and calibration software. The updated SCIAMACHY HDO/H2O product shows somewhat steeper latitudinal and temporal gradients and a steeper Rayleigh distillation curve, strengthening previous conclusions that current isotope-enabled general circulation models underestimate the variability in the near-surface HDO/H2O ratio.

Assessing visual disturbance conditions on the Custer National Forest

Assessment of soil disturbance on the Custer National Forest was conducted during two summers to determine if the U.S. Forest Service Forest Soil Disturbance Monitoring Protocol (FSDMP) was able to distinguish post-harvest soil conditions in a chronological sequence of sites harvested using different ground-based logging systems. Results from the first year of sampling suggested that the FSDMP point sampling method may not be sensitive enough to measure post-harvest disturbance in stands with low levels of disturbance. Therefore, a revised random transect method was used during the second sampling season to determine the actual extent of soil disturbance in these cutting units. Using combined data collected from both summers I detected statistically significant differences (p < 0.05) in fine fraction bulk density measurements between FSDMP disturbance classes across all sites. Disturbance class 3 (most severe) had the highest reported bulk density, which suggest that the FSDMP visual class estimates are defined adequately allowing for correlations to be made between visual disturbance and actual soil physical characteristics. Forest site productivity can be defined by its ability to retain carbon and convert it to above- and belowground biomass. However, forest management activities that alter basic site characteristics have the potential to alter productivity. Soil compaction is one critical management impact that is important to understand; compaction has been shown to impede the root growth potential of plants, reduce water infiltration rates increasing erosion potential, and alter plant available water and nutrients, depending on soil texture. A new method to assess ground cover, erosion, and other soil disturbances was recently published by the U.S. Forest Service, as the Forest Soil Disturbance Protocol (FSDMP). The FSDMP allows soil scientists to visually assign a disturbance class estimate (0 – none, 1, 2, 3 – severe) from field measures of consistently defined soil disturbance indicators (erosion, fire, rutting, compaction, and platy/massive/puddled structure) in small circular (15 cm) plots to compare soil quality properties pre- and post- harvest condition. Using this protocol we were able to determine that ground-based timber harvesting activities occurring on the Custer National Forest are not reaching the 15% maximum threshold for detrimental soil disturbance outlined by the Region 1 Soil Quality Standards.

SO2 emissions at Volcan de Colima, 2003-2007

Volcán de Colima has been continuously erupting since the onset of dome growth in 1998. This period of unrest has had 4 prominent periods; 1998-1999, 2003, 2004-2005, and the current dome growth that began in February of 2007. Each of these episodes was marked by lava extrusion forming a dome and lava flows, followed by explosions that destroyed the dome. The Correlation Spectrometer (COSPEC) was used to determine SO2 emission rates on 164 days from May 2003 to February 2007, using both stationary ground based scans and some flight traverses. Scans were separated into the categories of explosive degassing and passive, or background degassing. These scans show variation in the SO2 flow rate from below detection limit (~3 t/d depending on environmental conditions) during background, passive emissions to a peak of 2949 t/d (34 kilograms/second) during an explosion on 9 October, 2004. Both passive and explosive degassing increased when there was lava extrusion in 2004 and with the increased explosive activity in 2005. These two different processes of degassing wax with each other when activity increases and wane together as well, indicating a parallel cyclicity in the volcanic eruption and degassing rates, where the conduit partially seals (pressurizes) between explosions. Colima’s gas and eruptive behavior is compared to similar systems such as Santiaguito and Soufrière Hills, Montserrat. About 2/3 of Colima’s SO2 degassing, amounting to 1.3 x 105 tonnes in 3.74 yrs has come in short lived small (VEI=0-1) vertical explosions that occurred at the rate of 100-3000explosions/ month, and the remaining third has occured in continuous passive degassing. Colima emits sulfur at a rate equivalent to about 0.04 to 0.08 wt % S, similar to other andesitic convergent plate boundary volcanoes. There has been an explosive destruction of the dome in every cycle for that past 5 years, and it is assumed that the current dome which began growth in February, 2007 (just at the end of this study) will be destroyed. Higher emission rates seen in the quiescence of 2006 may have eased the pressure at the time, resulting in the slow effusion of the current dome and lack of explosivity.

Direction finding in the presence of a more realistic environment model

Direction-of-arrival (DOA) estimation is susceptible to errors introduced by the presence of real-ground and resonant size scatterers in the vicinity of the antenna array. To compensate for these errors pre-calibration and auto-calibration techniques are presented. The effects of real-ground constituent parameters on the mutual coupling (MC) of wire type antenna arrays for DOA estimation are investigated. This is accomplished by pre-calibration of the antenna array over the real-ground using the finite element method (FEM). The mutual impedance matrix is pre-estimated and used to remove the perturbations in the received terminal voltage. The unperturbed terminal voltage is incorporated in MUSIC algorithm to estimate DOAs. First, MC of quarter wave monopole antenna arrays is investigated. This work augments an existing MC compensation technique for ground-based antennas and proposes reduction in MC for antennas over finite ground as compared to the perfect ground. A factor of 4 decrease in both the real and imaginary parts of the MC is observed when considering a poor ground versus a perfectly conducting one for quarter wave monopoles in the receiving mode. A simulated result to show the compensation of errors direction of arrival (DOA) estimation with actual realization of the environment is also presented. Secondly, investigations for the effects on received MC of λ/2 dipole arrays placed near real-earth are carried out. As a rule of thumb, estimation of mutual coupling can be divided in two regions of antenna height that is very near ground 0