Qualitative and quantitative applications of LiDAR imagery in fluvial geomorphology

The potential for geomorphological mapping and quantitative calculations of light detection and ranging (LiDAR) data within fluvial geomorphology was studied for two river catchments within Belgium (Dijle and Amblève), which differ in physical settings and floodplain morphology. Two commercial, of‐the‐shelf LiDAR datasets with different specifications (horizontal resolution and vertical accuracy) were available for parts of the floodplains of both catchments. Real‐time kinematic (RTK) Global Positioning System (GPS) data were used as ground truth for error calculations. Qualitative analysis of LiDAR data allowed the identification of former channel patterns, levees, colluvial hillslope and fan deposits. These results were confirmed by field data, topographic surveys and historical maps. The pixel resolution proved to be an important factor in the identification of small landforms: only features with a width equal to or larger than LiDAR resolution can be detected. This poses limits on the usability of regionally available LiDAR data, which often have a horizontal resolution of several metres. The LiDAR data were also used in a quantitative analysis of channel dynamics. In the study area, the width of the Dijle River channel increased 3 m on average between 1969 and 2003. A sediment budget of channel processes for the period 1969–2003 indicated a total river bank erosion of 16·1 10³ m³ and a total within channel deposition of 7·1 10³ m³, resulting in a net river erosion of 9·0 10³ m³ or c. 0·4 Mg year^?1 per metre river length. Sequential LiDAR data can in theory be used to calculate vertical sedimentation rates, as long as there is control on the error of the reference levels used. Copyright © 2008 John Wiley and Sons, Ltd.     

The critical role of climate and saprolite weathering in landscape evolution

Landscapes evolve in response to external forces, such as tectonics and climate, that influence surface processes of erosion and weathering. Internal feedbacks between erosion and weathering also play an integral role in regulating the landscapes response. Our understanding of these internal and external feedbacks is limited to a handful of field‐based studies, only a few of which have explicitly examined saprolite weathering. Here, we report rates of erosion and weathering in saprolite and soil to quantify how climate influences denudation, by focusing on an elevation transect in the western Sierra Nevada Mountains, California. We use an adapted mass balance approach and couple soil‐production rates from the cosmogenic radionuclide (CRN) ¹⁰Be with zirconium concentrations in rock, saprolite and soil. Our approach includes deep saprolite weathering and suggests that previous studies may have underestimated denudation rates across similar landscapes. Along the studied climate gradient, chemical weathering rates peak at middle elevations (1200–2000 m), averaging 112·3 ± 9·7 t km^–2 y^–1 compared to high and low elevation sites (46·8 ± 5·2 t km^?2 y^?1). Measured weathering rates follow similar patterns with climate as those of predicted silica fluxes, modeled using an Arrhenius temperature relationship and a linear relationship between flux and precipitation. Furthermore, chemical weathering and erosion are tightly correlated across our sites, and physical erosion rates increase with both saprolite weathering rates and intensity. Unexpectedly, saprolite and soil weathering intensities are inversely related, such that more weathered saprolites are overlain by weakly weathered soils. These data quantify exciting links between climate, weathering and erosion, and together suggest that climate controls chemical weathering via temperature and moisture control on chemical reaction rates. Our results also suggest that saprolite weathering reduces bedrock coherence, leading to faster rates of soil transport that, in turn, decrease material residence times in the soil column and limit soil weathering. Copyright © 2009 John Wiley & Sons, Ltd.     

Rain properties controlling soil splash detachment

Laboratory experiments have been conducted to study the effects of various rain properties on sand detachment resulting from raindrop impact. Splash cups were exposed to simulated rainfall intensities ranging between 10 and 140 mm h⁻¹. The detached sand was collected and weighed whereas rain intensity, equivalent drop diameter and fall velocity of raindrops were measured with an optical spectro‐pluviometer (OSP). The properties of the simulated rain (i.e. median volume diameter and kinetic energy) were compared with those observed in natural conditions. Statistical analyses have been undertaken in order to evaluate which rain variable best predicts the mass of sand detached. Linear and non‐linear correlations between the mass of detached sediment and the product of drop size (d) by drop velocity (v), i.e. D^αV^β, with values of α varying between 1 and 6 and β between 0 and 3, have been computed. The results indicate that the coefficient of determination (R²) for α ranging between 3 and 5 and β lower or equal to 2 are satisfying. Although kinetic energy (D³V²) described splash detachment relatively well, the product of momentum by drop diameter (D⁴V) was slightly superior in describing splash detachment. Therefore, the momentum multiplied by the drop diameter is recommended as the best rain variable to describe splash detachment. Copyright © 2000 John Wiley & Sons, Ltd.     

Characterization of seasonal local recharge using electrical resistivity tomography and magnetic resonance sounding

A groundwater recharge process of heterogeneous hard rock aquifer in the Moole Hole experimental watershed, south India, is being studied to understand the groundwater flow behaviour. Significant seasonal variations in groundwater level are observed in boreholes located at the outlet area indicating that the recharge process is probably taking place below intermittent streams. In order to localize groundwater recharge zones and to optimize implementation of boreholes, a geophysical survey was carried out during and after the 2004 monsoon across the outlet zone. Magnetic resonance soundings (MRS) have been performed to characterize the aquifer and measure groundwater level depletion. The results of MRS are consistent with the observation in boreholes, but it suffers from degraded lateral resolution. A better resolution of the regolith/bedrock interface is achieved using electrical resistivity tomography (ERT). ERT results are confirmed by resistivity logging in the boreholes. ERT surveys have been carried out twice—before and during the monsoon—across the stream area. The major feature of recharge is revealed below the stream with a decrease by 80% of the calculated resistivity. The time‐lapse ERT also shows unexpected variations at a depth of 20 m below the slopes that could have been interpreted as a consequence of a deep seasonal water flow. However, in this area time‐lapse ERT does not match with borehole data. Numerical modelling shows that in the presence of a shallow water infiltration, an inversion artefact may take place thus limiting the reliability of time‐lapse ERT. A combination of ERT with MRS provides valuable information on structure and aquifer properties respectively, giving a clue for a conceptual model of the recharge process: infiltration takes place in the conductive fractured‐fissured part of the bedrock underlying the stream and clayey material present on both sides slows down its lateral dissipation. Copyright © 2007 John Wiley & Sons, Ltd.     

Widespread occurrence of aquifers currently undetectable with the MRS technique in the Grenville geological province, Canada

In the summer of 2003, a series of field tests were made over aquifers within or next to the Grenville geological province, Canada. Twelve sites distributed within three general areas of the province and nearby to its southern border were surveyed with the MRS (Magnetic Resonance Sounding) technique. Each site was characterized by saturated shallow water bearing layers of high effective porosity. The tests were planned to evaluate the MRS technology in such environments for groundwater resources quantification and aquifer/aquitard media characterization. A priori, the environmental conditions in Eastern Canada were assumed favorable for MRS work such as low to moderate ambient noise, earth's field value above 55,000 nT, targets consisting of shallow and high porosity aquifers. Despite such conditions, no clear MRS responses were observed. Following supplementary investigations, this is now attributed to the wide-scale occurrence of disseminated magnetite within the aquifers which shorten the signal decay time. Although the lack of MRS response over aquifers has been reported before, this was not the case for such widespread occurrence. This also supports the adaptation of the MRS technology to conditions similar to those observed in the Grenville province.     

Rotational bulge and one plume convection pattern: Influence on Martian true polar wander

Motion of the entire solid planet with respect to its spin axis have been proposed on Mars. This movement is known as True Polar Wander (TPW). According to the conservation of angular momentum with no external torque, on geological time scales the axis of maximum inertia of a planet is aligned with the rotation axis. Then rearrangement of masses within the mantle disturbs the planet's inertia and induces TPW. The convection pattern on Mars is possibly controlled by a sequence of single plumes originating from the core-mantle boundary. Using a homogeneous model of the martian mantle and modelling the plume as a sphere, we calculate the inertial tensor perturbations caused by the plume mass anomaly. We investigate the stabilizing influence of the remnant rotational bulge due to the lithosphere elasticity on these perturbations. It appears that, during early martian history, the elastic lithosphere was thin enough to allow its fractures under the inertia perturbations induced by a hot plume. Consequently, the lithosphere's behaviour became effectively viscoelastic and the plume could induce large TPW. We conclude that one plume convection pattern should have greatly influenced the rotation pole behaviour during early Mars history: around 4 Gyr ago, Mars already could have experienced two TPW events lasting possibly only a few million years each. We then compare our scenario with others already published in the literature.     

The effects of hydrothermal eruptions and a tectonic earthquake on a cycling crater lake (Inferno Crater Lake,Waimangu, New Zealand)

The lake level and temperature of Inferno Crater Lake, a hot lake in the Waimangu Geothermal Field, vary in quasi-regular cycles with a period of about 30 to 50 days. Between 1970 and 1991, several hydrothermal eruptions at other sites within the Waimangu system affected this cycling, especially the cycle period. This indicated that although the cycling was most evident in one small lake, it was sensitive to any changes in the surrounding geothermal system. After the Raupo Pond eruption, which produced the largest effects, the recovery of the temperature and cycle period was consistent with extra energy diffusing out of the previously proposed heat pipe system. The 1987 Edgecumbe earthquake, whose hypocenter was 50 km from Waimangu, also seems to have significantly affected the hydrothermal system, reducing the cycle period for at least a year afterwards.     

Merging active and passive data sets in traveltime tomography: The case study of Campi Flegrei caldera (Southern Italy)

We propose a strategy for merging both active and passive data sets in linearized tomographic inversion. We illustrate this in the reconstruction of 3D images of a complex volcanic structure, the Campi Flegrei caldera, located in the vicinity of the city of Naples, southern Italy. The caldera is occasionally the site of significant unrests characterized by large ground uplifts and seismicity. The P and S velocity models of the caldera structure are obtained by a tomographic inversion based on travel times recorded during two distinct experiments. The first data set is composed of 606 earthquakes recorded in 1984 and the second set is composed of recordings for 1528 shots produced during the SERAPIS experiment in 2001. The tomographic inversion is performed using an improved method based on an accurate finite‐difference traveltime computation and a simultaneous inversion of both velocity models and earthquake locations. In order to determine the adequate inversion parameters and relative data weighting factors, we perform massive synthetic simulations allowing one to merge the two types of data optimally. The proper merging provides high resolution velocity models, which allow one to reliably retrieve velocity anomalies over a large part of the tomography area. The obtained images confirm the presence of a high P velocity ring in the southern part of the bay of Pozzuoli and extends its trace inland as compared to previous results. This annular anomaly represents the buried trace of the rim of the Campi Flegrei caldera. Its shape at 1.5 km depth is in good agreement with the location of hydrothermalized lava inferred by gravimetric data modelling. The Vp/Vs model confirms the presence of two characteristic features. At about 1 km depth a very high Vp/Vs anomaly is observed below the town of Pozzuoli and is interpreted as due to the presence of rocks that contain fluids in the liquid phase. A low Vp/Vs body extending at about 3–4 km depth below a large part of the caldera is interpreted as the top of formations that are enriched in gas under supercritical conditions.     

Potential and limits of double-difference tomographic methods

Double-difference tomographic methods use directly accurate time delays computed between similar signals. Such methods are designed to image very heterogeneous media, such as volcanoes or fault zones. In seismological applications, similar signals are recorded at a given station from earthquakes sharing similar and close-by sources. In seismic exploration experiments, similar signals are often recorded at neighbouring receivers. After a brief presentation of the tomographic algorithm used, a seismological application is summarized. The potential and limits of double-difference tomographic methods are explored using various numerical experiments. They show that two effects are competing in double-difference tomography: (i) the degradation of the stability of the inversion due to the geometrical proximity of the rays used in the differentiation and (ii) the decrease in modelling error, which allows improving the stability of the inversion and using smaller quantities of a priori information when data are sufficiently accurate. The best resolution is obtained for an optimal value of the inter-source or inter-receiver distance. For optimal values of these distances and a priori information, tomography using traveltime differences provides significantly better resolved results than using traveltimes.     

Dynamics and fluxes of organic carbon and nitrogen in two Guiana Shield river basins impacted by deforestation and mining activities

Deforestation and mining activities have proven to be very damaging to rivers because these activities disturb the environmental characteristics of rivers. Thus, the concentrations of dissolved organic carbon (DOC), particulate organic carbon (POC), particulate nitrogen (PN), and Chlorophyll‐a (Chl‐a) were measured monthly during 2 hydrological years in the Maroni and Oyapock Rivers to assess the dynamics and fluxes of organic carbon and nitrogen in these 2 Guiana Shield basins, which have been strongly (Maroni) and weakly (Oyapock) impacted by deforestation and mining activities. The 2‐year time series show that DOC, POC, PN, and Chl‐a concentrations vary seasonally with discharge in both rivers, indicating a hydrologically dominated control. Temporal patterns of DOC, POC, and PN indicate that these variables show maximum concentrations in rising waters due to the yield of organic matter and nitrogen accumulated in soils, which are incorporated into the rivers during rainfall. However, the Chl‐a concentrations were at a maximum during low‐water stages. The C/N and C/Chl‐a ratios also showed a seasonal trend, with lower values during the low water periods due to an increase in algal biomass. During high water, the POC in both rivers is the result of terrestrial organic matter, whereas during low‐water autochthonous organic matter can reach up to 34% of the POC. The mean annual fluxes of TOC and PN were higher (4.56 × 10⁵ tonC year^?1 and 1.77 × 10⁴ tonN year^?1, respectively) in the Maroni River than those (1.84 × 10⁵ tonC year^?1 and 0.54 × 10⁴ tonN year^?1, respectively) in the Oyapock River. However, the specific fluxes of DOC, POC, and PN from both basins were nearly the same. Although gold mining activities are performed in both basins, there is no conclusive evidence regarding the impact of these activities on the dynamics of organic matter and particulate nitrogen in the Maroni and Oyapock Rivers.