Accurate measures of human effects on landscape processes require consideration of both the direct impacts from human activities and the indirect consequences of the interactions between humans and the landscape. This is particularly evident in systems experiencing regular natural disturbances such as in the mountainous areas of southwestern China, where the remaining population of giant pandas (Ailuropoda melanoleuca) is supported. Here the spatiotemporal patterns of human impacts, forests, and bamboo episodic die‐offs combine to determine the distribution of panda habitat. To study the complex interactions of humans and landscapes, we developed an integrated spatiotemporally explicit model of household activities, natural vegetation dynamics, and their impacts on panda habitat. Using this model we examined the direct consequences of local fuelwood collection and household creation on areas of critical giant panda habitat and the indirect impacts when coupled with vegetation dynamics. Through simulations, we found that over the next 30 years household impacts would result in the loss of up to 30% of the habitat relied on by pandas during past bamboo die‐offs. The accumulation and spatial distribution of household impacts would also have a considerable indirect influence on the spatial distribution of understory bamboo. While human impacts influence both bamboo die‐off and regeneration, over 19% of pre‐existing low‐elevation bamboo habitat may be lost following an episodic die‐off depending on the severity of the impacts and timing of the die‐offs. Our study showed not only the importance of the spatial distribution of direct household impacts on habitat, but also the far‐reaching effects of the indirect interactions between humans and the landscapes they are modifying.

Accurate digital mapping of soil organic carbon (SOC) is important in understanding the global carbon cycle and its implications in mitigating climate change. Visible and near-infrared hyperspectral imaging technology provides an alternative for mapping SOC efficiently and accurately, especially at regional and global scales. However, there is a lack of understanding of the impacts of spatial resolution of hyperspectral images and spatial autocorrelation of spectral information on the accuracy of SOC retrievals. In this study, the hyperspectral images (380–1700 nm) with a spatial resolution of 1 m were acquired by Headwall Micro-Hyperspec airborne sensors. Then, hyperspectral images were resampled into three different spatial resolutions of 10 m, 30 m, and 60 m by near neighbor (NN), bilinear interpolation (BI), and cubic convolution (CC) resampling methods. The geographically weighted regression (GWR) model was used to explore the role of spatial autocorrelation in predicting SOC contrast with the partial least squares regression (PLSR) model. Results showed that (1) the hyperspectral images can be used to predict SOC and the spatial autocorrelation can improve the prediction accuracy, as the ratio of performance to interquartile range (RPIQ) values of PLSR and GWR were 1.957 and 2.003; (2) The SOC prediction accuracy decreased with the degradation of spatial resolution, and the RPIQ values of PLSR were from 1.957 to 1.134, and of GWR were from 2.003 to 1.136; (3) Three resampling methods had a much weaker influence than spatial resolution on SOC predictions because the differences of RPIQ values of NN, BI, and CC resampling methods were 0.146, 0.175, and 0.025 in the spatial resolutions of 10 m, 30 m, and 60 m, respectively; (4) Finally, the Global Moran’s I and the Anselin Local Moran’s I proved the existence of the spatial autocorrelation in SOC maps. We hope that this study can offer valuable information for digital soil mapping by satellite hyperspectral images in the near future.

The rapid monitoring and accurate estimation of dynamic changes in soil organic carbon (SOC) can make great efforts in understanding the global carbon cycle. Traditional field survey is the main approach to obtain soil data and measure SOC content. However, the limited number of soil samples and the sampling cost hinder the quality of digital soil mapping. This research aims to explore the sensitive of sampling density in digital soil mapping, and then design a suitable soil sampling plan based on a series of sampling indices. Headwall hyperspectral images (400–1700 nm) were used to estimate the SOC map by partial least squares regression (PLSR) and PLSR kriging (PLSRK). Three traditional soil sampling methods (random, grid, and Latin hypercube sampling) with 10 classes of sampling densities (6.26, 2.79, 1.57, 1.01, 0.69, 0.53, 0.39, 0.30, 0.26, and 0.20 ha⁻¹) were designed. The R², root mean square error (RMSE) and ratio of standard deviation to RMSE (RPD) were used to evaluate the prediction accuracy in digital soil mapping by ordinary kriging. Three new indices, namely, the ratio of sampling efficiency to performance (RSEP), the density of soil samples index and the comprehensive evaluation index of prediction accuracy, were used to select a suitable soil sampling plan. Results showed that (1) the prediction accuracy of PLSRK (RPD = 2.00) was higher by approximately 11.73% than that of PLSR (RPD = 1.79), and the hyperspectral images provided an actual referential SOC map for the study of soil sampling; (2) the grid sampling plan performed better than the random and Latin hypercube sampling methods, and the quality of SOC map improves with the increase of the sampling density, and (3) the computer simulation and field verification indicated that RSEP is one feasible index in designing a suitable soil sampling plan.

The global coverage of Chlorophyll-a concentration (Chl-a) has been continuously available from ocean color satellite sensors since September 1997 and the Chl-a data (1997–2019) were used to produce a climatological dataset by averaging Chl-a values at same locations and same day of year. The constructed climatology can remarkably reduce the variability of satellite data and clearly exhibit the seasonal cycles, demonstrating that the growth and decay of phytoplankton recurs with similarly seasonal cycles year after year. As the shapes of time series of the climatology exhibit strong periodical change, we wonder whether the seasonality of Chl-a can be expressed by a mathematic equation. Our results show that sinusoid functions are suitable to describe cyclical variations of data in time series and patterns of the daily climatology can be matched by sine equations with parameters of mean, amplitude, phase, and frequency. Three types of sine equations were used to match the climatological Chl-a with Mean Relative Differences (MRD) of 7.1%, 4.5%, and 3.3%, respectively. The sine equation with four sinusoids can modulate the shapes of the fitted values to match various patterns of climatology with small MRD values (less than 5%) in about 90% of global oceans. The fitted values can reflect an overall pattern of seasonal cycles of Chl-a which can be taken as a time series of biomass baseline for describing the state of seasonal variations of phytoplankton. The amplitude images, the spatial patterns of seasonal variations of phytoplankton, can be used to identify the transition zone chlorophyll fronts. The timing of phytoplankton blooms is identified by the biggest peak of the fitted values and used to classify oceans as different bloom seasons, indicating that blooms occur in all four seasons with regional features. In global oceans within latitude domains (48°N–48°S), blooms occupy approximately half of the ocean (50.6%) during boreal winter (December–February) in the northern hemisphere and more than half (58.0%) during austral winter (June–August) in the southern hemisphere. Therefore, the sine equation can be used to match the daily Chl-a climatology and the fitted values can reflect the seasonal cycles of phytoplankton, which can be used to investigate the underlying phenological characteristics.

The severity of wildfires is increasing and has driven increases in nighttime fire activity. Enhanced capability to detect the active burning regions of wildfires at night could significantly improve the effectiveness of wildfire management operations. Potassium line emission in the NIR near 770 nm is a signature of active burning. We test the use of multi-band imaging from an aircraft at night to distinguish a wood-burning fire from artificial light sources. We find that a simple ratio of the signals in two broad bands, one including 770 nm, effectively discriminates the fire from artificial light sources. This offers the possibility of nighttime fire detection with high spatial resolution using silicon sensors sensitive in the NIR.

Harmful algal blooms (HABs) are one of the major environmental concerns, as they have various negative effects on public and environmental health, recreational services, and economics. HAB modeling is challenging due to inconsistent and insufficient data, as well as the nonlinear nature of algae formation data. However, it is crucial for attaining sustainable development goals related to clean water and sanitation. From this point of view, we employed the sparse identification nonlinear dynamics (SINDy) technique to model microcystin, an algal toxin, utilizing dissolved oxygen as a water quality metric and evaporation as a meteorological parameter. SINDy is a novel approach that combines a sparse regression and machine learning method to reconstruct the analytical representation of a dynamical system. The model results indicate that MAPE values of approximately 2% were achieved in three out of four lakes, while the MAPE value of the remaining lake is 11%. Moreover, a model-driven and web-based interactive tool was created to develop environmental education, raise public awareness on HAB events, and produce more effective solutions to HAB problems through what-if scenarios. This interactive and user-friendly web platform allows tracking the status of HABs in lakes and observing the impact of specific parameters on harmful algae formation.