Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatiotemporal patterns”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

Stable localized patterns in thin liquid films

We study a 2-D nonlinear evolution equation which describes the 3-D spatiotemporal behavior of the air-liquid interface of a thin liquid film lying on the underside of a cooled horizontal plate. We show that the Marangoni effect can stabilize the destabilizing effect of gravity (the Rayleigh-Taylor instability) allowing for the existence of stable localized axisymmetric solutions for a wide range of parameter values. Various properties of these structures are discussed.

Deissler, Robert J.↗

Stable localized patterns in thin liquid films

A two-dimensional nonlinear evolution equation is studied which describes the three-dimensional spatiotemporal behavior of the air-liquid interface of a thin liquid film lying on the underside of a cooled horizontal plate. It is shown that the equation has a Liapunov functional, and this fact is exploited to demonstrate that the Marangoni effect can stabilize the destabilizing effect of gravity (the Rayleigh-Taylor instability), allowing for the existence of stable localized axisymmetric solutions for a wide range of parameter values. Various properties of these structures are discussed.

Deissler, Robert J.↗

Vestibular convergence patterns in vestibular nuclei neurons of alert primates

Sensory signal convergence is a fundamental and important aspect of brain function. Such convergence may often involve complex multidimensional interactions as those proposed for the processing of otolith and semicircular canal (SCC) information for the detection of translational head movements and the effective discrimination from physically congruent gravity signals. In the present study, we have examined the responses of primate rostral vestibular nuclei (VN) neurons that do not exhibit any eye movement-related activity using 0.5-Hz translational and three-dimensional (3D) rotational motion. Three distinct neural populations were identified. Approximately one-fourth of the cells exclusively encoded rotational movements (canal-only neurons) and were unresponsive to translation. The canal-only central neurons encoded head rotation in SCC coordinates, exhibited little orthogonal canal convergence, and were characterized with significantly higher sensitivities to rotation as compared to primary SCC afferents. Another fourth of the neurons modulated their firing rates during translation (otolith-only cells). During rotations, these neurons only responded when the axis of rotation was earth-horizontal and the head was changing orientation relative to gravity. The remaining one-half of VN neurons were sensitive to both rotations and translations (otolith + canal neurons). Unlike primary otolith afferents, however, central neurons often exhibited significant spatiotemporal (noncosine) tuning properties and a wide variety of response dynamics to translation. To characterize the pattern of SCC inputs to otolith + canal neurons, their rotational maximum sensitivity vectors were computed using exclusively responses during earth-vertical axis rotations (EVA). Maximum sensitivity vectors were distributed throughout the 3D space, suggesting strong convergence from multiple SCCs. These neurons were also tested with earth-horizontal axis rotations (EHA), which would activate both vertical canals and otolith organs. However, the recorded responses could not be predicted from a linear combination of EVA rotational and translational responses. In contrast, one-third of the neurons responded similarly during EVA and EHA rotations, although a significant response modulation was present during translation. Thus this subpopulation of otolith + canal cells, which included neurons with either high- or low-pass dynamics to translation, appear to selectively ignore the component of otolith-selective activation that is due to changes in the orientation of the head relative to gravity. Thus contrary to primary otolith afferents and otolith-only central neurons that respond equivalently to tilts relative to gravity and translational movements, approximately one-third of the otolith + canal cells seem to encode a true estimate of the translational component of the imposed passive head and body movement.

Non-NASA Center↗

Abrupt Reduction in Shipping Emission as an Inadvertent Geoengineering Termination Shock Produces Substantial Radiative Warming

Human activities affect the Earth’s climate through modifying the composition of the atmosphere, which then creates radiative forcing that drives climate change. The warming effect of anthropogenic greenhouse gases has been partially balanced by the cooling effect of anthropogenic aerosols. In 2020, fuel regulations abruptly reduced the emission of sulfur dioxide from international shipping by about 80% and created an inadvertent geoengineering termination shock with global impact. Here we estimate the regulation leads to a radiative forcing of +0.2±0.11Wm-2 averaged over the global ocean.The amount of radiative forcing could lead to a doubling (or more) of the warming rate in the 2020s compared with the rate since 1980 with strong spatiotemporal heterogeneity.The warming effect is consistent with the recent observed strong warming in 2023 and expected to make the 2020s anomalously warm. The forcing represents 80% to the measured increase in planetary heat uptake since 2020. The radiative forcing also has strong hemispheric contrast, which has important implications for precipitation pattern changes. Our result suggests marine cloud brightening may be a viable geoengineering method in temporarily cooling the climate that has its unique challenges due to inherent spatiotemporal heterogeneity.

Tianle Yuan↗

Spatiotemporal Filtering Using Principal Component Analysis and Karhunen-Loeve Expansion Approaches for Regional GPS Network Analysis

Spatial filtering is an effective way to improve the precision of coordinate time series for regional GPS networks by reducing so-called common mode errors, thereby providing better resolution for detecting weak or transient deformation signals. The commonly used approach to regional filtering assumes that the common mode error is spatially uniform, which is a good approximation for networks of hundreds of kilometers extent, but breaks down as the spatial extent increases. A more rigorous approach should remove the assumption of spatially uniform distribution and let the data themselves reveal the spatial distribution of the common mode error. The principal component analysis (PCA) and the Karhunen-Loeve expansion (KLE) both decompose network time series into a set of temporally varying modes and their spatial responses. Therefore they provide a mathematical framework to perform spatiotemporal filtering.We apply the combination of PCA and KLE to daily station coordinate time series of the Southern California Integrated GPS Network (SCIGN) for the period 2000 to 2004. We demonstrate that spatially and temporally correlated common mode errors are the dominant error source in daily GPS solutions. The spatial characteristics of the common mode errors are close to uniform for all east, north, and vertical components, which implies a very long wavelength source for the common mode errors, compared to the spatial extent of the GPS network in southern California. Furthermore, the common mode errors exhibit temporally nonrandom patterns.

displacement↗

Spatially-Explicit Monitoring of Crop Photosynthetic Capacity Through the Use of Space-Based Chlorophyll Fluorescence Data

Plant functional traits such as photosynthetic capacity are critical parameters for terrestrial biosphere models. However, their spatial and temporal characteristics are still poorly represented. In this study, we used satellite observations of sun-induced fluorescence (SIF) to estimate top-of-canopy photosynthetic capacity (maximum carboxylation rate, V(sub cmax) at a reference temperature of 25 deg C) for crops, which was in turn utilized to simulate regional gross primary production (GPP). We first estimate the key parameter, V(sub cmax), in the widely-used FvCB photosynthesis model using field measurements of CO2 and water fluxes during 2007-2012 at seven crop eddy covariance flux sites over the US Corn Belt. The results showed that satellite far-red SIF retrievals have a stronger link to V(sub cmax) at the seasonal scale (R(squared)=0.70 for C4 and R(squared)=0.63 for C3 crop) as compared with widely-used vegetation indices. We calibrate an empirical model linking V(sub cmax) with SIF that was used to estimate spatially and temporally varying crop V(sub cmax) for the US Corn Belt region. The resulting V(sub cmax) maps are used together with meteorological data from MERRA reanalysis data and vegetation structural parameters derived from the satellite- based spectral reflectance data to constrain the Soil-Canopy Observation of Photosynthesis and Energy (SCOPE) balance model in order to estimate regional crop GPP. Our results show a substantial improvement in the seasonal and spatial patterns of cropland GPP when compared with crop yield inventory data. The evaluation with tall tower atmospheric CO2 measurements further supports our estimation of spatiotemporal V(sub cmax) from space-borne SIF. Considering that SIF has a direct link to photosynthetic activity, our findings highlight the potential to infer regional V(sub cmax) using remotely sensed SIF data and to use this information for a better quantification of regional cropland carbon cycles.

Drought↗

Vertical Structure of a Springtime Smoky and Humid Troposphere Over the Southeast Atlantic From Aircraft and Reanalysis

The springtime atmosphere over the southeast Atlantic Ocean (SEA) is subjected to a consistent layer of biomass burning (BB) smoke from widespread fires on the African continent. An elevated humidity signal is coincident with this layer, consistently proportional to the amount of smoke present. The combined humidity and BB aerosol has potentially significant radiative and dynamic impacts. Here, we use aircraft-based observations from the NASA ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) deployments in conjunction with reanalyses to characterize covariations in humidity and BB smoke across the SEA. The observed plume–vapor relationship, and its agreement with the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis version 5 (ERA5) and Copernicus Atmosphere Monitoring Service (CAMS) reanalysis, persists across all observations, although the magnitude of the relationship varies as the season progresses. Water vapor is well represented by the reanalyses, while CAMS tends to underestimate carbon monoxide especially under high BB. While CAMS aerosol optical depth (AOD) is generally overestimated relative to ORACLES AOD, the observations show a consistent relationship between carbon monoxide (CO) and aerosol extinction, demonstrating the utility of the CO tracer to understanding vertical aerosol distribution. We next use k-means clustering of the reanalyses to examine multi-year seasonal patterns and distributions. We identify canonical profile types of humidity and of CO, allowing us to characterize changes in vapor and BB atmospheric structures, and their impacts as they covary. While the humidity profiles show a range in both total water vapor concentration and in vertical structure, the CO profiles primarily vary in terms of maximum concentration, with similar vertical structures in each. The distribution of profile types varies spatiotemporally across the SEA region and through the season, ranging from largely one type in the northeast and southwest to more evenly distributed between multiple types where air masses meet in the middle of the SEA. These distributions follow patterns of transport from the humid, smoky source region (greatest influence in the northeast of the SEA) and the seasonal changes in both humidity and smoke (increasing and decreasing through the season, respectively). With this work, we establish a framework for a more complete analysis of the broader radiative and dynamical effects of humid aerosols over the SEA.

Kristina Pistone↗

Direction of Moving Plaids is Biased by Asymmetric Viewing Windows

Directionally selective V1 neurons are tuned to particular spatiotemporal frequencies and respond to local 1-D edge motion. At least some MT neurons however appear to respond to the actual velocity of moving 2D patterns (Movshon et al., EBR, 11:117, 1986). To better understand how the 1D local motion information available from V1 is integrated to derive a 2D velocity signal we investigated human perception of moving plaids, 2-D patterns composed of the sum of two 1-D sinusoidal gratings of different orientations. We measured the effect of the shape of the viewing window on the perceived direction of plaid motion. The plaids were spatially windowed by 2-D Gaussians with unequal standard deviations (sigma l, sigma 2). Four observers indicated perceived direction by adjusting a pointer. Direction errors were measured as a function of the difference between window orientation and true plaid direction (delta theta) for several grating spatial frequencies (SF = 0.3, 0.6, 1.2 c/d) and window aspect ratios (AR = sigma 1/sigma 2 = 1, 1.4, 2,4). Observers showed systematic errors in perceived direction (approx. 15 degrees for R = 4 and SF = 0.6 c/d) that peaked at delta theta approximately 40 degrees. The errors increased for increasing aspect ratio and decreased for increasing spatial frequencies (or number of cycles). These results show that despite the unambiguous motion of the plaids, within asymmetric windows, human can systematically misperceive plaid direction. These data constrain models of motion integration within extrastriate cortex and, in particular, are inconsistent with algorithms that use either the Intersection of Constraints rule or cross correlation to compute velocity.

Beutter, Brent Robert↗

Reduction of display artifacts by random sampling

The application of random-sampling techniques to remove visible artifacts (such as flicker, moire patterns, and paradoxical motion) introduced in TV-type displays by discrete sequential scanning is discussed and demonstrated. Sequential-scanning artifacts are described; the window of visibility defined in spatiotemporal frequency space by Watson and Ahumada (1982 and 1983) and Watson et al. (1983) is explained; the basic principles of random sampling are reviewed and illustrated by the case of the human retina; and it is proposed that the sampling artifacts can be replaced by random noise, which can then be shifted to frequency-space regions outside the window of visibility. Vertical sequential, single-random-sequence, and continuously renewed random-sequence plotting displays generating 128 points at update rates up to 130 Hz are applied to images of stationary and moving lines, and best results are obtained with the single random sequence for the stationary lines and with the renewed random sequence for the moving lines.

Ahumada, A. J., Jr.↗

Identifying and Analyzing Uncertainty Structures in the TRMM Microwave Imager Precipitation Product over Tropical Ocean Basins

Despite continuous improvements in microwave sensors and retrieval algorithms, our understanding of precipitation uncertainty is quite limited, due primarily to inconsistent findings in studies that compare satellite estimates to in situ observations over different parts of the world. This study seeks to characterize the temporal and spatial properties of uncertainty in the Tropical Rainfall Measuring Mission Microwave Imager surface rainfall product over tropical ocean basins. Two uncertainty analysis frameworks are introduced to qualitatively evaluate the properties of uncertainty under a hierarchy of spatiotemporal data resolutions. The first framework (i.e. 'climate method') demonstrates that, apart from random errors and regionally dependent biases, a large component of the overall precipitation uncertainty is manifested in cyclical patterns that are closely related to large-scale atmospheric modes of variability. By estimating the magnitudes of major uncertainty sources independently, the climate method is able to explain 45-88% of the monthly uncertainty variability. The percentage is largely resolution dependent (with the lowest percentage explained associated with a 1 deg x 1 deg spatial/1 month temporal resolution, and highest associated with a 3 deg x 3 deg spatial/3 month temporal resolution). The second framework (i.e. 'weather method') explains regional mean precipitation uncertainty as a summation of uncertainties associated with individual precipitation systems. By further assuming that self-similar recurring precipitation systems yield qualitatively comparable precipitation uncertainties, the weather method can consistently resolve about 50 % of the daily uncertainty variability, with only limited dependence on the regions of interest.

Climate↗

First Directly Retrieved Global Distribution of Tropospheric Column Ozone from GOME: Comparison with the GEOS-CHEM Model

We present the first directly retrieved global distribution of tropospheric column ozone from Global Ozone Monitoring Experiment (GOME) ultraviolet measurements during December 1996 to November 1997. The retrievals clearly show signals due to convection, biomass burning, stratospheric influence, pollution, and transport. They are capable of capturing the spatiotemporal evolution of tropospheric column ozone in response to regional or short time-scale events such as the 1997-1998 El Nino event and a 10-20 DU change within a few days. The global distribution of tropospheric column ozone displays the well-known wave-1 pattern in the tropics, nearly zonal bands of enhanced tropospheric column ozone of 36-48 DU at 20degS-30degS during the austral spring and at 25degN-45degN during the boreal spring and summer, low tropospheric column ozone of <30 DU uniformly distributed south of 35 S during all seasons, and relatively high tropospheric column ozone of >33 DU at some northern high-latitudes during the spring. Simulation from a chemical transport model corroborates most of the above structures, with small biases of <+/-5 DU and consistent seasonal cycles in most regions, especially in the southern hemisphere. However, significant positive biases of 5-20 DU occur in some northern tropical and subtropical regions such as the Middle East during summer. Comparison of GOME with monthly-averaged Measurement of Ozone and Water Vapor by Airbus in-service Aircraft (MOZAIC) tropospheric column ozone for these regions usually shows good consistency within 1 a standard deviations and retrieval uncertainties. Some biases can be accounted for by inadequate sensitivity to lower tropospheric ozone, the different spatiotemporal sampling and the spatiotemporal variations in tropospheric column ozone.

Liu, Xiong↗

Spatiotemporal Prediction of Fine Particulate Matter Using High-Resolution Satellite Images in the Southeastern US 2003-2011

Numerous studies have demonstrated that fine particulate matter (PM(sub 2.5), particles smaller than 2.5 micrometers in aerodynamic diameter) is associated with adverse health outcomes. The use of ground monitoring stations of PM(sub 2.5) to assess personal exposure, however, induces measurement error. Land-use regression provides spatially resolved predictions but land-use terms do not vary temporally. Meanwhile, the advent of satellite-retrieved aerosol optical depth (AOD) products have made possible to predict the spatial and temporal patterns of PM(sub 2.5) exposures. In this paper, we used AOD data with other PM(sub 2.5) variables, such as meteorological variables, land-use regression, and spatial smoothing to predict daily concentrations of PM(sub 2.5) at a 1 sq km resolution of the Southeastern United States including the seven states of Georgia, North Carolina, South Carolina, Alabama, Tennessee, Mississippi, and Florida for the years from 2003 to 2011. We divided the study area into three regions and applied separate mixed-effect models to calibrate AOD using ground PM(sub 2.5) measurements and other spatiotemporal predictors. Using 10-fold cross-validation, we obtained out of sample R2 values of 0.77, 0.81, and 0.70 with the square root of the mean squared prediction errors of 2.89, 2.51, and 2.82 cu micrograms for regions 1, 2, and 3, respectively. The slopes of the relationships between predicted PM2.5 and held out measurements were approximately 1 indicating no bias between the observed and modeled PM(sub 2.5) concentrations. Predictions can be used in epidemiological studies investigating the effects of both acute and chronic exposures to PM(sub 2.5). Our model results will also extend the existing studies on PM(sub 2.5) which have mostly focused on urban areas because of the paucity of monitors in rural areas.

PM2.5↗

Spatiotemporal Organization of Energy Release Events in the Quiet Solar Corona

Using data from the STEREO and SOHO spacecraft, we show that temporal organization of energy release events in the quiet solar corona is close to random, in contrast to the clustered behavior of flaring times in solar active regions. The locations of the quiet-Sun events follow the meso- and supergranulation pattern of the underling photosphere. Together with earlier reports of the scale-free event size statistics, our findings suggest that quiet solar regions responsible for bulk coronal heating operate in a driven self-organized critical state, possibly involving long-range Alfvenic interactions.

sun activity↗

Evaluation of Machine Learning and Deep Learning Algorithms for Fire Prediction in Southeast Asia

Vegetation fires are prevalent in South/Southeast Asian countries, making fire prediction crucial due to their potential environmental, economic, and social impacts. Accurate predictions of fires facilitate timely interventions, helping to mitigate uncontrolled fires that can lead to biodiversity loss and air quality issues. In this study, we utilize VIIRS satellite-derived fire data alongside six machine learning and deep learning models—Simple Persistence, Multi-Layer Perceptron (MLP), Convolutional Neural Network (CNN), Long Short-Term Memory (LSTM), CNN-LSTM, and ConvLSTM—to determine the most effective fire prediction model, using Root Mean Square Error (RMSE) as the metric. Our results indicate that the CNN model is the most reliable in regions with spatial dependencies, such as Brunei, Indonesia, Malaysia, the Philippines, Timor-Leste, and Thailand. Conversely, the ConvLSTM model excels in countries with complex spatiotemporal dynamics like Laos, Myanmar, and Vietnam. The CNN-LSTM hybrid model also performed well in Cambodia, suggesting a need for a balanced approach in areas requiring both spatial and temporal feature extraction. Furthermore, simpler models like Persistence and MLP showed limitations in capturing dynamic patterns and temporal dependencies. Our findings highlight the importance of evaluating models before implementing any decision support systems (DSS) in fire management. By tailoring models to specific regional fire data, we can enhance prediction accuracy and responsiveness, ultimately improving fire risk management in Southeast Asia and beyond.

Deep learning↗

Regional Analysis of the 2015–16 Lower Mekong River Basin Drought Using NASA Satellite Observations

Study region Lower Mekong River Basin (LMRB) Study focus Satellite remote sensing products are widely used for monitoring droughts. Using NASA satellite sensors of precipitation (Global Measurement Mission, GPM), soil moisture (Soil Moisture Active and Passive, SMAP), and terrestrial water storage (Gravity Recovery and Climate Experiment, GRACE), this study evaluates the historical drought in the LMRB during 2015–16. SMAP soil moisture was validated against in-situ soil moisture, and GPM precipitation and SMAP soil moisture were cross-validated with streamflow observations. The spatiotemporal dynamics of soil moisture were also examined in different ranges of catchment areas. In performing the analysis, we used lagged correlations between hydrological variables and the indices of the Standardized Precipitation Index (SPI) and Standardized Streamflow Index (SSI). New hydrological insights for the regions Spatio-temporal patterns of drought in 2015–16 were examined from the entire basin to small watersheds. A mismatch occurs when using GRACE data to study droughts in small watersheds (many of the small watersheds would be a fraction of the few 100 km2 spatial resolutions of GRACE pixel). In smaller watersheds, hydrological drought (SSI) was closely defined with SMAP soil moisture downscaled to 1 km rather than the meteorological drought index (SPI). By leveraging satellite-based observations across a range of spatial scales, this study highlights the utility of Earth observations in informing water resources and land management decisions at the regional scale.

Mekong↗

Quantifying Spatiotemporal Variability of Glacier Algal Blooms and the Impact on Surface Albedo in Southwestern Greenland

Albedo reduction due to light-absorbing impurities can substantially enhance ice sheet surface melt by increasing surface absorption of solar energy. Glacier algae have been suggested to play a critical role in darkening the ablation zone in southwestern Greenland. It was very recently found that the Sentinel-3 Ocean and Land Colour Instrument (OLCI) band ratio R709 nm∕R673 nm can characterize the spatial patterns of glacier algal blooms. However, Sentinel-3 was launched in 2016, and current data are only available over three melting seasons (2016–2019). Here, we demonstrate the capability of the MEdium Resolution Imaging Spectrometer (MERIS) for mapping glacier algae from space and extend the quantification of glacier algal blooms over southwestern Greenland back to the period 2004–2011. Several band ratio indices (MERIS chlorophyll a indices and the impurity index) were computed and compared with each other. The results indicate that the MERIS two-band ratio index (2BDA) R709 nm∕R665 nm is very effective in capturing the spatial distribution and temporal dynamics of glacier algal growth on bare ice in July and August. We analyzed the interannual (2004–2011) and summer (July–August) trends of algal distribution and found significant seasonal and interannual increases in glacier algae close to the Jakobshavn Isbrae Glacier and along the middle dark zone between the altitudes of 1200 and 1400 m. Using broadband albedo data from the Moderate Resolution Imaging Spectroradiometer (MODIS), we quantified the impact of glacier algal growth on bare ice albedo, finding a significant correlation between algal development and albedo reduction over algae-abundant areas. Our analysis indicates the strong potential for the satellite algal index to be used to reduce bare ice albedo biases in regional climate model simulations.

surface albedo↗

The Dawn of Geostationary Air Quality Monitoring: Case Studies from Seoul and Los Angeles

With the near-future launch of geostationary pollution monitoring satellite instruments over North America, East Asia, and Europe, the air quality community is preparing for an integrated global atmospheric composition observing system at unprecedented spatial and temporal resolutions. One of the ways that NASA has supported this community preparation is through demonstration of future space-borne capabilities using the Geostationary Trace gas and Aerosol Sensor Optimization (GeoTASO) airborne instrument. This paper integrates repeated high-resolution maps from GeoTASO, ground-based Pandora spectrometers, and low Earth orbit measurements from the Ozone Mapping and Profiler Suite (OMPS), for case studies over two metropolitan areas: Seoul, South Korea on June 9, 2016 and Los Angeles, California on June 27, 2017. This dataset provides a unique opportunity to illustrate how geostationary air quality monitoring platforms and ground-based remote sensing networks will close the current spatiotemporal observation gap. GeoTASO observes large differences in diurnal behavior between these urban areas, with NO2 accumulating within the Seoul Metropolitan Area through the day but NO2 peaking in the morning and decreasing throughout the afternoon in the Los Angeles Basin. In both areas, the earliest morning maps exhibit spatial patterns similar to emission source areas (e.g., urbanized valleys, roadways, major airports). These spatial patterns change later in the day due to boundary layer dynamics, horizontal transport, and chemistry. The nominal resolution of GeoTASO is finer than will be obtained from geostationary platforms, but when NO2 data over Los Angeles are up-scaled to the expected resolution of TEMPO, spatial features discussed are conserved. Pandora instruments installed in both metropolitan areas capture the diurnal patterns observed by GeoTASO, continuously and over longer time periods, and will play a critical role in validation of the next generation of satellite measurement. These case studies demonstrate that different regions can have diverse diurnal patterns and that day-to-day variability due to meteorology or anthropogenic patterns such as weekday/weekend variations in emissions is large. Low Earth orbit measurements, despite their inability to capture the diurnal patterns at fine spatial resolution, will be essential for intercalibrating the geostationary radiances and cross-validating the geostationary retrievals in an integrated global observing system.

Geostationary↗

Recovery: Fast and Slow—Vegetation Response During the 2012–2016 California Drought

The 2012–2016 California Drought severely impacted natural vegetation across a wide range of environmental gradient. Although several studies have reported an increase in plant water stress and mortality, the spatiotemporal variations of ecosystem productivity responses and the associated environmental and biological drivers remain unclear. Here, using Enhanced Vegetation Index from the Moderate resolution imaging spectrometer, we found that 45% of the natural ecosystems showed an abrupt change (breakpoint [BP]) in productivity during 2012–2016. There were three major contrasting temporal patterns of productivity responses: (i) a steady increase under higher temperature followed by a decline due to accumulated moisture depletion (high elevation forest) or temperature decrease (high elevation nonforest), (ii) gradual decline during the drought followed by a rapid recovery within 1 year after drought stress was partially relieved, and (iii) both a gradual decline and an abrupt decline. The magnitude of abrupt changes was negatively correlated (r = −0.80, p < 0.001) with initial gradual changes. Overall, changes during BP offset, on average, 57% of the preceding gradual responses. The spatial variability in ecosystem response patterns is driven by both environmental and biological factors. Particularly, for forests, positive BP was driven by increasing rainfall and decreasing temperature, while negative BP was mainly driven by the precipitation anomaly. By 2019, 33% of the natural vegetation have recovered to the level of EVI in 2010. Ecosystem responses to multiyear droughts can influence ecosystem dynamics in a complex pattern. Multiple ecohydrological factors should be considered to understand and predict the long-term drought impacts on ecosystems.

Xi Yang↗