3D numerical simulation of particle-particle collisions in saltation mode near stream beds

The importance of particle-particle collisions in sediment saltation in the bed-load layer is analyzed herein by means of numerical simulation. The particle saltation theoretical/numerical model follows a Lagrangian approach, and addresses the motion of sediment particles in an open channel flow described by a logarithmic velocity profile. The model is validated with experimental data obtained from the literature. In order to evaluate the importance of the phenomenon, simulations with and without particle-particle collisions were carried out. Results for two different sediment concentrations are presented, namely 0.13% and 2.33%. For each concentration of particles, three different flow intensities were considered, and trajectories of two different particle sizes, within the sand range were computed. Changes in particle rotation, particle velocity, and angle of trajectory before and after particle-particle collisions appear to be relatively important at lower shear stresses, whereas they decrease in significance with increasing flow intensities. Analyses of the evolution in time of the second order moment of particle location suggest that inter-particle collisions introduce transverse diffusion in saltating particles in the span-wise direction.     

Integrated provenance analysis of a convergent retroarc foreland system: U–Pb ages,heavy minerals,Nd isotopes,and sandstone compositions of the Middle Magdalena Valley basin,northern Andes,Colombia

Sediment provenance analysis remains a powerful method for testing hypotheses on the temporal and spatial evolution of uplifted source regions, but issues such as recycling, nonunique sources, and pre- and post-depositional modifications may complicate interpretation of results from individual provenance techniques. Convergent retroarc systems commonly contain sediment sources that are sufficiently diverse (continental magmatic arc, fold–thrust belt, and stable craton) to enable explicit provenance assessments. In this paper, we combine detrital zircon U–Pb geochronology, heavy mineral identification, Nd isotopic analyses, conventional sandstone petrography, and paleocurrent measurements to reconstruct the clastic provenance history of a long-lived sedimentary basin now exposed in an intermontane zone of the northern Andean hinterland of Colombia. The Middle Magdalena Valley basin, situated between the Central Cordillera and Eastern Cordillera, contains a 5–10 km-thick succession of Upper Cretaceous to Quaternary fill. The integrated techniques show a pronounced change in provenance during the Paleocene transition from the lower to upper Lisama Formation. We interpret this as a shift from an eastern cratonic source to a western Andean source composed of magmatic-arc rocks uplifted during initial shortening of the Central Cordillera. The appearance of detrital chloritoid and a shift to more negative ε_Nd(t=0) values in middle Eocene strata of the middle La Paz Formation are attributed to shortening-related exhumation of a continental basement block (La Cira–Infantas paleohigh), now buried, along the axis of the Magdalena Valley. The diverse provenance proxies also show distinct changes during middle to late Eocene deposition of the Esmeraldas Formation that likely reflect initial rock uplift and exhumation of the fold–thrust belt defining the Eastern Cordillera. Upsection, detrital zircon U–Pb ages and heavy mineral assemblages for Oligocene and younger clastic deposits indicate that the Mesozoic sedimentary cover of the Eastern Cordillera was recycled during continued Cenozoic shortening. Our multidisciplinary provenance study refines the tectonic history of the Colombian Andes and demonstrates that uncertainties related to sediment recycling, nonunique sources, source heterogeneity, and climate in interpreting provenance data can be minimized via an integrated approach.     

Integrated geophysics and soil gas profiles as a tool to characterize active faults: the Amer fault example (Pyrenees, NE Spain)

The combined use of geophysical and soil gas composition exploration methods allows to rapidly obtain at relative low cost information that might be related to seismic activity conditions. In this study, we carried out geochemical soil gas sampling (²²²Rn, ²²⁰Rn and CO₂), electrical resistivity tomography and seismic refraction profiles in two selected zones near the town of Amer in the Spanish Pyrenees, where the presence of recent fractures is evident in the field. Data analysis clearly reveals anomalous values for each gas at specific positions along the electrical imaging transects. Geomorphologic and hydrogeologic data and the integration of geophysical data and soil gas measurements indicate that: (1) endogene gases radon (²²²Rn) and carbon dioxide (CO₂) are released from the meta-sedimentary basement rocks across the main fractured zones with higher permeability values, while lower Cenozoic detrital sedimentary formations act as an impervious boundary; (2) sites with highest radon concentrations (52?kBq?m^?3) coincide with the zones in the Amer fault showing more recent geomorphic evidence of activity, and more specifically with those areas covered by thinner surficial formations; (3) the lowest ²²²Rn values (0.2?C0.4?kBq?m^?3) were recorded just on the master active fault plane. This pattern could be explained by a dilution effect resulting from high rates of soil CO₂ efflux (267?g?m^?2?day^?1); (4) soil thoron (²²⁰Rn) activity is maximum (143?kBq?m^?3) in areas with high surficial fracturing; (5) groundwater pumping may cause important distortions in the natural flow dynamics and in the measured concentrations of gases. The agreement between the different data (geochemical, geophysical, and hydrogeological) and field observations (geology and geomorphology) leads us to propose a preliminary tectonic-gravitational model for the study area.     

Evaluating restoration of man‐made slopes: a threshold approach balancing vegetation and rill erosion

The management of reclaimed slopes derived from industrial and civil activities (e.g. surface mining and road construction) requires the development of practical stability analysis approaches that integrate the processes and mechanisms that rule the dynamics of these ubiquitous emerging ecosystems. This work describes a new modelling approach focused on stability analysis of water‐limited reclaimed slopes, where interactive relationships between rill erosion and vegetation regulate ecosystem stability. Our framework reproduces two main groups of possible trends along the temporal evolution of reclaimed slopes: successful trends, characterized by widespread vegetation development and the effective control of rill erosion processes; and gullying trends, characterized by the progressive loss of vegetation and a sharp logistic increase in erosion rates. Furthermore, this analytical approach allows the determination of threshold values for the state variables (i.e. vegetation cover and rill erosion) that drive the system's stability, facilitating the identification of critical situations that require specific human intervention (e.g. revegetation or, in very problematic cases, revegetation combined with rill network destruction) to ensure the long‐term sustainability of the restored ecosystem. The application of our threshold analysis framework in Mediterranean‐dry reclaimed slopes derived from surface coal mining (the Teruel coalfield in central‐eastern Spain) showed a good field‐based performance. Therefore, we believe that this model is a valuable contribution for the management of water‐limited reclaimed systems, including those associated with rill erosion, as it provides a tool for the evaluation of restoration success and can play an important role in decision‐making during ecosystem restoration in severely disturbed landscapes. Copyright © 2011 John Wiley & Sons, Ltd.     

Two hundred years of environmental change in Picos de Europa National Park inferred from sediments of Lago Enol, northern Iberia

We present a study of two short sediment cores recovered from Lago Enol, in the Picos de Europa National Park, Cantabrian Mountains, northern Iberia. We inferred past climate conditions and anthropogenic impacts using geochemical and biological (pollen and diatoms) variables in the dated sequences, in conjunction with temperature and precipitation data collected since 1871 at meteorological stations in the region. The record provides evidence of environmental changes during the last 200 years. At the end of the Little Ice Age (~1800?C1875 AD) the region was characterized by an open landscape. Long-term use of the area for mixed livestock grazing in the mountains, and cultivation of rye during the nineteenth century, contributed to the expansion of grassland at the expense of forest. Warmer temperatures since the end of the nineteenth century are inferred from a change in diatom assemblages and development of the local forest. Socioeconomic transformation during the twentieth century, such as livestock changes related to dairy specialization, planting of non-native trees, mining activities, and management of the national park since its creation in 1918, caused profound changes in the catchment and in the lake ecology. The last several decades (~1970?C2007 AD) of the Lago Enol sediment record are strikingly different from previous periods, indicating lower runoff and increasing lake productivity, particularly since AD 2000. Today, the large number of tourists who visit the area cause substantial impacts on this ecosystem.     

Geochemical processes in a Mediterranean Lake: a high-resolution study of the last 4,000 years in Zoñar Lake,southern Spain

High-resolution geochemical analysis of a 6-m-long sediment core from Zoñar Lake, southern Spain, provides a detailed characterization of major changes in lake and watershed processes during the last 4,000 years. Geochemical variables were used as paleolimnological indicators and complement Zoñar Lakes’s paleoenvironmental reconstruction based on sedimentological and biological proxies, which define periods of increasing allochthonous input to the lake and periods of dominant autochthonous sedimentation. Chemical ratios identify periods of endogenic carbonate formation (higher Ca/Al, Sr/Al and Ba/Al ratios), evaporite precipitation (higher S/Al, Sr/Al ratios), and anoxic conditions (higher Mo/Al, U/Th ratios and Eu anomaly). Higher productivity is marked by elevated organic carbon content and carbonate precipitation (Mg/Ca). Hydrological reconstruction for Zoñar Lake was based on sedimentological, mineralogical and biological proxies, and shows that lower lake levels are characterized by Sr-rich sediments (a brackish lake with aragonite) and S-rich sediments (a saline lake with gypsum), while higher lake levels are characterized by sediments enriched in elements associated with alumino-silicates (Al, K, Ti, Fe, trace and rare earth elements), reflecting fresher conditions. Geochemical indicators also mark periods of higher detrital input to the lake related to human activity in the watershed: (1) during the Iberian Roman Humid Period (650 BC–AD 300), around the onset of the Little Ice Age (AD 1400), during the relatively drier Post-Roman and Middle Ages (AD 800–1400), and over the last 50 years, due to mechanized farming practices. Heavy metal enrichment in the sediments (Cu and Ni) suggests intensification of human activities during the Iberian Roman Period, and the use of fertilizers during the last 50 years.     

Interhemispheric Linkage of Paleoclimate During the Last Glaciation

Combined glacial geologic and palynologic data from the southern Lake District, Seno Reloncaví, and Isla Grande de Chiloé in middle latitudes (40°35’–42°25’S) of the Southern Hemisphere Andes suggest (1) that full-glacial or near-full-glacial climate conditions persisted from about 29,400 to 14,550 ¹⁴C yr BP in late Llanquihue time, (2) that within this late Llanquihue interval mean summer temperature was depressed 6°–8°C compared to modern values during major glacier advances into the outer moraine belt at 29,400, 26,760, 22,295–22,570, and 14,550–14,805 ¹⁴C yr BP , (3) that summer temperature depression was as great during early Llanquihue as during late Llanquihue time, (4) that climate deteriorated from warmer conditions during the early part to colder conditions during the later part of middle Llanquihue time, (5) that superimposed on long-term climate deterioration are Gramineae peaks on Isla Grande de Chiloé that represent cooling at 44,520–47,110 ¹⁴C yr BP (T-11), 32,105–35,764 ¹⁴C yr BP (T-9), 24,895–26,019 ¹⁴C yr BP (T-7), 21,430–22,774 ¹⁴C yr BP (T-5), and 13,040–15,200 ¹⁴C yr BP (T-3), (6) that the initial phase of the glacial/interglacial transition of the last termination involved at least two major steps, one beginning at 14,600 ¹⁴C yr BP and another at 12,700–13,000 ¹⁴ C yr BP , and (7) that a late-glacial climate reversal of ≥2–3° C set in close to 12,200 ¹⁴C yr BP , after an interval of near-interglacial warmth, and continued into Younger Dryas time. The late-glacial climate signal from the southern Chilean Lake District ties into that from proglacial Lago Mascardi in the nearby Argentine Andes, which shows rapid ice recession peaking at 12,400 ¹⁴C yr BP , followed by a reversal of trend that culminated in Younger-Dryas-age glacier readvance at 11,400–10,200 ¹⁴C yr BP . Many full- and late-glacial climate shifts in the southern Lake District match those from New Zealand at nearly the same Southern Hemisphere middle latitudes. At the last glacial maximum (LGM), snowline lowering relative to present-day values was nearly the same in the Southern Alps (875 m) and the Chilean Andes (1000 m). Particularly noteworthy are the new Younger-Dryas-age exposure dates of the Lake Misery moraines in Arthur's Pass in the Southern Alps. Moreover, pollen records from the Waikato lowlands on North Island show that a major vegetation shift at close to 14,700 ¹⁴C yr BP marked the beginning of the last glacial/interglacial transition (Newnham et al. 1989). The synchronous and nearly uniform lowering of snowlines in Southern Hemisphere middle-latitude mountains compared with Northern Hemisphere values suggests global cooling of about the same magnitude in both hemispheres at the LGM. When compared with paleoclimate records from the North Atlantic region, the middle-latitude Southern Hemisphere terrestrial data imply interhemispheric symmetry of the structure and timing of the last glacial/interglacial transition. In both regions atmospheric warming pulses are implicated near the beginning of Oldest Dryas time (～14,600 ¹⁴C yr BP) and near the Oldest Dryas/Bölling transition (～12,700–13,000 ¹⁴ C yr BP ). The second of these warming pulses was coincident with resumption of North Atlantic thermohaline circulation similar to that of the modern mode, with strong formation of Lower North Atlantic Deep Water in the Nordic Seas. In both regions, the maximum Bölling-age warmth was achieved at 12,200–12,500 ¹⁴ C yr BP , and was followed by a reversal in climate trend. In the North Atlantic region, and possibly in middle latitudes of the Southern Hemisphere, this reversal culminated in a Younger-Dryas-age cold pulse. Although changes in ocean circulation can redistribute heat between the hemispheres, they cannot alone account either for the synchronous planetary cooling of the LGM or for the synchronous interhemispheric warming steps of the abrupt glacial-to-interglacial transition. Instead, the dominant interhemispheric climate linkage must feature a global atmospheric signal. The most likely source of this signal is a change in the greenhouse content of the atmosphere. We speculate that the Oldest Dryas warming pulse originated from an increase in atmospheric water-vapor production by half-precession forcing in the tropics. The major thermohaline switch near the Oldest Dryas/Bölling transition then couldhave triggered another increase in tropical water-vapor production to near-interglacial values.     

Intercomparison of coastal area morphodynamic models

Five different coastal area morphodynamic models have been set up to run on the same offshore breakwater layout and an intercomparison carried out on the hydrodynamic and morphodynamic output produced by each scheme. In addition, the predicted morphodynamics was checked against available laboratory and field data.It is concluded that the models are capable of producing realistic estimates for the dominant morphodynamic features associated with offshore breakwaters. Coupling of the wave, current and sediment transport components of each scheme is shown to yield bathymetry which attains a state of equilibrium, unlike models which are based on the initial transport field only.