DomiRank Centrality reveals structural fragility of complex networks via node dominance.

Determining the key elements of interconnected infrastructure and complex systems is paramount to ensure system functionality and integrity. This work quantifies the dominance of the networks' nodes in their respective neighborhoods, introducing a centrality metric, DomiRank, that integrates local and global topological information via a tunable parameter. We present an analytical formula and an efficient parallelizable algorithm for DomiRank centrality, making it applicable to massive networks. From the networks' structure and function perspective, nodes with high values of DomiRank highlight fragile neighborhoods whose integrity and functionality are highly dependent on those dominant nodes. Underscoring this relation between dominance and fragility, we show that DomiRank systematically outperforms other centrality metrics in generating targeted attacks that effectively compromise network structure and disrupt its functionality for synthetic and real-world topologies. Moreover, we show that DomiRank-based attacks inflict more enduring damage in the network, hindering its ability to rebound and, thus, impairing system resilience. DomiRank centrality capitalizes on the competition mechanism embedded in its definition to expose the fragility of networks, paving the way to design strategies to mitigate vulnerability and enhance the resilience of critical infrastructures.

Climate-driven changes in the predictability of seasonal precipitation.

Climate-driven changes in precipitation amounts and their seasonal variability are expected in many continental-scale regions during the remainder of the 21st century. However, much less is known about future changes in the predictability of seasonal precipitation, an important earth system property relevant for climate adaptation. Here, on the basis of CMIP6 models that capture the present-day teleconnections between seasonal precipitation and previous-season sea surface temperature (SST), we show that climate change is expected to alter the SST-precipitation relationships and thus our ability to predict seasonal precipitation by 2100. Specifically, in the tropics, seasonal precipitation predictability from SSTs is projected to increase throughout the year, except the northern Amazonia during boreal winter. Concurrently, in the extra-tropics predictability is likely to increase in central Asia during boreal spring and winter. The altered predictability, together with enhanced interannual variability of seasonal precipitation, poses new opportunities and challenges for regional water management.

Historical spatiotemporal changes in fire danger potential across biomes.

This study 1) identifies the seasons and biomes that exhibit significant (1980-2019) changes in fire danger potential, as quantified by the Canadian Fire Weather Index (FWI); 2) explores what types of fire behavior potentials may be contributing to changes in fire danger potential, as quantified by the United States Energy Release Component (ERC) and the Ignition Component (IC); 3) provides spatiotemporal insight on how fire danger potential and fire behavior potential are responding in relation to changes in seasonal precipitation totals and seasonal mean air temperature across biomes. Time series of these fire potentials, as well as seasonal mean temperature, and seasonal precipitation totals are generated using data from the 0.25° ECMWF spatial resolution Reanalysis 5th Generation (ERA5) and the Climatic Research Unit gridded Time Series (CRU TS). The Mann-Kendall test is then applied to identify significant spatiotemporal trends across each biome. Results indicate that the September-November season (SON) exhibits the greatest rate of increase in fire danger potential, followed by the June-August season (JJA), December, January-February season (DJF), and March-May season (MAM), and this is predominant over the Tropical and Subtropical Moist Broadleaf Forest Biome, as well as all vegetation types of the temperate biomes. Similarly, the temperate biomes experience the greatest rate of increase in fire intensity potential and ignition potential, but prevalent during the DJF and MAM seasons. Furthermore, there is a significant positive correlation between fire danger potential and seasonal mean air temperature during JJA in the Northern Hemisphere for the temperate biomes in North America and Europe, as well as the Tropical and Subtropical biomes in Africa. Our analysis provides quantitative insight as to how fire danger potential and fire behavior potential have been responding to changes in seasonal mean temperature and seasonal precipitation totals across different ecoregions around the world.

Hotspots of Predictability: Identifying Regions of High Precipitation Predictability at Seasonal Timescales From Limited Time Series Observations.

Precipitation prediction at seasonal timescales is important for planning and management of water resources as well as preparedness for hazards such as floods, droughts and wildfires. Quantifying predictability is quite challenging as a consequence of a large number of potential drivers, varying antecedent conditions, and small sample size of high-quality observations available at seasonal timescales, that in turn, increases prediction uncertainty and the risk of model overfitting. Here, we introduce a generalized probabilistic framework to account for these issues and assess predictability under uncertainty. We focus on prediction of winter (Nov-Mar) precipitation across the contiguous United States, using sea surface temperature-derived indices (averaged in Aug-Oct) as predictors. In our analysis we identify "predictability hotspots," which we define as regions where precipitation is inherently more predictable. Our framework estimates the entire predictive distribution of precipitation using copulas and quantifies prediction uncertainties, while employing principal component analysis for dimensionality reduction and a cross validation technique to avoid overfitting. We also evaluate how predictability changes across different quantiles of the precipitation distribution (dry, normal, wet amounts) using a multi-category 3 × 3 contingency table. Our results indicate that well-defined predictability hotspots occur in the Southwest and Southeast. Moreover, extreme dry and wet conditions are shown to be relatively more predictable compared to normal conditions. Our study may help with water resources management in several subregions of the United States and can be used to assess the fidelity of earth system models in successfully representing teleconnections and predictability.

California wildfire spread derived using VIIRS satellite observations and an object-based tracking system.

Changing wildfire regimes in the western US and other fire-prone regions pose considerable risks to human health and ecosystem function. However, our understanding of wildfire behavior is still limited by a lack of data products that systematically quantify fire spread, behavior and impacts. Here we develop a novel object-based system for tracking the progression of individual fires using 375 m Visible Infrared Imaging Radiometer Suite active fire detections. At each half-daily time step, fire pixels are clustered according to their spatial proximity, and are either appended to an existing active fire object or are assigned to a new object. This automatic system allows us to update the attributes of each fire event, delineate the fire perimeter, and identify the active fire front shortly after satellite data acquisition. Using this system, we mapped the history of California fires during 2012-2020. Our approach and data stream may be useful for calibration and evaluation of fire spread models, estimation of near-real-time wildfire emissions, and as means for prescribing initial conditions in fire forecast models.

Critical Tokunaga model for river networks.

The hierarchical organization and self-similarity in river basins have been topics of extensive research in hydrology and geomorphology starting with the pioneering work of Horton in 1945. Despite significant theoretical and applied advances, however, the mathematical origin of and relation among Horton's laws for different stream attributes remain unsettled. Here we capitalize on a recently developed theory of random self-similar trees to elucidate the origin of Horton's laws, Hack's laws, basin fractal dimensions, power-law distributions of link attributes, and power-law relations between distinct attributes. In particular, we introduce a one-parametric family of self-similar critical Tokunaga trees that includes the celebrated Shreve's random topology model and extends to trees that approximate the observed river networks with realistic exponents. The results offer tools to increase our understanding of landscape organization under different hydroclimatic forcings, and to extend scaling relationships useful for hydrologic prediction to resolutions higher than those observed.

Probabilistic Evaluation of Drought in CMIP6 Simulations.

As droughts have widespread social and ecological impacts, it is critical to develop long-term adaptation and mitigation strategies to reduce drought vulnerability. Climate models are important in quantifying drought changes. Here, we assess the ability of 285 CMIP6 historical simulations, from 17 models, to reproduce drought duration and severity in three observational data sets using the Standardized Precipitation Index (SPI). We used summary statistics beyond the mean and standard deviation, and devised a novel probabilistic framework, based on the Hellinger distance, to quantify the difference between observed and simulated drought characteristics. Results show that many simulations have less than ± 10 % error in reproducing the observed drought summary statistics. The hypothesis that simulations and observations are described by the same distribution cannot be rejected for more than 80 % of the grids based on our H distance framework. No single model stood out as demonstrating consistently better performance over large regions of the globe. The variance in drought statistics among the simulations is higher in the tropics compared to other latitudinal zones. Though the models capture the characteristics of dry spells well, there is considerable bias in low precipitation values. Good model performance in terms of SPI does not imply good performance in simulating low precipitation. Our study emphasizes the need to probabilistically evaluate climate model simulations in order to both pinpoint model weaknesses and identify a subset of best-performing models that are useful for impact assessments.

Integrated assessment modeling reveals near-channel management as cost-effective to improve water quality in agricultural watersheds.

Despite decades of policy that strives to reduce nutrient and sediment export from agricultural fields, surface water quality in intensively managed agricultural landscapes remains highly degraded. Recent analyses show that current conservation efforts are not sufficient to reverse widespread water degradation in Midwestern agricultural systems. Intensifying row crop agriculture and increasing climate pressure require a more integrated approach to water quality management that addresses diverse sources of nutrients and sediment and off-field mitigation actions. We used multiobjective optimization analysis and integrated three biophysical models to evaluate the cost-effectiveness of alternative portfolios of watershed management practices at achieving nitrate and suspended sediment reduction goals in an agricultural basin of the Upper Midwestern United States. Integrating watershed-scale models enabled the inclusion of near-channel management alongside more typical field management and thus directly the comparison of cost-effectiveness across portfolios. The optimization analysis revealed that fluvial wetlands (i.e., wide, slow-flowing, vegetated water bodies within the riverine corridor) are the single-most cost-effective management action to reduce both nitrate and sediment loads and will be essential for meeting moderate to aggressive water quality targets. Although highly cost-effective, wetland construction was costly compared to other practices, and it was not selected in portfolios at low investment levels. Wetland performance was sensitive to placement, emphasizing the importance of watershed scale planning to realize potential benefits of wetland restorations. We conclude that extensive interagency cooperation and coordination at a watershed scale is required to achieve substantial, economically viable improvements in water quality under intensive row crop agricultural production.

Underestimated MJO variability in CMIP6 models.

The Madden-Julian Oscillation (MJO) is the leading mode of intra-seasonal climate variability, having profound impacts on a wide range of weather and climate phenomena. Here, we use a wavelet-based spectral Principal Component Analysis (wsPCA) to evaluate the skill of 20 state-of-the-art CMIP6 models in capturing the magnitude and dynamics of the MJO. By construction, wsPCA has the ability to focus on desired frequencies and capture each propagative physical mode with one principal component (PC). We show that the MJO contribution to the total intra-seasonal climate variability is substantially underestimated in most CMIP6 models. The joint distribution of the modulus and angular frequency of the wavelet PC series associated with MJO is used to rank models relatively to the observations through the Wasserstein distance. Using Hovmöller phase-longitude diagrams, we also show that precipitation variability associated with MJO is underestimated in most CMIP6 models for the Amazonia, Southwest Africa, and Maritime Continent.

Zonally contrasting shifts of the tropical rainbelt in response to climate change.

Future changes in the position of the intertropical convergence zone (ITCZ; a narrow band of heavy precipitation in the tropics) with climate change could affect the livelihood and food security of billions of people. Although models predict a future narrowing of the ITCZ, uncertainties remain large regarding its future position, with most past work focusing on zonal-mean shifts. Here we use projections from 27 state-of-the-art (CMIP6) climate models and document a robust zonally-varying ITCZ response to the SSP3-7.0 scenario by 2100, with a northward shift over eastern Africa and the Indian Ocean, and a southward shift in the eastern Pacific and Atlantic Oceans. The zonally-varying response is consistent with changes in the divergent atmospheric energy transport, and sector-mean shifts of the energy flux equator. Our analysis provides insight about mechanisms influencing the future position of the tropical rainbelt, and may allow for more robust projections of climate change impacts.

Rotated Spectral Principal Component Analysis (rsPCA) for Identifying Dynamical Modes of Variability in Climate Systems.

Spectral PCA (sPCA), in contrast to classical PCA, offers the advantage of identifying organized spatiotemporal patterns within specific frequency bands and extracting dynamical modes. However, the unavoidable trade-off between frequency resolution and robustness of the PCs leads to high sensitivity to noise and overfitting, which limits the interpretation of the sPCA results. We propose herein a simple nonparametric implementation of sPCA using the continuous analytic Morlet wavelet as a robust estimator of the cross-spectral matrices with good frequency resolution. To improve the interpretability of the results, especially when several modes of similar amplitude exist within the same frequency band, we propose a rotation of the complex-valued eigenvectors to optimize their spatial regularity (smoothness). The developed method, called rotated spectral PCA (rsPCA), is tested on synthetic data simulating propagating waves and shows impressive performance even with high levels of noise in the data. Applied to global historical geopotential height (GPH) and sea surface temperature (SST) daily time series, the method accurately captures patterns of atmospheric Rossby waves at high frequencies (3-60-day periods) in both GPH and SST and El Niño-Southern Oscillation (ENSO) at low frequencies (2-7-yr periodicity) in SST. At high frequencies the rsPCA successfully unmixes the identified waves, revealing spatially coherent patterns with robust propagation dynamics.

Climate Signatures on Lake And Wetland Size Distributions in Arctic Deltas.

Understanding how thermokarst lakes on arctic river deltas will respond to rapid warming is critical for projecting how carbon storage and fluxes will change in those vulnerable environments. Yet, this understanding is currently limited partly due to the complexity of disentangling significant interannual variability from the longer-term surface water signatures on the landscape, using the short summertime window of optical spaceborne observations. Here, we rigorously separate perennial lakes from ephemeral wetlands on 12 arctic deltas and report distinct size distributions and climate trends for the two waterbodies. Namely, we find a lognormal distribution for lakes and a power-law distribution for wetlands, consistent with a simple proportionate growth model and inundated topography, respectively. Furthermore, while no trend with temperature is found for wetlands, a statistically significant decreasing trend of mean lake size with warmer temperatures is found, attributed to colder deltas having deeper and thicker permafrost preserving larger lakes.

From turbulence to landscapes: Logarithmic mean profiles in bounded complex systems.

We show that similarly to the logarithmic mean-velocity profile in wall-bounded turbulence, the landscape topography presents an intermediate region with a logarithmic mean-elevation profile. Such profiles are present in complex topographies with channel branching and fractal river networks resulting from model simulation, controlled laboratory experiments, and natural landscapes. Dimensional and self-similarity arguments are used to corroborate this finding. We also tested the presence of logarithmic profiles in discrete, minimalist models of networks obtained from optimality principles (optimal channel networks) and directed percolation. The emergence of self-similar scaling appears as a robust outcome in dynamically different, but spatially bounded, complex systems, as a dimensional consequence of length-scale independence.

Forecasting Global Fire Emissions on Subseasonal to Seasonal (S2S) Time Scales.

Fire emissions of gases and aerosols alter atmospheric composition and have substantial impacts on climate, ecosystem function, and human health. Warming climate and human expansion in fire-prone landscapes exacerbate fire impacts and call for more effective management tools. Here we developed a global fire forecasting system that predicts monthly emissions using past fire data and climate variables for lead times of 1 to 6 months. Using monthly fire emissions from the Global Fire Emissions Database (GFED) as the prediction target, we fit a statistical time series model, the Autoregressive Integrated Moving Average model with eXogenous variables (ARIMAX), in over 1,300 different fire regions. Optimized parameters were then used to forecast future emissions. The forecast system took into account information about region-specific seasonality, long-term trends, recent fire observations, and climate drivers representing both large-scale climate variability and local fire weather. We cross-validated the forecast skill of the system with different combinations of predictors and forecast lead times. The reference model, which combined endogenous and exogenous predictors with a 1 month forecast lead time, explained 52% of the variability in the global fire emissions anomaly, considerably exceeding the performance of a reference model that assumed persistent emissions during the forecast period. The system also successfully resolved detailed spatial patterns of fire emissions anomalies in regions with significant fire activity. This study bridges the gap between the efforts of near-real-time fire forecasts and seasonal fire outlooks and represents a step toward establishing an operational global fire, smoke, and carbon cycle forecasting system.

Channel Network Control on Seasonal Lake Area Dynamics in Arctic Deltas.

The abundant lakes dotting arctic deltas are hotspots of methane emissions and biogeochemical activity, but seasonal variability in lake extents introduces uncertainty in estimates of lacustrine carbon emissions, typically performed at annual or longer time scales. To characterize variability in lake extents, we analyzed summertime lake area loss (i.e., shrinkage) on two deltas over the past 20 years, using Landsat-derived water masks. We find that monthly shrinkage rates have a pronounced structured variability around the channel network with the shrinkage rate systematically decreasing farther away from the channels. This pattern of shrinkage is predominantly attributed to a deeper active layer enhancing near-surface connectivity and storage and greater vegetation density closer to the channels leading to increased evapotranspiration rates. This shrinkage signal, easily extracted from remote sensing observations, may offer the means to constrain estimates of lacustrine methane emissions and to develop process-based estimates of depth to permafrost on arctic deltas.

Graph-Guided Regularized Regression of Pacific Ocean Climate Variables to Increase Predictive Skill of Southwestern U.S. Winter Precipitation.

Understanding the physical drivers of seasonal hydroclimatic variability and improving predictive skill remains a challenge with important socioeconomic and environmental implications for many regions around the world. Physics-based deterministic models show limited ability to predict precipitation as the lead time increases, due to imperfect representation of physical processes and incomplete knowledge of initial conditions. Similarly, statistical methods drawing upon established climate teleconnections have low prediction skill due to the complex nature of the climate system. Recently, promising data-driven approaches have been proposed, but they often suffer from overparameterization and overfitting due to the short observational record, and they often do not account for spatiotemporal dependencies among covariates (i.e., predictors such as sea surface temperatures). This study addresses these challenges via a predictive model based on a graph-guided regularizer that simultaneously promotes similarity of predictive weights for highly correlated covariates and enforces sparsity in the covariate domain. This approach both decreases the effective dimensionality of the problem and identifies the most predictive features without specifying them a priori. We use large ensemble simulations from a climate model to construct this regularizer, reducing the structural uncertainty in the estimation. We apply the learned model to predict winter precipitation in the southwestern United States using sea surface temperatures over the entire Pacific basin, and demonstrate its superiority compared to other regularization approaches and statistical models informed by known teleconnections. Our results highlight the potential to combine optimally the space-time structure of predictor variables learned from climate models with new graph-based regularizers to improve seasonal prediction.

Beyond the Pixel: Using Patterns and Multiscale Spatial Information to Improve the Retrieval of Precipitation from Spaceborne Passive Microwave Imagers.

The quantitative estimation of precipitation from orbiting passive microwave imagers has been performed for more than 30 years. The development of retrieval methods consists of establishing physical or statistical relationships between the brightness temperatures (TBs) measured at frequencies between 5 and 200 GHz and precipitation. Until now, these relationships have essentially been established at the "pixel" level, associating the average precipitation rate inside a predefined area (the pixel) to the collocated multispectral radiometric measurement. This approach considers each pixel as an independent realization of a process and ignores the fact that precipitation is a dynamic variable with rich multiscale spatial and temporal organization. Here we propose to look beyond the pixel values of the TBs and show that useful information for precipitation retrieval can be derived from the variations of the observed TBs in a spatial neighborhood around the pixel of interest. We also show that considering neighboring information allows us to better handle the complex observation geometry of conical-scanning microwave imagers, involving frequency-dependent beamwidths, overlapping fields of view, and large Earth incidence angles. Using spatial convolution filters, we compute "nonlocal" radiometric parameters sensitive to spatial patterns and scale-dependent structures of the TB fields, which are the "geometric signatures" of specific precipitation structures such as convective cells. We demonstrate that using nonlocal radiometric parameters to enrich the spectral information associated to each pixel allows for reduced retrieval uncertainty (reduction of 6%-11% of the mean absolute retrieval error) in a simple k-nearest neighbors retrieval scheme.