Land motions are significantly widespread in the Yellow River delta (YRD). There is, however, a lack of understanding of the delta-wide comprehensive deformation mode and its dynamic mechanism, especially triggered by groundwater extraction. This paper adopts an integrated analysis of multidisciplinary data of image geodesy, geophysics, geology and hydrogeology to provide insights into Earth surface displacement patterns and dynamics in the YRD. Delta-scale land motions were measured for the first time using L-band ALOS images processed using multi-temporal InSAR, illustrating multiple obvious surface sinking regions and a maximum annual subsidence velocity of up to 130 mm. Then, the InSAR-constrained distributed point source model with optimal kernel parameters, a smoothness factor of 10 and a model grid size of 300 m was established and confirmed to be rational, reliable and accurate for modeling analysis over the YRD. Remarkable horizontal surface displacements, moving towards and converging on a sinking center, were recovered by means of modeling and measured using InSAR, with a maximum rate of up to 60 mm per year, which can trigger significant disasters, such as ground fissures and building damage. In addition, the annual total water storage variation at the delta scale, the most meaningful outcome, can be calculated and reaches a total of approximately 12,010 × 10³ m³ in Guangrao city, efficiently filling the gap of GRACE and in situ investigations for delta-wide aquifer monitoring. Finally, a comparative analysis of time series InSAR measurements, modeling outcomes, and fault and groundwater data was conducted, and the strong agreement demonstrates that faults control aquifer distribution and hence the spatial distribution of groundwater-withdrawal-related regional land subsidence. Moreover, the obvious asymmetric displacements, demonstrating a northeasterly displacement trend, further reveal that faults control aquifer distribution and Earth surface deformation. These findings are useful for understanding the land motion patterns and dynamics, helping to sustainably manage groundwater and control disasters in the YRD and elsewhere worldwide.

Although recent research suggests that the morphology and crustal structure of slow–ultraslow-spreading ridges are mainly controlled by melt supply, there is a lack of quantitative understanding of the effect of systematic changes in melt supply on the seafloor spreading state of mid-ocean ridges. In this study, we used bathymetry, free-air gravity anomaly, and sediment thickness data to calculate the residual bathymetry, mantle Bouguer gravity and crustal thickness of the Reykjanes Ridge. According to the gradient of changes in crustal thickness and residual bathymetry along the axis, the influence of melt supply on the spreading state of the Reykjanes Ridge can be divided into three zones: ultra-strong effect zone (0–160 km), strong effect zone (160–610 km), and weak effect zone (610–930 km). In the ultra-strong effect zone, excess melt supply and a higher melting degree result in a strong upwelling and large melt eruption. The change in relative position between the Reykjanes Ridge and the Iceland hotspot results in the spreading state of the Reykjanes Ridge transforming from asymmetric spreading to symmetric spreading. In the strong effect zone, the decrease in melt supply and melting degree weakens the mantle upwelling and enhances the viscosity of the dehydrated mantle layer. Sufficient viscosity of the dehydrated mantle layer forces asymmetric asthenosphere rise along the sloping boundary of the lithosphere, resulting in symmetric spreading. In the weak effect zone, the pattern of magma upwelling becomes a focused magma supply pattern similar to that of the slow–ultraslow-spreading of the mid-ocean ridge, and tectonics dominate the spreading process. The asymmetry of this weak effect zone may be due to the concentration of tectonic and magmatic activity on one flank of the ridge.