On the variability of the flow along the Meso-American Barrier Reef system: a numerical model study of the influence of the Caribbean current and eddies

A high resolution (3–8 km grid), 3D numerical ocean model of the West Caribbean Sea (WCS) is used to investigate the variability and the forcing of flows near the Meso-American Barrier Reef System (MBRS) which runs along the coasts of Mexico, Belize, Guatemala and Honduras. Mesoscale variations in velocity and temperature along the reef were found in seasonal model simulations and in observations; these variations are associated with meandering of the Caribbean current (CC) and the propagation of Caribbean eddies. Diagnostic calculations and a simple assimilation technique are combined to infer the dynamically adjusted flow associated with particular eddies. The results demonstrate that when a cyclonic eddy (negative sea surface height anomaly (SSHA)) is found near the MBRS the CC shifts offshore, the cyclonic circulation in the Gulf of Honduras (GOH) intensifies, and a strong southward flow results along the reef. However, when an anticyclonic eddy (positive SSHA) is found near the reef, the CC moves onshore and the flow is predominantly westward across the reef. The model results help to explain how drifters are able to propagate in a direction opposite to the mean circulation when eddies cause a reversal of the coastal circulation. The effect of including the Meso-American Lagoon west of the Belize Reef in the model topography was also investigated, to show the importance of having accurate coastal topography in determining the variations of transports across the MBRS. The variations found in transports across the MBRS (on seasonal and mesoscale time scales) may have important consequences for biological activities along the reef such as spawning aggregations; better understanding the nature of these variations will help ongoing efforts in coral reef conservation and maintaining the health of the ecosystem in the region.     

Construction of detrital mineral populations: insights from mixing of U–Pb zircon ages in Himalayan rivers

Fission‐track, U–Pb and Pb–Pb analyses of detrital heavy mineral populations in depositional basins and modern river sediments are widely used to infer the exhumational history of mountain belts. However, relatively few studies address the underlying assumption that detrital mineral populations provide an accurate representation of their entire source region. Implicit in this assumption is the idea that all units have equal potential to contribute heavy minerals in proportion to their exposure area in the source region. In reality, the detrital mineral population may be biased by variable concentrations of minerals in bedrock and differential erosion rates within the source region. This study evaluates the relative importance of these two variables by using mixing of U–Pb zircon ages to trace zircon populations from source units, through the fluvial system, and into the foreland. The first part of the study focuses on the Marsyandi drainage in central Nepal, using tributaries that drain single formations to define the U–Pb age distributions of individual units and using trunk river samples to evaluate the relative contributions from each lithology. Observed mixing proportions are compared with proportions predicted by a simple model incorporating lithologic exposure area and zircon concentration. The relative erosion rates that account for the discrepancy between the observed and predicted mixing proportions are then modelled and compared with independent erosional proxies. The study also compares U–Pb age distributions from four adjacent drainages spanning ～250 km along the Himalayan front using the Kolmogorov–Smirnov statistic and statistical estimates of the proportion of zircon derived from each upstream lithology. Results show that, along this broad swath of rugged mountains, the U–Pb age distributions are remarkably similar, thereby allowing data from more localized sources to be extrapolated along strike.     

Rivers turned to rock: Late Quaternary alluvial induration influencing the behaviour and morphology of an anabranching river in the Australian monsoon tropics

Late Quaternary alluvial induration has greatly influenced contemporary channel morphology on the anabranching Gilbert River in the monsoon tropics of the Gulf of Carpentaria. The Gilbert, one of a number of rivers in this region, has contributed to an extensive system of coalescing low-gradient and partly indurated riverine plains. Extensive channel sands were deposited by enhanced flow conditions during marine oxygen isotope (OI) Stage 5. Subsequent flow declined, probably associated with increased aridity, however, enhanced runoff recurred again in OI Stages 4–3 (65–50 ka). Aridity then capped these plains with 4–7 m of mud. A widespread network of sandy distributary channels was incised into this muddy surface from sometime after the Last Glacial Maximum (LGM) to the mid Holocene during a fluvial episode more active than the present but less so than those of OI Stages 5 and 3. This network is still partly active but with channel avulsion and abandonment now occurring largely proximal to the main Gilbert flow path.A tropical climate and reactive catchment lithology have enhanced chemical weathering and lithification of alluvium along the river resulting in the formation of small rapids, waterfalls and inset gorges, features characteristic more of bedrock than alluvial systems. Thermoluminescence (TL) and comparative optically stimulated luminescence (OSL) ages of the sediments are presented along with U/Th ages of pedogenic calcrete and Fe/Mn oxyhydroxide/ oxide accumulations. They show that calcrete precipitated during the Late Quaternary at times similar to those that favoured ferricrete formation, possibly because of an alternating wet–dry climate. Intense chemical alteration of the alluvium leading to induration appears to have prevailed for much of the Late Quaternary but, probably due to exceptional dryness, not during the LGM. The result has been restricted channel migration and a reduced capacity for the channel to adjust and accommodate sudden changes in bedload. Consequent avulsions have caused local stream powers to increase by an order of magnitude, inducing knickpoint erosion, local incision and the sudden influx of additional bedload that has triggered further avulsions. The Gilbert River, while less energetic than its Pleistocene ancestors, is clearly an avulsive system, and emphasizes the importance in some tropical rivers of alluvial induration for reinforcing the banks, generating nickpoints, reworking sediment and thereby developing and maintaining an indurated and anabranching river style.     

The Okavango giant mafic dyke swarm (NE Botswana): its structural significance within the Karoo Large Igneous Province

The structural organization of a giant mafic dyke swarm, the Okavango complex, in the northern Karoo Large Igneous Province (LIP) of NE Botswana is detailed. This N110°E-oriented dyke swarm extends for 1500 km with a maximum width of 100 km through Archaean basement terranes and Permo-Jurassic sedimentary sequences. The cornerstone of the study is the quantitative analysis of N>170 (exposed) and N>420 (detected by ground magnetics) dykes evidenced on a ca. 80-km-long section lying in crystalline host-rocks, at high-angle to the densest zone of the swarm (Shashe area). Individual dykes are generally sub-vertical and parallel to the entire swarm. Statistical analysis of width data indicates anomalous dyke frequency (few data <5.0 m) and mean dyke thickness (high value of 17 m) with respect to values classically obtained from other giant swarms. Variations of mean dyke thicknesses from 17 (N110°E swarm) to 27 m (adjoining and coeval N70°E giant swarm) are assigned to the conditions hosting fracture networks dilated as either shear or pure extensional structures, respectively, in response to an inferred NNW–SSE extension. Both fracture patterns are regarded as inherited brittle basement fabrics associated with a previous (Proterozoic) dyking event. The Okavango N110°E dyke swarm is thus a polyphase intrusive system in which total dilation caused by Karoo dykes (estimated frequency of 87%) is 12.2% (6315 m of cumulative dyke width) throughout the 52-km-long projected Shashe section. Assuming that Karoo mafic dyke swarms in NE Botswana follow inherited Proterozoic fractures, as similarly applied for most of the nearly synchronous giant dyke complexes converging towards the Nuanetsi area, leads us to consider that the resulting triple junction-like dyke/fracture pattern is not a definitive proof for a deep mantle plume in the Karoo LIP.     

How a wet tropical rainforest copes with repeated volcanic destruction

The Holocene Period for the province of West New Britain, Papua New Guinea, is characterised by periodic catastrophic volcanism. The region is mantled in dense wet tropical rainforest, and has been occupied by people since the Pleistocene. Analyses of peat from two nearby sites within a lowland rainforest environment provide us with a macro-level landscape account of the periodic destruction and recovery of the coastal forests during seven periods of volcanic activity in the latter part (2900 yr ago to present) of the Holocene. Radiocarbon dating shows the very close correlation of the peat and tephra layers at both sites, yet the pollen analysis reveals different vegetation communities. These initial results allow us to begin identifying the processes of recovery, and to recognise different ecological pressures placed on vegetation at these neighbouring sites. Evidence of hydrological changes are observed beginning with a marine incursion recorded at Garu Site 3 1360 ¹⁴C yr B.P. The distinct differences in the vegetation re-establishment and community regeneration rates suggest the greater level of disturbance at Garu Site 1 could be related to the depth of the ashfall, although the proximity of a known human settlement may also be a contributing factor. Of note, palynologically, we found that the fern spore flora is particularly rich and believe it will be useful for ecological interpretation.     

Precessional forcing of tropical vegetation carbon storage

Since the Mid Pleistocene Revolution, which occurred about one million years ago, global temperatures have fluctuated with a quasi‐periodicity of ca. 100 ka. The pattern of past change in the extent of woodlands, and therefore by inference vegetation carbon storage, has been demonstrated to have a strong positive link with this global temperature change at high and mid latitudes. However, understanding of climate systems and ecosystem function indicates that the pattern of woodland change at low latitudes may follow a fundamentally different pattern. We present output from the intermediate complexity model GENIE‐1, comprising a single transient simulation over the last 800 ka and a 174‐member ensemble of 130 ka transient simulations over the last glacial cycle. These simulations suggest that while vegetation carbon storage in mid–high northern latitudes robustly follows the characteristic ca. 100 ka cycle, this signal is not a robust feature of tropical vegetation, which is subject to stronger direct forcing by the precessional (21 ka) orbital cycle (albeit with a highly uncertain response). We conclude that the correlation of palaeoenvironmental records from low latitudes with global temperature change must be done with caution. Copyright © 2011 John Wiley & Sons, Ltd.     

Integration of aquifer geology,groundwater flow and arsenic distribution in deltaic aquifers – A unifying concept

Groundwater arsenic (As) presents a public health risk of great magnitude in densely populated Asian delta regions, most acutely in the Bengal Basin (West Bengal, India and Bangladesh). Research has focused on the sources, mobilisation, and heterogeneity of groundwater As, but a consistent explanation of As distribution from local to basin scale remains elusive. We show for the Bengal Aquifer System that the numerous, discontinuous silt‐clay layers together with surface topography impose a hierarchical pattern of groundwater flow, which constrains As penetration into the aquifer and controls its redistribution towards discharge zones, where it is re‐sequestered to solid phases. This is particularly so for the discrete periods of As release to groundwater in the shallow subsurface associated with sea level high‐stand conditions of Quaternary inter‐glacial periods. We propose a hypothesis concerning groundwater flow ( S ilt‐clay layers I mpose H ierarchical groundwater flow patterns constraining A rsenic progression [SIHA]), which links consensus views on the As source and history of sedimentation in the basin to the variety of spatial and depth distributions of groundwater As reported in the literature. SIHA reconciles apparent inconsistencies between independent, in some cases contrasting, field observations. We infer that lithological and topographic controls on groundwater flow, inherent to SIHA, apply more generally to deltaic aquifers elsewhere. The analysis suggests that groundwater As may persist in the aquifers of Asian deltas over thousands of years, but in certain regions, particularly at deeper levels, As will not exceed low background concentrations unless groundwater flow systems are short‐circuited by excessive pumping.     

An analytical approach to ascertain saturation‐excess versus infiltration‐excess overland flow in urban and reference landscapes

Uncontrolled overland flow drives flooding, erosion, and contaminant transport, with the severity of these outcomes often amplified in urban areas. In pervious media such as urban soils, overland flow is initiated via either infiltration‐excess (where precipitation rate exceeds infiltration capacity) or saturation‐excess (when precipitation volume exceeds soil profile storage) mechanisms. These processes call for different management strategies, making it important for municipalities to discern between them. In this study, we derived a generalized one‐dimensional model that distinguishes between infiltration‐excess overland flow (IEOF) and saturation‐excess overland flow (SEOF) using Green–Ampt infiltration concepts. Next, we applied this model to estimate overland flow generation from pervious areas in 11 U.S. cities. We used rainfall forcing that represented low‐ and high‐intensity events and compared responses among measured urban versus predevelopment reference soil hydraulic properties. The derivation showed that the propensity for IEOF versus SEOF is related to the equivalence between two nondimensional ratios: (a) precipitation rate to depth‐weighted hydraulic conductivity and (b) depth of soil profile restrictive layer to soil capillary potential. Across all cities, reference soil profiles were associated with greater IEOF for the high‐intensity set of storms, and urbanized soil profiles tended towards production of SEOF during the lower intensity set of storms. Urban soils produced more cumulative overland flow as a fraction of cumulative precipitation than did reference soils, particularly under conditions associated with SEOF. These results will assist cities in identifying the type and extent of interventions needed to manage storm water produced from pervious areas.