Upstream perturbation and floodplain formation effects on chute-cutoff-dominated meandering river pattern and dynamics

A sustained dynamic inflow perturbation and bar–floodplain conversion are considered crucial to dynamic meandering. Past experiments, one-dimensional modelling and linear theory have demonstrated that the initiation and persistence of dynamic meandering require a periodic transverse motion of the inflow. However, it remains unknown whether the period of the inflow perturbation affects self-formed meander dynamics. Here, we numerically study the effect of the inflow perturbation period on the development and meander dynamics of a chute-cutoff-dominated river, which requires two-dimensional modelling with vegetation forming floodplain on bars. We extended the morphodynamic model Nays2D with growth and mortality rules of vegetation to allow for meandering. We tested the effect of a transversely migrating inflow boundary by varying the perturbation period between runs over an order of magnitude around typical modelled meander periods. Following the cutoff cascade after initial meander formation from a straight channel, all runs with sufficient vegetation show series of growing meanders terminated by chute cutoffs. This generates an intricate channel belt topography with point bar complexes truncated by chutes, oxbow lakes, and scroll-bar-related vegetation age patterns. The sinuosity, braiding index and meander period, which emerge from the inherent biomorphological feedback loops, are unrelated to the inflow perturbation period, although the spin-up to dynamic equilibrium takes a longer time and distance for weak and absent inflow perturbations. This explains why, in previous experimental studies, dynamic meandering was only accomplished with a sustained upstream perturbation in flumes that were short relative to the meander wavelength. Our modelling of self-formed meander patterns is evidence that scroll-bar-dominated and chute-cutoff-dominated meanders develop from downstream convecting instabilities. This insight extends to many more fluvial, estuarine and coastal systems in morphological models and experiments, which require sustained dynamic perturbations to form complex patterns and develop natural dynamics. © 2019 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.     

Carbon isotope evidence for CAM photosynthesis in the Mesozoic

Crassulacean acid metabolism (CAM) photosynthesis, a flexible adaptation to water‐stress, has not been shown to occur before the Pleistocene, despite suggestions that CAM first evolved in the Mesozoic or earlier. Here we report on multiple (≥7) positive organic carbon isotope excursions (δ¹³C_org > ?20 ‰ ) in Late Triassic and Early Jurassic (c. 200 Ma) sedimentary rocks from an arid terrestrial ecosystem in South Africa. These excursions are interpreted as evidence of the episodic dominance of CAM plants in response to heightened aridity in the heartland of Gondwana, and hint at climatic instability in a greenhouse world.     

Towards a conjoint-based, context-dependent model of task allocation in activity settings: Some numerical experiments

 None of the currently developed activity-based models of transport demand explicitly models task allocation among household members. To fill this gap, the present paper suggests to complement activity-based models of activity scheduling with a context-dependent model of task allocation. That is, it is assumed that the allocation of tasks within households is partly based on such contextual variables as the amount of time a member has to spend on mandatory activities and car availability. In particular, the paper advocates a conjoint-based approach, based on an assignment task as opposed to the traditional ranking, rating or choice response formats. By definition, an assignment task involves a combinatorial explosion of choice alternatives, implying that additional operational decisions to estimate the context-dependent model are required. This study presents the results of various numerical experiments, conducted to better understand the impacts of those decisions on the degree of bias in the parameter estimates of the choice model. The results of these simulations indicate that under particular assumptions, the development and estimation of a conjoint-based, context-dependent model of task allocation within households is feasible. Received: March 1999 / Accepted: June 2001     

Application of a mixing-ratios based formulation to model mixing-driven dissolution experiments

We address the question of how one can combine theoretical and numerical modeling approaches with limited measurements from laboratory flow cell experiments to realistically quantify salient features of complex mixing-driven multicomponent reactive transport problems in porous media. Flow cells are commonly used to examine processes affecting reactive transport through porous media, under controlled conditions. An advantage of flow cells is their suitability for relatively fast and reliable experiments, although measuring spatial distributions of a state variable within the cell is often difficult. In general, fluid is sampled only at the flow cell outlet, and concentration measurements are usually interpreted in terms of integrated reaction rates. In reactive transport problems, however, the spatial distribution of the reaction rates within the cell might be more important than the bulk integrated value. Recent advances in theoretical and numerical modeling of complex reactive transport problems [De Simoni M, Carrera J, Sanchez-Vila X, Guadagnini A. A procedure for the solution of multicomponent reactive transport problems. Water Resour Res 2005;41:W11410. doi: 10.1029/2005WR004056, De Simoni M, Sanchez-Vila X, Carrera J, Saaltink MW. A mixing ratios-based formulation for multicomponent reactive transport. Water Resour Res 2007;43:W07419. doi: 10.1029/2006WR005256] result in a methodology conducive to a simple exact expression for the space–time distribution of reaction rates in the presence of homogeneous or heterogeneous reactions in chemical equilibrium. The key points of the methodology are that a general reactive transport problem, involving a relatively high number of chemical species, can be formulated in terms of a set of decoupled partial differential equations, and the amount of reactants evolving into products depends on the rate at which solutions mix. The main objective of the current study is to show how this methodology can be used in conjunction with laboratory experiments to properly describe the key processes that occur in a complex, geochemically-active system under chemical equilibrium conditions. We model three CaCO₃ dissolution experiments reported in Singurindy et al. [Singurindy O, Berkowitz B, Lowell RP. Carbonate dissolution and precipitation in coastal environments: Laboratory analysis and theoretical consideration. Water Resour Res 2004;40:W04401. doi: 10.1029/2003WR002651, Singurindy O, Berkowitz B, Lowell RP. Correction to Carbonate dissolution and precipitation in coastal environments: laboratory analysis and theoretical consideration. Water Resour Res 2005;41:W11701. doi: 10.1029/2005WR004433], in which saltwater and freshwater were mixed in different proportions. The integrated reaction rate within the cell estimated from the experiments are modeled independently by means of (a) a state-of-the-art reactive transport code, and (b) the uncoupled methodology of [12, 13], both of which use dispersivity as a single, adjustable parameter. The good agreement between the results from both methodologies demonstrates the feasibility of using simple solutions to design and analyze laboratory experiments involving complex geochemical problems.     

Petrographic,fluid inclusion and isotopic study of the Dikulushi Cu–Ag deposit,Katanga (D.R.C.): implications for exploration

The Dikulushi Cu–Ag vein-type deposit is located on the Kundelungu Plateau, in the southeastern part of the Democratic Republic of Congo (D.R.C.). The Kundelungu Plateau is situated to the north of the Lufilian Arc that hosts the world-class stratiform Cu–Co deposits of the Central African Copperbelt. A combined petrographic, fluid inclusion and stable isotope study revealed that the mineralisation at Dikulushi developed during two spatially and temporally distinct mineralising episodes. An early Cu–Pb–Zn–Fe mineralisation took place during the Lufilian Orogeny in a zone of crosscutting EW- and NE-oriented faults and consists of a sequence of sulphides that precipitated from moderate-temperature, saline H₂O–NaCl–CaCl₂-rich fluids. These fluids interacted extensively with the country rocks. Sulphur was probably derived from thermochemical reduction of Neoproterozoic seawater sulphate. Undeformed, post-orogenic Cu–Ag mineralisation remobilised the upper part of the Cu–Pb–Zn–Fe mineralisation in an oxidising environment along reactivated and newly formed NE-oriented faults in the eastern part of the deposit. This mineralisation is dominated by massive Ag-rich chalcocite that precipitated from low-temperature H₂O–NaCl–KCl fluids, generated by mixing of moderate- and low-saline fluids. The same evolution in mineralisation assemblages and types of mineralising fluids is observed in three other Cu deposits on the Kundelungu Plateau. Therefore, the recognition of two distinct types of (vein-type) mineralisation in the study area has a profound impact on the exploration in the Kundelungu Plateau region. The identification of a Cu–Ag type mineralisation at the surface could imply the presence of a Cu–Pb–Zn–Fe mineralisation at depth.     

Climate entropy budget of the HadCM3 atmosphere–ocean general circulation model and of FAMOUS, its low-resolution version

The entropy budget is calculated of the coupled atmosphere–ocean general circulation model HadCM3. Estimates of the different entropy sources and sinks of the climate system are obtained directly from the diabatic heating terms, and an approximate estimate of the planetary entropy production is also provided. The rate of material entropy production of the climate system is found to be ～50 mW m^?2 K^?1, a value intermediate in the range 30–70 mW m^?2 K^?1 previously reported from different models. The largest part of this is due to sensible and latent heat transport (～38 mW m^?2 K^?1). Another 13 mW m^?2 K^?1 is due to dissipation of kinetic energy in the atmosphere by friction and Reynolds stresses. Numerical entropy production in the atmosphere dynamical core is found to be about 0.7 mW m^?2 K^?1. The material entropy production within the ocean due to turbulent mixing is ～1 mW m^?2 K^?1, a very small contribution to the material entropy production of the climate system. The rate of change of entropy of the model climate system is about 1 mW m^?2 K^?1 or less, which is comparable with the typical size of the fluctuations of the entropy sources due to interannual variability, and a more accurate closure of the budget than achieved by previous analyses. Results are similar for FAMOUS, which has a lower spatial resolution but similar formulation to HadCM3, while more substantial differences are found with respect to other models, suggesting that the formulation of the model has an important influence on the climate entropy budget. Since this is the first diagnosis of the entropy budget in a climate model of the type and complexity used for projection of twenty-first century climate change, it would be valuable if similar analyses were carried out for other such models.     

Sm–Nd isotopic mapping of lithospheric growth and stabilization in the eastern Kaapvaal craton

Whole-rock Sm–Nd isotope systematics of 79 Archean granitoids from the eastern Kaapvaal craton, southern Africa, are used to delineate lithospheric boundaries and to constrain the timescale of crustal growth, assembly and geochemical differentiation c. 3.66–2.70 Ga. Offsets in ε_Nd values for 3.2–3.3 Ga granitoids across the Barberton greenstone belt (BGB) are consistent with existing models for c. 3.23 Ga accretion of newly formed lithosphere north of the BGB onto pre-existing c. 3.66 Ga lithosphere south of the BGB along a doubly verging subduction margin. The Nd isotopic signature of c. 3.3–3.2 Ga magmatic rocks show that significant crustal growth occurred during subduction–accretion. After c. 3.2 Ga, however, the Nd signature of intrusive rocks c. 3.1 and 2.7 Ga is dominated by intracrustal recycling rather than by new additions from the mantle, signalling cratonic stability.     

Radiative transfer in disc galaxies – I. A comparison of four methods to solve the transfer equation in plane-parallel geometry

Accurate photometric and kinematic modelling of disc galaxies requires the inclusion of radiative transfer models. Because of the complexity of the radiative transfer equation (RTE), sophisticated techniques are required. Various techniques have been employed for the attenuation in disc galaxies, but a quantitative comparison of them is difficult, because of the differing assumptions, approximations and accuracy requirements that are adopted in the literature. In this paper, we present an unbiased comparison of four methods to solve the RTE, in terms of accuracy, efficiency and flexibility. We apply them all to one problem that can serve as a first approximation of large portions of disc galaxies: a one-dimensional plane-parallel geometry, with both absorption and multiple scattering taken into account, with arbitrary vertical distributions of stars and dust and an arbitrary angular redistribution of the scattering. We find that the spherical harmonics method is by far the most efficient way to solve the RTE, whereas both Monte Carlo simulations and the iteration method, which are straightforward to extend to more complex geometries, have a cost that is about 170 times larger.     

Gliding and overthrust nappe tectonics in the Barberton Greenstone Belt

Regional recumbent folds, inverted stratigraphy, nappes and olistostromes are described from the southern part of the 3.3–3.5 Ga Barberton Greenstone Belt. Overthrusting of thin rigid silicified slabs with minimum dimensions of 25 km² and up to 500 m in thickness, occurred over minimum distances of 86 km. More ductile and coherent units were overfolded up to at least 2 km during their emplacement. The glide planes on which these nappes travelled were zones of high fluid pressures related to hydrothermal fluid circulation patterns, driven by heat sources from igneous intrusions. The upwelling areas of the geothermal convection cells were sites of mud-pools and hydrothermal vents which may mark the trailing edges (pull aparts) of the overthrust units. Progressive silica and carbonate precipitation due to decreasing temperatures, within the zones of fluid migration distant from the areas of high heat flow, probably acted as built-in braking systems below the travelling slabs. Active sedimentation and metasomatism during this tectonism indicates a protracted history for the evolution of the greenstone belt. The recognition of nappe and overthrust tectonics in the Barberton Belt, processes which may have been commonplace in Archaean terrains, necessitates a re-evaluation of the stratigraphy of this belt.