Numerical study of turbulent flow over complex aeolian dune fields: the White Sands National Monument

Mesoscale structure of the atmospheric boundary layer across a natural roughness transition

The structure and intensity of turbulence in the atmospheric boundary layer (ABL) drive fluxes of sediment, contaminants, heat, moisture, and CO[Formula: see text] at the Earth's surface. Where ABL flows encounter changes in roughness-such as cities, wind farms, forest canopies, and landforms-a new mesoscopic flow scale is introduced: the internal boundary layer (IBL), which represents a near-bed region of transient flow adjustment that develops over kilometers. Measurement of this new mesoscopic scale lies outside present observational capabilities of ABL flows, and simplified models fail to capture the sensitive dependence of turbulence on roughness geometry. Here, we use large-eddy simulations, run over high-resolution topographic data and validated against field observations, to examine the structure of the ABL across a natural roughness transition: the emergent sand dunes at White Sands National Park. We observe that development of the IBL is triggered by the abrupt transition from smooth playa surface to dunes; however, continuous changes in the size and spacing of dunes over several kilometers influence the downwind patterns of boundary stress and near-bed turbulence. Coherent flow structures grow and merge over the entire [Formula: see text]10 km distance of the dune field and modulate the influence of large-scale atmospheric turbulence on the bed. Simulated boundary stresses in the developing IBL counter existing expectations and explain the observed downwind decrease in dune migration, demonstrating a mesoscale coupling between flow and form that governs landscape dynamics. More broadly, our findings demonstrate the importance of resolving both turbulence and realistic roughness for understanding fluid-boundary interactions in environmental flows.

Stability of the Interaction between Two Sand Dunes in an Idealized Laboratory Experiment

Sand dunes, which arise spontaneously due to the dynamical interplay between a sedimentary interface and a fluid flow, are one of the most famous examples of emergence in a geological system. The large scale organization of a dune field is believed to be controlled by pairwise (either remote or direct) dune-dune interactions. Recent studies have shown that remote long-range feedback is closely related to the turbulent wake structure forming downstream of a dune. Here, we study the stability of an idealized two-dune system arising as a consequence of such remote, wake-induced interactions. The system is realized in a subaqueous quasi-2D laboratory experiment and the results are compared with a qualitative dynamical systems model. Despite its simplicity, the system exhibits rich dynamical behavior. In particular, we show that, depending on the parameter regime, the dune-dune feedback can either stabilize or destabilize the symmetric dune configuration, and we demonstrate the existence of an asymmetric attracting state coupling dunes of different sizes.

Wake Induced Long Range Repulsion of Aqueous Dunes

Sand dunes rarely occur in isolation, but usually form vast dune fields. The large scale dynamics of these fields is hitherto poorly understood, not least due to the lack of longtime observations. Theoretical models usually abstract dunes in a field as self-propelled autonomous agents, exchanging mass, either remotely or as a consequence of collisions. In contrast to the spirit of these models, here we present experimental evidence that aqueous dunes interact over large distances without the necessity of exchanging mass. Interactions are mediated by turbulent structures forming in the wake of a dune, and lead to dune-dune repulsion, which can prevent collisions. We conjecture that a similar mechanism may be present in wind driven dunes, potentially explaining the observed robust stability of dune fields in different environments.

Measuring Hyperscale Topographic Anisotropy as a Continuous Landscape Property

Several landforms are known to exhibit topographic anisotropy, defined as a directional inequality in elevation. The quantitative analysis of topographic anisotropy has largely focused on measurements taken from specific landforms, ignoring the surrounding landscape. Recent research has made progress in measuring topographic anisotropy as a distributed field in natural landscapes. However, current methods are computationally inefficient, as they require specialized hardware and computing environments, or have a limited selection of scales that undermines the feasibility and quality of multiscale analyses by introducing bias. By necessity, current methods operate with a limited set of scales, rather than the full distribution of possible landscapes. Therefore, we present a method for measuring topographic anisotropy in the landscape that has the computational efficiency required for hyperscale analysis by using the integral image filtering approach to compute oriented local topographic position (LTP) measurements, coupled with a root-mean-square deviation (RMSD) model that compares directional samples to an omnidirectional sample. Two tools were developed: One to output a scale signature for a single cell, and the other to output a raster containing the maximum anisotropy value across a range of scales. The performances of both algorithms were tested using two data sets containing repetitive, similarly sized and oriented anisotropic landforms, including a dune field and a drumlin field. The results demonstrated that the method presented has the robustness and sensitivity to identify complex hyperscale anisotropy such as nested features (e.g., a drumlin located within a valley).