Techniques for Lagrangian modelling of dispersion in geophysical flows

Basic concepts and definitions relative to Lagrangian Particle Dispersion Models (LPDMs)for the description of turbulent dispersion are introduced. The study focusses on LPDMs that use as input, for the large scale motion, fields produced by Eulerian models, with the small scale motions described by Lagrangian Stochastic Models (LSMs). The data of two different dynamical model have been used: a Large Eddy Simulation (LES) and a General Circulation Model (GCM). After reviewing the small scale closure adopted by the Eulerian model, the development and implementation of appropriate LSMs is outlined. The basic requirement of every LPDM used in this work is its fullfillment of the Well Mixed Condition (WMC). For the dispersion description in the GCM domain, a stochastic model of Markov order 0, consistent with the eddy-viscosity closure of the dynamical model, is implemented. A LSM of Markov order 1, more suitable for shorter timescales, has been implemented for the description of the unresolved motion of the LES fields. Different assumptions on the small scale correlation time are made. Tests of the LSM on GCM fields suggest that the use of an interpolation algorithm able to maintain an analytical consistency between the diffusion coefficient and its derivative is mandatory if the model has to satisfy the WMC. Also a dynamical time step selection scheme based on the diffusion coefficient shape is introduced, and the criteria for the integration step selection are discussed. Absolute and relative dispersion experiments are made with various unresolved motion settings for the LSM on LES data, and the results are compared with laboratory data. The study shows that the unresolved turbulence parameterization has a negligible influence on the absolute dispersion, while it affects the contribution of the relative dispersion and meandering to absolute dispersion, as well as the Lagrangian correlation.

Integrated Geophysical Investigation of the St. James Fault Complex: A Case Study

We noninvasively detected the characteristics and location of a regional fault in an area of poor bedrock exposure complicated by karst weathering features in the subsurface. Because this regional fault is associated with sinkhole formation, its location is important for hazard avoidance. The bedrock lithologies on either side of the fault trace are similar; hence, we chose an approach that capitalized on the complementary strengths of very low frequency (VLF) electromagnetic, resistivity, and gravity methods. VLF proved most useful as a first-order reconnaissance tool, allowing us to define a narrow target area for further geophysical exploration. Fault-related epikarst was delineated using resistivity. Ultimately, a high-resolution gravity survey and subsequent inverse modeling using the results of the resistivity survey helped to further constrain the location and approximate orientation of the fault. The combined results indicated that the location of the fault trace needed to be adjusted 53 m south of the current published location and was consistent with a north-dipping thrust fault. Additionally, a gravity low south of the fault trace agreed with the location of conductive material from the resistivity and VLF surveys. We interpreted these anomalies to represent enhanced epikarst in the fault footwall. We clearly found that a staged approach involving a progression of methods beginning with a reconnaissance VLF survey, followed by high-resolution gravity and electrical resistivity surveys, can be used to characterize a fault and fault-related karst in an area of poor bedrock surface exposure.

Geophysical Analysis of Seasonal Montandon Gravel Ridge Water Table Fluctuation and Moisture Gradient Variation due to Storm Events

The purpose of this research project is to continue exploring the Montandon Long-Term Hydrologic Research Site(LTHR) by using multiple geophysical methods to obtain more accurate and precise information regarding subsurface hydrologic properties of a local gravel ridge,which are important to both the health of surrounding ecosystems and local agriculture. Through using non-invasive geophysical methods such as seismic refraction, Direct Current resistivity and ground penetrating radar (GPR) instead of invasive methods such as boreholedrilling which displace sediment and may alter water flow, data collection is less likely to bias the data itself. In addition to imaging the gravel ridge subsurface, another important researchpurpose is to observe how both water table elevation and the moisture gradient (moisture content of the unsaturated zone) change over a seasonal time period and directly after storm events. The combination of three types of data collection allows the strengths of each method combine together and provide a relatively strongly supported conclusions compared to previous research. Precipitation and geophysical data suggest that an overall increase in precipitation during the summer months causes a sharp decrease in subsurface resistivity within the unsaturated zone. GPR velocity data indicate significant immediate increase in moisture content within the shallow vadose zone (< 1m), suggesting that rain water was infiltrating into the shallow subsurface. Furthermore, the combination of resistivity and GPR results suggest that the decreased resistivity within the shallow layers is due to increased ion content within groundwater. This is unexpected as rainwater is assumed to have a DC resistivity value of 3.33*105 ohm-m. These results may suggest that ions within the sediment must beincorporated into the infiltrating water.

Monitoring Stormwater Redistribution Into Low Resistivity Soils Using Noninvasive Geophysical Techniques

Water held in the unsaturated zone is important for agriculture and construction and is replenished by infiltrating rainwater. Monitoring the soil water content of clay soils using ground-penetrating radar (GPR) has not been researched, as clay soils cause attenuation of GPR signal. In this study, GPR common-midpoint soundings (CMPs) are used in the clayey soils of the Miller Run floodplain to monitor changes in the soil water content (SWC) before and after rainfall events. GPR accomplishes this task because increases in water content will increase the dielectric constant of the subsurface material, and decrease the velocity of the GPR wave. Using an empirical relationship between dielectric constant and SWC, the Topp relation, we are able to calculate a SWC from these velocity measurements. Non-invasive electromagnetics, resistivity, and seismic were performed, and from these surveys, the layering at the field site was delineated. EM characterized the horizontal variation of the soil, allowing us to target the most clay rich area. At the CMP location, resistivity indicates the vertical structure of the subsurface consists of a 40 cm thick layer with a resistivity of 100 ohm*m. Between 40 cm and 1.5 m is a layer with a resistivity of 40 ohm*m. The thickness estimates were confirmed with invasive auger and trenching methods away from the CMP location. GPR CMPs were collected relative to a July 2013 and September 2013 storm. The velocity observations from the CMPs had a precision of +/- 0.001 m/ns as assessed by repeat analysis. In the case of both storms, the GPR data showed the expected relationship between the rainstorms and calculated SWC, with the SWC increasing sharply after the rainstorm and decreasing as time passed. We compared these data to auger core samples collected at the same time as the CMPs were taken, and the volumetric analysis of the cores confirmed the trend seen in the GPR, with SWC values between 3 and 5 percent lower than the GPR estimates. Our data shows that we can, with good precision, monitor changes in the SWC of conductive soils in response to rainfall events, despite the attenuation induced by the clay.

Using Micro-Gravity Techniques to Map Alluvium Thickness and Pleistocene Location of the West Branch of the Susquehanna River Near Muncy, Pennsylvania

Laurentide glaciation during the early Pleistocene (~970 ka) dammed the southeast-flowing West Branch of the Susquehanna River (WBSR), scouring bedrock and creating 100-km-long glacial Lake Lesley near the Great Bend at Muncy, Pennsylvania (Ramage et al., 1998). Local drill logs and well data indicate that subsequent paleo-outwash floods and modern fluvial processes have deposited as much as 30 meters of alluvium in this area, but little is known about the valley fill architecture and the bedrock-alluvium interface. By gaining a greater understanding of the bedrock-alluvium interface the project will not only supplement existing depth to bedrock information, but also provide information pertinent to the evolution of the Muncy Valley landscape. This project determined if variations in the thickness of the valley fill were detectable using micro-gravity techniques to map the bedrock-alluvium interface. The gravity method was deemed appropriate due to scale of the study area (~30 km2), ease of operation by a single person, and the available geophysical equipment. A LaCoste and Romberg Gravitron unit was used to collect gravitational field readings at 49 locations over 5 transects across the Muncy Creek and Susquehanna River valleys (approximately 30 km2), with at least two gravity base stations per transect. Precise latitude, longitude and ground surface elevation at each location were measured using an OPUS corrected Trimble RTK-GPS unit. Base stations were chosen based on ease of access due to the necessity of repeat measurements. Gravity measurement locations were selected and marked to provide easy access and repeat measurements. The gravimeter was returned to a base station within every two hours and a looping procedure was used to determine drift and maximize confidence in the gravity measurements. A two-minute calibration reading at each station was used to minimize any tares in the data. The Gravitron digitally recorded finite impulse response filtered gravity measurements every 20 seconds at each station. A measurement period of 15 minutes was used for each base station occupation and a minimum of 5 minutes at all other locations. Longer or multiple measurements were utilized at some sites if drift or other externalities (i.e. train or truck traffic) were effecting readings. Average, median, standard deviation and 95% confidence interval were calculated for each station. Tidal, drift, latitude, free-air, Bouguer and terrain corrections were then applied. The results show that the gravitational field decreases as alluvium thickness increases across the axes of the Susquehanna River and Muncy Creek valleys. However, the location of the gravity low does not correspond with the present-day location of the West Branch of the Susquehanna River (WBSR), suggesting that the WBSR may have been constrained along Bald Eagle Mountain by a glacial lobe originating from the Muncy Creek Valley to the northeast. Using a 3-D inversion model, the topography of the bedrock-alluvium interface was determined over the extent of the study area using a density contrast of -0.8 g/cm3. Our results are consistent with the bedrock geometry of the area, and provide a low-cost, non-invasive and efficient method for exploring the subsurface and for supplementing existing well data.

Geophysical Investigation of Subsurface Three-Dimensional Architecture of Icy Debris Fans With Ground Penetrating Radar in the Wrangell Mountains, Alaska

Icy debris fans have are newly-described landforms (Kochel and Trop, 2008 and 2012) as landforms developed immediately after deglaciation on Earth and similar features have been observed on Mars. Subsurface characteristics of Icy debris fans have not been previously investigated. Ground penetrating radar (GPR) was used to non-invasively investigate the subsurface characteristics of icy debris fans near McCarthy, Alaska, USA. The three fans investigated in Alaska are the East, West, and Middle fans (Kochel and Trop, 2008 and 2012) which below the Nabesna ice cap and on top of the McCarthy Creek Glacier. Icy debris fans in general are a largely unexplored suite of paraglacial landforms and processes in alpine regions. Recent field studies focused on direct observations and depositional processes. Their results showed that the fan's composition is primarily influenced by the type and frequency of depositional processes that supply the fan. Photographic studies show that the East Fan receives far more ice and snow avalanches whereas the Middle and West Fans receive fewer mass wasting events but more clastic debris is deposited on the Middle and West fan from rock falls and icy debris flows. GPR profiles and Wide-angle reflection and refraction (WARR) surveys consisting of both, common mid-point (CMP), and common shot-point (CSP) surveys investigated the subsurface geometry of the fans and the McCarthy Creek Glacier. All GPR surveys were collected in July of 2013 with 100MHz bi-static antennas. Four axial profiles and three cross-fan profiles were done on the West and Middle fans as well as the McCarthy Creek Glacier in order to investigate the relationship between the three features. GPR profiles yielded reflectors that were continuous for 10+ m and hyperbolic reflections in the subsurface. The depth to these reflections in the subsurface requires knowledge of the velocity of the subsurface. To find the velocity of the subsurface eight WARR surveys collected on the fans and on the McCarthy Creek glacier to provide information on variability of subsurface velocities. The profiles of the Middle and West fan have more reflections in their profiles compared to profiles done on the McCarthy Creek Glacier. Based on the WARR surveys, we interpret the lower energy return in the glacier to be caused by two reasons. 1) The increased attenuation due to wet ice versus drier ice and on the fan with GPR velocities >0.15m/ns. 2) Lack of interfaces in the glacier compared to those in the fans which are inferred to be produced by the alternating layers of stratified ice and lithic-rich layers. The GPR profiles on the West and Middle Fans show the shallow subsurface being dominated by lenticular reflections interpreted to be consistent with the shape of surficial deposits. The West Fan is distinguished from the Middle Fan by the nature of its reflections patterns and thicknesses of reflection packages that clearly shows the Middle fan with a greater thickness. The changes in subsurface reflections between the Middle and West Fans as well as the McCarthy Creek Glacier are thought to reflect the type and frequency of depositional processes and surrounding bedrock and talus slopes.

Numerical and granulometric approaches to geophysical granular flows

Volcanic ash clouds can be fed by an upward-directed eruption column (Plinian column) or by elutriation from extensive pyroclastic-flows (coignimbrite cloud). For large-scale eruptions, there is considerable uncertainty about which mechanism is dominant. Here we analyze in a novel way a comprehensive grainsize database for pyroclastic deposits. We demonstrate that the Mount Pinatubo climactic eruption deposits were substantially derived from coignimbrite clouds, and not only by a Plinian cloud as generally thought. Coignimbrite ash-fall deposits are much richer in breathable <10 m ash (5–25 wt%) than pure Plinian ash at most distances from the source volcano. We also show that coignimbrite ash clouds, as at Pinatubo, are expected to be more water rich than Plinian clouds, leading to removal of more HCl prior to stratospheric injection, thereby reducing their atmospheric impact.

Seasonal Geophysical Monitoring of Biogenic Gases in a Northern Peatland: Implications for Temporal and Spatial Variability in Free Phase Gas Production Rates

A set of high resolution surface ground penetrating radar (GPR) surveys, combined with elevation rod ( to monitor surface deformation) and gas flux measurements, were used to investigate in situ biogenic gas dynamics within a northern peatland (Caribou Bog, Maine). Gas production rates were directly estimated from the time series of GPR measurements. Spatial variability in gas production was also investigated by comparing two sites with different geological and ecological attributes, showing differences and/or similarities depending on season. One site characterized by thick highly humified peat deposits (5-6 m), wooded heath vegetation and open pools showed large ebullition events during the summer season, with estimated emissions (based on an assumed range of CH(4) concentration) between 100 and 172 g CH(4) m(-2) during a single event. The other site characterized by thinner less humified peat deposits (2-3 m) and shrub vegetation showed much smaller ebullition events during the same season (between 13 and 23 g CH(4) m(-2)). A consistent period of free-phase gas (FPG) accumulation during the fall and winter, enhanced by the frozen surficial peat acting as a confining layer, was followed by a decrease in FPG after the snow/ice melt that released estimated fluxes between 100 and 200 g CH(4) m(-2) from both sites. Estimated FPG production rates during periods of biogenic gas accumulation ranged between 0.22 and 2.00 g CH(4) m(3) d(-1) and reflected strong seasonal and spatial variability associated with differences in temperature, peat soil properties, and/or depositional attributes (e. g., stratigraphy). Periods of decreased atmospheric pressure coincided with short-period increases in biogenic gas flux, including a very rapid decrease in FPG content associated with an ebullition event that released an estimated 39 and 67 g CH(4) m(-2) in less than 3.5 hours. These results provide insights into the spatial and seasonal variability in production and emission of biogenic gases from northern peatlands.

Geophysical and Hydrological Evaluation of Two Bog Complexes in a Northern Peatland: Implications for the Distribution of Biogenic Gases at the Basin Scale

Ground penetrating radar (GPR) was used to determine peat basin geometry and the spatial distribution of free-phase biogenic gasses in two separate units of a northern peatland (Central and Southern Unit of Caribou Bog, Maine). The Central Unit is characterized by a deep basin structure (15 m maximum depth) and a raised (eccentric) bog topographic profile (up to 2 m topographic variation). Here numerous regions of electromagnetic (EM) wave scattering are considered diagnostic of the presence of extensive free-phase biogenic gas. In contrast, the Southern Unit is shallower (8 m maximum depth) and has a slightly convex upwards bog profile (less than 1 m topographic variation), and areas of EM wave scattering are notably absent. The biogenic gas zones interpreted from GPR in the Central Unit are associated with: (1) wooded heath vegetation at the surface, (2) open pools at the surface, (3) high water table elevations near the center of the basin, and (4) a region of overpressure (at approximately 5 m depth) immediately below the zone of free-phase gas accumulation. The latter suggests (1) a transient pressure head associated with low hydraulic conductivity resulting from the biogenic gasses themselves or confining layers in the peat that restrict both gas release and groundwater flow and/or (2) overpressure in the peat column as a result of the gas buildup itself. In contrast, the Southern Unit, where zones of EM scattering are absent, is characterized by: (1) predominantly shrub vegetation, (2) a lack of open pools, (3) only minor variations (less than 1 m) in water table elevation throughout the entire unit; and (4) generally upward groundwater flow throughout the basin. The results illustrate the nonuniformity of free-phase biogenic gas distribution at the peat basin scale and provide insights into the processes and controls associated with CH4 and CO2 accumulation in peatlands.

Geophysical measurements on sediment cores from North Pond, Mid Atlantic Ridge

Die Porosität ist einer der wichtigsten gesteinsphysikalischen Parameter von Sedimenten. Daher werden in dieser Arbeit die Abhängigkeiten zwischen dem zentralen gesteinsphysikalischen Parameter Porosität und den Parametern Wärmeleitfähigkeit und spezifischer elektrischer Widerstand an Sedimenten von North Pond, eines Sedimentbeckens an der westlichen Flanke des Mittelatlantischen Rückens bei ca. 23°N, qualitativ und empirisch untersucht. Die Messung des spezifischen Widerstands wurde mittels eines miniaturisierten Wenner-Arrays mit einem Zentimeter Elektrodenabstand noch an Bord der FS Maria S. Merian auf der Reise MSM11/1 vom 10. Februar bis 13. März 2009 an Sedimentkernhälften durchgeführt, ebenso wie die Messung der Wärmeleitfähigkeit mit einer Nadelsonde. Die Messung der Gamma-Dichte an einem Multi Sensor Core Logger erfolgte im Rahmen dieser Arbeit und wurde wie die Messung der Porosität an Proben im Marum in Bremen durchgeführt. Die zugrundeliegende und zu untersuchende Beziehung zwischen spezifischem Widerstand und Porosität ist Archie's Law, die untersuchte Beziehung zwischen Wärmeleitfähigkeit und Porosität ist das Modell des geometrischen Mittels der Sedimentkomponenten. Die physikalische Beziehung zwischen der Porosität und den beiden Parametern spezifischer Widerstand und Wärmeleitfähigkeit konnte an den Sedimenten von North Pond bestätigt werden.

Geophysical investigations in the European Alps

Joint interpretation of magnetotelluric and geomagnetic depth sounding data in the western European Alps offer new insights into the conductivity structure of the Earth's crust and mantle. This first large scale electromagnetic study in the Alps covers a cross-section from Germany to northern Italy and shows the importance of the alpine mountain chain as an interrupter of continuous conductors. Poor data quality due to the highly crystalline underground is overcome by Remote Reference and Robust Processing techniques. 3d-forward-modelling reveals on the one hand interrupted dipping crustal conductors with maximum conductance of 4960 S and on the other hand a lithosphere thickening up to 208 km beneath the central western Alps. Graphite networks arising from Paleozoic sedimentary deposits are considered to be accountable for the occurrence of high conductivity and the distribution pattern of crustal conductors. The influence of huge sedimentary molasse basins on the electromagnetic data is suggested to be minor compared with the influence of crustal conductors. In conclusion, electromagnetic results can be attributed to the geological, tectonic and palaeogeographical background. Dipping direction (S-SE) and maximum angle (10.1°) of the northern crustal conductor reveal the main thrusting conditions beneath the Helvetic Alps whereas the existence of a crustal conductor in the Briançonnais supports theses about its palaeographic belonging to the Iberian Peninsula.

Grain size distribution and clay mineralogy of upper Miocene sediments of ODP Hole 107-652A

The upper Miocene sedimentary sequence of Site 652, located on the lower continental margin of eastern Sardinia, was cored and logged during Ocean Drilling Program (ODP) Leg 107. Geophysical and geochemical logs from the interval 170-365 m below seafloor (mbsf), as well as various core measurements (CaCO3, grain size, X-ray diffraction), provide a mineralogical-geochemical picture that is interpreted in the framework of the climatic and tectonic evolution of the western Tyrrhenian. The results indicate the presence of short- and long-term mineralogical variations. Short-term variations are represented by calcium-carbonate fluctuations in which the amount of CaCO3 is correlated to the grain size of the sediments; coarser sediments are associated with high carbonate content and abundant detrital material. Long-term variation corresponds to a gross grain-size change in the upper part of the sequence, where predominantly fine-grained sediments may indicate a gradual deepening of the lacustrine basin towards the Pliocene. Regional climatic changes and rift-related tectonism are possible causes of this variability in the sedimentation patterns. The clay association is characterized by chlorite, illite, and smectite as dominant minerals, as well as mixed-layers clays, kaolinite, and palygorskite. Chlorite, mixed-layers clays, and illite increase at the expense of smectite below the pebble zone (335 mbsf). This is indicative of diagenetic processes related to the high geothermal gradient and to the chemistry of the evaporative pore waters, rather than to changes in the depositional environment.