Search

This paper presents the results of the Ensemble Riemannian Data Assimilation for relatively high-dimensional nonlinear dynamical systems, focusing on the chaotic Lorenz-96 model and a two-layer quasi-geostrophic (QG) model of atmospheric circulation. The analysis state in this approach is inferred from a joint distribution that optimally couples the background probability distribution and the likelihood function, enabling formal treatment of systematic biases without any Gaussian assumptions. Despite the risk of the curse of dimensionality in the computation of the coupling distribution, comparisons with the classic implementation of the particle filter and the stochastic ensemble Kalman filter demonstrate that with the same ensemble size, the presented methodology could improve the predictability of dynamical systems. In particular, under systematic errors, the root mean squared error of the analysis state can be reduced by 20 % (30 %) in Lorenz-96 (QG) model.

Bedload particle hops are defined as successive motions of a particle from start to stop, characterizing one of the most fundamental processes describing bedload sediment transport in rivers. Although two transport regimes have been recently identified for short- and long-hops, respectively (Wu et al., Water Resour Res, 2020), there still lacks a theory explaining how the mean hop distance-travel time scaling may extend to cover the phenomenology of bedload particle motions. Here we propose a velocity-variation based formulation, and for the first time, we obtain analytical solution for the mean hop distance-travel time relation valid for the entire range of travel times, which agrees well with the measured data (Wu et al., J Fluid Mech, 2021). Regarding travel times, we identify three distinct regimes in terms of different scaling exponents: respectively as ~1.5 for an initial regime and ~5/3 for a transition regime, which define the short-hops; and 1 for the so-called Taylor dispersion regime defining long-hops. The corresponding probability density function of the hop distance is also analytically obtained and experimentally verified.

The channelization cascade observed in terrestrial landscapes describes the progressive formation of large channels from smaller ones starting from diffusion-dominated hillslopes. This behavior is reminiscent of other non-equilibrium complex systems, particularly fluids turbulence, where larger vortices break down into smaller ones until viscous dissipation dominates. Based on this analogy, we show that topographic surfaces emerging between parallel zero-elevation boundaries present a logarithmic scaling in the mean-elevation profile, which resembles the well-known logarithmic velocity profile in wall-bounded turbulence. Within this region of elevation fluctuation, the power spectrum exhibits a power-law decay resembling the Kolmogorov -5/3 scaling of turbulence. We also demonstrate that similar scaling behaviors emerge in surfaces from a laboratory experiment, natural basins, and constructed following optimality principles. In general, we show that the steady-state solutions of the governing equations of landscape evolution are the stationary surfaces of a functional defined as the average domain elevation. Depending on the exponent of the specific drainage area in the erosion term (m), the steady-state surfaces are local minimum (m1) of the average domain elevation.