Erodibility of soft freshwater sediments in Markermeer: the role of bioturbation by meiobenthic fauna

The Markermeer is a large and shallow man-made freshwater lake in the Netherlands, characterized by its high turbidity. As part of a study aiming to mitigate this high turbidity, we studied the water–bed exchange processes of the lake’s muddy bed. The upper centimeter’s–decimeter’s of the lake bed sediments mainly consists of soft anoxic mud. Recent measurements have proved the existence of a thin oxic layer on top of this soft anoxic mud. This oxic layer, which is much easier to be eroded than the anoxic mud, is believed to be related with Markermeer’s high-turbidity levels. Our hypothesis is that the thin oxic layer develops from the anoxic mud, enhanced by bioturbation. Actually, we will demonstrate that it is the bioturbated state of the bed that increases its erodability, and not the oxidation state of the sediments. In particular, we will refer to bioturbation caused by meiobenthic fauna. The objective of this study is therefore to determine the influence of the development of the thin oxic layer on the water–bed exchange processes, as well as to establish the role of bioturbation on those processes. This is done by quantifying the erosion rate as a function of bed shear stresses, and at different stages of the development of the oxic layer. Our experiments show that bioturbation increases the rate at which Markermeer sediments are eroded by almost an order of magnitude. The short-term fine sediment dynamics in Markermeer are found to be driven by the complex and highly dynamic interactions between physics, chemistry, and biology. Finally, the long-term fine sediment dynamics are driven by the erosion of the historical deposits in the lake’s bed, which is only possible after bioturbation, and which leads to an increase of the stock of sediments in the lake’s muddy bed.     

Tectono‐sedimentary evolution of the western Corinth rift (Central Greece)

The Corinth rift (Greece) is one of the world's most active rifts. The early Plio‐Pleistocene rift is preserved in the northern Peloponnese peninsula, south of the active Corinth rift. Although chronostratigraphic resolution is limited, new structural, stratigraphic and sedimentological data for an area >400 km² record early rift evolution in three phases separated by distinct episodes of extension rate acceleration and northward fault migration associated with major erosion. Minimum total N–S extension is estimated at 6.4–7.7 km. The earliest asymmetrical, broad rift accommodated slow extension (0.6–1 mm a^?1) over >3 Myrs and closed to the west. North‐dipping faults with throws of 1000–2200 m defined narrow blocks (4–7 km) with little footwall relief. A N‐NE flowing antecedent river system infilled significant inherited relief (Lower group). In the earliest Pleistocene, significant fluvial incision coincided with a 15 km northward rift margin migration. Extension rates increased to 2–2.5 mm a^?1. The antecedent rivers then built giant Gilbert‐type fan deltas (Middle group) north into a deepening lacustrine/marine basin. N‐dipping, basin margin faults accommodated throws <1500 m. Delta architecture records initiation, growth and death of this fault system over ca. 800 ka. In the Middle Pleistocene, the rift margin again migrated 5 km north. Extension rate increased to 3.4–4.8 mm a^?1. This transition may correspond to an unconformity in offshore lithostratigraphy. Middle group deltas were uplifted and incised as new hangingwall deltas built into the Gulf (Upper group). A final increase to present‐day extension rates (11–16 mm a^?1) probably occurred in the Holocene. Fault and fault block dimensions did not change significantly with time suggesting control by crustal rheological layering. Extension rate acceleration may be due to strain softening or to regional tectonic factors.     

Sea‐level and ocean‐current control on carbonate‐platform growth,Maldives, Indian Ocean

Multichannel high‐resolution seismic and multibeam data were acquired from the Maldives‐isolated carbonate platform in the Indian Ocean for a detailed characterization of the Neogene bank architecture of this edifice. The goal of the research is to decipher the controlling factors of platform evolution, with a special emphasis on sea‐level changes and changes of the oceanic currents. The stacking pattern of Lower to Middle Miocene depositional sequences, with an evolution of a ramp geometry to a flat‐topped platform, reflects variations of accommodation, which here are proposed to be primarily governed by fluctuations of relative sea level. Easterly currents during this stage of bank growth controlled an asymmetric east‐directed progradation of the bank edge. During the late middle Miocene, this system was replaced by a twofold configuration of bank development. Bank growth continued synchronously with partial bank demise and associated sediment‐drift deposition. This turnover is attributed to the onset and/or intensification of the Indian monsoon and related upwelling and occurrence of currents, locally changing environmental conditions and impinging upon the carbonate system. Mega spill over lobes, shaped by reversing currents, formed as large‐scale prograding complexes, which have previously been interpreted as deposits formed during a forced regression. On a regional scale, a complex carbonate‐platform growth can occur, with a coexistence of bank‐margin progradation and aggradation, as well as partial drowning. It is further shown that a downward shift of clinoforms and offlapping geometries in carbonate platforms are not necessarily indicative for a sea‐level driven forced regression. Findings are expected to be applicable to other examples of Cenozoic platforms in the Indo‐Pacific region.     

Mapping the Land Cover of Mexico Using AVHRR Time‐Series Data Sets

Abstract

An important methodological and analytical requirement for analyzing spatial relationships between regional habitats and species distributions in Mexico is the development of standard methods for mapping the country's land cover/land use formations. This necessarily involves the use of global data such as that produced by the Advanced Very High Resolution Radiometer (AVHRR). We created a nine‐band time‐series composite image from AVHRR Normalized Difference Vegetation Index (NDVI) bi‐weekly data. Each band represented the maximum NDVI for a particular month of either 1992 or 1993. We carried out a supervised classification approach, using the latest comprehensive land cover/vegetation map created by the Mexican National Institute of Geography (INEGI) as reference data. Training areas for 26 land cover/vegetation types were selected and digitized on the computer's screen by overlaying the INEGI vector coverage on the NDVI image. To obtain specific spectral responses for each vegetation type, as determined by its characteristic phenology and geographic location, the statistics of the spectral signatures were subjected to a cluster analysis. A total of 104 classes distributed among the 26 land cover types were used to perform the classification. Elevation data were used to direct classification output for pine‐oak and coastal vegetation types. The overall correspondence value of the classification proposed in this paper was 54%; however, for main vegetation formations correspondence values were higher (60‐80%). In order to obtain refinements in the proposed classification we recommend further analysis of the signature statistics and adding topographic data into the classification algorithm.     

The 2011 submarine volcanic eruption in El Hierro (Canary Islands)

Forty years after the Teneguía Volcano (La Palma, 1971), a submarine eruption took place off the town of La Restinga, south of El Hierro, the smallest and youngest island of the Canarian Archipelago. Precursors allowed an early detection of the event and its approximate location, suggesting it was submarine. Uncertainties derived from insufficient scientific information available to the authorities during the eruption, leading to disproportionate civil protection measures, which had an impact on the island's economy—based primarily on tourism—while residents experienced extra fear and distress.     

Soil–Water Characteristic Curves Obtained for a Colluvial and Lateritic Soil Profile Considering the Macro and Micro Porosity

This study aims to determine and to model the relationship of matric suction versus water content, known as soil?Cwater characteristic curve (SWCC), for a tropical soil profile at the southeastern Brazil. This soil profile consists of a colluvial, lateritic silty clay, with thickness of about 6.5?m. The filter paper and pressure plate techniques were used to determine the SWCC. Specimens were trimmed from the undisturbed soil samples collected along soil profile depth and submitted to drying paths. Impregnated thin-layer plates and a petrographic microscope were used to examine the structure and mineralogical composition of the soil samples. Mercury intrusion porosimetry tests were performed on some soil samples to obtain the frequency histogram of the pores. SWCC with bimodal pore-size distribution were obtained by the filter paper technique and SWCC with unimodal pore-size distribution were obtained by the pressure plate technique. The SWCC showed values of air-entry ranging from 1 to 3?kPa, depending on the technique used, resembling soils with macroporosity, as the lateritic soils. Combining both techniques, the SWCC could be fitted by a model that takes into consideration soils with macro and microporosity.     

Exponential time integrators for stochastic partial differential equations in 3D reservoir simulation

The transport of chemically reactive solutes (e.g. surfactants, CO₂ or dissolved minerals) is of fundamental importance to a wide range of applications in oil and gas reservoirs such as enhanced oil recovery and mineral scale formation. In this work, we investigate exponential time integrators, in conjunction with an upwind weighted finite volume discretisation in space, for the efficient and accurate simulation of advection–dispersion processes including non-linear chemical reactions in highly heterogeneous 3D oil reservoirs. We model sub-grid fluctuations in transport velocities and uncertainty in the reaction term by writing the advection–dispersion–reaction equation as a stochastic partial differential equation with multiplicative noise. The exponential integrators are based on the variation of constants solution and solve the linear system exactly. While this is at the expense of computing the exponential of the stiff matrix representing the finite volume discretisation, the use of real Léja point or the Krylov subspace technique to approximate the exponential makes these methods competitive compared to standard finite difference-based time integrators. For the deterministic system, we investigate two exponential time integrators, the second-order accurate exponential Euler midpoint (EEM) scheme and exponential time differencing of order one (ETD1). All our numerical examples demonstrate that our methods can compete in terms of efficiency and accuracy compared with standard first-order semi-implicit time integrators when solving (stochastic) partial differential equations that model mixing and chemical reactions in 3D heterogeneous porous media. Our results suggest that exponential time integrators such as the ETD1 and EEM schemes could be applied to typical 3D reservoir models comprising tens to hundreds of thousands unknowns.     

EVA: GPS-based extended velocity and acceleration determination

In this work, a new GPS carrier phase-based velocity and acceleration determination method is presented that extends the effective range of previous techniques. The method is named ‘EVA’, and may find applications in fields such as airborne gravimetry when rough terrain or water bodies make difficult or impractical to set up nearby GPS reference receivers. The EVA method is similar to methods such as Kennedy (Precise acceleration determination from carrier phase measurements. In: Proceedings of the 15th international technical meeting of the satellite division of the Institute of Navigation. ION GPS 2002, Portland pp 962–972, 2002b) since it uses L1 carrier phase observables for velocity and acceleration determination. However, it introduces a wide network of stations and it is independent of precise clock information because it estimates satellite clock drifts and drift rates ‘on-the-fly’, requiring only orbit data of sufficient quality. Moreover, with EVA the solution rate is only limited by data rate, and not by the available precise satellite clocks data rate. The results obtained are more robust for long baselines than the results obtained with the reference Kennedy method. An advantage of being independent of precise clock information is that, beside IGS Final products, also the Rapid, Ultra-Rapid (observed) and Ultra-Rapid (predicted) products may be used. Moreover, the EVA technique may also use the undifferenced ionosphere-free carrier phase combination (LC), overcoming baseline limitations in cases where ionosphere gradients may be an issue and very low biases are required. During the development of this work, some problems were found in the velocity estimation process of the Kennedy method. The sources of the problems were identified, and an improved version of the Kennedy method was used for this research work. An experiment was performed using a light aircraft flying over the Pyrenees, showing that both EVA and the improved Kennedy methods are able to cope with the dynamics of mountainous flight. A RTK-derived solution was also generated, and when comparing the three methods to a known zero-velocity reference the results yielded similar performance. The EVA and the improved-Kennedy methods outperformed the RTK solutions, and the EVA method provided the best results in this experiment. Finally, both the improved version of the Kennedy method and the EVA method were applied to a network in equatorial South America with baselines of more than 1,770 km, and during local noon. Under this tough scenario, the EVA method showed a clear advantage for all components of velocity and acceleration, yielding better and more robust results.