A fundamental long‐standing question regarding Mars history is whether the flat and low‐lying northern plains ever hosted an ocean. The best opportunity to solve this problem is provided by stratigraphic observations of sedimentary deposits onlapping the crustal dichotomy. Here, we use high‐resolution imagery and topography to analyze a branching network of inverted channel and channel lobe deposits in the Aeolis Dorsa region, just north of the dichotomy boundary. Observations of stacked, cross‐cutting channel bodies and stratal geometries indicate that these landforms represent exhumed distributary channel deposits. Observations of depositional trunk feeder channel bodies, a lack of evidence for past topographic confinement, channel avulsions at similar elevations, and the presence of a strong break in dip slope between topset and foreset beds suggest that this distributary system was most likely a delta, rather than an alluvial fan or submarine fan. Sediment transport calculations using both measured and derived channel geometries indicate a minimum delta deposition time on the order of 400 years. The location of this delta within a thick and widespread clastic wedge abutting the crustal dichotomy boundary, unconfined by any observable craters, suggests a standing body of water potentially 10⁵ km²in extent or greater and is spatially consistent with hypotheses for a northern ocean.

Floodplain and terrace features can provide information about current and past fluvial processes, including channel response to varying discharge and sediment flux, sediment storage, and the climatic or tectonic history of a catchment. Previous methods of identifying floodplain and terraces from digital elevation models (DEMs) tend to be semi-automated, requiring the input of independent datasets or manual editing by the user. In this study we present a new method of identifying floodplain and terrace features based on two thresholds: local gradient, and elevation compared to the nearest channel. These thresholds are calculated statistically from the DEM using quantile–quantile plots and do not need to be set manually for each landscape in question. We test our method against field-mapped floodplain initiation points, published flood hazard maps, and digitised terrace surfaces from seven field sites from the US and one field site from the UK. For each site, we use high-resolution DEMs derived from light detection and ranging (lidar) where available, as well as coarser resolution national datasets to test the sensitivity of our method to grid resolution. We find that our method is successful in extracting floodplain and terrace features compared to the field-mapped data from the range of landscapes and grid resolutions tested. The method is most accurate in areas where there is a contrast in slope and elevation between the feature of interest and the surrounding landscape, such as confined valley settings. Our method provides a new tool for rapidly and objectively identifying floodplain and terrace features on a landscape scale, with applications including flood risk mapping, reconstruction of landscape evolution, and quantification of sediment storage and routing.

The mass of carbon stored as organic matter in terrestrial systems is sufficiently large to play an important role in the global biogeochemical cycling of CO₂ and O₂. Field measurements of radiocarbon-depleted particulate organic carbon (POC) in rivers suggest that terrestrial organic matter persists in surface environments over millennial (or greater) timescales, but the exact mechanisms behind these long storage times remain poorly understood. To address this knowledge gap, we developed a numerical model for the radiocarbon content of riverine POC that accounts for both the duration of sediment storage in river deposits and the effects of POC cycling. We specifically target rivers because sediment transport influences the maximum amount of time organic matter can persist in the terrestrial realm and river catchment areas are large relative to the spatial scale of variability in biogeochemical processes. Our results show that rivers preferentially erode young deposits, which, at steady state, requires that the oldest river deposits are stored for longer than expected for a well-mixed sedimentary reservoir. This geometric relationship can be described by an exponentially tempered power-law distribution of sediment storage durations, which allows for significant aging of biospheric POC. While OC cycling partially limits the effects of sediment storage, the consistency between our model predictions and a compilation of field data highlights the important role of storage in setting the radiocarbon content of riverine POC. The results of this study imply that the controls on the terrestrial OC cycle are not limited to the factors that affect rates of primary productivity and respiration but also include the dynamics of terrestrial sedimentary systems.

Bedrock river valleys are fundamental components of many landscapes, and their morphologies—from slot canyons with incised meanders to wide valleys with strath terraces—may record environmental history. Several formation mechanisms for particular valley types have been proposed that involve changes in climatic and tectonic forcing, but the uniqueness of valley evolution pathways and the long‐term stability of valley morphology under constant forcing are unknown and are not predicted in existing numerical models for vertically incising rivers. Because rivers often migrate more rapidly through alluvium than through bedrock, we explore the hypothesis that the distribution of bank materials strongly influences river meandering kinematics and can explain the diversity of bedrock river valley morphology. Simulations using a numerical model of river meandering with vector‐based bank‐material tracking indicate that channel lateral erosion rate in sediment and bedrock, vertical erosion rate, and initial alluvial‐belt width explain first‐order differences in bedrock valley type; that bedrock‐bound channels can evolve under steady forcing from alluvial states; and that weak bedrock and low vertical incision rates favor wide, shallow valleys, while resistant bedrock and high vertical incision rates favor narrow, deep valleys. During vertical incision, sustained planation of the valley floor is favored when bedrock boundaries restrict channel migration to a zone of thin sediment fill. The inherent unsteadiness of river meandering in space and time is enhanced by evolving spatial contrasts in bank strength between sediment and bedrock and can account for several valley features—including strath terraces and underfit valleys—commonly ascribed to external drivers.