ORIGINAL RESEARCH article

Front. Earth Sci., 25 May 2021

Sec. Environmental Informatics and Remote Sensing

Volume 9 - 2021 | https://doi.org/10.3389/feart.2021.669440

Blowout Morphometrics and Mass Balances

  • 1. Department of Earth Sciences, University College, London, United Kingdom

  • 2. Earth and Planetary Sciences, Birkbeck, University of London, London, United Kingdom

  • 3. Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou, China

  • 4. Department of Geography, University of Sheffield, Sheffield, United Kingdom

Abstract

The Gonghe Basin on the Qinghai-Tibet Plateau has a cold, arid climate and has suffered severe land degradation. Climate change as well as anthropogenic activities including overgrazing have resulted in widespread blowout development and the formation of some of Earth’s largest blowouts. The blowouts are part of an aeolian dominated landscape that passes from deflation zone to grass covered plain, and then through blowouts of increasing size and complexity to transverse barchanoid dunes that are migrating into the valley of the Yellow River. A combination of structure-from-motion (SfM) optical drone mapping, ground-penetrating radar (GPR) and soil pits are used to investigate blowout scour hollows and depositional lobes. Comparisons of the volumes of sediment removed from the scour hollows with the volumes of sediment deposited within adjacent lobes varies between sites. The lobe volume is invariably less than the volume of the scour hollow. This can, in part, be attributed to aeolian reworking of the lobe, distributing sand further downwind and uplifting of dust. However, much of the difference in volumes between the scour and lobe can be attributed to the measurement technique, particularly where GPR was employed to calculate lobe volumes. The wavelength of the GPR limits its ability to resolve thin layers of sand resulting in an underestimate of the deposited sand at the margins of a lobe where the sand thickness is equal to, or less than, the wavelength of the GPR. For thin sand layers, beneath the resolution of the GPR, soil pits suggest a closer match between the volume of sand eroded from the scour and the volume of the lobe, albeit with large measurement uncertainty. We put forth two hypotheses to explain the spatio-temporal evolution of the blowout dune field. The downwind increase in blowout dune size could either reflect a downwind propagation of aeolian instability; or it could result from an upwind propagation of the instability, which started at the highest points in the landscape and has subsequently migrated in a northwesterly direction, towards lower elevations. Recent optically stimulated luminescence dating appear to support the latter hypothesis.

1 Introduction

Blowouts are aeolian features consisting of an erosional depression and an associated downwind depositional lobe or apron. They play a critical role in diagnosing landscape changes, acting as an initial source of sand in reactivating dune fields as well as supplying sediment to downwind features (). Blowouts are found in coastal, semi-arid, sandy grassland, and desert landscapes (; ) and are found across the world, including Europe (van Boxel et al., 1997; ; ; Smyth et al., 2013; ), North America (; ; ), Africa (), and Asia (Sun et al., 2016; ; ). Blowout morphology varies, with scour hollows being classified as saucer, bowl, or trough-shaped and the eroded sediment deposited immediately downwind as either a large lobe or a thinly spread layer of sediment called an apron (Smyth et al., 2020).

Blowouts develop when a stabilised dune’s surface is disturbed, exposing the underlying sand to wind erosion. A positive feedback mechanism is then initiated with the reduction in surface roughness facilitating an increase in wind velocity, allowing for further aeolian deflation to occur (). Additional erosion and subsequent destabilisation of the surface, including any anchorage vegetation, then leads to further reductions in surface roughness, and an increasing rate of deflation. Over time, a small area of exposed sand may erode to form a large hollow, supplying sediment to the downwind environment (). Once initiated, the positive feedback mechanism is self-sustaining and can continue even if conditions are no longer conducive to blowout initialisation (). Subsequent morphologic development can be restricted by physical characteristics such as the size of the original stabilised dune, a layer of calcrete, or an armoured surface (). Seasonal changes, such as the height of the water table, surface moisture levels, or the magnitude and direction of prevailing winds, can also limit blowout development (; ; ). Blowout stabilization requires primary succession of pioneer flora to re-establish vegetation cover, mitigating blowout development by increasing surface roughness and decreasing bed shear stress (Schwarz et al., 2018).

The initial exposure of underlying sediment to wind erosion can occur through natural processes as well as anthropogenic activity. Natural factors which lead to blowout formation include coastal wave erosion, storm events, vegetation cover change, and increased wind speeds (; Smyth et al., 2013; ). Biotic factors associated with the grazing pressures from wild fauna, such as rabbits and bison, have also been found to influence blowout formation (van Boxel et al., 1997; ; ). Natural processes may occur seasonally, with reduced efficacy in months with more precipitation or a weaker wind regime (; ). In addition, anthropogenic activities can cause blowout initiation through direct and indirect actions. Off-road vehicle disturbance, pathway trampling, and poor land-use practices can remove anchorage vegetation facilitating blowout genesis (; ; ; ). Pastoral farming poses a dual threat of increased grazing pressures as well as high-intensity trampling involved with moving livestock between paddocks (; van der Hagen et al., 2020). Indirectly, infrastructure projects and urbanisation may also divert wind patterns, changing the spatiality of sand transport mechanisms and influencing the formation of deflation landforms (; Wernette et al., 2020). Abiotic and biotic factors rarely work independently (). As such, blowout initiation may be exacerbated by a number of processes at any given study location.

Measuring up to 40 m deep and 1,000 m long (), the blowouts in Gonghe Basin, Qinghai-Tibet Plateau (QTP) distinguish themselves from previously studied, smaller features in other parts of the world (). In comparison to their coastal counterparts, the mega-blowouts of the QTP defy conventional models of blowout development, exceeding a size in which coastal blowouts would tend towards stabilisation (). Yet the scarcity of detailed morphological investigations into blowouts on the QTP means that much remains unknown about these mega features. Of particular interest is the difference in volume between the depositional lobes and erosional scour hollows. An inequality in their mass balance would quantify the amount of sediment lost to the surrounding area. This is either as a fine layer extending beyond the depositional lobe or as atmospherically entrained dust.

In this paper we present morphometrics and a mass balance for a range of blowouts on the QTP. An innovative methodology which combines ground-penetrating radar (GPR) surveys, Structure-from-Motion (SfM) modeling from optical drone footage, and geospatial analysis is used to determine a mass balance between the erosional hollow and depositional lobe of three blowouts. A further two compound mega-blowouts are also assessed. Such quantifications are vital to address both the lack of detailed morphometrics for the full range of blowout features and the lack of parametrization in sediment reactivation processes. This research focuses on the Gonghe Basin on the QTP, but the concepts and methodology can be applied elsewhere.

2 Study Site

The Gonghe Basin on the QTP in China is an important region for agriculture, pastoral farming, and power production. The semi-arid steppe is comprised of grassland, farmland, and semi-shifting dunes (Zhang et al., 2003; Yu and Jia, 2014). In recent years, the cold alpine desert region of the Gonghe Basin has experienced severe land degradation from anthropogenic activities and local climatic warming (Wang et al., 2017). Aeolian sand reactivation, predominantly from wind erosion, has substantially eroded cultivated farmlands (; Zhang et al., 2003) and proceeds to threaten downwind infrastructure projects from the risk of sand encroachment. Of particular note, the Longyangxia Dam and Solar Power Park are located downwind from such eroded farmlands in Gonghe Basin. Attempts have been made to reduce soil degradation by enclosing grasslands with fencing; however, unprotected areas remain overgrazed and the basin continues to be cultivated (Zhang et al., 2003). The threat of further wind erosion is therefore still present. Despite this, the morphology, distributions, and impacts of wind erosion and blowout development in the Gonghe Basin remains a largely understudied phenomena ().

The mega-blowouts of the QTP are some of the largest on Earth () and are themselves part of a wider sand flow system from the Gonghe Basin to the Yellow River valley. The Gonghe Basin (Figure 1) covers an area of 13,800 km2 () and is a northeast to southwest aligned endorheic basin flanked by the Qinghai Nanshan mountains to the north and the Wahongshan Mountains to the south. There is a remarkable coincidence between the orientation of the valley, the mountains on either side, and the winds. Winds blowing along the axis of the basin have caused widespread erosion (Zhang et al., 2003; ). Sand deflated from the valley has accumulated at the downwind end of the valley forming a gently inclined plain that rises around 450 m over 50 km, a gradient of about 1:100, before descending very steeply into the valley of the Yellow River that is around 620 m deep. Across much of the plain the sand-covered slope is vegetated and described as an “alpine meadow” (Yu et al., 2017). The topography includes northeast to southwest elongated ridges of sand that are possibly vegetated linear dunes.

FIGURE 1

The Gonghe Basin experiences average annual temperatures of 2.4—4.1°C and an annual average precipitation of 314.3—414.8 mm (). Approximately 80% of this rainfall occurs between June and September (Yu and Jia, 2014; Yu et al., 2017). Vegetation cover during this wet season increases to a maximum of 20% but remains below 5% during the rest of the year (Wang et al., 2018). The cold alpine desert region in the northeastern QTP is an unusual aeolian environment due to its low temperature and low air pressure (Lu et al., 2005). Sitting at 3750 m above sea level and at the intersection of the westerlies wind belt, Asian, and Indian monsoon, the fragile ecosystem is extremely sensitive to global climate change (Qiang et al., 2016; Luo et al., 2019b). Strong northerly and northwest winds of up to 40 ms−1 occur from February to May (Qiang et al., 2013) with average annual speeds ranging from 2.1–2.7 ms−1 (Yan and Shi, 2004; Yu and Jia, 2014).

The Longyangxia Power Park is located within the Basin and is comprised of a hydroelectric dam and a solar power farm. It was the largest photovoltaic power station in the world until 2017 () and generates an average of 5.94 GW h hydroelectric power and 0.48 GW h photovoltaic power a year (). Aeolian sediment transport in the basin poses three main threats to the Longyangxia Power Park, including the burial of solar panels, obscuration of solar panels by dust, and the in-filling of the Longyangxia Dam Reservoir.

Five blowouts were surveyed in the Gonghe Basin (Table 1), located approximately 50 km south-west of Qabqa, Gonghe County, Qinghai Province (Figure 1). Blowouts A, B, and C are located within grazed pastures. Vehicular disturbance and animal activity in the area contribute heavily to blowout initiation. Access to the gated pastures intensifies trampling pressures at egresses and channels vehicles through corridors between paddocks. Additionally, small depressions in the landscape are used by cattle, either for shelter of dust bathing, causing similar fauna pressures on the landscape to those described in Canadian prairies by . Downwind, the size of the blowouts increases until the sand cover increases to form a field of transverse and barchanoid dunes that migrate into the valley of the Yellow River. The larger, more complex blowouts D and E are close to the interface between the blowouts and the active sand dunes and appear to be developed within vegetation-stabilised parabolic dunes. Blowouts B, D, and E were previously studied by , who looked at their development and growth from 1967 to 2015 with historical satellite imagery, and monitored the expansion of their erosional hollows with annual surveys from 2015 to 2018 as well as drone images (). Here we present updated morphometrics for these blowouts using finer-resolution DEM reconstructions and imaging of the internal structures within depositional lobes. The five blowouts exemplify the range of features found in the Gonghe Basin and are discussed in order of ascending complexity.

TABLE 1

BlowoutLatitude (DD)Longitude (DD)Hollow depth (m)Elevation (m)
A35.9512100.30345.123095
B35.9641100.27608.303101
C35.9514100.29854.813096
D35.8967100.28973329
E35.8600100.28013235

Location of the studied blowout features (centroid of the erosional hollow).

3 Methodology

This study employs an innovative combination of surface and subsurface imaging. SfM optical drone mapping was combined with GPR surveys to provide a high-resolution, non-invasive approach that leaves the fragile landscape undisturbed. We extend the technical work of to provide geomorphometry and geospatial analysis of the Gonghe Basin using a consistent suit of SfM data. The fieldwork for this study was conducted in July 2018.

Two Phantom four drones were used to survey the full extent of the blowouts. Images were first collected at nadir angles with additional images then captured inside the scour hollows to ensure lateral coverage. A Digital Elevation Model (DEM) and ortho-mosaic were created for each blowout feature using the established SfM approach within Agisoft Photoscan Professional Edition (Scarelli et al., 2016; ). Firstly, the optical drone images were aligned and a dense point cloud calculated. Outlying points within the dense clouds were removed manually and a mesh built from the cleaned dataset. This mesh was converted into a DEM and an orthorectified true colour composite texture of the optical drone images was created. Ground control points were recorded using a differential GPS, allowing the SfM data to be precisely georeferenced. Flight details and DEM reconstruction metrics are listed in the Supplementary Material (Supplementary Table S4). The three-dimensional models of the study areas were then used to determine the erosional hollow boundary and lobe extent for each blowout feature where it could be discerned. The elevation of the erosional hollow boundary was sampled at 0.2 m intervals and an original surface interpolated using ordinary kriging. This surface was differenced from the DEM to determine the volume of the erosional hollow.

GPR surveys have been successfully conducted to investigate the internal structure of aeolian features (; ), but relatively few studies have combined the technique with remote sensing to analyse blowout features (Neal and Roberts, 2001; ; ). The low conductivity of sand and large-scale structures contained within aeolian features lend themselves well to GPR surveys () which use the transmission of high-frequency electromagnetic fields to image subsurface structures. The GPR profiles within this paper were collected using 200 MHz antennas, a step size of 0.2 m, and an antenna separation of 1 m. CMP surveys indicated a velocity of 0.08 mns−1, that of damp sand. This corresponds well with field observations as short rainfall events occurred frequently throughout the duration of the fieldwork. Similar to the drone mapping, a differential GPS was used to precisely co-locate the GPR survey transects. Topographic corrections were reconstructed from the SfM-determined DEMs. Where a depositional lobe was evident, its depth was traced and extracted at the same resolution as the step size (0.2 m). Ordinary kriging was then used to interpolate a basal surface of the lobe, which was subtracted from the DEM to determine lobe volume.

The uncertainty of the volumetric calculations was determined by the GPR parameters. The signal had a wavelength of 0.4 m, and thus a resolution one quarter that of the wavelength, 0.1 m. The volumetric uncertainty (m3) was therefore calculated by multiplying the area of the lobe footprint (m2) by ±0.05 m.

4 Ground-Penetrating Radar Image Analysis

Interpretation of the GPR profiles follows the radar stratigraphy and radar facies methodology (; Neal and Roberts, 2001; ), where reflections on the GPR profile are interpreted as stratigraphic horizons and reflection terminations mark breaks in deposition or truncation due to erosion. The top of the soil horizon at the base of the lobe forms a continuous high-amplitude sub-horizontal reflection at Blowout A (Figures 2, 3).

FIGURE 2

FIGURE 3

At Blowout B the basal reflection has a very low-angle dip downwind parallel to the blowout lobe and a relatively steep slope to the west perpendicular to the lobe (Figures 4, 5), because the lobe runs along the flank of a vegetated linear dune. The GPR profile along the axis of the lobe at Blowout A shows inclined reflections that dip in the downwind direction (Figure 2). These are interpreted as former positions of the lobe surface and record the extension of the lobe downwind. Changes in the dip of the reflections suggest that the downwind margin of the lobe steepened and then flattened out. The dip of the reflections shows that the downwind margin did not reach the angle of repose. Truncation of the reflections within the lobe indicate periods of erosion, when the lobe was reshaped by changes in the wind. Reflections on the GPR profiles across the lobe are dominated by convex reflections that are interpreted to come from former positions of the lobe surface (Figure 3).

FIGURE 4

FIGURE 5

The convex morphology mirrors the lobe surface at the time of survey. However, the convex reflections are not symmetrical; some are truncated and others are only preserved on the east side of the lobe. Truncation of the reflections marks erosion, most notably at the upwind end of the lobe where a central scour has eroded into the upwind end of the lobe. Truncation of reflections within the lobe are interpreted as reactivation surfaces where the lobe was reshaped by a change in the wind. The asymmetric development of the lobe with local advance towards the west indicates a change in wind direction.

The presence of downlap and truncation surfaces within the lobes indicates reshaping of the lobe surface due to changes in the wind strength and direction similar to reactivation surfaces identified in sand dunes () as well as GPR profiles across sand dunes (e.g., ). The abundance of truncation and downlap terminations within the lobe strata indicates extensive reworking of the lobe, with episodes of erosion as well as periods of deposition.

5 Morphometric Analysis

5.1 Blowout A

Blowout A (Figure 6) exemplifies the standard blowout model with one active erosional hollow and a subsequent clearly defined downwind depositional lobe. The blowout, including the deflation hollow and depositional lobe, measures 170 m in length, 33 m in width, and has a total footprint of 4,015 m2. The hollow accounts for just under half of the total feature length (79 m). The blowout is oriented to the south-east yet deposited sediment is evident on the south-west edge of the hollow. This is likely caused by a slight deviation from a perfect elongated concave shape, allowing eroded material to be deposited at an angle when exiting the hollow. The maximum width of the hollow is similar to that of the lobe, reaching 33 m and 31 m, respectively.

FIGURE 6

A GPR survey was conducted to determine the volume of the depositional lobe. A single longitudinal transect bisected the depositional lobe from NW to SE (Figure 2). Eight transects were then made from SW to NE at 10 m intervals perpendicular to the longitudinal transect (Figure 3). The continuous lines at the top of the transects in Figures 2, 3 denote the direct signals between the transmitter and the receiver. The first signal is the airwave, the fastest direct signal that travels from the transmitter to the receiver directly through the air at the speed of light. The second signal represents the groundwave, a signal that travels directly to the receiver through the near-surface without reflecting off of any soil horizons at a lower velocity. These direct arrivals effectively obscure the top 0.5 m of strata, limiting the use of the GPR to detect the thickness of layers of sand where the thickness is less than 0.5 m. The NW-SE bisect shows inclined reflections dipping in the downwind direction as well as a reflection from the base of the lobe (Figure 2). The continuous sub-horizontal reflection at the base of the lobe marks the contact between the lobe and the underlying soil horizon and is coloured blue on Figure 2. Inclined reflections within the lobe increase in dip downwind until it reaches an angle of around 30 degrees and then decreases again (Figure 2) The downwind dipping reflections are evidence of progradation, with high-angle reflections indicating the angle of repose and development of a slipface on the downwind end of the lobe. Truncation of the inclined reflections at the upwind end of the lobe indicates erosion and is attributed to aeolian reworking of the lobe deposits. Bounding surfaces between sets of inclined reflections are attributed to aeolian reworking of the lobe, possibly during storm events. In conjunction with the sub-horizontal reflections in the perpendicular transects, the lobe shows that tabular cross-strata have formed at the leeward slope. The GPR profiles across the lobe show a central scour at the upwind end that is eroding through the lobe deposits and the underlying soil, truncating the basal blue reflection (Figure 2). Reflections within the cross-section profiles show convex-upwards reflections as well as low-angle inclined reflections (Figure 3). The convex-upward reflections are attributed to the dome-like morphology of the lobe recording aggradation on the lobe surface. This is most likely wind ripple laminae that are beneath the resolution of the GPR. The low-angle inclined reflections overlay bounding surfaces and are interpreted to indicate changes in the accretion direction of the lobe, with the lobe building out towards the east, most likely a response to changes in wind direction.

The DEM indicates that the scour hollow is 3 m deeper than the lobe is thick. However, the lobe extends over a larger area, 2,133 m2 compared to 1,882 m2 of the hollow. The erosional hollow has a volume of 5,806 m3 and the lobe has a volume of 2,230 110 m3. The volumes calculated for the depositional lobe appear to account for only 38% of the material eroded from the hollow, suggesting that a significant amount of sediment (3,570 m3) is lost to the surrounding area, either as a fine layer of sediment or entrained in the atmosphere as dust. It is likely that most of this volume can be attributed to downwind distribution of sand and dispersion into the surrounding grassland in layers that are too thin to be imaged using a 200 MHz GPR system.

5.2 Blowout B

Blowout B (Figure 7) differs from Blowout A (Figure 6) in that it is eroding into the flank of an elongated ridge, possibly a vegetated linear dune. The blowout lobe is not readily apparent in the DEM or in the orthorectified true colour composite. Instead, an apron of sediment is evident, measuring 633 m in length and 73 m in width at its maximum.

FIGURE 7

A partially broken barbed wire fence precluded the first 30 m of a bisecting GPR transect from being surveyed. Despite this, a 84 m NW-SE bisect was completed (Figure 4). Five SW-NE transects were made perpendicular to the bisect at approximately 20 m intervals (Figure 5). In contrast to Blowout A, the transects show a much more complex history. Low-angle reflections and truncated erosion surfaces indicate several periods of reconstruction and erosion. Reworking of the lobe sediment has produced clear scour surfaces and saddle morphology is evident. The increased internal complexity of the lobe could be indicative of maturity, or result from the slope exposing the lobe to a more complex wind regime.

When extracting the basal surface from the GPR transects, the resulting lobe depth is surprisingly symmetrical despite the topographic gradient. Analysis based solely on the DEM without knowledge of the subsurface would therefore overestimated the plume depth. Interpolating the basal surface covered by GPR surveys produces a volume estimate of 6,310 260 m3. To estimate the total lobe volume, the first 10 m of the baseline was extrapolated to cover the area underneath the barbed wire. This produced an estimated volume of 8,700 m3. The maximum volume, calculated by extending the base of the lobe to the identified hollow-lobe transition area, is 9,920 m3. The lobe therefore accounts for roughly half the volume of the hollow (19,900 1,900 m3). The remaining 50–56% of the eroded material may exist as a thin layer extending downwind which is too thin for the GPR to image. Alternatively, the “missing” sand may have escaped the depositional lobe by saltation towards and into the Longyangxia reservoir, or could have been released into the atmosphere as a fine dust.

5.3 Blowout C

The apron of Blowout C (Figure 8) is easy to discern in the orthorectified true colour composite, but is harder to delimit in the DEM due to a topographic incline. The apron is more than twice the maximum width of the hollow (58 m and 28 m, respectively), and three times that of its length (257 m and 87 m, respectively). Additionally, the hollow has a footprint only a tenth that of the apron (1,270 m2 and 13,000 m2, respectively). Five GPR surveys were conducted across the depositional plume, comprised of two 150 m longitudinal transects and three 50 m latitudinal transects. However, the apron was too thin to be resolved in the GPR surveys. Instead, 88 small pits were dug into the surface at random locations to measure the soil horizon depth (Supplementary Table S3). The depth of sand at each pit was measured in centimetres using an extending steel ruler. Pit locations were recorded using DGPS. These soil horizon depths underwent ordinary kriging to interpolate a basal surface (Figure 8D). The maximum depth of the horizon was measured at 0.46 m, and its volume was calculated at 2,290 m3, albeit with a substantial kriging interpolation uncertainty of 75%. This volume of the depositional lobe is only 980 m3 less than the erosional hollow’s volume (3,270 64 m3). The relatively small difference between the volume of the topsoil and the hollow may be indicative of lobe deflation with sand distributed into the surrounding grassland and lost as dust. However it would be premature to apply this conclusion to the other blowouts given the large uncertainty for the volume estimate of the soil horizon.

FIGURE 8

5.4 Blowout D

Blowout D shows a complex mega-blowout feature (Figure 9). The boundary between hollow and lobe is hard to discern with a shallow gradient existing throughout its 1,220 m length and 470 m span. The feature covers an area of 446,000 m2, an order of magnitude larger than Blowouts A–C. With this larger magnitude comes increasing complexity. Multiple secondary hollows are evident on the south-eastern edge of the feature. Some of these hollows are in the process of merging into each other while also simultaneously eroding the side of the primary hollow perpendicular to the direction of the prevailing wind. Juxtaposed to the simplistic model of Blowout A, which has one erosional hollow and one depositional lobe/apron, Blowout D is simultaneously eroding at multiple points along its boundary. Smaller erosional hollows appear on the edge of the larger primary hollow, feeding sediment into the larger system.

FIGURE 9

The sides of the erosional hollow are oversteepened (Figure 9E), with sediment along the scour sides approaching or sitting at the angle of repose. Resistant soils and plant roots flanking the scour sides are maintaining these slopes.

The topography within the hollow also differs from that of Blowouts A–C, with areas of sediment transport evident inside the erosional part of the feature. The mega-blowout appears to facilitate a wind regime that is reworking the sediment below the level of the hollow edge. Further downwind, ripple and proto-dune movement are evident. The hollow-lobe transition is elongated, with a gradual incline over which sand is transported out of the hollow. This complex system is not witnessed in the smaller-scale blowouts, where eroded sediment is immediately transported and deposited at the downwind margin of the hollow.

5.5 Blowout E

Blowout E (Figure 10) is another complex mega-blowout. A primary erosional hollow 180,110 m2 in size has deposited sediment over an incline to the east while simultaneously infilling historic and vegetated hollows to the south. This primary deposit covers an area of 238,000 m2. Multiple smaller hollows are evident further south of the primary deposit, and a secondary apron is present. Vegetation is evident both within the second largest hollow and around the edges of the secondary apron.

FIGURE 10

Two possible interpretations of the mega-blowout’s history of sediment reactivation may be drawn from our analysis. Firstly, the secondary erosional hollow could be mature and inactive, slowly becoming infilled by sediment produced in the active erosional hollow to the north-west. Secondly, the lobe could be in active transport through a larger hollow. The similar depths between the primary and secondary erosional hollows suggest the feature may be a single large hollow. A limit, such as a perched water table, may have been reached with sediment now passing through the feature in a similar way to Blowout D.

6 Discussion

suggest that the erosion in a blowout will equal the volume of sediment in the depositional apron due to the conservation of mass. However, we are not aware of any other studies of blowout deflation hollow/depositional lobe volumes with which to compare our results. As a consequence, we can only speculate that the differences observed in this study are not unique but should apply more widely. assume that the apron is 100% efficient in trapping wind-blown sand. By analogy with sand dunes it has been shown that dunes with slipfaces are not very efficient sand traps, trapping only a small percentage of the saltation flux, and dunes without slipfaces even less so (). The clear implication is that sand derived from deflation hollows can accumulate locally within depositional lobes, but equally it will be distributed more widely downwind. This combination of local sand accumulation and downwind reworking can account for most of the variations in morphology observed in the blowout lobes. The other significant influence is the wind direction, which can affect the direction of sediment transport and resulting lobe development.

Truncation of reflections on the GPR profiles provide evidence for erosion and indicates breaks in deposition, as well as reworking of the sand on the surface of the blowout lobe. Changes in the accretion direction of reflections indicates changes in the lobe morphology that is attributed to changes in wind direction. GPR also images the base of the lobes and this has been used to calculate the volume of sediment within two lobes. While the results of the GPR survey in imaging the lobe strata and the base of the lobe are positive, the calculation of lobe volumes using GPR was limited by the inability of the 200 MHz GPR system to image sand layers less than 0.5 m thickness, the section that is covered by arrival of the direct air and ground waves. Higher-frequency GPR systems e.g. 500 MHz antennas with a shorter wavelength would have helped in this respect. The use of soil pits to determine lobe thickness and volume shows a much closer agreement between the volume of the scour hollow and the volume of the depositional lobe at location C. Nonetheless, the kriging interpolation of these point measurements is associated with large statistical uncertainty (Table 2).

TABLE 2

BlowoutABC
Volume of hollow (m3)5806±10019900±19003270±64
Volume of depositional lobe (m3)2230±1108700–9920570–(4000)
Volume deficit (%)62±253±630 +52/−30

Volume deficit measurement of the three small blowouts.

It would be advantageous for future work to further constrain depositional volumes. There exists, however, two difficulties in this. Firstly, an invasive field campaign risks disturbing the fragile environment and contributing to further land degradation. Secondly, the attribution of widely spread fine material to a particular blowout becomes problematic further downwind. Delimiting the input of upwind features to the amount of sediment in the topsoil surrounding a blowout would require a comprehensive monitoring campaign and/or foreknowledge of sediment emissivity. We have shown that a large amount of material is lost to the surrounding downwind environment, and thus features further downwind are unlikely to be acting in a closed system. The proximity of other blowouts and their level of activity would therefore need to be carefully considered when constraining depositional volumes further. This issue notwithstanding, repeated drone surveys could be conducted to determine lobe accretion and scour deflation through time. While sandy environments pose challenges to photogrammetric methods, studies quantifying sand deposition on beaches using UAV surveys have produced DEMs with centimetre and sub-centimetre vertical accuracies (; , ). Similar applications would be useful to determine volumetric changes through time. However, these studies would be limited to quantifications of topographic changes as sub-surface imaging is still required to identify the lobe base.

6.1 Elevation and Scour Depth

The scour depth of blowouts in the coastal zone is limited by the presence of the water table and associated capillary zone that restricts sand entrainment, as well as the presence of coarse-grained beach deposits beneath coastal dunes. As a consequence the scour depth in coastal dunes is restricted because wind cannot deflate sand beneath the water table, which is close to mean sea level at the coast, nor can it entrain coarse-grained beach deposits that are commonly raised slightly above the high tide line in a backshore setting. In contrast, the blowouts of the QTP lack a shallow water table and the deflation depth is limited instead by the thickness of aeolian sand and depth to underlying strata. In the Gonghe Basin the sand thickness and blowout dimensions appear to increase in the downwind direction until they join the active sand dunes. We speculate that the increased sand thickness away from the deflation zone in the Gonghe Basin is a legacy of a downwind decrease in transport capacity as wind flow expanded away from the topographic confines of the basin. Drought conditions during the Little Ice Age are believed to have resulted in the reactivation of the blowouts (Luo et al., 2019b).

6.2 Upwind Propagation of Instability and Blowouts

The upwind decrease in blowout dimensions, away from the active dunes, and towards the vegetated plain suggests an upwind propagation of the blowout system, propagating upwind, away from the active dunes. This pattern can be viewed as a systematic downwind propagation of blowout morphology, from small blowouts (localities A, B, and C, this study), to larger mega-blowouts downwind (localities D and E, this study) and sites described by Luo et al. (2019b). This pattern suggests a downwind propagation of instability with blowout dimensions increasing downwind. However, an alternative hypothesis is that the sizes of the blowouts relate to their age. In which case the larger blowouts which are downwind started earlier than the smaller blowouts upwind. This explanation suggests the opposite, that the instability started at the downwind end and propagated upwind even though individual blowouts tend to propagate downwind. This hypothesis appears to be supported by OSL ages (; Luo et al., 2019b). One possible reason for what appears at first to be a counter-intuitive explanation for the blowout pattern is that the downwind areas which are at a higher elevation (Figure 1) were closer to the threshold for blowout development and it is only as anthropogenic impacts have increased in the grazed pastures at lower elevations that blowouts have developed upwind, suggesting a bio-geomorphic control on the wider pattern of climatically driven blowout development.

7 Conclusion

High resolution topographic surveys of blowouts including erosional hollows and adjacent depositional lobes were collected using SfM optical drone mapping and DGPS to create digital elevation models (DEMs) of the blowouts. In addition, ground-penetrating radar (GPR) profiles were collected along and across the depositional lobes at three of the smaller and less complex blowouts. Volumes of the erosional hollows and their adjacent depositional lobes are calculated from the DEMs, GPR data, and soil pits. The results show that although the areas of the hollows and lobes are similar, the volume of the lobes is much less than the volume of the adjacent hollows with the lobe volume amounting to around 50% of the hollow. Part of the difference in volumes can be attributed to the method used to calculate the lobe volumes. The wavelength of the GPR prevents measurement of sand layers that are thinner than the section covered by the first arrivals of the direct groundwaves and airwaves. Despite this, the results indicate that blowout lobes are not very efficient sand traps, capturing only 50% of the sand eroded from the hollow. The remaining 50% is probably distributed further downwind, as a thin sand layer captured by local vegetation, with an undetermined fraction transported further downwind as saltating sand or a “dust” plume during storm events.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author contributions

Drone footage collected by RB and TM with support from TB, WL, and PV. GIS field data collection from PV, TB, and WL. GPR surveys conducted by ED and CB, with additional support from TB. GPR interpretations by CB. Models created by TB and RB. Morphometrics and data analysis from TB. Manuscript written by TB, with substantial contributions from CB and edits from PV.

Funding

TB is supported by the Natural Environment Research Council (NERC) through Grant No. NE/L002485/1. This reseach was funded by Grant No. 41771015 of the Chinese National Natural Science Foundation, awarded to WL, and a Thesiger Oman fellowship (THES 02/18) of the Royal Geographical Society, awarded to PV.

Acknowledgments

We would like to thank Mark Bateman, Mei Shao, Songnian Luo, Wanli Luo, and Binhua Xu for field assistance, and the farmers of Gonghe for granting us access to their land. We would also like to thank LG-R and the anonymous reviewers for their critical and constructive comments.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.2021.669440/full#supplementary-material

References

Summary

Keywords

structure from motion, Qinghai-Tibet Plateau, ground-penetrating radar, aeolian geomorphology, geospatial reconstruction, blowouts

Citation

Baird TR, Bristow CS, Luo W, Du E, Bryant RG, Mitchell TM and Vermeesch P (2021) Blowout Morphometrics and Mass Balances. Front. Earth Sci. 9:669440. doi: 10.3389/feart.2021.669440

Received

18 February 2021

Accepted

04 May 2021

Published

25 May 2021

Volume

9 - 2021

Edited by

Zhiwei Xu, Nanjing University, China

Reviewed by

Zhengyao Liu, Shaanxi Normal University, China

Yang Bai, Huizhou University, China

Leví García-Romero, University of Las Palmas de Gran Canaria, Spain

Updates

Copyright

*Correspondence: Thomas R Baird,

This article was submitted to Environmental Informatics and Remote Sensing, a section of the journal Frontiers in Earth Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics