From a river valley to a ria: Evolution of an incised valley (Ría de Ferrol,north‐west Spain) since the Last Glacial Maximum

The evolution of incised valleys is an important area of research due to the invaluable data it provides on sea‐level variations and depositional environments. In this article the sedimentary evolution of the Ría de Ferrol (north‐west Spain) from the Last Glacial Maximum to the present is reconstructed using a multidisciplinary approach, combining seismic and sedimentary facies, and supported by radiocarbon data and geochemical proxies to distinguish the elements of sedimentary architecture within the ria infill. The main objectives are: (i) to analyse the ria environment as a type of incised valley to evaluate the response of the system to the different drivers; (ii) to investigate the major controlling factors; and (iii) to explore the differentiation between rias and estuaries. As a consequence of the sea‐level rise subsequent to the Last Glacial Maximum (ca 20 kyr bp ), an extensive basin, drained by a braided palaeoriver, evolved into a tide‐dominated estuary and finally into a ria environment. Late Pleistocene and Holocene high‐frequency sea‐level variations were major factors that modulated the type of depositional environments and their evolution. Another major modulating factor was the antecedent morphology of the ria, with a rock‐incised narrow channel in the middle of the basin (the Ferrol Strait), which influenced the evolution of the ria as it became flooded during Holocene transgression. The strait acted as a rock‐bounded ‘tidal inlet’ enhancing the tidal erosion and deposition at both ends, i.e. with an ebb‐tidal delta in the outer sector and tidal sandbanks in the inner sector. The final step in the evolution of the incised valley into the modern‐defined ria system was driven by the last relative sea‐level rise (after 4 kyr bp ) when the river mouths retreated landward and a single palaeoriver was converted into minor rivers and streams with scattered mouths in an extensive coastal area.     

Topographical controls on small‐volume pyroclastic flows

Emplacement of small‐volume (<0·1 km³) pyroclastic flows is significantly influenced by topography. The Arico ignimbrite on Tenerife (Canary Islands) is a characteristic small‐volume pyroclastic flow deposit emplaced on high relief topography. The pyroclastic flow flowed down pre‐existing valleys on the southern slopes of the island. In proximal areas deep (up to 100 m) valleys acted as efficient conduits for the pyroclastic flow, which was mostly channelled; in this particular area the ignimbrite corresponds to a homogeneous, moderately welded deposit, consisting of flattened pumices in an abundant ashy matrix with a relatively low lithic fragment content. In intermediate zones significant changes occur in the steepness of the slope and, although still channelled, here the pyroclastic flow was influenced by hydraulic jumps. In this area, two different units can be clearly distinguished in the ignimbrite: the lower unit is composed of a lithic‐rich ground‐layer deposit that formed at the turbulent, highly concentrated head of the flow; the upper unit consists of a well welded pumice‐rich deposit that occasionally reveals a basal layer formed by shearing with the lower part. This division into two units is maintained as far as distal areas near the present‐day coastline, where the slope is very gentle or null and the ignimbrite is not channelled. The ground layer is not found in distal areas. The ignimbrite here only consists of the upper unit, which is occasionally repeated due to a surging process provoked by the lower flow speed, as the pyroclastic flow spread out of the channelled zone. A theoretical model on how topography controlled the deposition of the Arico ignimbrite is derived by interpreting the observed lithological and sedimentological variations in terms of changes in topography and bedrock morphology. This new model is of general applicability and will help to explain other deposits of similar characteristics.     

ANALYSIS OF PILES USED FOR SLOPE STABILIZATION

Piles used for the stabilization of slopes have to be adequately designed to resist the induced lateral loads due to the movement of the unstable slope. In this paper, a numerical method is presented for the analysis of this problem. In this approach, the piles are modelled using beam finite elements. The soil response at the individual piles is modelled using the modulus of subgrade reaction and pile–soil–pile interaction considered using the theory of elasticity. Two case histories, one for single pile and the other for pile group, are analysed which show that the numerical model can predict the general characteristics of the piles reasonably well. The study suggests that the design of the piles based on the computed response from single pile analysis, ignoring group effects, may be unduly conservative.     

Guandishan Granitoids of the Paleoproterozoic Lüliang Metamorphic Complex in the Trans-North China Orogen: SHRIMP Zircon Ages, Petrogenesis and Tectonic Implications

Abstract: The Paleoproterozoic Lüliang Metamorphic Complex (PLMC) is situated in the middle segment of the western margin of the Trans-North China Orogen (TNCO), North China Craton (NCC). As the most important lithological assemblages in the southern part of the PLMC, Guandishan granitoids consist of early gneissic tonalities, granodiorites and gneissic monzogranites, and younger gneissic to massive monzogranites. Petrochemical features reveal that the early gneissic tonalities and granodiorites belong to the medium-K calc-alkaline series; the early gneissic monzogranites are transitional from high-K calc-alkaline to the shoshonite series; the younger gneissic to massive monzogranites belong to the high-k calc-alkaline series, and all rocks are characterized by right-declined REE patterns and negative Nb, Ta, Sr, P, and Ti anomalies in the primitive mantle normalized spidergrams. SHRIMP zircon U–Pb isotopic dating reveals that the early gneissic tonalities and granodiorites formed at ~2.17 Ga, the early gneissic monzogranites at ~2.06 Ga, and the younger gneissic to massive monzogranites at ~1.84 Ga. Sm–Nd isotopic data show that the early gneissic tonalities and granodiorites have εNd(t) values of +0.48 to ?3.19 with Nd-depleted mantle model ages (TDM) of 2.76–2.47 Ga, and early gneissic monzogranites have εNd(t) values of ?0.53 to ?2.51 with TDM of 2.61–2.43 Ga, and the younger gneissic monzogranites have εNd(t) values of ?6.41 to ?2.78 with a TDM of 2.69–2.52 Ga.These geochemical and isotopic data indicate that the early gneissic tonalities, granodiorites, and monzogranites were derived from the partial melting of metamorphosed basaltic and pelitic rocks, respectively, in a continental arc setting. The younger gneissic to massive monzogranites were derived by partial melting of metamorphosed greywackes within the continental crust. Combined with previously regional data, we suggest that the Paleoproterozoic granitoid magmatism in the Guandishan granitoids of the PLMC may provide the best geological signature for the complete spectrum of Paleoproterozoic geodynamic processes in the Trans-North China Orogen from oceanic subduction, through collisional orogenesis, to post-orogenic extension and uplift.     

East Asia winter monsoon variations on a millennial time‐scale before the last glacial–interglacial cycle

Grain size and magnetic susceptibility measurements on samples from a typical loess–palaeosol sequence on the central Chinese Loess Plateau are used to reconstruct the Pleistocene East Asian monsoon climate. The coarse‐grained fraction, i.e. the weight percentage > 30 μm of the bulk grain‐size distribution, is used as a sensitive proxy index of the East Asia winter monsoon strength. On the basis of an absolute time‐scale, time‐series variations of this proxy show that winter monsoon strengths varied on millennial time‐scales during the periods 145–165, 240–280, 320–350, 390–440, 600–640, 860–890, 900–930 and 1330–1400 kyr BP. The wavelength of these climatic oscillations varied between 1.89 and 4.0 kyr, as is shown by spectral analysis using the multitaper method. Although numerical simulation experiments show that high frequencies also can arise from measurement errors in the grain‐size analysis, the frequencies prove to be sufficiently stable when the spectral analysis is repeated with a different number of tapers. For the time being, we do not correlate these climatic oscillations with palaeoclimatic records in the North Atlantic deep‐sea sediments because both time‐scales need to be further improved. Our data, however, certainly demonstrate that millennial‐scale East Asian winter monsoon variations in the last 1.4 million years can be detected from terrestrial loess records. Copyright © 1999 John Wiley & Sons, Ltd.     

交叉裂隙岩体裂纹扩展试验及混合有限-离散元数值模拟研究

对预制交叉裂隙岩石试样进行裂纹扩展试验,研究了裂纹萌生、扩展、聚合过程,分析了主裂隙和轴向载荷夹角及主次裂隙夹角对裂纹起裂应力和聚合应力的影响,并用混合有限元-离散元程序,即图形处理器并行化的3D Y-HFDEM代码对试验进行了仿真计算,实现了岩石破坏从连续介质向非连续介质的过渡,对裂纹类及损伤破坏模式进行了识别,捕捉到了试验中难以发现的现象。研究表明：随主裂隙与轴向载荷夹角增加,裂纹聚合区的拉伸裂纹数量增加;裂纹起裂和聚合应力与主裂隙与轴向载荷夹角成正比;主次裂隙夹角增加,岩石的破坏模式由拉伸破坏转为剪切破坏,交叉裂隙加剧岩石破碎程度;主裂隙尖端萌生扩展的拉伸-剪切混合裂缝引起的破坏在岩石破坏中占主导地位,是导致岩体失去承载能力的主控裂纹;混合有限元-离散元仿真软件GPGPU并行化的3D Y-HFDEM IDE在岩石裂纹扩展研究中具有优势,可以捕捉实验室难以发现的损伤断裂类型,可以作为岩石裂纹扩展研究的有力工具。     

Controls on low-pressure anatexis

Low-pressure anatexis, whereby rocks melt in place after passing through the andalusite stability field, develops under more restricted conditions than does low-pressure metamorphism. Our thermal modelling and review of published work indicate that the following mechanisms, operating alone, may induce anatexis in typical pelitic rocks without inducing wholesale melting in the lower crust: (i) magmatic advection by pervasive flow; (ii) crustal-scale detachment faulting; and (iii) the presence of a high heat-producing layer. Of these, only magmatic advection by pervasive flow and crustal-scale detachment faulting have been shown quantitatively to provide sufficient heat to cause widespread melting. Combinations of the above mechanisms with pluton-scale magmatic advection, shear heating, removal of the lithospheric mantle, or with each other provide additional means of developing suitable high temperatures at shallow crustal levels to generate low-pressure anatexis.     

Climate change impacts on groundwater recharge- uncertainty, shortcomings, and the way forward?

An integrated approach to assessing the regional impacts of climate and socio-economic change on groundwater recharge is described from East Anglia, UK. Many factors affect future groundwater recharge including changed precipitation and temperature regimes, coastal flooding, urbanization, woodland establishment, and changes in cropping and rotations.Important sources of uncertainty and shortcomings in recharge estimation are discussed in the light of the results. The uncertainty in, and importance of, socio-economic scenarios in exploring the consequences of unknown future changes are highlighted. Changes to soil properties are occurring over a range of time scales, such that the soils of the future may not have the same infiltration properties as existing soils. The potential implications involved in assuming unchanging soil properties are described.To focus on the direct impacts of climate change is to neglect the potentially important role of policy, societal values and economic processes in shaping the landscape above aquifers. If the likely consequences of future changes of groundwater recharge, resulting from both climate and socio-economic change, are to be assessed, hydrogeologists must increasingly work with researchers from other disciplines, such as socio-economists, agricultural modellers and soil scientists.

A way to synchronize models with seismic faults for earthquake forecasting: Insights from a simple stochastic model

Numerical models are starting to be used for determining the future behaviour of seismic faults and fault networks. Their final goal would be to forecast future large earthquakes. In order to use them for this task, it is necessary to synchronize each model with the current status of the actual fault or fault network it simulates (just as, for example, meteorologists synchronize their models with the atmosphere by incorporating current atmospheric data in them). However, lithospheric dynamics is largely unobservable: important parameters cannot (or can rarely) be measured in Nature. Earthquakes, though, provide indirect but measurable clues of the stress and strain status in the lithosphere, which should be helpful for the synchronization of the models.The rupture area is one of the measurable parameters of earthquakes. Here we explore how it can be used to at least synchronize fault models between themselves and forecast synthetic earthquakes. Our purpose here is to forecast synthetic earthquakes in a simple but stochastic (random) fault model. By imposing the rupture area of the synthetic earthquakes of this model on other models, the latter become partially synchronized with the first one. We use these partially synchronized models to successfully forecast most of the largest earthquakes generated by the first model. This forecasting strategy outperforms others that only take into account the earthquake series. Our results suggest that probably a good way to synchronize more detailed models with real faults is to force them to reproduce the sequence of previous earthquake ruptures on the faults. This hypothesis could be tested in the future with more detailed models and actual seismic data.     

The vulnerability of karst springs—a case study of the Hubelj spring (SW Slovenia)

The hydraulic behaviour of the karst aquifer in the Hubelj spring catchment area (SW Slovenia) was studied by using an indirect research method based on natural tracers. The variations of natural tracers (in precipitation and in groundwater) during the storm event made possible the separation of the Hubelj spring storm hydrograph by the three-component separation technique. The results produced information on the aquifer recharge, storage and discharge processes, as well as on the mechanisms that affected them. They verified the so-called epikarst hypothesis presuming that an important part of a karst aquifer recharge reaches rapidly and intensively from the epikarst zone. It was demonstrated that epikarst water could occupy up to 50% of the spring discharge during precipitation events. This phenomenon could have important consequences on protection and management of the problems of karst aquifers, including engineering problems in karst areas. With this respect the results could give way to new engineering ideas.