Semi-analytical solutions of groundwater flow in multi-zone (patchy) wedge-shaped aquifers

Alluvial fans are potential sites of potable groundwater in many parts of the world. Characteristics of alluvial fans sediments are changed radially from high energy coarse-grained deposition near the apex to low energy fine-grained deposition downstream so that patchy wedge-shaped aquifers with radial heterogeneity are formed. The hydraulic parameters of the aquifers (e.g. hydraulic conductivity and specific storage) change in the same fashion. Analytical or semi-analytical solutions of the flow in wedge-shaped aquifers are available for homogeneous cases. In this paper we derive semi-analytical solutions of groundwater flow to a well in multi-zone wedge-shaped aquifers. Solutions are provided for three wedge boundary configurations namely: constant head–constant head wedge, constant head–barrier wedge and barrier–barrier wedge. Derivation involves the use of integral transforms methods. The effect of heterogeneity ratios of zones on the response of the aquifer is examined. The results are presented in form of drawdown and drawdown derivative type curves. Heterogeneity has a significant effect on over all response of the pumped aquifer. Solutions help understanding the behavior of heterogeneous multi-zone aquifers for sustainable development of the groundwater resources in alluvial fans.     

Patterns and dynamics of river–aquifer exchange with variably-saturated flow using a fully-coupled model

The parallel physically-based surface–subsurface model PARFLOW was used to investigate the spatial patterns and temporal dynamics of river–aquifer exchange in a heterogeneous alluvial river–aquifer system with deep water table. Aquifer heterogeneity at two scales was incorporated into the model. The architecture of the alluvial hydrofacies was represented based on conditioned geostatistical indicator simulations. Subscale variability of hydraulic conductivities (K) within hydrofacies bodies was created with a parallel Gaussian simulation. The effects of subscale heterogeneity were investigated in a Monte Carlo framework. Dynamics and patterns of river–aquifer exchange were simulated for a 30-day flow event. Simulation results show the rapid formation of saturated connections between the river channel and the deep water table at preferential flow zones that are characterized by high conductivity hydrofacies. Where the river intersects low conductivity hydrofacies shallow perched saturated zones immediately below the river form, but seepage to the deep water table remains unsaturated and seepage rates are low. Preferential flow zones, although only taking up around 50% of the river channel, account for more than 98% of total seepage. Groundwater recharge is most efficiently realized through these zones. Subscale variability of K_sat slightly increased seepage volumes, but did not change the general seepage patterns (preferential flow zones versus perched zones). Overall it is concluded that typical alluvial heterogeneity (hydrofacies architecture) is an important control of river–aquifer exchange in rivers overlying deep water tables. Simulated patterns and dynamics are in line with field observations and results from previous modeling studies using simpler models. Alluvial heterogeneity results in distinct patterns and dynamics of river–aquifer exchange with implications for groundwater recharge and the management of riparian zones (e.g. river channel-floodplain connectivity via saturated zones).     

Application of kinematic wave theory for predicting flash flood hazards on coupled alluvial fan–piedmont plain landforms

A rainstorm that caused a severe flash flood on the piedmont plain at the toe positions of two alluvial fans located to the west of the Organ Mountains in Dona Ana County, New Mexico, USA, is analysed. The space–time distributions of rainfall are evaluated from the Next Generation Weather Radar (NEXRAD) and overland flow is modelled as kinematic wave. The spatial distribution of rainfall shows a topographic control. The greatest rainfall depth, duration, and intensity occurred at the higher elevation mountain slopes and decreased with decreasing elevation from the alluvial fans to the piedmont plain. The alluvial fan–piedmont plain system is modelled by coupling divergent and rectangular overland flow planes. Explicit finite difference approximations, hybridized with the analytical method of characteristics, are made to the kinematic wave equations to account for the spatial and temporal distribution of the rainfall and variable boundary conditions. Simulation results indicate that sheet‐flow floodwater elevations rise (1) in a nonlinear fashion from the apex to toe positions of the alluvial fans, and (2) near linearly from the toe positions of the alluvial fans onto the piedmont plains with the formation of kinematic shocks near the middle to the upstream end of the plane at times between the initiation of the rainstorm and the time of concentration of the plane. Thus, the maximum flooding occurs at the middle or upstream sections of the piedmont plains regardless of the pattern of space–time variability of rainfall. These results are in agreement with observed geomorphologic features suggesting that piedmont plains are naturally flood‐prone areas. This case study demonstrates that flood hazards on piedmont plains can exceed those on alluvial fans. The models presented in this study suggest that the flood hazard zones on coupled alluvial fan–piedmont plain landforms should be delineated transverse to the flow directions, as opposed to the flood hazard zones with boundaries in the longitudinal direction of the axis of an alluvial fan. Copyright © 2003 John Wiley & Sons, Ltd.     

Topographic predictors of susceptibility to alluvial fan flooding,Southern Apennines

The flooding susceptibility of alluvial fans in the Southern Apennines has long been neglected. To partly address this oversight, we focus on the region of Campania which contains highly urbanized piedmont areas particularly vulnerable to flooding. Our findings are based on stratigraphic analysis of the fans and morphometric analysis of the basin‐fan systems. Using geomorphological analysis we recognized active alluvial fans while stratigraphic analysis together with statistical analysis of the morphometric variables was used to classify the fans in terms of the transport process involved. The results indicate that in the geological context examined, the best discrimination between debris flow (Df) and water flood (Wf) processes is achieved by means of two related variables, one for the basin (feeder channel inclination, Cg) and one for the fan (fan length, Fl). The probability that an unclassified fan belongs to group Wf is computed by applying a logistic function in which a P value exceeding 0.5 indicates that a basin/fan system belongs to group Wf. This important result led to the classification of the entire basin/fan system data. As regards process intensity, debris flow‐dominated fans are susceptible to the occurrence of flows with high viscosity and hence subject to more severe events than water flood‐dominated fans. Bearing this in mind, the data gathered in this study allow us to detect where alluvial fan flooding might occur and give information on the different degrees of susceptibility at a regional scale. Regrettably, urban development in recent decades has failed to take the presence of such alluvial fans into account due to the long recurrence time (50–100 years) between floods. This paper outlines the distribution of such susceptibility scenarios throughout the region, thereby constituting an initial step to implementing alluvial fan flooding control and mitigation. Copyright © 2012 John Wiley & Sons, Ltd.     

The coupling status of alluvial fans and debris cones: a review and synthesis

Alluvial fans and debris cones link two zones of the fluvial system (e.g. hillslope gully systems to stream channels; mountain catchment sediment source areas to main river systems or to sedimentary basins) and therefore have important coupling or buffering roles. These roles may be both functional and preservational. The functional role includes debris‐cone coupling, which controls sediment supply from hillslope gully systems to stream channels, influencing channel morphology. Coupling through larger alluvial fans, expressed by fanhead trenching, causes a distal shift in sedimentation zones, or when expressed by through‐fan trenching, causes complete sediment by‐pass. The preservational role stems from the fact that fans and cones are temporary sediment storage zones, and may preserve a record of source–area environmental change more sensitively than would sediments preserved further downsystem. Fan coupling mechanisms include distally‐induced coupling (basal scour, ‘toe cutting’, marginal incision) and proximally‐induced coupling (fanhead and midfan trenching). These mechanisms lead initially to partial coupling, either extending the immediate sediment source area to the stream system or shifting the focus of sedimentation distally. Complete coupling involves transmission of sediment from the feeder catchment through the fan environment into the downstream drainage or a sedimentary basin. The implications of coupling relate to downstream channel response, fan morphology, sedimentation patterns and vertical sedimentary sequences. Temporal and spatial scales of coupling are related, and with increasing scales the dominant controls shift from storm events to land cover to climatic and base‐level change and ultimately to the relationships between tectonics and accommodation space. Finally, future research challenges are identified. Modern dating techniques and sophisticated analysis of remotely sensed data can greatly improve our understanding of fan dynamics, and should lead to better cross‐scale integration between short‐term process‐based approaches and long‐term sedimentological applications, while maintaining high quality field‐based observations. Copyright © 2011 John Wiley & Sons, Ltd.     

Geometry of alluvial fans: Effect of discharge and sediment size

The slope of an alluvial fan increases with increasing debris size and sediment concentration in the flow, and decreases with increasing discharge. Laboratory studies suggest that the discharge which controls this slope, or dominant discharge, is that which is equalled or exceeded one quarter to one third of the time that flow occurs on the fan. In contrast, the dominant discharge in perennial alluvial rivers is equalled or exceeded only about 5 per cent of the time that flow occurs in the river. The dominant discharge on fans increases with increasing debris size, reflecting the importance of threshold stress. The slope of some natural and most laboratory alluvial fans is steepest on the flanks and gentlest along the axis. Consideration of the momentum of water debouching onto a fan at its apex suggests that the difference in slope between axis and flank should be greatest on steep fans composed of relatively non-cohesive materials because on such fans higher discharges tend to flow down the axis, whereas lower discharges can be turned to course down the flanks. On fans with gentle slopes or composed of more cohesive material the higher discharges can also be turned toward the flanks, so on such fans the difference in slope between the axis and flank is less pronounced. Field and laboratory observations support this interpretation. Because deposition at any one time on an alluvial fan is localized, some areas aggrade while others remain at a fixed elevation. This process is treated as a Markov process with the probability of diversion from an area of active deposition into an adjacent lower area increasing as the height of the active area above the mean or ‘ideal’ surface increases. Analysis of data from laboratory and natural fans suggests that the amplitude of such surface irregularities is greater on fans composed of coarser material. The data on natural fans also suggest an increase in amplitude of the irregularities with increasing fan area.     

Modelling alluvial fan morphology

A mathematical model which estimates the scale-independent sediment surface profile of alluvial fans has been developed. This model utilizes a diffusive sediment transport model and an unsteady, radial flow, conservation relationship. These equations are approximately solved assuming a quasi-steady-state closure with appropriate modelling assumptions for two end member fan types: (1) fans where most of the fan surface is depositionally active (denoted here as ‘homogeneous’) and (2) fans characterized by channelling and sediment sorting processes. The fundamental result for these two fan types is a dimensionless sediment profile relationship which approximates most fan surfaces. The model suggests that the overall dimensionless morphology of alluvial fans is governed more by fundamental diffusion principles in sediment deposition than by individual environmental or basin characteristics. Additionally, this work potentially can be extended to model temporal variation in fan development. Preliminary comparison with alluvial fan profiles is reasonable, indicating that this model provides useful qualitative and quantitative information relating to alluvial fan process and morphology. Copyright © 1999 John Wiley & Sons, Ltd.     

Hyperconcentrated flows on a forested alluvial fan of eastern Canada: geomorphic characteristics,return period,and triggering scenarios

The assessment of the dominant flow type on alluvial fans usually refers to two categories: debris‐flow fans (i.e. sediment gravity flows) and fluvial fans (i.e. fluid gravity flows). Here we report the results of combined morphometric, stratigraphic and sedimentological approaches which suggest that hyperconcentrated flows, a transitional process rheologically distinct from debris flows and floods and sometimes referred to as debris floods, mud floods, or transitional debris flows, are the dominant fan building process in eastern Canada. These flows produce transitional facies between those of debris flows which consist of a cohesive matrix‐supported diamicton, and those of river flows which display more distinct stratification. The size of the blocks in the channels and the abrasion scars at the base of several trees attest to the high transport capacity of these flows. The fan channels are routed according to various obstacles comprised primarily of woody debris that impede sediment transit. However, these conditions of sediment storage are combined with readily available sediment due to the friable nature of the local lithology. Tree‐ring analysis allowed the reconstruction of eight hydrogeomorphic events which are characterized by a return period of 9.25 years for the period 1934–2008, although most of the analyzed events occurred after 1970. Historical weather data analysis indicates that they were related to rare hydrometeorological events at regional and local scales. This evidence led to the elaboration of weather scenarios likely responsible for triggering flows on the fan. According to these scenarios, two distinct hydrologic regimes emerge: the torrential rainfall regime and the nival regime related to snowmelt processes. Hydrogeomorphic processes occurring in a cold‐temperate climate, and particularly on small forested alluvial fans of north‐eastern North America, should receive more attention from land managers given the hazard they represent, as well as because of their sensitivity to various meteorological parameters. Copyright © 2014 John Wiley & Sons, Ltd.     

Controls on modern alluvial fan processes in the central Alps,northern Italy

Alluvial fan development in Alpine areas is often affected by catastrophic sedimentary processes associated with extreme ?oods events, causing serious risks for people living on the fans. Hazard assessment in these areas depends on proper identi?cation of the dominant sedimentary processes on the fans. Data from a set of 209 alluvial fans from the central Alps of Italy are presented in this paper and analysed with the help of various statistical techniques (linear regression, principal components analysis, cluster analysis, discriminant analysis and logistic regression). First, we used modern sedimentary facies and historical records (?ood events since 15th century), to distinguish between the two dominant sedimentary processes on alluvial fans: debris ?ows and stream?ows. Then, in order to analyse the main controls on past and present fan processes, 36 morphological, geological and land‐use variables were analysed. As with observations for arid‐environment fans, catchment morphology is the most in?uential factor in the study area, whereas geology and land use are minor controls. The role of climatic change and landsliding within the catchments also seems to be very important and is discussed. Statistical techniques also help in differentiating groups of alluvial fans by sets of controlling factors, including stage and type of evolution. Finally, by using discriminant analysis and logistic regression, we classi?ed alluvial fans according to the dominant sedimentary process, with a success rate ranging between 75 and 92 per cent. Copyright © 2004 John Wiley & Sons, Ltd.     

THE ROLE OF ALLUVIAL FANS IN THE MOUNTAIN FLUVIAL SYSTEMS OF SOUTHEAST SPAIN: IMPLICATIONS OF CLIMATIC CHANGE

The mountain fluvial systems of southeast Spain involve sediment supply from steep mountain slopes into headwater channels. Alluvial fans often occur where these headwater channels emerge from the mountain areas, and may influence the connectivity of the sediment transport system from the mountain source areas to the main lowland drainages. Critical in this role is whether the alluvial fans are aggrading or dissecting, and whether there is a break or continuity in the channel through the fan environment. Previous work has identified some of the factors influencing the behaviour of the alluvial fans in southeast Spain. This paper deals with the mountain front alluvial fans in the semi-arid areas of Murcia and Almeria provinces. It attempts, by mapping the location of alluvial fans, then their classification into aggrading or dissecting fans, to identify the extent to which the mountain fluvial systems are buffered by aggrading alluvial fans or exhibit channel continuity through the mountain front environment. It further considers the implications of climatically induced changes between aggradational and dissectional behaviour on alluvial fans.     

Depositional process and alluvial fan-drainage basin morphometric relationships near banff,Alberta, Canada

Two distinct types of alluvial fans occur in the Bow River Valley, Alberta, Canada: fluvially dominated and debris flow dominated. Large, gently sloping fans dominated by fluvial processes are associated with large and less rugged drainage basins, and small rugged basins have produced small, steep fans dominated by debris flow processes. Quantitative analysis demonstrates that strong fan-basin morphometric relationships occur despite a short fan history. Statistical analysis of fan area-basin area relationships indicate that debris flow fan areas do not increase in size as quickly as contributing basins. The relationship of fluvial fan area to basin area is not statistically significant. However, this relationship is probably affected by fan erosion. Examination of fan slope to basin ruggedness relationships indicates that fan slope increases more rapidly than basin ruggedness for both fan types. This is likely related to non-linear discharge and sediment size effects on fluvial fans, and reworking of larger fan surfaces by fluvial processes on debris flow fans.     

Catchment lithology as a major control on alluvial megafan development,Kohrud Mountain range,central Iran

The relative importance of tectonics, climate, base level and source lithology as primary factors on alluvial‐fan evolution, fan morphology and sedimentary style remain in question. This study examines the role of catchment lithology on development and evolution of alluvial megafans (>30 km in length), along the flanks of the Kohrud Mountain range, NE Esfahan, central Iran. These fans toe out at axial basin river and playa‐fringe sediments towards the centre of basin and tectonics, climatic change and base‐level fluctuations, were consistent for their development. They formed in a tectonically active basin, under arid to semiarid climate and a long term (Plio‐Pleistocene to Recent) change from wetter to drier conditions. The key differences between two of these fans, Soh and Zefreh fans, along the west and south flanks of this mountain range, is that their catchments are underlain by dissimilar bedrock types. The source‐area lithologies of the Soh and Zefreh fans are in sedimentary and igneous terrains, respectively, and these fans developed their geometry mainly in response to different weathering intensities of their catchment bedrock lithologies. Fan surface mapping (based on 1/50000 topographic maps, satellite images, and fieldwork), reveals that the geomorphic evolution of these fans differs in that the relatively large‐scale incision and through trenching of the Soh fan is absent in the Zefreh fan. Whereas the limited sediment supply of the Soh fan has resulted in a deep incised channel, the Zefreh fan has remained aggradational with little or no trenching into proximal to medial fan surface due to its catchment bedrock geology, composed mainly by physically weathered volcaniclastic lithology and characterized by high sediment supply for delivery during episodic flash floods. Sediment supply, which is mainly a function of climate and source lithology, is a dominant driver behind the development of fan sequences in alluvial megafans. Copyright © 2012 John Wiley & Sons, Ltd.     

EXTRACTING FEATURES OF ALLUVIAL FAN AND DISCUSSING LANDFORMS EVOLUTION BASED ON HIGH-RESOLUTION TOPOGRAPHY DATA: TAKING ALLUVIAL FAN OF LAOHUSHAN ALONG HAIYUAN FAULT ZONE AS AN INSTANCE

Range-front alluvial fan deposition in arid and semiarid environments records vast amounts of climatic and tectonic information. Differentiating and characterizing alluvial fan morphology is an important part in Quaternary alluvial fan research. Traditional method such as field observations is a most important part of deciphering and mapping the alluvial fan. Large-scale automatically mapping of alluvial fan stratigraphy before traditional field observations could provide guidance for mapping alluvial fan morphology, thus improving subsequent field work efficiency. In this research, high-resolution topographic data were used to quantify relief and roughness of alluvial fan within the Laohushan. These data suggest that mean surface roughness plotted against the size of the moving window is characterized by an initial increase in surface roughness with increased window size, but it shows no longer increase as a function of windows size. These data also suggest that alluvial fans in this study site smooth out with time until a threshold is crossed where roughness increases at greater wavelength with age as a result of surface runoff and headward tributary incision into the oldest surfaces which suggests the evolution process of alluvial fan. Researchers usually differentiate alluvial morphology by mapping characteristics of fan surface in the field by describing surface clast size, rock varnish accumulation, and desert pavement development and analysis of aerial photographs or satellite imagery. Recently, the emergence of high-resolution topographic data has renewed interest in the quantitative characterization of alluvial and colluvium landforms. Surface morphology that fan surface initially tends to become smoother with increasing age due to the formation of desert pavement and the degradation of bar-and-swale topography and subsequently, landforms become more dissected due to tectonics and climatic change induced increased erosion and channelization of the surface with time is widely used to distinguish alluvial fan types. Those characteristics would reflect various kinds of morphology metrics extracted from high-resolution topographic data. In the arid and semiarid regions of northwestern China, plenty of alluvial fans are preserved completely for lack of artificial reforming, and there exists sparse surface vegetation. In the meantime, range-front alluvial fan displaced by a number of active faults formed a series of dislocated landforms with different offsets which is a major reference mark in fault activity research. In this research, six map units(Qf₆-Qf₁), youngest to oldest, were observed in the study area by mapping performed by identifying geomorphic features in the field that are spatially discernible using hill-shade and digital orthophoto map. Alluvial fan relief and roughness were computed across multiple observation scales(2m×2m to 100m×100m)based on the topographic parameters of altitude difference and standard deviation of slope, curvature and aspect. In this research, mean relief keeps increasing with increased window size while mean surface roughness is characterized by a rapid increase over wavelengths of 6~15m, representing the typical length scale of bar-and-swale topography. At longer wavelengths, surface roughness values increase by only minor amounts, suggesting the topographic saturation length is 6~15m for those fan surfaces in which saturation length of standard deviation of curvature is less than 8m. Box and whisker plot of surface roughness averaged over 8m² for each alluvial fan unit in the study area suggests that the pattern of surfaces smoothing out with age and then starting to become rougher again as age increases further beyond Qf₄ or Qf₃ unit. The younger alluvial fan is characterized by prominent bar-and-swale while the older alluvial fan is characterized by tributaries headward incision. Cumulative frequency distributions of relief and surface roughness in Figure 8 are determined in an 8m by 8m moving window for the comparison of six alluvial fan units in the northeast piedmont of Laohushan. From these distributions we know that Qf₆ and Qf₁ reflect the prominent relief which is related to bar-and-swale and tributaries headward incision respectively, while Qf₄ and Qf₃ reflect the moderate relief which is related to subdued topography. Surface roughness, in addition to facilitating the characterization of individual fan units, lends insight to alluvial landform development. We summarize an alluvial landform evolutionary scheme which evolves four stages depending on characteristics of alluvial fan morphology development and features of relief and roughness. The initial stage in this study site is defined as the active alluvial fan channels with bars of coarse cobbles and boulders and swales consisting of finer-grained pebbles and sand which could be reflected by high mean relief and mean roughness values. As time goes, bar-and-swale topography is still present, but an immature pavement, composed of finer grained clasts, has started to form. In the third stage, the bar-and-swale topography on the fan surface is subdued, yet still observable, with clasts ranging from pebbles to cobbles in size and there exists obvious headward tributary incision. Eventually, tributary channels form from erosion by surface runoff. Headward incision of these tributaries wears down the steep walls of channels that are incised through the stable, planar surface, transforming the oldest alluvial landforms into convex hillslopes, leaving only small remnants of the planar surface intact. Those evolutionary character suggests that alluvial fans in this area smooth out with time, however, relief or roughness would be translated to increase at greater wavelength with age until a threshold is crossed. This research suggests that relief and roughness calculated from high-resolution topographic data of this study site could reflect alluvial fan morphology development and provide constraint data to differentiate alluvial fan unit.     

USING UAV PHOTOGRAMMETRY TECHNOLOGY TO EXTRACT INFORMATION OF TECTONIC ACTIVITY OF COMPLEX ALLUVIAL FAN——A CASE STUDY OF AN ALLUVIAL FAN IN THE SOUTHERN MARGIN OF BARKOL BASIN

Alluvial fans that are in the process of development always show complex geomorphic features due to natural modification. Accordingly, analyzing these fans whether to be influenced by tectonic deformation is one of the technique difficulties in active tectonic studies. Complex alluvial fans are the focus of the study of active tectonics such as fracture mapping and activity behavior analysis, for they have often retained important structural information. Traditional measurement methods, such as satellite remote sensing, RTK GPS and Lidar, are difficult to meet the demand for the study of micro tectonic deformation because of the reason of accuracy or cost performance. The recent UAV photogrammetry technology, due to its many advantages such as low cost, high resolution, and efficiency of exporting DEM and DOM data, has been widely used in three-dimensional modeling, ground mapping and other fields. In the quantitative study of active tectonics, this technology fills up the deficiency in the research of the micro structure of the traditional measurement. Through detailed field investigations and paleoseismic trenching, we further used this technology to obtain the topographic data of a complex alluvial fan located at the southern marginal fault of Barkol Basin, Xinjiang. Pointing at the alluvial fans that are in the process of development, and on the basis of topographic analysis and image processing for DEM, we take the research method of secondary partitions of the geomorphic surface and cut the alluvial fans longitudinally according to the difference of its age. Through the establishment of profile cluster within each partition, separate analysis and data contrast with the adjacent partitions, we acquired the tectonic activity information during the development of alluvial fan. The tectonic vertical deformation of this alluvial fan is about 2.5m.     

Accounting for Aquifer Heterogeneity from Geological Data to Management Tools

A nested workflow of multiple‐point geostatistics (MPG) and sequential Gaussian simulation (SGS) was tested on a study area of 6 km² located about 20 km northwest of Quebec City, Canada. In order to assess its geological and hydrogeological parameter heterogeneity and to provide tools to evaluate uncertainties in aquifer management, direct and indirect field measurements are used as inputs in the geostatistical simulations to reproduce large and small‐scale heterogeneities. To do so, the lithological information is first associated to equivalent hydrogeological facies (hydrofacies) according to hydraulic properties measured at several wells. Then, heterogeneous hydrofacies (HF) realizations are generated using a prior geological model as training image (TI) with the MPG algorithm. The hydraulic conductivity (K) heterogeneity modeling within each HF is finally computed using SGS algorithm. Different K models are integrated in a finite‐element hydrogeological model to calculate multiple transport simulations. Different scenarios exhibit variations in mass transport path and dispersion associated with the large‐ and small‐scale heterogeneity respectively. Three‐dimensional maps showing the probability of overpassing different thresholds are presented as examples of management tools.     

Applied Mathematics in EM Studies with Special Emphasis on an Uncertainty Quantification and 3-D Integral Equation Modelling

Despite impressive progress in the development and application of electromagnetic (EM) deterministic inverse schemes to map the 3-D distribution of electrical conductivity within the Earth, there is one question which remains poorly addressed—uncertainty quantification of the recovered conductivity models. Apparently, only an inversion based on a statistical approach provides a systematic framework to quantify such uncertainties. The Metropolis–Hastings (M–H) algorithm is the most popular technique for sampling the posterior probability distribution that describes the solution of the statistical inverse problem. However, all statistical inverse schemes require an enormous amount of forward simulations and thus appear to be extremely demanding computationally, if not prohibitive, if a 3-D set up is invoked. This urges development of fast and scalable 3-D modelling codes which can run large-scale 3-D models of practical interest for fractions of a second on high-performance multi-core platforms. But, even with these codes, the challenge for M–H methods is to construct proposal functions that simultaneously provide a good approximation of the target density function while being inexpensive to be sampled. In this paper we address both of these issues. First we introduce a variant of the M–H method which uses information about the local gradient and Hessian of the penalty function. This, in particular, allows us to exploit adjoint-based machinery that has been instrumental for the fast solution of deterministic inverse problems. We explain why this modification of M–H significantly accelerates sampling of the posterior probability distribution. In addition we show how Hessian handling (inverse, square root) can be made practicable by a low-rank approximation using the Lanczos algorithm. Ultimately we discuss uncertainty analysis based on stochastic inversion results. In addition, we demonstrate how this analysis can be performed within a deterministic approach. In the second part, we summarize modern trends in the development of efficient 3-D EM forward modelling schemes with special emphasis on recent advances in the integral equation approach.     

Geomorphological evolution and sediment stratigraphy of numerically simulated alluvial fans

When a sediment laden river reaches a flat basin area the coarse fraction of their sediment load is deposited in a cone shaped structure called an alluvial fan. In this article we used the State Space Soil Production and Assessment Model (SSSPAM) coupled landform–soilscape evolution model to simulate the development of alluvial fans in two- and three-dimensional landforms. In SSSPAM the physical processes of erosion and armouring, soil weathering and sediment deposition were modelled using state-space matrices, in both two and three dimensions. The results of the two-dimensional fan showed that the fan grew vertically and laterally keeping a concave up long profile. It also showed a downstream fining of the sediments along the fan profile. Both of these observations are in agreement with available literature concerning natural and experimental fan formations. Simulations with the three-dimensional landform produced a fan with a semicircular shape with concave up long profiles and concave down cross profiles which is typical for fans found in nature and ones developed in laboratory conditions. During the simulation the main channel which brings sediment to the fan structure changed its position constantly leading to the semicircular shape of the fan. This behaviour is similar to the autogenic process of ‘fanhead trenching’ which is the major mechanism of sediment redistribution while the fan is developing. The three-dimensional fan simulation also exhibited the downstream fining of sediments from the fan apex to the peripheries. Further, the simulated fan also developed complex internal sediment stratification which is modelled by SSSPAM. Currently such complex sediment stratification is thought to be a result of allogenic processes. However, this simulation shows that, such complex internal sediment structures can develop through autogenic processes as well. © 2020 John Wiley & Sons, Ltd.     

Mechanisms for avulsion on alluvial fans: Insights from high-frequency topographic data

Avulsion is a key process in building alluvial fans, but it is also a formidable natural hazard. Based on laboratory experiments monitored with novel high-frequency photogrammetry, we present a new model for avulsion on widely graded gravel fans. Previous experimental studies of alluvial fans have suggested that avulsion occurs in a periodic autogenic cycle, that is thought to be mediated by the gradient of the fan and fan-channel. However, those studies measured gradients at low spatial or temporal resolutions, which capture temporally or spatially averaged topographic evolution. Here, we present high-resolution (1 mm), high-frequency (1-minute) topographic data and orthophotos from an alluvial fan experiment. Avulsions in the experiment were rapid and, in contrast to some previous experimental studies, avulsion occurrence was aperiodic. Moreover, we found little evidence of the back-filling observed at coarser temporal and spatial resolutions. Our observations suggest that avulsion is disproportionately affected by sediment accumulation in the channel, particularly around larger, less mobile grains. Such in-channel deposition can cause channel shifting that interrupts the autogenic avulsion cycle, so that avulsions are aperiodic and their timing is more difficult to predict.     

济阳坳陷第三系隐蔽藏储层预测配套技术

利用沉积学、层序地层学理论分析和研究沉积相特征,得出探区内发育的储层主要有冲积扇、河流、三角洲、浊积扇和滩坝五大沉积体系,提出了上第三系”水流控砂”、下第三系”坡槽控砂”的沉积模式,得到储层按构造带分布的规律.隐蔽性油气藏以”断导”模式成藏,与储层发育相匹配形成四个大的隐蔽性油气藏群.研究地震储层描述方法的适应性和敏感性,利用正、反演类比方法研究地震反射特征,实现了地震相和沉积相的转换.剖析典型的油气藏类型,分析成功的经验,总结失败的教训,得到了针对河流、冲积扇、三角洲、滩坝和浊积扇的描述流程和配套技术系列.     

Erosion of fixed dunes in the Sahel,Central Niger

During the 1974 rainy season gullies of the order of 150–300 m long were active on the flanks of fixed Pleistocene dunes in the vicinity of Janjari, central Niger. The gullies terminate on small (0.5–1 × 10⁴m²) alluvial fans, where sand deposition had occurred below the intersection point. Gully activity was probably the result of a particularly severe storm. Comparison between fixed dunes subject to differing climatic conditions suggests that gullying and alluvial fan formation may be characteristic of a restricted morphoclimatic zone.