Although sea-level highstands are typically associated with sediment-starved continental shelves, high sea level does not hinder major river floods. Turbidity currents generated by plunging of sediment-laden rivers at the fluvial-marine interface, known as hyperpycnal flows, allow for cross-shelf transport of suspended sand beyond the coastline. Hyper pycnal flows in southern California have deposited six subaqueous fans on the shelf of the northern Santa Barbara Channel in the Holocene. Using eight cores and nine grab samples, we describe the deposits, age, and stratigraphic architecture of two fans in the Santa Barbara Channel. Fan lobes have up to 3 m of relief and are composed of multiple hyperpycnite beds ~5 cm to 40 cm thick. Deposit architecture and geometry suggest the hyperpycnal flows became positively buoyant and lifted off the seabed, resulting in well-sorted, structureless, elongate sand lobes. Contrary to conventional sequence stratigraphic models, the presence of these features on the continental shelf suggests that active-margin shelves may locally develop high-quality reservoir sand bodies during sea-level highstands, and that such shelves need not be solely the site of sediment bypass. These deposits may provide a Quaternary analogue to many well-sorted sand bodies in the rock record that are interpreted as turbidites but lack typical Bouma-type features.

Deltaic river networks naturally reorganize as interconnected channels move to redistribute water, sediment, and nutrients across the delta plain. Network change is documented in decades of satellite imagery and laboratory experiments, but our ability to measure and understand channel movements is limited: existing methods are difficult to employ efficiently and struggle to distinguish between gradual movements (channel migration) and abrupt shifts in river course (channel avulsions). Here, we present a method to extract channel migration from plan‐view imagery using particle image velocimetry (PIV). Although originally designed to track particles moving in a fluid, PIV can be adapted to track channels moving on the delta surface, based on input estimates of channel width, migration timescale, and maps of the wet‐dry interface. Results for a delta experiment show that PIV‐derived vector fields accurately capture channel‐bank movements, as compared to manually drawn maps and an independent image‐registration technique. Unlike other methods, PIV targets the process of channel migration, excluding changes associated with channel avulsions and overbank flow. PIV‐derived migration rates from the experiment span an order of magnitude and are reduced under lower sediment supply and during sea‐level rise, supporting recent models. Together, results indicate that PIV offers a fast and reliable way to measure channel migration in river networks, that channel migration rates under non‐cohesive conditions can displace channels a distance comparable to their width in the time needed to aggrade ∼10% of the channel depth, and that migration direction is ∼60% orthogonal to mean flow direction and ∼40% flow‐parallel overall.

Understanding subsurface structure and groundwater flow in deltaic aquifers is essential for evaluating the vulnerability of groundwater resources in delta systems. Deltaic aquifers contain coarse‐grained paleochannels that preserve a record of former surface river channels as well as fine‐grained floodplain deposits. The distribution of these deposits and how they are interconnected control groundwater flow and contaminant transport. In this work, we link depositional environments of deltaic aquifers to stratigraphic (static) and flow and transport (dynamic) connectivity metrics. Numerical models of deltaic stratigraphy were generated using a reduced‐complexity numerical model (DeltaRCM) with different input sand fractions (ISF) and rates of sea‐level rise (SLR). The groundwater flow and advective transport behavior of these deltas were simulated using MODFLOW and MODPATH. By comparing the static and dynamic metrics calculated from these numerical models, we show that groundwater behavior can be predicted by particular aspects of the subsurface architecture, and that horizontal and vertical connectivity display different characteristics. We also evaluate relationships between connectivity metrics and two environmental controls on delta evolution: ISF and SLR rate. The results show that geologic setting strongly influences both static and dynamic connectivity in different directions. These results provide insights into quantitatively differentiated subsurface hydraulic behavior between deltas formed under different external forcing (ISF and SLR rate) and they are a potential link in using information from delta surface networks and depositional history to predict vulnerability to aquifer contamination.

Hyperpycnal gravity currents rapidly transport sediment across shore from rivers to the continental shelf and deep sea. Although these geophysical processes are important sediment dispersal mechanisms, few distinct geomorphic features on the continental shelf can be attributed to hyperpycnal flows. Here we provide evidence of large depositional features derived from hyperpycnal plumes on the continental shelf of the northern Santa Barbara Channel, California, from the combination of new sonar, lidar, and seismic reflection data. These data reveal lobate fans directly offshore of the mouths of several watersheds known to produce hyperpycnal concentrations of suspended sediment. The fans occur on an upwardly concave section of the shelf where slopes decrease from 0.04 to 0.01, and the location of these fans is consistent with wave‐ and auto‐suspending sediment gravity current theories. Thus, we provide the first documentation that the morphology of sediment deposits on the continental shelf can be dictated by river‐generated hyperpycnal flows.

River deltas are densely populated regions of the world with vulnerable groundwater reserves. Contamination of these groundwater aquifers via saline water intrusion and pollutant transport is a growing threat due to both anthropogenic and climate changes. The arrangement and composition of subsurface sediment is known to have a significant impact on aquifer contamination; however, developing accurate depictions of the subsurface is challenging. In this work, we explore the relationship between surface and subsurface properties and identify the metrics most sensitive to different forcing conditions. To do so, we simulate river delta evolution with the rule‐based numerical model, DeltaRCM, and test the influence of input sand fraction and steady sea level rise (SLR) on delta evolution. From the model outputs, we measure a variety of surface and subsurface metrics chosen based on their applicability to imagery and modeling results. The Kullback‐Leibler (KL) divergence is then used to quantitatively gauge which metrics are most indicative of the imposed forcings. Both qualitative observations and the KL divergence analysis suggest that estimates of subsurface connectivity can be constrained using surface information. In particular, more variable shoreline roughness values and higher surface wetted fraction values correspond to increased subsurface connectivity. These findings complement traditional methods of estimating subsurface structure in river‐dominated delta systems and represent a step toward the identification of a direct link between surface observations and subsurface form.

We present a geometric, sediment mass-balance model for the interaction of axial and transverse alluvial systems in a subsiding basin. By comparing the model result with a flume experiment that employed a simplified half-graben tectonic geometry with axial and transverse sediment sources, we quantify rates of axial-transverse erosional sediment mixing. In the experiment, the lateral migration rate of the axial-transverse boundaries due to the sediment mixing scales with sediment supplies delivered by transverse drainages, but not with water (or sediment) discharge from the axial channel or with tectonic tilting rate. Using an empirical lateral erosion rate, the model shows how sediment supply partitioning among the axial, hanging-wall, and footwall drainages controls the width and the location of the axial-channel belt. Comparing the modeling results with field cases demonstrates that transverse sediment fluxes could slow the axial-channel migration or even reverse the movement against the tectonic forcing.

Hyperpycnal currents are river‐derived turbidity currents capable of transporting significant volumes of sediment from the shoreline onto the shelf and potentially further to deep ocean basins. However, their capacity to deposit sand bodies on the continental shelf is poorly understood. Shelf hyperpycnites remain an overlooked depositional element in source to sink systems, primarily due to their limited recognition in the rock record. Recent discoveries of modern shelf hyperpycnites, and previous work describing hyperpycnites deposited in slope or deep‐water settings, provide a valuable framework for understanding and recognizing shelf hyperpycnites in the rock record. This article describes well‐sorted lobate sand bodies on the continental shelf of the Neuquén Basin, Argentina, interpreted to have been deposited by hyperpycnal currents. These hyperpycnites of the Jurassic Lajas Formation are characterized by well‐sorted, medium‐grained, parallel‐laminated sandstones with hundreds of metre extensive, decimetre thick beds encased by organic‐rich, thinly laminated sandstone and siltstone. These deposits represent slightly obliquely‐migrating sand lobes fed by small rivers and deposited on the continental shelf. Hyperpycnites of the Lajas Formation highlight several unique characteristics of hyperpycnal deposits, including their distinctively thick horizontal laminae attributed to pulsing of the hyperpycnal currents, the extraction of coarse gravel due to low flow competence, and the extraction of mud due to lofting of light interstitial fluid. Recognition of shelf hyperpycnites in the Lajas Formation of the Neuquén Basin allows for a broader understanding of shelf processes and adds to the developing facies models of hyperpycnites. Recognizing and understanding the geometry and internal architecture of shelf hyperpycnites will improve current understanding of sediment transfer from rivers to deeper water, will improve palaeoenvironmental interpretations of sediment gravity‐flow deposits, and has implications for modelling potentially high‐quality hydrocarbon reservoirs.

Characterization of groundwater aquifers and hydrocarbon reservoirs requires an understanding of the distribution and connectivity of subsurface sandbodies. In deltaic environments, distributary channel networks serve as the primary conduits for water and sediment. Once these networks are buried and translated into the subsurface, the coarse-grained channel fills serve as primary conduits for subsurface fluids such as water, oil or gas. The temporal evolution of channels on the surface therefore plays a first-order role in the 3D permeability and connectivity of subsurface networks. Land surface imagery is more broadly available than topographic or subsurface data, and time-series imagery of river networks can hold useful information for constraining the shallow subsurface.However, these reconstructions require an understanding of the degree to which channel bathymetry and river kinematics affect connectivity of subsurface sandbodies. Here, we present a novel method for building synthetic cross sections using overhead images of an experimental delta. We use principal components analysis to extract river networks from surface imagery, then couple this with an inverse-CDF method to estimate channel bathymetry, to generate a time-series of synthetic delta topography. This synthetic topography is then transformed, accounting for deposition and subsidence, to produce synthetic stratigraphy that differentiates coarse-grained channel fill from overbank and offshore deposition. We find that large-scale subsurface architecture is relatively insensitive to details of channel bathymetry, but instead is primarily controlled by channel location and kinematics. We analyse the connectivity of sand bodies and the geometries of barriers to flow and find that periods of rapid sea-level rise have more variability in sand body connectivity. We also find that barrier width decreases downstream during all sea-level phases. Our method generates synthetic stratigraphy that closely resembles the large-scale architecture and 2-dimensional connectivity of the real stratigraphy built during the experiment it was based