A prediction–correction solver is presented here for rapid simulation of free-surface flows in dendritic and looped river networks. Rather than solving a large global algebraic system over the entire domain of river networks, the model solves subsystems for subdomains of branches in two steps: prediction and correction. With the help of the prediction–correction method (PCM), the partial linearization technique, the semi-implicit method, and the Eulerian–Lagrangian method (ELM), the model only needs to solve tridiagonal linear systems for branches and is free of any iteration. The new model was tested using a hypothetical looped river network with regular cross sections, the Three Gorges Reservoir (TGR) dendritic river network, and the Jing-south looped river system (with seasonally flooding branches). In the first test, a time-step sensitivity study was conducted and the model was revealed to produce accurate simulations at large time steps when the condition for application of the PCM to river networks was satisfied. In the TGR test, the PCM model provided almost the same histories of water levels and discharges as those simulated by the HEC-RAS model. In the Jing-south test, the mean absolute error in simulated water levels was 0.07–0.24 m, and the relative error in simulated cross-section water flux was 0.5–4.9% compared with field data; the conservation error was generally 2 × 10⁻⁴ to 3 × 10⁻⁴. The PCM model was revealed to be 2–4 times as fast as a reported model, which solves local nonlinear subsystems using two-layer iterations, and 1.2–1.4 times as fast as the HEC-RAS. Using a time step of 1200 s, it took the sequential code 26.8 and 23.1 s to complete a simulation of a one-year unsteady flow process, respectively, in the TGR river networks (with 588 cells) and the Jing-south river system (with 662 cells).

A semi–implicit numerical model with a linear solver is proposed for the free-surface and pressurized mixed flows in hydraulic systems. It solves the two flow regimes within a unified formulation, and is much simpler than existing similar models for mixed flows. Using a local linearization and an Eulerian–Lagrangian method, the new model only needs to solve a tridiagonal linear system (arising from velocity-pressure coupling) and is free of iterations. The model is tested using various types of mixed flows, where the simulation results agree with analytical solutions, experiment data and the results reported by former researchers. Sensitivity studies of grid scales and time steps are both performed, where a common grid scale provides grid-independent results and a common time step provides time-step-independent results. Moreover, the model is revealed to achieve stable and accurate simulations at large time steps for which the CFL is greater than 1. In simulations of a challenging case (mixed flows characterized by frequent flow-regime conversions and a closed pipe with wide-top cross-sections), an artificial slot (A-slot) technique is proposed to cope with possible instabilities related to the discontinuous main-diagonal coefficients of the linear system. In this test, a slot-width sensitivity study is also performed, and the suitable slot-width ratio (ε) for the linear solver is suggested to be 0.05–0.1.

Because of special morphologies and complex runoff–tide interactions, the landward floodtide flows in Yangtze Estuary are observed to spill over from the North to the South Branches, carrying a lot of sediment. To quantitatively clarify the spillover problem, a two-dimensional numerical model using a high-resolution channel-refined unstructured grid is developed for the entire Yangtze Estuary from Datong to river mouths (620 km) and part of the East Sea. The developed model ensures a good description of the river-coast-ocean coupling, the irregular boundaries, and local river regimes in the Yangtze Estuary. In tests, the simulated histories of the tidal level, depth-averaged velocity, and sediment concentration agree well with field data. The spillover of sediment in the Yangtze Estuary is studied using the condition of a spring and a neap tide in dry seasons. For a representative cross-section in the upper reach of the North Branch (QLG), the difference of the cross-sectional sediment flux (CSSF) between floodtide and ebbtide durations is 43.85–11.26 × 10⁴ t/day, accounting for 37.5–34.9% of the landward floodtide CSSF. The mechanics of sediment spillover in Yangtze Estuary are clarified in terms of a successive process comprising the source, transport, and drainage of the spillover sediment.

Responses of irregular rotational tidal flows to an outlet regulation (the Guyuan Sand (GYS) regulation) in the three-level branching Yangtze Estuary are studied by a high-resolution numerical model and theoretical analysis. The project is launched around GYS at the outlet of the North Branch of the Yangtze Estuary. The tidal flows around GYS are rotational and become irregular under the influences of the runoff-tide interactions, rapidly varying topographies and complex solid boundaries in coastal areas. Three designs of GYS regulation were studied, including various diversion dikes and new outlets of different widths. The regulation disturbs the irregular rotational flows around GYS, and further changes the estuarine tidal processes and the water exchange between different branches of the branching Yangtze Estuary. It was interesting to find that additional current and additional storage are formed along the North Branch when a southward outlet and the clockwise rotational flow met around GYS. This special phenomenon is named “guide effect” in this study. The guide effect, together with common resist effect (arising from the narrowed outlet channel), reshapes the estuarine tidal processes. Based on the simulation result and a theoretical analysis, response mechanics of irregular rotational tidal flows to the outlet regulation in complex branching estuaries are quantitatively studied.

Dam construction often changes downstream fluvial processes by reducing sediment supply. Taking Tiaoguan reach and Laijiapu reach of the lower Jingjiang Reach downstream of the Three Gorges Project as examples, three-dimensional flow velocity, sediment, and bed elevation were observed in the two bends for investigating the impact of flow structure and sediment transport of different discharges on sharp bend morphology. Results indicated that the flow structure and sediment transport process in curved channels depended upon the flow stages, which affected the patterns of erosion and deposition along the point bars and concave banks. Flow separation and development of secondary flow were depended on the shapes of point bars and flow depths nearby, and the strength of secondary flow increased with flow discharge. The high flow discharges, which had high sediment carrying capacity and stream power, provided the main driving force for erosion on upstream point bar, thus the type and duration of floods were crucial factors in the morphological evolution of meandering bends. The reduction of sediment supply should be responsible for erosion on the point bars, causing the flow to migrate toward the convex banks. In meandering rivers with reduced sediment supply, retreats (push inward) of inner (convex) banks dominated advances (pull inward) of outer (concave) banks. In addition, the formation and development of concave-bank bars might relate particularly to meander curvature. This study is expected to constitute a reference for bank protection and river management in meandering bends downstream of reservoirs.

Due to the inherent limitations in remote sensing image quality, seasonal variations, and radiometric inconsistencies, river extraction based on remote sensing image classification often results in omissions. These challenges are particularly pronounced in the detection of narrow and complex river networks, where fine river features are frequently underrepresented, leading to fragmented and discontinuous water body extraction. To address these issues and enhance both the completeness and accuracy of fine river identification, this study proposes an advanced fine river extraction and optimization method. Firstly, a linear river feature enhancement algorithm for preliminary optimization is introduced, which combines Frangi filtering with an improved GA-OTSU segmentation technique. By thoroughly analyzing the global features of high-resolution remote sensing images, Frangi filtering is employed to enhance the river linear characteristics. Subsequently, the improved GA-OTSU thresholding algorithm is applied for feature segmentation, yielding the initial results. In the next stage, to preserve the original river topology and ensure stripe continuity, a river skeleton refinement algorithm is utilized to retain critical skeletal information about the river networks. Following this, river endpoints are identified using a connectivity domain labeling algorithm, and the bounding rectangles of potential disconnected regions are delineated. To address discontinuities, river endpoints are shifted and reconnected based on structural similarity index (SSIM) metrics, effectively bridging gaps in the river network. Finally, nonlinear water optimization combined K-means clustering segmentation, topology and spectral inspection, and small-area removal are designed to supplement some missed water bodies and remove some non-water bodies. Experimental results demonstrate that the proposed method significantly improves the regularization and completeness of river extraction, particularly in cases of fine, narrow, and discontinuous river features. The approach ensures more reliable and consistent river delineation, making the extracted results more robust and applicable for practical hydrological and environmental analyses.