Groundwater in the wetlands of the Okavango Delta, Botswana, and its contribution to the structure and function of the ecosystem

The Okavango Delta of northern Botswana is a large (40,000 km²) alluvial fan located at the terminus of the Okavango River. The river discharges about 10 km³ of water onto the fan each year, augmented by about 6 km³ of rainfall, which sustains about 2500 km² of permanent wetland and up to 8000 km² of seasonal wetland. Interaction between this surface water and the groundwater strongly influences the structure and function of the wetland ecosystem. The climate is semi-arid, and only 2% of the water leaves as surface flow and probably very little as groundwater flow. The bulk of the water is lost to the atmosphere. The Okavango River also delivers about 170,000 tonnes of bedload sediment and about 360,000 tonnes of solutes to the Delta each year, most of which are deposited on the fan. Bedload is deposited in the proximal, permanent wetland, whilst much of the solute load is deposited in the seasonal wetland. Notwithstanding the high evapotranspirational loss, saline surface water is rare. Between 80 and 90% of the seasonal flood water infiltrates the ground, recharging the groundwater beneath the flood plains and the many islands on the flood plains. The remainder is lost by evaporation. This groundwater reservoir is transpired into the atmosphere by both aquatic vegetation on the flood plains and terrestrial vegetation on the islands, and the water table is steadily lowered following passage of the seasonal flood. Trees, which are almost exclusively confined to islands, are particularly important, as they lower the water table beneath islands relative to the surrounding wetlands. There is therefore a net flow of groundwater towards islands. Accumulation of dissolved salts in this groundwater leads to precipitation of solutes (mainly of silica and calcite) in the soils beneath island fringes and the islands grow by vertical expansion. Islands are thus an expression of the chemical sedimentation taking place on the fan. Sodium bicarbonate accumulates in the groundwater beneath island centres, and this impacts on the vegetation, leading ultimately to barren island interiors. Dense saline brine thus produced subsides under density-driven flow. This cycling of seasonal flood water through the groundwater reservoir thus plays a key role in creating and maintaining the biological and habitat diversity of the wetland, and inhibits the formation of saline surface water.     

Reassessing southern African silcrete geochemistry: implications for silcrete origin and sourcing of silcrete artefacts

A synthesis of the geochemistry of silcretes and their host sediments in the Kalahari Desert and Cape coastal zone, using isocon comparisons, shows that silcretes in the two regions are very different. Kalahari Desert silcretes outcrop along drainage-lines and within pans, and formed by groundwater silicification of near-surface Kalahari Group sands. Silicification was approximately isovolumetric. Few elements were lost; silicon (Si) and potassium (K) were gained as microquartz precipitated in the sediment porosity and glauconite formed in the sub-oxic groundwater conditions. The low titanium (Ti) content reflects the composition of the host sands. Additional elements in the Kalahari Desert silcretes were supplied in river water and derived from weathering of silicates in basement rocks. Evaporation under an arid climate produced high-pH groundwater that mobilized and precipitated Si; this process is still occurring. In the Cape coastal zone, pedogenic silcretes cap hills and plateaus, overlying deeply weathered argillaceous bedrock. Silicification resulted from intensive weathering that destroyed the bedrock silicates, almost completely removing most elements and causing a substantial volume decrease. Some of the silica released formed a microcrystalline quartz matrix, and most Ti precipitated as anatase, so the Cape silcretes contain relatively high Ti levels. The intense weathering that formed the Cape silcretes could have occurred in the Eocene, during and after the Palaeocene-Eocene Thermal Maximum, when more acidic rainfall and high temperatures resulted in intensified silicate weathering worldwide. This could have been responsible for widespread formation of pedogenic silcretes elsewhere in Africa and around the globe. Trace element sourcing of silcrete artefacts to particular outcrops has most potential in the Cape, where differences between separate bedrock areas are reflected in the silcrete composition. In the Kalahari Desert, gains of some elements can override compositional differences of the parent material, and sourcing should be based on elements that show the least change during silicification. © 2020 John Wiley & Sons, Ltd.     

Dust provenance and its role as a potential fertilizing agent for the Okavango Delta,Botswana

Dust plays a globally important role in supplying biologically essential elements to landscapes underlain by nutrient-poor substrates. Here we show that dust may play a significant role in sustaining productivity in the vast wetlands of the Okavango Delta in southern Africa, one of the world's richest biodiversity hotspots. Dust accumulates preferentially on tree-covered islands in the seasonal swamps of the Delta, creating pockets of fine-grained, nutrient-rich material within the semi-arid landscape of the Kalahari Desert. Strontium and neodymium isotopes reveal that this dust likely originates predominantly from the Makgadikgadi salt pans, located 300 km away, and contributes 10–80% of the fine-grained material present in Okavango island soils. Surface material sourced from the Makgadikgadi Pans contains relatively high amounts of bioavailable phosphorus and iron, potentially influencing Okavango Delta biological productivity. We propose that long-term ecosystem productivity and nutrient availability in the Okavango may be strongly mediated by regional dust inputs. Understanding the influence of dust deposition on nutrient loads and biogeochemical cycling is thus critical for predicting the response of the Okavango Delta to future changes in climate. We suggest that dust inputs may play a significant role in the supply of nutrients to other large, global wetland systems located in dryland environments. © 2020 John Wiley & Sons, Ltd     

Controls on the transition from meandering to straight channels in the wetlands of the Okavango Delta,Botswana

Previous studies of alluvial rivers have shown that channel patterns form a continuum controlled by interactions among factors such as gradient, discharge, sediment size, and bank strength. Data from channels in the permanent wetlands of the Okavango Delta add to these ?ndings by focusing on pattern transitions in channels with banks formed by sedges and grasses that are rooted in peat and underlain by unconsolidated sediment. Channels are well de?ned, and transport ?ne–medium sand as bedload between the vertical, vegetation‐lined banks. Water depths, velocities, grain sizes, and bankline vegetation do not vary signi?cantly or systematically downstream, but the permeable banks allow water to leak from the channels, contributing to an overall downstream decrease in discharge and width. In addition, as the Okavango River ?ows from the <12 km wide ‘Panhandle’ and splits into distributaries in the broader ‘Fan’, valley gradient steepens by c. 60 per cent. These downstream changes result in channel pattern adjustments. In the Panhandle, the Okavango River is a relatively wide (c. 30–100 m), actively meandering, sinuous channel (P > 2·0), but further downstream in the Fan, the narrower (<40 m) distributaries follow laterally stable, less sinuous (‘straight’) courses (P < 1·75). Some channel pattern discrimination diagrams based on simple indices of gradient, discharge, sediment size or stream power are inadequate for analysing the meandering–straight transition in the Okavango but Parker's (1976) approach, based on ratios of depth–width and slope–Froude number, accurately characterizes the transition. Our ?eld observations, combined with the results from previous experimental studies, suggest that in relatively wide channels (w/d > 10), thalweg meandering results in scour of the unconsolidated sediment at the bank base, leading to undermining and collapse of the vegetation, and to slow meander migration. However, as channels narrow downstream (w/d < 10) with discharge losses, proportionally increasing sidewall drag exerted by bankline vegetation suppresses thalweg meandering and bank scour, and channels follow stable, less sinuous courses. Copyright © 2004 John Wiley & Sons, Ltd.     

Forecasting the spatial extent of the annual flood in the Okavango delta, Botswana

The pristine Okavango Delta wetland of northern Botswana is potentially under threat due to water abstraction from its tributaries. We have developed a statistical model which makes it possible to predict the extent of wetland loss which will arise from water abstraction. The model also permits prediction of the maximum area of flooding, and its spatial distribution, three months in advance of the flood maximum. The model was calibrated using maximum areas of seasonal inundation extracted from satellite imagery covering the period 1985–2000, which were correlated with rainfall and total flood discharge. A technique was developed to translate the modelled flood area into a flood map. The methodology can predict maximum area of flooding and its distribution with better than 90% accuracy. An important, although relatively minor, source of error in the spatial distribution of the flood arises from a secular change in flood distribution in the distal Delta which has taken place over the last 15 years. Reconstruction of flooding history back to 1934 suggests that the Delta may be subject to a quasi 80 year climatic oscillation. If this oscillation continues, the extent of flooding will increase in the coming decades.     

The role of biota in the initiation and growth of islands on the floodplain of the Okavango alluvial fan,Botswana

A group of islands of varying size on the floodplain of the Okavango alluvial fan, were studied to establish the processes which lead to the initiation and growth of islands. It was found that islands are initiated by the mound-building activities of the termite Macrotermes michaelseni. These termites import fine grained materials to use as a mortar for the construction of epigeal mounds. Their activities create a topographic feature, raised above the level of seasonal flooding, and also change the physical properties and nutrient status of the mound soil. Shrubs and trees are able to colonize these mounds, which results in increased transpiration. As a result, precipitation of calcite and silica from the shallow ground water occurs preferentially beneath the mounds, resulting in vertical and especially lateral growth, causing island expansion. © 1998 John Wiley & Sons, Ltd.     

The biogeomorphological life cycle of poplars during the fluvial biogeomorphological succession: a special focus on Populus nigra L.

Riverine ecosystems are recurrently rejuvenated during destructive flood events and vegetation succession starts again. Poplars (i.e. species from Populus genera) respond to hydrogeomorphological constraints, but, in turn, also influence these processes. Thus, poplar development on bare mineral substrates is not exclusively a one‐way vegetative process. Reciprocal interactions and adjustments between poplar species and sediment dynamics during their life cycle lead to the emergence of biogeomorphological entities within the fluvial corridor, such as vegetated islands, benches and floodplains. Based on a review of geomorphological, biological and ecological literature, we have identified and described the co‐constructing processes between riparian poplars and their fluvial environment. We have explored the possibility that the modification of the hydrogeomorphological environment exerted, in particular, by the European black poplar (Populus nigra L.), increases its fitness and thus results in positive niche construction. We focus on the fundamental phases of dispersal, recruitment and establishment until sexual maturity of P. nigra by describing the hierarchy of interactions and the pattern of feedbacks between biotic and abiotic components. We explicitly relate the biological life cycle of P. nigra to the fluvial biogeomorphic succession model by referring to the ‘biogeomorphological life cycle’ of P. nigra. Finally, we propose new research perspectives based on this theoretical framework. Copyright © 2014 John Wiley & Sons, Ltd.     

Fluvial chronology and valley floor evolution of the upper bowmont valley,borders region,Scotland

Two major fluvial terrace surfaces, and four less significant surfaces, are identified by aerial photographic interpretation, field mapping and levelling on one c. 2 km long reach of the upper Bowmont Water, an upland stream draining the northern Cheviot Hills, Scotland. The oldest terrace surface remains undated. Cartographic, radiocarbon and palynological dating suggest that later terraces formed very recently, within the last 250 years, with the most prominent terrace fill having aggraded in the 18th century. Incision below this terrace surface is dated to approximately the end of the 18th century, correlated with channel trenching at other sites in the region. Detailed palaeoecological and documentary evidence is used to examine whether climatic or land-use changes might have instigated aggradation, and a link with increased precipitation and flooding during part of the ‘Little Ice Age’ is suggested.     

Understanding mass fluvial erosion along a bank profile: using PEEP technology for quantifying retreat lengths and identifying event timing

This study provides fundamental examination of mass fluvial erosion along a stream bank by identifying event timing, quantifying retreat lengths, and providing ranges of incipient shear stress for hydraulically driven erosion. Mass fluvial erosion is defined here as the detachment of thin soil layers or conglomerates from the bank face under higher hydraulic shear stresses relative to surface fluvial erosion, or the entrainment of individual grains or aggregates under lower hydraulic shear stresses. We explore the relationship between the two regimes in a representative, US Midwestern stream with semi‐cohesive bank soils, namely Clear Creek, IA. Photo‐Electronic Erosion Pins (PEEPs) provide, for the first time, in situ measurements of mass fluvial erosion retreat lengths during a season. The PEEPs were installed at identical locations where surface fluvial erosion measurements exist for identifying the transition point between the two regimes. This transition is postulated to occur when the applied shear stress surpasses a second threshold, namely the critical shear stress for mass fluvial erosion. We hypothesize that the regimes are intricately related and surface fluvial erosion can facilitate mass fluvial erosion. Selective entrainment of unbound/exposed, mostly silt‐sized particles at low shear stresses over sand‐sized sediment can armor the bank surface, limiting the removal of the underlying soil. The armoring here is enhanced by cementation from the presence of optimal levels of sand and clay. Select studies show that fluvial erosion strength can increase several‐fold when appropriate amounts of sand and clay are mixed and cement together. Hence, soil layers or conglomerates are entrained with higher flows. The critical shear stress for mass fluvial erosion was found to be an order of magnitude higher than that of surface fluvial erosion, and proceeded with higher (approximately 2–4 times) erodibility. The results were well represented by a mechanistic detachment model that captures the two regimes. Copyright © 2017 John Wiley & Sons, Ltd.     

Where will widening occur in an outwash braidplain? A new approach to detecting controls on fluvial lateral erosion in a glacierized catchment (north-western Spitsbergen,Svalbard)

The dynamics of fluvial system evolution depend on fluvial processes and their driving forces associated with climatic variations, which affect changes in the morphology of river channels and floodplains. Neither channel slope and morphology, nor the properties of fluvial sediment have previously been considered as determinants of active braidplain widening on outwash plains (formed from valley/alpine glaciers and confined by pre-existing topography) in the High Arctic region and in the forefields of retreating glaciers. Factors determining widening of braidplain activity of the Waldemar River outwash (north-western Spitsbergen, Svalbard) were analysed on the basis of geomorphological, sedimentological, glaciological and meteorological research, and indicate significant multiple correlations between meltwater discharge, precipitation, braidplain width, morphology of the braided channel and the textural features of braidplain deposits. The capability of multivariate adaptive regression splines to detect these relationships was described, and threshold values were identified. The results indicate that the rate of active braidplain widening is proportional to the meltwater outflow in proglacial rivers, which can decrease despite a growing rate of glacial ablation. A proposed model enables us to predict zones of braidplain prone to widening activity in the High Arctic (humid) outwash fans and plains as well as fans developed in arid intermountain basins. The necessary conditions to activate braidplain widening processes were (1) spatial changes in the outwash feeding system due to glacier terminus retreat and (2) crossing the thresholds in passive factors (S and d₅₀) controlling lateral erosion intensity. As a result, the braidplain reached a new dynamic equilibrium, in which high-magnitude–low-frequency extreme meltwater discharges were of particular importance in terms of braidplain dynamics and are the dominant controls on the pattern of distributary channels.     

Salt transport on islands in the Okavango Delta: Numerical investigations

This study uses a numerical model to investigate the groundwater flow and salt transport mechanisms below islands in the Okavango Delta. Continuous evapotranspiration on the islands results in accumulation of solutes and the formation of a saline boundary layer, which may eventually become unstable. A novel Lagrangian method is employed in this study and compared to other numerical methods. The numerical results support the geophysical observations of density fingering on Thata Island. However, the process is slow and it takes some hundreds of years until density fingering is triggered. The results are sensitive to changes of the hydraulic gradient and the evapotranspiration rate. Small changes may lead to different plume developments. Results further demonstrate that density effects may be entirely overridden by lateral flow on islands embedded in a sufficiently high regional hydraulic gradient.     

Duricrust development and valley evolution: Process–landform links in the kalahari

Duricrusts are an important landscape component of the Kalahari region of central southern Africa. Their exposures within the dry valleys (mekgacha) of the Kalahari provide some of the most widespread surface outcrops of the terrestrial Jurassic to Holocene Kalahari Group sediments. Exposures have been extensively used in the construction of lithostrati-graphic sequences, on the assumption that valley systems have incised their courses through a pre-existing duricrust sequence. Recent work, however, has identified the role of groundwater erosion processes in valley development, which may have influenced duricrust formation. Studies of duricrusts from boreholes drilled within two mekgacha show that duricrust type is intrinsically related to the presence of a valley. Analyses of calcretes and silcretes in a series of profiles and thin sections from the Letlhakeng area of Botswana also indicate extensive alteration and diagenesis in association with former higher water tables. Sedimentary sequences within duricrust host materials can be identified but there is no evidence for correlation of duricrust cements between exposures. Profile studies from the Auob Valley in Namibia, however, suggest that this valley has incised through a sequence of duricrusts. Caution is advised in future attempts to correlate duricrust types on the basis of valley exposures, with the recommendation that where such exposures are used in a lithostratigraphic context, only duricrust host material characteristics and not cementing materials should be considered.     

The interpretation of visible band imagery of turbid shallow seas in terms of the distribution of suspended particulates

Visible band imagery from the Nimbus 7 Coastal Zone Colour Scanner (CZCS) and other satellite scanners indicates the occurrence of well-defined and persistent patterns of fine sediment concentration in the Irish Sea. The principal features of these distributions are related to variations in the intensity of tidal stirring with sharp transitions in the sub-surface reflectance occurring at frontal boundaries where there are abrupt changes in the vertical mixing regime. Observations of the suspended sediment concentration, both direct and by optical beam transmittance meter, have been used to demonstrate the correspondence between seston variations and changes in sub-surface reflectance as observed by satellite scanners.Radiometric measurements of sub-surface reflectance R are used to investigate the dependence of R (550 nm) on the concentration of suspended particulates Σ. The slope of the Σ- R relation is found to increase with the concentration of organic seston in a systematic way, a trend which is also apparent in an independent data set from the mouth of the Rhoˆne in the Mediterranean.These results may provide the basis for improved quantitative determination of seston from remotely sensed data if the mean level of organic seston concentration is known, or can be estimated from the spectral ratios in the reflectance data.     

The micro‐topography of the wetlands of the Okavango Delta,Botswana

The surface of the 40 000 km² Okavango alluvial fan is remarkably smooth, and almost everywhere lies within two to three metres of a perfectly smooth theoretical surface. Deviations from this perfect surface give rise to islands in the Okavango wetlands. This micro‐topography was mapped by assigning empirical elevations to remotely sensed vegetation community classes, based on the observation that vegetation is very sensitive to small, local differences in elevation. Even though empirical, the method produces fairly accurate results. The technique allows estimation of depths of inundation and therefore will be applicable even when high resolution radar altimetry becomes available. The micro‐topography has arisen as a result of clastic sedimentation in distributary channels, which produces local relief of less than two metres, and more importantly as a result of chemical precipitation in island soils, which produces similar local relief. The micro‐topography is, therefore, an expression of the non‐random sedimentation taking place on the fan. Volume calculations of islands extracted from the micro‐topography, combined with estimates of current sediment in?ux, suggest that the land surface of the wetland may only be a few tens of thousands of years old. Constant switching of water distribution, driven by local aggradation, has distributed sediment widely. Mass balance calculations suggest that over a period of c. 150 000 years all of the fan would at one time or other have been inundated, and thus subject to sedimentation. Coalescing of islands over time results in net aggradation of the fan surface. The amount of vertical aggradation on islands and in channels is restricted by the water depth. Restricted vertical relief, in turn, maximizes the distribution of water, limiting its average depth. Aggradation in the permanent swamps occurs predominantly by clastic sedimentation. Rates of aggradation here are very similar to those in the seasonal swamps, maintaining the overall gradient, possibly because of the operation of a feedback loop between the two. The limited amount of local aggradation arising from both clastic and chemical sedimentation, combined with constant changes in water distribution, has resulted in a near‐perfect conical surface over the fan. In addition to providing information on sedimentary processes, the micro‐topography has several useful hydrological applications. Copyright © 2004 John Wiley & Sons, Ltd.     

How do streamflow generation mechanisms affect watershed hypsometry?

Several studies have shown that the dominant streamflow generation mechanism in a river basin can leave distinct geomorphological signatures in basin topography. In particular, it has been suggested previously that basins generated by groundwater discharge tend to have a larger hypsometric integral than surface runoff basins because fluvial erosion is more focused in the valleys where groundwater discharge tends to occur. In this analysis, we aim to clarify this relationship by developing an alternative method to quantify the effects of streamflow generation mechanisms on basin hypsometry and by using a numerical model that can generate streamflow by different processes to evaluate the sensitivity of the results to the hydrological and geomorphological properties of the basin. The model results suggest that the hypsometric characteristics that are usually associated with groundwater discharge basins, such as a larger hypsometric integral, occur primarily when drainage networks are still advancing in the watershed. During later stages of development, an additional factor such as lithological controls or a distinct geomorphological process would be needed to preserve these features. The model results also show that the hypsometric effects are stronger when the parameters of the fluvial erosion process promote the influence of small discharge rates. Copyright © 2007 John Wiley & Sons, Ltd.     

Field application of close‐range digital photogrammetry (CRDP) for grain‐scale fluvial morphology studies

In situ measurement of grain‐scale fluvial morphology is important for studies on grain roughness, sediment transport and the interactions between animals and the geomorphology, topics relevant to many river practitioners. Close‐range digital photogrammetry (CRDP) and terrestrial laser scanning (TLS) are the two most common techniques to obtain high‐resolution digital elevation models (DEMs) from fluvial surfaces. However, field application of topography remote sensing at the grain scale is presently hindered mainly by the tedious workflow challenges that one needs to overcome to obtain high‐accuracy elevation data. A recommended approach for CRDP to collect high‐resolution and high‐accuracy DEMs has been developed for gravel‐bed flume studies. The present paper investigates the deployment of the laboratory technique on three exposed gravel bars in a natural river environment. In contrast to other approaches, having the calibration carried out in the laboratory removes the need for independently surveyed ground‐control targets, and makes for an efficient and effective data collection in the field. Optimization of the gravel‐bed imagery helps DEM collection, without being impacted by variable lighting conditions. The benefit of a light‐weight three‐dimensional printed gravel‐bed model for DEM quality assessment is shown, and confirms the reliability of grain roughness data measured with CRDP. Imagery and DEM analysis evidences sedimentological contrasts between gravel bars within the reach. The analysis of the surface elevations shows the effect variable grain‐size and sediment sorting have on the surface roughness. By plotting the two‐dimensional structure functions and surface slopes and aspects we identify different grain arrangements and surface structures. The calculation of the inclination index allows determining the surface‐forming flow direction(s). We show that progress in topography remote sensing is important to extend our knowledge on fluvial morphology processes at the grain scale, and how a technique customized for use by fluvial geomorphologists in the field benefits this progress. Copyright © 2016 John Wiley & Sons, Ltd.     

Dispersal of Mississippi and Atchafalaya sediment on the Texas-Louisiana shelf: Model estimates for the year 1993

A three-dimensional coupled hydrodynamic-sediment transport model for the Texas-Louisiana continental shelf was developed using the Regional Ocean Modeling System (ROMS) and used to represent fluvial sediment transport and deposition for the year 1993. The model included water and sediment discharge from the Mississippi River and Atchafalaya Bay, seabed resuspension, and suspended transport by currents. Input wave properties were provided by the Simulating WAves Nearshore (SWAN) model so that ROMS could estimate wave-driven bed stresses, critical to shallow-water sediment suspension. The model used temporally variable but spatially uniform winds, spatially variable seabed grain size distributions, and six sediment tracers from rivers and seabed.At the end of the year 1993, much of the modeled fluvial sediment accumulation was localized with deposition focused near sediment sources. Mississippi sediment remained within 20-40 km of the Mississippi Delta. Most Atchafalaya sediment remained landward of the 10-m isobath in the inner-most shelf south of Atchafalaya Bay. Atchafalaya sediment displayed an elongated westward dispersal pattern toward the Chenier Plain, reflecting the importance of wave resuspension and perennially westward depth-averaged currents in the shallow waters (<10 m). Due to relatively high settling velocities assumed for sediment from the Mississippi River as well as the shallowness of the shelf south of Atchafalaya Bay, most sediment traveled only a short distance before initial deposition. Little fluvial sediment could be transported into the vicinity of the “Dead Zone” (low-oxygen area) within a seasonal-annual timeframe. Near the Mississippi Delta and Atchafalaya Bay, alongshore sediment-transport fluxes always exceeded cross-shore fluxes. Estimated cumulative sediment fluxes next to Atchafalaya Bay were episodic and “stepwise-like” compared to the relatively gradual transport around the Mississippi Delta. During a large storm in March 1993, strong winds helped vertically mix the water column over the entire shelf (up to 100-m isobath), and wave shear stress dominated total bed stress. During fair-weather conditions in May 1993, however, the freshwater plumes spread onto a stratified water column, and combined wave-current shear stress only exceeded the threshold for suspending sediment in the inner-most part of the shelf.