The timing and isotopic character of regional hydrothermal alteration and associated epigenetic mineralization in the western sector of the Kaapvaal Craton (South Africa)

The Kaapvaal Craton of South Africa comprises an Archaean core of ≈3.5 Ga lithospheric and crustal rocks surrounded by younger accreted terrains of ≈3.0–2.7 and ≈2.1–1.9 Ga. The craton is covered by relatively undeformed 3.0–2.4 Ga supracrustal rocks, which show the effects of thermal and hydrothermal interaction. Part of this activity is manifested by a large number of epigenetic Pb–Zn (±Ag, Au, Cu, F) deposits in the cover rocks of the Kaapvaal Craton. These include small volcanic and breccia hosted deposits in mafic and felsic volcanic rocks of the 2.7 Ga Ventersdorp Supergroup and the Mississippi Valley-type (MVT) deposits in the carbonates of the Transvaal Supergroup.MVT mineralization at the Pering (and other Zn–Pb deposits) is hosted in fracture-generated N–S breccia bodies in the Paleoproterozoic carbonate succession of the western Kaapvaal Craton. The fluids carrying the metals were focused in vertical bodies within the fracture zones (FZ), the metals and the sulphur being carried together and precipitated in organic-rich sectors of the basin. Two small Pb–Zn deposits within mafic rocks of the Ventersdorp Supergroup, stratigraphically below the basin-hosted MVTs on the southwestern part of the Kaapvaal Craton have secondary chlorite which is extremely Rb-rich, associated with the mineralization. This chlorite and the associated altered basaltic host rocks give a Rb–Sr date of ≈1.98 Ga, and the associated galena Pb isotope data plot on the same array as those of other Pb–Zn deposits, the radiogenic intercept giving a date of ≈2.0 Ga. We interpret these data to indicate a craton-wide epigenetic fluid-infiltration event, which exploited the Maquassie Quartz Porphyry (MQP) as the aquifer and metal source.Sr isotopic results for the ore-zone gangue minerals show highly radiogenic ⁸⁷Sr/⁸⁶Sr ratios (>0.710) which support earlier models that the origin of radiogenic Sr isotopic composition in the calcite cements is the felsic tuffs (MQP) of the Ventersdorp Supergroup occurring at deeper levels within the basin. Relationships between δ¹⁸O and δ¹³C performed on carbonate cements within the aquifers are complex: the range in δ¹³C for some of the cements represents a mixture from two sources and with a progression from heavy carbon in the host to somewhat lighter carbon in the cements. Similarly, the lighter δ¹⁸O values have a narrow range indicative of rapid exchanges between hydrous fluid and rock.     

Modelling the influence of spatially varying hydrodynamics on the cross-sectional stability of double inlet systems

The cross-sectional stability of double inlet systems is investigated using an exploratory model that combines Escoffier’s stability concept for the evolution of the inlet’s cross-sectional area with a two-dimensional, depth-averaged (2DH) hydrodynamic model for tidal flow. The model geometry consists of four rectangular compartments, each with a uniform depth, associated with the ocean, tidal inlets and basin. The water motion, forced by an incoming Kelvin wave at the ocean’s open boundary and satisfying the linear shallow water equations on the f -plane with linearised bottom friction, is in each compartment written as a superposition of eigenmodes, i.e. Kelvin and Poincaré waves. A collocation method is employed to satisfy boundary and matching conditions. The analysis of resulting equilibrium configurations is done using flow diagrams. Model results show that internally generated spatial variations in the water motion are essential for the existence of stable equilibria with two inlets open. In the hydrodynamic model used in the paper, both radiation damping into the ocean and basin depth effects result in these necessary spatial variations. Coriolis effects trigger an asymmetry in the stable equilibrium cross-sectional areas of the inlets. Furthermore, square basin geometries generally correspond to significantly larger equilibrium values of the inlet cross-sections. These model outcomes result from a competition between a destabilising (caused by inlet bottom friction) and a stabilising mechanism (caused by spatially varying local pressure gradients over the inlets).     

Transverse structure of tidal and residual flow and sediment concentration in estuaries

An analytical and a numerical model are used to understand the response of velocity and sediment distributions over Gaussian-shaped estuarine cross-sections to changes in tidal forcing and water depth. The estuaries considered here are characterized by strong mixing and a relatively weak along-channel density gradient. It is also examined under what conditions the fast, two-dimensional analytical flow model yields results that agree with those obtained with the more complex three-dimensional numerical model. The analytical model reproduces and explains the main velocity and sediment characteristics in large parts of the parameter space considered (average tidal velocity amplitude, 0.1–1 m s^− 1 and maximum water depth, 10–60 m). Its skills are lower for along-channel residual flows if nonlinearities are moderate to high (strong tides in deep estuaries) and for transverse flows and residual sediment concentrations if the Ekman number is small (weak tides in deep estuaries). An important new aspect of the analytical model is the incorporation of tidal variations in the across-channel density gradient, causing a double circulation pattern in the transverse flow during slack tides. The gradient also leads to a new tidally rectified residual flow component via net advection of along-channel tidal momentum by the density-induced transverse tidal flow. The component features landward currents in the channel and seaward currents over the slopes and is particularly effective in deeper water. It acts jointly with components induced by horizontal density differences, Coriolis-induced tidal rectification and Stokes discharge, resulting in different along-channel residual flow regimes. The residual across-channel density gradient is crucial for the residual transverse circulation and for the residual sediment concentration. The clockwise density-induced circulation traps sediment in the fresher water over the left slope (looking up-estuary in the northern hemisphere). Model results are largely consistent with available field data of well-mixed estuaries.     

The use of scenarios as the basis for combined assessment of climate change mitigation and adaptation

Scenarios are used to explore the consequences of different adaptation and mitigation strategies under uncertainty. In this paper, two scenarios are used to explore developments with (1) no mitigation leading to an increase of global mean temperature of 4 °C by 2100 and (2) an ambitious mitigation strategy leading to 2 °C increase by 2100. For the second scenario, uncertainties in the climate system imply that a global mean temperature increase of 3 °C or more cannot be ruled out. Our analysis shows that, in many cases, adaptation and mitigation are not trade-offs but supplements. For example, the number of people exposed to increased water resource stress due to climate change can be substantially reduced in the mitigation scenario, but adaptation will still be required for the remaining large numbers of people exposed to increased stress. Another example is sea level rise, for which, from a global and purely monetary perspective, adaptation (up to 2100) seems more effective than mitigation. From the perspective of poorer and small island countries, however, stringent mitigation is necessary to keep risks at manageable levels. For agriculture, only a scenario based on a combination of adaptation and mitigation is able to avoid serious climate change impacts.     

Three-Dimensional Structure of Optimal Nonlinear Excitation for Decadal Variability of the Thermohaline Circulation

The decadal variability of the North Atlantic thermohaline circulation（THC） is investigated within a three-dimensional ocean circulation model using the conditional nonlinear optimal perturbation method. The results show that the optimal initial perturbations of temperature and salinity exciting the strongest decadal THC variations have similar structures： the perturbations are mainly in the northwestern basin at a depth ranging from 1500 to 3000 m. These temperature and salinity perturbations act as the optimal precursors for future modifications of the THC, highlighting the importance of observations in the northwestern basin to monitor the variations of temperature and salinity at depth. The decadal THC variation in the nonlinear model initialized by the optimal salinity perturbations is much stronger than that caused by the optimal temperature perturbations, indicating that salinity variations might play a relatively important role in exciting the decadal THC variability. Moreover, the decadal THC variations in the tangent linear and nonlinear models show remarkably different characteristics, suggesting the importance of nonlinear processes in the decadal variability of the THC.     

Late Variscan exhumation histories of the southern Rhenohercynian Zone and western Mid-German Crystalline Rise: results from thermal modeling

Thermal modeling techniques constrained by published petrological and thermo-chronometric data were applied to examine late orogenic burial and exhumation at a Variscan suture zone in Central Europe. The suture separates the southern Rhenohercynian zone from the Mid-German Crystalline Rise and traces the former site of a small oceanic basin. Closure of this basin during Variscan subduction and subsequent collision of continental units were responsible for different tectono-metamorphic evolutions in the suture's footwall and hanging wall. Relative convergence rates between the southern Rhenohercynian zone and western Mid-German Crystalline Rise can be inferred from the pressure-temperature-time evolution of the Northern Phyllite Zone. During Late Viséan-Early Namurian times, horizontal thrusting velocities were at least 20 mm/a. Thermal modeling suggests that exhumation of the Mid-German Crystalline Rise occurred temporarily at rates of more than 3 mm/a. Such rapid exhumation cannot be produced by erosion only, but requires a substantial contribution of extensional strain. Exhumation by upper crustal extension occurred contemporaneously with convergence and is explained by continuous underplating of crustal slices and thrusting along faults with ramp-flat geometry. Finally, implications for the tectono-metamorphic history of the study area and the thermal state of the crust during late Variscan exhumation are discussed.