The rapid development of urbanization and population growth in China has posed a major threat to the green sustainable development of the ecological environment. However, the impact of urbanization on the eco-environmental quality (EEQ) in China remains to be developed. Understanding their interactive coupling mechanism is of great significance to achieve the urban sustainable development goals. By using multi-source remote sensing data and the coupling coordination degree model (CCDM), we intended to answer the question “What are the temporal and spatial characteristics of urbanization and EEQ in China on the pixel scale during 2000–2013, and what is the coupling mechanism between the urbanization and the EEQ?”. To answer these questions, we explored the coupling mechanism between urbanization and the EEQ in China with a combined mathematical and graphics model. The results show that the urbanization and the coupling coordination degree (CCD) of the whole region continually increased from 2000 to 2013, especially in the three major urban agglomerations, with a spatial distribution pattern that was “high in the east and low in the west”. Most importantly, from 2000 to 2013, the CCD type of cities in China gradually evolved from uncoordinated cities to coordinated cities. Additionally, the decisive factor affecting the CCD from 2000 to 2013 was the development of urbanization, and the degree at which urbanization had an impact on CCD was about 8.4 times larger than that of the EEQ. At the same time, the rapid urbanization that has occurred in some areas has led to a significant decline in the EEQ, thus indicating that China needs to increase its protection of the ecological environment while pursuing social and economic development in the future. This study makes up for the deficiencies in the existing literature and investigates the long-term coupling of the EEQ and urbanization in China, thereby providing a new research perspective for the sustainable development of China and even the world in the future.

Parameterization schemes (bulk formulae) are widely used to estimate all-sky surface downward longwave radiation (SDLR) due to the simple, readily available inputs and acceptable accuracy from local to regional scales. Seven widely used bulk formulae are evaluated using the ground measurements collected from 44 globally distributed flux measurement sites of five networks. The Bayesian model averaging (BMA) method is introduced to integrate multiple bulk formulae to obtain an estimate of cloudy-sky SDLR for the first time. The second multiple linear regression model of Carmona et al. (2014) performs the best, with BIAS, RMSE, and R² of zero, 20.13 W·m⁻² and 0.87, respectively. The BMA method can achieve balanced results that are close to the accuracy of the second multiple linear regression model of Carmona et al. (2014) and better than the average accuracy of seven bulk formulae, with BIAS, RMSE, and R² of −1.08 W·m⁻², 21.99 W·m⁻² and 0.87, respectively. In addition, the bulk formula of Crawford and Duchon (1999) is preferred if there is insufficient data to calibrate the bulk formulae because it does not need local calibration and has an acceptable accuracy, with BIAS, RMSE, and R² of 0.96 W·m⁻², 26.58 W·m⁻² and 0.82, respectively. The effects of climate type, land cover type, and surface elevation are also investigated to fully assess the applicability of each bulk formula and BMA. In general, there is no cloudy-sky bulk parametrization scheme that can be successfully applied everywhere.

Anomalies displaying significant deviations between terrain elevation products acquired from spaceborne full-waveform LiDAR and reference elevations are frequently observed in assessment studies. While the predominant focus is on “normal” data, recognizing anomalies within datasets obtained from the Geoscience Laser Altimeter System (GLAS) and the Global Ecosystem Dynamics Investigation (GEDI) is essential for a comprehensive understanding of widely used spaceborne full-waveform data, which not only facilitates optimal data utilization but also enhances the exploration of potential applications. Nevertheless, our comprehension of anomalies remains limited as they have received scant specific attention. Diverging from prevalent practices of directly eliminating outliers, we conducted a targeted exploration of anomalies in forested areas using both transmitted and return waveforms from the GLAS and the GEDI in conjunction with airborne LiDAR point cloud data. We unveiled that elevation anomalies stem not from the transmitted pulses or product algorithms, but rather from scattering sources. We further observed similarities between the GLAS and the GEDI despite their considerable disparities in sensor parameters, with the waveforms characterized by a low signal-to-noise ratio and a near exponential decay in return energy; specifically, return signals of anomalies originated from clouds rather than the land surface. This discovery underscores the potential of deriving cloud-top height from spaceborne full-waveform LiDAR missions, particularly the GEDI, suggesting promising prospects for applying GEDI data in atmospheric science—an area that has received scant attention thus far. To mitigate the impact of abnormal return waveforms on diverse land surface studies, we strongly recommend incorporating spaceborne LiDAR-offered terrain elevation in data filtering by establishing an elevation-difference threshold against a reference elevation. This is especially vital for studies concerning forest parameters due to potential cloud interference, yet a consensus has not been reached within the community.