After the formation of the bend ice cover, the ice thickness of the bend is not uniformly distributed, and an open-water area is usually formed downstream of the bend. The spatial and temporal variation of the ice thickness in seven cross sections was determined via Unmanned Aerial Vehicle Ground Penetrating Radar (UAV-GPR) technology and traditional borehole measurements. The plane morphology change of the open water was observed by Sentinel-2. The results show that the average dielectric permittivity of GPR was 3.231, 3.249, and 3.317 on three surveys (5 January 2022, 16 February 2022, and 25 February 2022) of the Yellow River ice growing period, respectively. The average ice thickness of the three surveys was 0.402 m, 0.509 m, and 0.633 m, respectively. The ice thickness of the concave bank was larger than that of the convex bank. The plane morphology of the open water first shrinks rapidly longitudinally and then shrinks slowly transversely. The vertical boundary of the open water was composed of two arcs, in which the slope of Arc I (close to the water surface) was steeper than that of Arc II, and the hazardous distance of the open-water boundary was 10.3 m. The increased flow mostly affected the slope change of Arc I. Finally, we discuss the variation of hummocky ice and flat ice in GPR images and the physical factors affecting GPR detection accuracy, as well as the ice-thickness variation of concave and convex banks in relation to channel curvature.

Unfrozen free and non-free water between ice crystals in flat and hummock ice in the Yellow River exists as water films with varying contents based on ice temperature. These contents can affect the radar wave velocity of the ice despite its theoretical dependence on the crystal structure and ice body components. The unfrozen water content in ice depends on the ice temperature, which is controlled by the air temperature, solar radiation, and ice thickness. Winter air temperature and radar-detected ice thickness data observed at the Shisifenzi bend in the Yellow River from 2020 to 2021 were analyzed. The unfrozen water content in the ice was the primary factor influencing the accuracy of flat ice thickness detection. The heat flux at the ice–water interface in the Yellow River was determined. The evolution of ice thickness and temperature were simulated using a one-dimensional (1D) ice thermodynamic model forced by the local weather station data (i.e., air temperature, solar radiation, wind speed, and cloud cover). On this basis, the measured ice thickness data of 13 drill holes were combined to calculate 1251 thermodynamically simulated ice thicknesses consistent with the ice thickness detection time of the radar; therefore, statistical relationships regarding the influence of air temperature and the combined action of air temperature and ice thickness on the radar wave velocity in granular and columnar ice during air temperature increases and decreases were determined. Finally, the statistical relationship between the combined influence of air temperature and ice thickness on radar wave velocity was selected as a parameterization scheme to dynamically correct the radar wave velocity of flat ice. To enhance the radar detection accuracy for flat ice thickness, the radar wave velocity of ice was parameterized as a function. Given the presence of unfrozen frazil ice and accumulated broken ice blocks in the Yellow River, radar is suggested to detect the thickness of different types of ice in future research.

Acquisition of continuous drift ice characteristic parameters such as ice size, shape, concentration, and drift velocity are of great importance for understanding river freezing and thawing processes. This study acquired hourly oblique images captured by a shore-based camera in the winter of 2021–2022 on the Yellow River, China. The pixel point scale method for correcting oblique images is provided. The 61 lines were measured at the calibration site and the absolute error between the measured value and the calculated value was in the range of 0.009–0.850 m, with a mean error of 0.236 m. After the correction of oblique images, the true equivalent diameter of drift ice during the freezing period ranged from 0.52–13.10 m with a mean size of 3.36 m, which was larger than that of 2.30 m during the thawing period which ranged from 0.20–12.54 m. It was found that the size of drift ice increased with time during the freezing period and decreased with time during the thawing period. The fractal dimension and roundness were used to represent drift ice shape. The fractal dimension ranged from 1.0–1.3 and the roundness ranged from 0.1–1.0. A Gaussian distribution was used to estimate drift ice size and shape distributions. There is a nonlinear relationship between ice concentration and drift velocity, which can be well expressed by the logistic function. In the future, drift ice parameters for more years and hydrometeorological data for the same time need to be accumulated, which helps to analyze the freezing and thawing patterns of river ice.