Regression tree models have been widely used for remote sensing-based ecosystem mapping. Improper use of the sample data (model training and testing data) may cause overfitting and underfitting effects in the model. The goal of this study is to develop an optimal sampling data usage strategy for any dataset and identify an appropriate number of rules in the regression tree model that will improve its accuracy and robustness. Landsat 8 data and Moderate-Resolution Imaging Spectroradiometer-scaled Normalized Difference Vegetation Index (NDVI) were used to develop regression tree models. A Python procedure was designed to generate random replications of model parameter options across a range of model development data sizes and rule number constraints. The mean absolute difference (MAD) between the predicted and actual NDVI (scaled NDVI, value from 0–200) and its variability across the different randomized replications were calculated to assess the accuracy and stability of the models. In our case study, a six-rule regression tree model developed from 80% of the sample data had the lowest MAD (MAD_training = 2.5 and MAD_testing = 2.4), which was suggested as the optimal model. This study demonstrates how the training data and rule number selections impact model accuracy and provides important guidance for future remote-sensing-based ecosystem modeling.

The U.S. Geological Survey’s Land Change Monitoring, Assessment, and Projection (LCMAP) initiative involves detecting changes in land cover, use, and condition with the goal of producing land change information to improve the understanding of the Earth system and provide insights on the impacts of land surface change on society. The change detection method ingests all available high-quality data from the Landsat archive in a time series approach to identify the timing and location of land surface change. Annual thematic land cover maps are then produced by classifying time series models. In this paper, we describe the optimization of the classification method used to derive the thematic land cover product. We investigated the influences of auxiliary data, sample size, and training from different sources such as the U.S. Geological Survey’s Land Cover Trends project and National Land Cover Database (NLCD 2001 and NLCD 2011). The results were evaluated and validated based on independent data from the training dataset. We found that refining the auxiliary data effectively reduced artifacts in the thematic land cover map that are related to data availability. We improved the classification accuracy and stability considerably by using a total of 20 million training pixels with a minimum of 600,000 and a maximum of 8 million training pixels per class within geographic windows consisting of nine Analysis Ready Data tiles (450 by 450 km²). Comparisons revealed that the NLCD 2001 training data delivered the best classification accuracy. Compared to the original LCMAP classification strategy used for early evaluation (e.g., Trends training data, 20,000 samples), the optimized classification strategy improved the annual land cover map accuracy by an average of 10%.

This paper presents the methodology and results of two ecological-based net ecosystem production (NEP) regression tree models capable of up scaling measurements made at various flux tower sites throughout the U.S. Great Plains. Separate grassland and cropland NEP regression tree models were trained using various remote sensing data and other biogeophysical data, along with 15 flux towers contributing to the grassland model and 15 flux towers for the cropland model. The models yielded weekly mean daily grassland and cropland NEP maps of the U.S. Great Plains at 250 m resolution for 2000–2008. The grassland and cropland NEP maps were spatially summarized and statistically compared. The results of this study indicate that grassland and cropland ecosystems generally performed as weak net carbon (C) sinks, absorbing more C from the atmosphere than they released from 2000 to 2008. Grasslands demonstrated higher carbon sink potential (139 g C·m⁻²·year⁻¹) than non-irrigated croplands. A closer look into the weekly time series reveals the C fluctuation through time and space for each land cover type.

Rapid urbanization is reshaping thermal environments worldwide, with the strongest impacts occurring at the interface between urban and non-urban areas. Impervious surfaces, as key indicators of urban expansion, are critical for monitoring urban growth and assessing surface urban heat island (SUHI) effects. Land use and land cover change (LULCC) provides an essential link between urban dynamics and their environmental and societal consequences. Here, we integrated the U.S. Geological Survey (USGS) Climate Global Issues (CGI) Land Cover Product with Landsat thermal time-series to investigate SUHI evolution in two contrasting metropolitan regions: Wuhan, China, and Brasília, Brazil. Using data spanning 1986–2023, we analyzed the relationships between land cover, Landsat-based land surface temperature (LST), and SUHI intensity, and identified persistent thermal hotspots. Results demonstrate that the land cover data utilized increases the accuracy of impervious surface mapping along urban–rural gradients. Average SUHI intensities were 3.4 °C in Wuhan and 3.3 °C in Brasília, with statistically significant warming trends of 0.04 °C/year and 0.01 °C/year, respectively. Maximum temperature proved to be a robust indicator of SUHI intensification, capturing long-term upward trends. Our findings highlight the important role of urban land cover dynamics in shaping temporal SUHI variability and hotspot emergence. This prototype framework demonstrates the scientific and policy value of combining long-term land cover monitoring information with satellite thermal monitoring to quantify and track SUHI at city scale, supporting sustainable urban planning and climate adaptation strategies.