Numerical modelling of desiccation cracking

The ability to model and predict the formation of desiccation cracks is potentially beneficial in many applications such as clay liner design, earth dam construction, and crop science, etc. However, most studies have focused on statistical analysis of crack patterns and qualitative study of contributing factors to crack development rather than prediction. Because it is exceedingly difficult to capture the nonlinear processes during desiccation in analytical modelling, most such models handle crack formation without considering variation of material properties with time, and are unattractive to use in realistic modelling. The data obtained from laboratory experiments on clay soil desiccating in moulds were used as a basis to develop a more refined model of desiccation cracking. In this study, the properties, such as matric suction, stiffness and tensile strength of soil, and base adhesion, could be expressed approximately as functions of moisture content. The initial conditions and the development of suction due to desiccation and the varying material properties were inputted to UDEC, a distinct element code, using its internal programming language FISH. The model was able to capture some essential physical aspects of crack evolution in soil contained in moulds with varying lengths, heights, and materials of construction. Extension of this methodology is potentially beneficial not only for modelling desiccation cracking in clay, but also in other systems with evolving material properties such as concrete structures and road pavements. Copyright © 2010 John Wiley & Sons, Ltd.     

Reaction softening by dissolution–precipitation creep in a retrograde greenschist facies ductile shear zone,New Hampshire,USA

We describe strain localization by a mixed process of reaction and microstructural softening in a lower greenschist facies ductile fault zone that transposes and replaces middle to upper amphibolite facies fabrics and mineral assemblages in the host schist of the Littleton Formation near Claremont, New Hampshire. Here, Na‐poor muscovite and chlorite progressively replace first staurolite, then garnet, and finally biotite porphyroblasts as the core of the fault zone is approached. Across the transect, higher grade fabric‐forming Na‐rich muscovite is also progressively replaced by fabric‐forming Na‐poor muscovite. The mineralogy of the new phyllonitic fault‐rock produced is dominated by Na‐poor muscovite and chlorite together with late albite porphyroblasts. The replacement of the amphibolite facies porphyroblasts by muscovite and chlorite is pseudomorphic in some samples and shows that the chemical metastability of the porphyroblasts is sufficient to drive replacement. In contrast, element mapping shows that fabric‐forming Na‐rich muscovite is selectively replaced at high‐strain microstructural sites, indicating that strain energy played an important role in activating the dissolution of the compositionally metastable muscovite. The replacement of strong, high‐grade porphyroblasts by weaker Na‐poor muscovite and chlorite constitutes reaction softening. The crystallization of parallel and contiguous mica in the retrograde foliation at the expense of the earlier and locally crenulated Na‐rich muscovite‐defined foliation destroys not only the metastable high‐grade mineralogy, but also its stronger geometry. This process constitutes both reaction and microstructural softening. The deformation mechanism here was thus one of dissolution–precipitation creep, activated at considerably lower stresses than might be predicted in quartzofeldspathic rocks at the same lower greenschist facies conditions.     

Unlocking the applications of automated mineral analysis

Identifying and quantifying the relative abundance of minerals is a fundamental part of many aspects of both pure and applied geology. Historically, quantitative mineralogy could be achieved using optical microscopy and point counting. This is a slow and operator dependent process, and practically impossible to achieve in, for example, very fine grained samples. Over the last decade a range of automated mineralogy technologies have arisen from the global mining industry and are being increasingly used in other branches of geology. These technologies, based on scanning electron microscopy with linked energy dispersive spectrometers, have the potential to revolutionise how we quantify mineralogy. In addition, during measurement, the sample textures are also captured, providing a wealth of valuable data for the geologist. In this article we review the current state of automated mineralogy and highlight the many areas of application for this technology.     

Centennial scale benthic foraminiferal record of late Holocene oceanographic variability in Disko Bugt,West Greenland

A new centennial scale benthic foraminiferal record of late Holocene climate variability and oceanographic changes off West Greenland (Disko Bugt) highlights substantial subsurface water mass changes (e.g. temperature and salinity) of the West Greenland Current (WGC) over the past 3.6 ka BP. Benthic foraminifera reveal a long-term late Holocene cooling trend, which may be attributed to increased advection of cold, low-salinity water masses derived from the East Greenland Current (EGC). Cooling becomes most pronounced from c. 1.7 ka BP onwards. At this point the calcareous Atlantic benthic foraminiferal fauna decrease significantly and is replaced by an agglutinated Arctic fauna. Superimposed on this cooling trend, centennial scale variability in the WGC reveals a marked cold phase at c. 2.5 ka BP, which may correspond to the 2.7 ka BP cooling-event recorded in marine and terrestrial archives elsewhere in the North Atlantic region. A warm phase recognized at c. 1.8 ka BP is likely to correspond to the ‘Roman Warm Period’ and represents the warmest bottom water conditions. During the time period of the ‘Medieval Climate Anomaly’ we observe only a slight warming of the WGC. A progressively more dominant cold water contribution from the EGC on the WGC is documented by the prominent rise in abundance of agglutinated Arctic water species from 0.9 ka BP onwards. This cooling event culminates at c. 0.3 ka BP and represents the coldest episode of the ‘Little Ice Age’.Gradually increased influence of cold, low-salinity water masses derived from the EGC may be linked to enhanced advection of Polar and Arctic water by the EGC. These changes are possibly associated with a reported shift in the large-scale North Atlantic Oscillation atmospheric circulation pattern towards a more frequent negative North Atlantic Oscillation mode during the late Holocene.     

High-resolution fluvial records of Holocene environmental changes in the Sahel: the Yamé River at Ounjougou (Mali,West Africa)

The Yamé river, in the Bandiagara Plateau, Dogon Country, Mali, is characterised by extensive alluvial sedimentary records, particularly in the 1 km long Ounjougou reach where Holocene floodplain pockets are inset in the Pleistocene formations. These alluvial records have been investigated via geomorphologic fieldwork and sedimentologic and micromorphologic analyses and are supported by 79 radiocarbon dates. The alluvial deposits of the valley floor correspond to a vertical accretion of 3–10 m. The reconstruction of fluvial style changes provides evidence of four main aggradation periods. From 11,500 to 8760 cal. BP, the alluvial architecture and grain-size parameters indicate a wandering river. This period included phases of pulsed high-energy floods and avulsion related to a northward shift of the summer monsoon to around 14°N after 11,500 cal. BP. From 7800 to 5300 cal. BP, a swampy floodplain environment with standing water pools within a Sudanian savanna/woodland mosaic corresponds to the culmination of the Holocene humid period. From 3800 cal. BP onwards, rhythmic sedimentation attests to an increase in the duration and/or intensity of the dry season, giving a precise date for the local termination of the Holocene Optimum period. During the last two millennia and for the first time during the Holocene, the alluvial formations are progressively restricted whereas the colluvial deposits increase, indicating strong soil erosion and redeposition within the watershed related to an increase in human impact. Four major periods are characterised by incision (I1: ante 11,500, I2: 8760–7800; I3: 6790–6500 cal. BP; I4; 2400–1700 cal. BP) pointing to dramatic changes in fluvial style. They result from high-energy flood flows during dry spells and confirm the capacity of the floodplain pocket in the upstream reach of the Sahelian belt to record rapid Holocene climatic change.     

Steady‐state analytical solutions of flow and deformation coupling due to a point sink within a finite fluid‐saturated poroelastic layer

An exact steady‐state closed‐form solution is presented for coupled flow and deformation of an axisymmetric isotropic homogeneous fluid‐saturated poroelastic layer with a finite radius due to a point sink. The hydromechanical behavior of the poroelastic layer is governed by Biot's consolidation theory. Boundary conditions on the lateral surface are specifically chosen to match the appropriate finite Hankel transforms and simplify the transforms of the governing equations. Ordinary differential equations in the transformed domain are solved, and then the analytical solutions in the physical space for the pore pressure and the displacements are finally obtained by using finite Hankel inversions. The analytical solutions at some special locations such as the top and bottom surfaces, lateral surface, and the symmetrical axis are given and analyzed. And a case study for the consolidation of a water‐saturated soft clay layer due to pumping is conducted. The analytical solution is verified against the finite element solution. Meanwhile, an analysis of coupled hydromechanical behavior is carried out herein. The presented analytical solution is an exact solution to the practical poroelastic problem within an axisymmetric finite layer. It can provide us a better understanding of the poroelastic behavior of the finite layer due to fluid extraction. Besides, it can be applied to calibrate numerical schemes of axisymmetric poroelasticity within finite domains. Copyright © 2017 John Wiley & Sons, Ltd.     

Seasonality of particle-associated trace element fluxes in the deep northeast Atlantic Ocean

In the framework of the German contribution to the Joint Global Ocean Flux Study (JGOFS), deep-water fluxes of particle-associated trace elements were measured in the northeast Atlantic Ocean. The sinking particles were collected almost continuously from 1992 to 1996 at three time-series stations, L1 (33°N/22°W), L2 (47°N/20°W), and L3 (54°N/21°W), using sediment traps. The focus of the present study is the temporal variability of the particle-associated elemental fluxes of Al, Ca, Cd, Co, Cu, Fe, Mn, Ni, P, Pb, Ti, V, and Zn at a depth of 2000 m.A clear seasonality of the fluxes that persisted for several years was documented for the southernmost station (L1) at stable oligotrophic conditions in the area of the North Atlantic Subtropical Gyre East (NASTE). At L2 and L3, an episodic nature of the elemental fluxes was determined. Mesoscale eddies are known to frequently cause temporal and spatial variability in the flux of biogenic components in that area. These events modified the simple seasonal pattern controlled by the annual cycle at L2, in the North Atlantic Drift Region (NADR), and at L3, which was influenced by the Atlantic Arctic province (ARCT). All stations were characterized by an additional episodic lithogenic atmospheric supply reaching the deep sea.The integrated annual fluxes during the multi-year study revealed similar flux magnitudes for lithogenic elements (Al, Co, Fe, Ti, and V) at L2 and L3 and roughly twofold fluxes at L1. Biogenic elements (Cd, P, and Zn) showed the opposite trend, i.e., two to fourfold higher values at L2 and L3 than at L1. For Mn, Ni, and Cu, the spatial differences were smaller, perhaps because of the intermediate behavior, between lithogenic and biogenic, of these elements. Similarly, among the three study sites, there were no noticeable differences in the total annual flux of Pb.The respective lithogenic fractions of the deep-sea fluxes of Cd, Co, Cu, Mn, Ni, V, and Zn were subtracted based on the amount of Al, with the average composition of the continental crust as reference. This procedure allowed estimation of the labile trace element fraction (TE_exc) of the particles, i.e., TE taken up or scavenged during particle production and sedimentation. The ratios of TE_exc/P clearly demonstrated an enrichment of TE over labile P from biogenic surface material to the deep sea for Zn (factor 4–6), Mn (12–27), Ni (3–5), and Cu (9–25); an intermediate status for Co (0.5–2.2); and depletion for Cd vs. P (0.2–0.4). Surprisingly, the recycling behavior of excess Co was found to be similar to that of P. Hence, Co_exc behaved like a biogenic element; this is in contrast to total Co, which is dominated by the refractory lithogenic fraction.Moreover, it is argued that these excess elemental fluxes caused a loss of the dissolved elements in upper waters, since their transport reaches the deep-sea waters at 2000 m, a depth far below of deep-winter mixing and upwelling. The annual amount of excess TE exported from surface waters was estimated to be 1.3×10⁹ mol Zn y^?1, 4.4×10⁹ mol Mn y^?1, 4.9×10⁸ mol Ni y^?1, 2.2×10⁷ mol Cd y^?1, 7.4×10⁸ mol Cu y^?1, and 2.7×10⁷ mol Co y^?1 for the whole North Atlantic Ocean. Important primary sources that could replenish these losses are the aeolian and fluvial supply processes.