Investigating Peatland Stratigraphy and Hydrogeology Using Integrated Electrical Geophysics

Hydrology has been suggested as the mechanism controlling vegetation and related surficial pore-water chemistry in large peatlands. Peatland hydrology influences the carbon dynamics within these large carbon reservoirs and will influence their response to global warming. A geophysical survey was completed in Caribou Bog, a large peatland in Maine, to evaluate peatland stratigraphy and hydrology. Geophysical measurements were integrated with direct measurements of peat stratigraphy from probing, fluid chemistry, and vegetation patterns in the peatland. Consistent with previous field studies, ground-penetrating radar (GPR) was an excellent method for delineating peatland stratigraphy. Prominent reflectors from the peat-lake sediment and lake sediment-mineral soil contacts were precisely recorded up to 8 m deep. Two-dimensional resistivity and induced polarization imaging were used to investigate stratigraphy beneath the mineral soil, beyond the range of GPR. We observe that the peat is chargeable, and that IP imaging is an alternative method for defining peat thickness. The chargeability of peat is attributed to the high surface-charge density on partially decomposed organic matter. The electrical conductivity imaging resolved glaciomarine sediment thickness (a confining layer) and its variability across the basin. Comparison of the bulk conductivity images with peatland vegetation revealed a correlation between confining layer thickness and dominant vegetation type, suggesting that stratigraphy exerts a control on hydrogeology and vegetation distribution within this peatland. Terrain conductivity measured with a Geonics EM31 meter correlated with confining glaciomarine sediment thickness and was an effective method for estimating variability in glaciomarine sediment thickness over approximately 18 km(2). Our understanding of the hydrogeology, stratigraphy, and controls on vegetation growth in this peatland was much enhanced from the geophysical study.

Spatial Variability in Biogenic Gas Accumulations in Peat Soils Is Revealed By Ground Penetrating Radar (GPR)

We performed surface and borehole ground penetrating radar (GPR) tests, together with moisture probe measurements and direct gas sampling to detect areas of biogenic gas accumulation in a northern peatland. The main findings are: (1) shadow zones (signal scattering) observed in surface GPR correlate with areas of elevated CH4 and CO2 concentration; (2) high velocities in zero offset profiles and lower water content inferred from moisture probes correlate with surface GPR shadow zones; (3) zero offset profiles depict depth variable gas accumulation from 0-10% by volume; (4) strong reflectors may represent confining layers restricting upward gas migration. Our results have implications for defining the spatial distribution, volume and movement of biogenic gas in peatlands at multiple scales.

Radiometric Method for Emissivity Retrieval in High Reflective Materials

High reflective materials in the microwave region play a very important role in the realization of antenna reflectors for a broad range of applications, including radiometry. These reflectors have a characteristic emissivity which needs to be characterized accurately in order to perform a correct radiometric calibration of the instrument. Such a characterization can be performed by using open resonators, waveguide cavities or by radiometric measurements. The latter consists of comparative radiometric observations of absorbers, reference mirrors and the sample under test, or using the cold sky radiation as a direct reference source. While the first two mentioned techniques are suitable for the characterization of metal plates and mirrors, the latter has the advantages to be also applicable to soft materials. This paper describes how, through this radiometric techniques, it is possible to characterize the emissivity of the sample relative to a reference mirror and how to characterize the absolute emissivity of the latter by performing measurements at different incident angles. The results presented in this paper are based on our investigations on emissivity of a multilayer insulation material (MLI) for space mission, at the frequencies of 22 and 90 GHz.

A Compensating Anastigmatic Submillimeter Array Imaging System for STEAMR

In this paper, we present a novel technique for the removal of astigmatism in submillimeter-wave optical systems through employment of a specific combination of so-called astigmatic off-axis reflectors. This technique treats an orthogonally astigmatic beam using skew Gaussian beam analysis, from which an anastigmatic imaging network is derived. The resultant beam is considered truly stigmatic, with all Gaussian beam parameters in the orthogonal directions being matched. This is thus considered an improvement over previous techniques wherein a beam corrected for astigmatism has only the orthogonal beam amplitude radii matched, with phase shift and phase radius of curvature not considered. This technique is computationally efficient, negating the requirement for computationally intensive numerical analysis of shaped reflector surfaces. The required optical surfaces are also relatively simple to implement compared to such numerically optimized shaped surfaces. This technique is implemented in this work as part of the complete optics train for the STEAMR antenna. The STEAMR instrument is envisaged as a mutli-beam limb sounding instrument operating at submillimeter wavelengths. The antenna optics arrangement for this instrument uses multiple off-axis reflectors to control the incident radiation and couple them to their corresponding receiver feeds. An anastigmatic imaging network is successfully implemented into an optical model of this antenna, and the resultant design ensures optimal imaging of the beams to the corresponding feed horns. This example also addresses the challenges of imaging in multi-beam antenna systems.

Homogeneous reprocessing of GPS, GLONASS and SLR observations

The International GNSS Service (IGS) provides operational products for the GPS and GLONASS constellation. Homogeneously processed time series of parameters from the IGS are only available for GPS. Reprocessed GLONASS series are provided only by individual Analysis Centers (i. e. CODE and ESA), making it difficult to fully include the GLONASS system into a rigorous GNSS analysis. In view of the increasing number of active GLONASS satellites and a steadily growing number of GPS+GLONASS-tracking stations available over the past few years, Technische Universität Dresden, Technische Universität München, Universität Bern and Eidgenössische Technische Hochschule Zürich performed a combined reprocessing of GPS and GLONASS observations. Also, SLR observations to GPS and GLONASS are included in this reprocessing effort. Here, we show only SLR results from a GNSS orbit validation. In total, 18 years of data (1994–2011) have been processed from altogether 340 GNSS and 70 SLR stations. The use of GLONASS observations in addition to GPS has no impact on the estimated linear terrestrial reference frame parameters. However, daily station positions show an RMS reduction of 0.3 mm on average for the height component when additional GLONASS observations can be used for the time series determination. Analyzing satellite orbit overlaps, the rigorous combination of GPS and GLONASS neither improves nor degrades the GPS orbit precision. For GLONASS, however, the quality of the microwave-derived GLONASS orbits improves due to the combination. These findings are confirmed using independent SLR observations for a GNSS orbit validation. In comparison to previous studies, mean SLR biases for satellites GPS-35 and GPS-36 could be reduced in magnitude from −35 and −38 mm to −12 and −13 mm, respectively. Our results show that remaining SLR biases depend on the satellite type and the use of coated or uncoated retro-reflectors. For Earth rotation parameters, the increasing number of GLONASS satellites and tracking stations over the past few years leads to differences between GPS-only and GPS+GLONASS combined solutions which are most pronounced in the pole rate estimates with maximum 0.2 mas/day in magnitude. At the same time, the difference between GLONASS-only and combined solutions decreases. Derived GNSS orbits are used to estimate combined GPS+GLONASS satellite clocks, with first results presented in this paper. Phase observation residuals from a precise point positioning are at the level of 2 mm and particularly reveal poorly modeled yaw maneuver periods.

Optical Tolerance Analysis of the Multi-Beam Limb Viewing Instrument STEAMR

In this paper, we report on an optical tolerance analysis of the submillimeter atmospheric multi-beam limb sounder, STEAMR. Physical optics and ray-tracing methods were used to quantify and separate errors in beam pointing and distortion due to reflector misalignment and primary reflector surface deformations. Simulations were performed concurrently with the manufacturing of a multi-beam demonstrator of the relay optical system which shapes and images the beams to their corresponding receiver feed horns. Results from Monte Carlo simulations show that the inserts used for reflector mounting should be positioned with an overall accuracy better than 100 μm (~ 1/10 wavelength). Analyses of primary reflector surface deformations show that a deviation of magnitude 100 μm can be tolerable before deployment, whereas the corresponding variations should be less than 30 μm during operation. The most sensitive optical elements in terms of misalignments are found near the focal plane. This localized sensitivity is attributed to the off-axis nature of the beams at this location. Post-assembly mechanical measurements of the reflectors in the demonstrator show that alignment better than 50 μm could be obtained.

Satellite laser ranging to GPS and GLONASS

Satellite laser ranging (SLR) to the satellites of the global navigation satellite systems (GNSS) provides substantial and valuable information about the accuracy and quality of GNSS orbits and allows for the SLR-GNSS co-location in space. In the framework of the NAVSTAR-SLR experiment two GPS satellites of Block-IIA were equipped with laser retroreflector arrays (LRAs), whereas all satellites of the GLONASS system are equipped with LRAs in an operational mode. We summarize the outcome of the NAVSTAR-SLR experiment by processing 20 years of SLR observations to GPS and 12 years of SLR observations to GLONASS satellites using the reprocessed microwave orbits provided by the center for orbit determination in Europe (CODE). The dependency of the SLR residuals on the size, shape, and number of corner cubes in LRAs is studied. We show that the mean SLR residuals and the RMS of residuals depend on the coating of the LRAs and the block or type of GNSS satellites. The SLR mean residuals are also a function of the equipment used at SLR stations including the single-photon and multi-photon detection modes. We also show that the SLR observations to GNSS satellites are important to validate GNSS orbits and to assess deficiencies in the solar radiation pressure models. We found that the satellite signature effect, which is defined as a spread of optical pulse signals due to reflection from multiple reflectors, causes the variations of mean SLR residuals of up to 15 mm between the observations at nadir angles of 0∘ and 14∘. in case of multi-photon SLR stations. For single-photon SLR stations this effect does not exceed 1 mm. When using the new empirical CODE orbit model (ECOM), the SLR mean residual falls into the range 0.1–1.8 mm for high-performing single-photon SLR stations observing GLONASS-M satellites with uncoated corner cubes. For best-performing multi-photon stations the mean SLR residuals are between −12.2 and −25.6 mm due to the satellite signature effect.

Facies investigations on sediment core CRP-3 from the Ross Sea, Antarctica

The Cenozoic Victoria Land Basin (VLB) stratigraphic section penetrated by CRP-3 is mostly of Early Oligocene age. It contains an array of lithofacies comprising fine-grained mudrocks, interlaminated and interbedded mudrocks/sandstones, mud-rich and mud-poor sandstones, conglomerates and diamctites that are together interpreted as the products of shallow marine to possibly non-marine environments of deposition, affected by the periodic advance and retreat of tidewater glaciers. This lithofacies assemblage can be readily rationalised using the facies scheme designed originally for CRP-2/2A, and published previously. The uppermost 330 metres below sea floor (mbsf) shows a cyclical arrangement of lithofacies also similar to that recognised throughout CRP-2/2A, and interpreted to reflect cyclical variations in relative sea-level driven by ice volume fluctuations ('Motif A'). Between 330 and 480 mbsf, a series of less clearly cyclical units, generally fining-upward but nonetheless incorporating a significant subset of the facies assemblage, has been identified and noted in the Initial Report as 'Motif B' Below 480 mbsf, the section is arranged into a repetitive succession of fining-upward units, each of which comprises dolerite clast conglomerate at the base passing upward into relatively thick intervals of sandstones. The cycles present down 480 mbsf are defined as sequences, each interpreted to record cyclical variation of relative sea-level. The thickness distribution of sequences in CRP-3 provides some insights into the geological variables controlling sediment accumulation in the Early Oligocene section. The uppermost part of the section in CRP-3 comprises two or three thick, complete sequences that show a broadly symmetrical arrangement of lithofacies (similar to Sequences 9-11 in CRP-2/2A). This suggests a period of relatively rapid tectonic subsidence, which allowed preservation of the complete facies cycle. Below Sequence 3, however, is a considerable interval of thin, incomplete and erosionally truncated sequences (4-23), which incorporates both the remainder of Motif A sequences and all Motif B sequences recognised. The thinner and more truncated sequences suggest sediment accumulation under conditions of reduced accommodation, and given the lack of evidence for glacial conditions (see Powell et al., this volume) tends to argue for a period of reduced tectonic subsidence. The section below 480 mbsf consists of a series of fining-upward, conglomerate to sandstone intervals which cannot be readily interpreted in terms of relative sea-level change. A relatively mudrock-rich interval above the basal conglomerate/breccia (782-762 mbsf) may record initial flooding of the basin during early rift subsidence. The lithostratigraphy summarised above has been linked to seismic reflection data using depth conversion techniques (Henrys et al., this volume). The three uppermost reflectors ('o', 'p' and 'q') correlate to the package of thick sequences 1-3, and several deeper reflectors can also be correlated to sequence boundaries. The package of thick Sequences 1-3 shows a sheet-like cross-sectional geometry on seismic reflection lines, unlike the similar package recognised in CRP-2/2A.

Age models and stable isotope ratios from the Mid-Miocene

Mid-Miocene pelagic sedimentary sections can be correlated using intermediate and high resolution oxygen and carbon isotopic records of benthic foraminifera. Precision of a few tens of thousands of years is readily achievable at sites with high sedimentation rates, for example, Deep Sea Drilling Project sites 289 and 574. The mid-Miocene carbon isotope records are characterized by an interval of high d13C values between 17 and 13.5 Ma (the Monterey Excursion of Vincent and Berger 1985) upon which are superimposed a series of periodic or quasi-periodic fluctuations in d13C values. These fluctuations have a period of approximately 440 kyr, suggestive of the 413 kyr cycle predicted by Milankovitch theory. Vincent and Berger proposed that the Monterey Excursion was the result of increased organic carbon burial in continental margins sediments. The increased d13C values (called 13C maxima) superimposed on the generally high mid-Miocene signal coincide with increases in d18O values suggesting that periods of cooling and/or ice buildup were associated with exceptionally rapid burial of organic carbon and lowered atmospheric CO2 levels. It is likely that during the Monterey Excursion the ocean/atmosphere system became progressively more sensitive to small changes in insolation, ultimately leading to major cooling of deep water and expansion of continental ice. We have assigned an absolute chronology, based on biostratigraphic and magneto-biostratigraphic datum levels, to the isotope stratigraphy and have used that chronology to correlate unconformities, seismic reflectors, carbonate minima, and dissolution intervals. Intervals of sediment containing 13C maxima are usually better preserved than the overlying and underlying sediments, indicating that the d13C values of TCO2 in deep water and the corrosiveness of seawater are inversely correlated. This again suggests that the 13C maxima were associated with rapid burial of organic carbon and reduced levels of atmospheric CO2. The absolute chronology we have assigned to the isotopic record indicates that the major mid-Miocene deepwater cooling/ice volume expansion took 2 m.y. and was not abrupt as had been reported previously. The cooling appears abrupt at many sites because the interval is characterized by a number of dissolution intervals. The cooling was not monotonic, and the 2 m.y. interval included an episode of especially rapid cooling as well as a brief return to warmer conditions before the final phase of the cooling period. The increase in d18O values of benthic foraminifera between 14.9 and 12.9 Ma was greatest at deeper water sites and at sites closest to Antarctica. The data suggest that the d18O value of seawater increased by no more than about 1.1 per mil during this interval and that the remainder of the change in benthic d18O values resulted from cooling in Antarctic regions of deepwater formation. Equatorial planktonic foraminifera from sites 237 and 289 exhibit a series of 0.4 per mil steplike increases in d13C values. Only one of these increases in planktonic d13C is correlated with any of the features in the mid-Miocene benthic carbon isotope record.

Turbidite sedimentology and history at DSDP Site 89-585

At Site 585 in the East Mariana Basin, a 900-m section of Aptian-Albian to Recent sediments was recovered. The upper 590 m are pelagic components (carbonate, siliceous, and clay); small-scale graded sequences and laminations are common. The underlying sediments are volcaniclastic sandstones with a large proportion of shallow-water carbonate debris; sedimentary structures including complete Bouma sequences, cross-laminae, and scouring are common. These structures indicate that the entire section was deposited by turbidity currents. The change in lithology upward in the section reflects the evolution of the surrounding seamounts, from their growth stages during the middle of the Cretaceous to the later subsidence phases. Several black layers containing pyritized organic debris and associated turbidite structures were cored near the Cenomanian/Turonian boundary; this material has been transported from the flanks of the seamounts where it was deposited within a shallow anoxic zone. Seismic data extends the stratigraphy across the entire Basin, showing the reflectors onlapping the seamounts, and indicating at least 1200 m of sediment at Site 585. The crust is placed at 6900 m after correcting for sediment loading, and the subsidence curve indicates that the Basin has been deeper than 5500 m since before the Aptian.

(Table 3) Sedimentation rates of stage 1 and last glacial maximum from different Holes during cruises of RRS Discovery in 1984-1985 and 1987-1988, and RRS James Clark Ross in 1992-1993

Within the Scotia Sea, the axis of the Antarctic Circumpolar Current (ACC) is geographically confined, and sediments therefore contain a record of palaeo-flow speed uncomplicated by ACC axis migration. We outline Holocene and Last Glacial Maximum (LGM) current-controlled sedimentation using data from 3.5-kHz profiles, cores and current meter moorings. Geophysical surveys show areas of erosion and deposition controlled by Neogene basement topography. Deposition occurs in mounded sediment drifts or flatter areas, where 500-1000 m of sediment overlies acoustic basement. 3.5-kHz profiles show parallel, continuous sub-bottom reflectors with highest sedimentation rates in the centre of the drifts, and reflectors converging towards marginal zones of non-deposition. Locally, on the flanks of continental blocks (e.g. South Georgia), downslope processes are dominant. The absence of mudwaves on the sediment drifts may result from the unsteadiness of ACC flow. A core transect from the ACC axis south to the boundary with the Weddell Gyre shows a southward decrease in biogenic content, controlled by the Polar Front and the spring sea-ice edge. Both these features lay farther north at LGM. The cores have been dated by relative abundance of the radiolarian Cycladophora davisiana, and by changes in the biogenic Ba content, a palaeoproductivity indicator. Sedimentation rates range from 3 to 17 cm/ka. The grain size of Holocene sediments shows a coarsening trend from south to north, consistent with strongest bottom-current flow near the ACC axis, though interpretation is complicated by the presence of biogenic grains. Year-long current meter records indicate mean speeds from 7 cm/s in the south to 12 cm/s in the north, with benthic storm frequency increasing northwards. LGM sediments are predominantly terrigenous and show a clearer northward-coarsening trend, with well-sorted silts in the northern Scotia Sea. Assuming a constant terrigenous source, this implies stronger ACC flow at the LGM, contrasting with weaker Weddell Gyre flow deduced from earlier work.

Chemical and Sr isotopic compositions of 3 step leachates and residues of bulk samples, zeolites, and smectites from Lower Miocene sediments at ODP Site 166-1007

Massive clinoptilolite authigenesis was observed at about 1105 meters below sea floor (mbsf) in lower Miocene wellcompacted carbonate periplatform sediments from the Great Bahama Bank [Ocean Drilling Program, ODP Leg 166, Site 1007]. The diagenetic assemblage comprises abundant zeolite crystallized within foraminifer tests and sedimentary matrix, as well as Mg smectites. In carbonate-rich deposits, the formation of the zeolite requires a supply of silica. Thus, the objective of the study is to determine the origin of the silica supply, its diagenetic evolution, and consequently the related implications on interpretation of the sedimentary record, in terms of local or global paleoceanographic change. For lack of evidence for any volcaniclastic input or traces of Si-enriched deep fluids circulation, an in situ biogenic source of silica is validated by isotopic data and chemical modeling for the formation of such secondary minerals in shallow-water carbonate sequences. Geochemical and strontium isotopic data clearly establish the marine signature of the diagenetic zeolite, as well as its contemporaneous formation with the carbonate deposition (Sr model ages of 19.6-23.2 Ma). The test of saturation for the pore fluids specifies the equilibrium state of the present mineralogical assemblage. Seawater-rock modeling specifies that clinoptilolite precipitates from the dissolution of biogenic silica, which reacts with clay minerals. The amount of silica (opal-A) involved in the reaction has to be significant enough, at least 10 wt.%, to account for the observed content of clinoptilolite occurring at the most zeolite-rich level. Modeling also shows that the observed amount of clinoptilolite (~19%) reflects an in situ and short-term reaction due to the high reactivity of primary biogenic silica (opal-A) until its complete depletion. The episodic occurrence of these well-lithified zeolite-rich levels is consistent with the occurrence of seismic reflectors, particularly the P2 seismic sequence boundary located at 1115 mbsf depth and dated as 23.2 Ma. The age range of most zeolitic sedimentary levels (biostratigraphic ages of 21.5-22 Ma) correlates well with that of the early Miocene glaciation Mi-1 and Mi-1a global events. Thus, the clinoptilolite occurrence in the shallow carbonate platform environment far from volcanogenic supply, or in other sensitive marine areas, is potentially a significant new proxy for paleoproductivity and oceanic global events, such as the Miocene events, which are usually recognized in deep-sea pelagic sediments and high latitude deposits.

Physical properties and elastic constants of ODP Holes 137-504B and 140-504B upper crustal rocks

Seismic velocities have been measured at confining pressures of 100 MPa and 600 MPa for sheeted dike samples recovered during Ocean Drilling Program Legs 137 and 140. The compressional- and shear-wave velocities show an increase with depth at Hole 504B, which is in sharp contrast to the atmospheric pressure velocity measurements performed as part of the shipboard analyses. Rocks exposed to different types of alteration and fracture patterns show distinct changes in their physical properties. The seismic reflectors observed on the vertical seismic profile (VSP) experiment performed during Leg 111 may have been caused by low velocity zones resulting from alteration. The amount of fracturing and hydrothermal alteration in several zones also may have contributed to the acoustic impedance contrast necessary to produce the E5 reflector. Poisson's ratios calculated from laboratory velocity measurements show several low values at depths ranging from 1600 mbsf to 2000 mbsf, which tends to follow similar trends obtained from previous oceanic refraction experiments. A comparison of physical properties between samples recovered from Hole 504B and ophiolite studies in the Bay of Islands and Oman shows a good correlation with the Bay of Islands but significant differences from the measurements performed in the Oman complex.