Large-scale remote sensing analysis reveals an increasing coupling of grassland vitality to atmospheric water demand.

Grasslands provide important ecosystem services to society, including biodiversity, water security, erosion control, and forage production. Grasslands are also vulnerable to droughts, rendering their future vitality under climate change uncertain. Yet, the grassland response to drought is not well understood, especially for heterogeneous Central European grasslands. We here fill this gap by quantifying the spatiotemporal sensitivity of grasslands to drought using a novel remote sensing dataset from Landsat/Sentinel-2 paired with climate re-analysis data. Specifically, we quantified annual grassland vitality at fine spatial scale and national extent (Germany) from 1985 to 2021. We analyzed grassland sensitivity to drought by testing for statistically robust links between grassland vitality and common drought indices. We furthermore explored the spatiotemporal variability of drought sensitivity for 12 grassland habitat types given their different biotic and abiotic features. Grassland vitality maps revealed a large-scale reduction of grassland vitality during past droughts. The unprecedented drought of 2018-2019 stood out as the largest multi-year vitality decline since the mid-1980s. Grassland vitality was consistently coupled to drought (R2 = .09-.22) with Vapor Pressure Deficit explaining vitality best. This suggests that high atmospheric water demand, as observed during recent compounding drought and heatwave events, has major impacts on grassland vitality in Central Europe. We found a significant increase in drought sensitivity over time with highest sensitivities detected in periods of extremely high atmospheric water demand, suggesting that drought impacts on grasslands are becoming more severe with ongoing climate change. The spatial variability of grassland drought sensitivity was linked to different habitat types, with declining sensitivity from dry and mesic to wet habitats. Our study provides the first large-scale, long-term, and spatially explicit evidence of increasing drought sensitivities of Central European grasslands. With rising compound droughts and heatwaves under climate change, large-scale grassland vitality loss, as in 2018-2019, will thus become more likely in the future.

Unveiling patterns in human dominated landscapes through mapping the mass of US built structures.

Built structures increasingly dominate the Earth's landscapes; their surging mass is currently overtaking global biomass. We here assess built structures in the conterminous US by quantifying the mass of 14 stock-building materials in eight building types and nine types of mobility infrastructures. Our high-resolution maps reveal that built structures have become 2.6 times heavier than all plant biomass across the country and that most inhabited areas are mass-dominated by buildings or infrastructure. We analyze determinants of the material intensity and show that densely built settlements have substantially lower per-capita material stocks, while highest intensities are found in sparsely populated regions due to ubiquitous infrastructures. Out-migration aggravates already high intensities in rural areas as people leave while built structures remain - highlighting that quantifying the distribution of built-up mass at high resolution is an essential contribution to understanding the biophysical basis of societies, and to inform strategies to design more resource-efficient settlements and a sustainable circular economy.

The global wildland-urban interface.

The wildland-urban interface (WUI) is where buildings and wildland vegetation meet or intermingle1,2. It is where human-environmental conflicts and risks can be concentrated, including the loss of houses and lives to wildfire, habitat loss and fragmentation and the spread of zoonotic diseases3. However, a global analysis of the WUI has been lacking. Here, we present a global map of the 2020 WUI at 10 m resolution using a globally consistent and validated approach based on remote sensing-derived datasets of building area4 and wildland vegetation5. We show that the WUI is a global phenomenon, identify many previously undocumented WUI hotspots and highlight the wide range of population density, land cover types and biomass levels in different parts of the global WUI. The WUI covers only 4.7% of the land surface but is home to nearly half its population (3.5 billion). The WUI is especially widespread in Europe (15% of the land area) and the temperate broadleaf and mixed forests biome (18%). Of all people living near 2003-2020 wildfires (0.4 billion), two thirds have their home in the WUI, most of them in Africa (150 million). Given that wildfire activity is predicted to increase because of climate change in many regions6, there is a need to understand housing growth and vegetation patterns as drivers of WUI change.

Beyond "greening" and "browning": Trends in grassland ground cover fractions across Eurasia that account for spatial and temporal autocorrelation.

Grassland ecosystems cover up to 40% of the global land area and provide many ecosystem services directly supporting the livelihoods of over 1 billion people. Monitoring long-term changes in grasslands is crucial for food security, biodiversity conservation, achieving Land Degradation Neutrality goals, and modeling the global carbon budget. Although long-term grassland monitoring using remote sensing is extensive, it is typically based on a single vegetation index and does not account for temporal and spatial autocorrelation, which means that some trends are falsely identified while others are missed. Our goal was to analyze trends in grasslands in Eurasia, the largest continuous grassland ecosystems on Earth. To do so, we calculated Cumulative Endmember Fractions (annual sums of monthly ground cover fractions) derived from MODIS 2002-2020 time series, and applied a new statistical approach PARTS that explicitly accounts for temporal and spatial autocorrelation in trends. We examined trends in green vegetation, non-photosynthetic vegetation, and soil ground cover fractions considering their independent change trajectories and relations among fractions over time. We derived temporally uncorrelated pixel-based trend maps and statistically tested whether observed trends could be explained by elevation, land cover, SPEI3, climate, country, and their combinations, all while accounting for spatial autocorrelation. We found no statistical evidence for a decrease in vegetation cover in grasslands in Eurasia. Instead, there was a significant map-level increase in non-photosynthetic vegetation across the region and local increases in green vegetation with a concomitant decrease in soil fraction. Independent environmental variables affected trends significantly, but effects varied by region. Overall, our analyses show in a statistically robust manner that Eurasian grasslands have changed considerably over the past two decades. Our approach enhances remote sensing-based monitoring of trends in grasslands so that underlying processes can be discerned.

High-resolution data and maps of material stock, population, and employment in Austria from 1985 to 2018.

High-resolution maps of material stocks in buildings and infrastructures are of key importance for studies of societal resource use (social metabolism, circular economy, secondary resource potentials) as well as for transport studies and land system science. So far, such maps were only available for specific years but not in time series. Even for single years, data covering entire countries with high resolution, or using remote-sensing data are rare. Instead, they often have local extent (e.g., [1]), are lower resolution (e.g., [2]), or are based on other geospatial data (e.g., [3]). We here present data on the material stocks in three types of buildings (commercial and industrial, single- and multifamily houses) and three types of infrastructures (roads, railways, other infrastructures) for a 33-year time series for Austria at a spatial resolution of 30 m. The article also presents data on population and employment in Austria for the same time period, at the same spatial resolution. Data were derived with the same method applied in a recent study for Germany [4].

Ten facts about land systems for sustainability.

Land use is central to addressing sustainability issues, including biodiversity conservation, climate change, food security, poverty alleviation, and sustainable energy. In this paper, we synthesize knowledge accumulated in land system science, the integrated study of terrestrial social-ecological systems, into 10 hard truths that have strong, general, empirical support. These facts help to explain the challenges of achieving sustainability in land use and thus also point toward solutions. The 10 facts are as follows: 1) Meanings and values of land are socially constructed and contested; 2) land systems exhibit complex behaviors with abrupt, hard-to-predict changes; 3) irreversible changes and path dependence are common features of land systems; 4) some land uses have a small footprint but very large impacts; 5) drivers and impacts of land-use change are globally interconnected and spill over to distant locations; 6) humanity lives on a used planet where all land provides benefits to societies; 7) land-use change usually entails trade-offs between different benefits-"win-wins" are thus rare; 8) land tenure and land-use claims are often unclear, overlapping, and contested; 9) the benefits and burdens from land are unequally distributed; and 10) land users have multiple, sometimes conflicting, ideas of what social and environmental justice entails. The facts have implications for governance, but do not provide fixed answers. Instead they constitute a set of core principles which can guide scientists, policy makers, and practitioners toward meeting sustainability challenges in land use.

The Collaborative Research Center FONDA.

Today's scientific data analysis very often requires complex Data Analysis Workflows (DAWs) executed over distributed computational infrastructures, e.g., clusters. Much research effort is devoted to the tuning and performance optimization of specific workflows for specific clusters. However, an arguably even more important problem for accelerating research is the reduction of development, adaptation, and maintenance times of DAWs. We describe the design and setup of the Collaborative Research Center (CRC) 1404 "FONDA -- Foundations of Workflows for Large-Scale Scientific Data Analysis", in which roughly 50 researchers jointly investigate new technologies, algorithms, and models to increase the portability, adaptability, and dependability of DAWs executed over distributed infrastructures. We describe the motivation behind our project, explain its underlying core concepts, introduce FONDA's internal structure, and sketch our vision for the future of workflow-based scientific data analysis. We also describe some lessons learned during the "making of" a CRC in Computer Science with strong interdisciplinary components, with the aim to foster similar endeavors.

National-scale mapping of building height using Sentinel-1 and Sentinel-2 time series.

Urban areas and their vertical characteristics have a manifold and far-reaching impact on our environment. However, openly accessible information at high spatial resolution is still missing at large for complete countries or regions. In this study, we combined Sentinel-1A/B and Sentinel-2A/B time series to map building heights for entire Germany on a 10 m grid resolving built-up structures in rural and urban contexts. We utilized information from the spectral/polarization, temporal and spatial dimensions by combining band-wise temporal aggregation statistics with morphological metrics. We trained machine learning regression models with highly accurate building height information from several 3D building models. The novelty of this method lies in the very fine resolution yet large spatial extent to which it can be applied, as well as in the use of building shadows in optical imagery. Results indicate that both radar-only and optical-only models can be used to predict building height, but the synergistic combination of both data sources leads to superior results. When testing the model against independent datasets, very consistent performance was achieved (frequency-weighted RMSE of 2.9 m to 3.5 m), which suggests that the prediction of the most frequently occurring buildings was robust. The average building height varies considerably across Germany with lower buildings in Eastern and South-Eastern Germany and taller ones along the highly urbanized areas in Western Germany. We emphasize the straightforward applicability of this approach on the national scale. It mostly relies on freely available satellite imagery and open source software, which potentially permit frequent update cycles and cost-effective mapping that may be relevant for a plethora of different applications, e.g. physical analysis of structural features or mapping society's resource usage.

Gridded population mapping for Germany based on building density, height and type from Earth Observation data using census disaggregation and bottom-up estimates.

Gridded population data is widely used to map fine scale population patterns and dynamics to understand associated human-environmental processes for global change research, disaster risk assessment and other domains. This study mapped gridded population across Germany using weighting layers from building density, building height (both from previous studies) and building type datasets, all created from freely available, temporally and globally consistent Copernicus Sentinel-1 and Sentinel-2 data. We first produced and validated a nation-wide dataset of predominant residential and non-residential building types. We then examined the impact of different weighting layers from density, type and height on top-down dasymetric mapping quality across scales. We finally performed a nation-wide bottom-up population estimate based on the three datasets. We found that integrating building types into dasymetric mapping is helpful at fine scale, as population is not redistributed to non-residential areas. Building density improved the overall quality of population estimates at all scales compared to using a binary building layer. Most importantly, we found that the combined use of density and height, i.e. volume, considerably increased mapping quality in general and with regard to regional discrepancy by largely eliminating systematic underestimation in dense agglomerations and overestimation in rural areas. We also found that building density, type and volume, together with living floor area per capita, are suitable to produce accurate large-area bottom-up population estimates.

High-Resolution Maps of Material Stocks in Buildings and Infrastructures in Austria and Germany.

The dynamics of societal material stocks such as buildings and infrastructures and their spatial patterns drive surging resource use and emissions. Two main types of data are currently used to map stocks, night-time lights (NTL) from Earth-observing (EO) satellites and cadastral information. We present an alternative approach for broad-scale material stock mapping based on freely available high-resolution EO imagery and OpenStreetMap data. Maps of built-up surface area, building height, and building types were derived from optical Sentinel-2 and radar Sentinel-1 satellite data to map patterns of material stocks for Austria and Germany. Using material intensity factors, we calculated the mass of different types of buildings and infrastructures, distinguishing eight types of materials, at 10 m spatial resolution. The total mass of buildings and infrastructures in 2018 amounted to ∼5 Gt in Austria and ∼38 Gt in Germany (AT: ∼540 t/cap, DE: ∼450 t/cap). Cross-checks with independent data sources at various scales suggested that the method may yield more complete results than other data sources but could not rule out possible overestimations. The method yields thematic differentiations not possible with NTL, avoids the use of costly cadastral data, and is suitable for mapping larger areas and tracing trends over time.

Mapping urban-rural gradients of settlements and vegetation at national scale using Sentinel-2 spectral-temporal metrics and regression-based unmixing with synthetic training data.

The increasing impact of humans on land and ongoing global population growth requires an improved understanding of land cover (LC) and land use (LU) processes related to settlements. The heterogeneity of built-up areas and infrastructures as well as the importance of not only mapping, but also characterizing anthropogenic structures suggests using a sub-pixel mapping approach for analysing related LC from space. We implement a regression-based unmixing approach for mapping built-up surfaces and infrastructure, woody vegetation and non-woody vegetation for all of Germany and Austria at 10 m resolution to demonstrate the potential of sub-pixel mapping. We map LC fractions for one point in time, using all available Sentinel-2 data from 2017 and 2018 (<70% cloud cover). We combine the concept of synthetically mixed training data with statistical aggregations from spectral-temporal metrics (STM) derived from Sentinel-2 reflectance time series. We specifically examine how STM can be used for creating synthetically mixed training data. STM are known to facilitate large area mapping by being largely independent of image acquisition dates and inherently incorporate phenological information. Vegetation is an important part of settlements and time series information supports its mapping. Synthetically mixed training data facilitates a streamlined training by using pure reference spectra to generate artificial mixtures as input to regression modelling of LC fractions in mixed pixels. We here show how combining both offers great potential for wall-to-wall LC fraction mapping. We further investigate the positive effect of STM on map results by comparing the performance of different subsets of STM combinations. Our results indicate that many STM combinations containing spectral variability and vegetation indices provide suitable input to creating synthetic training data for regression-based fraction mapping. Results for built-up surfaces and infrastructure (MAE 0.13/RMSE 0.18 at 20 m resolution), woody vegetation (0.18, 0.22) and non-woody vegetation (0.14, 0.19) are highly consistent across Germany and Austria. Only a few surface types were not accurately predicted in our nation-wide mapping. Further research is required to optimize mapping of temporally invariant bare soil and rock surfaces that show spectral similarity to built-up surfaces and infrastructure. The proposed methodology combines benefits of both regression-based modelling with synthetically mixed training data and STM, and thus facilitates mapping of LC fractions on a national scale and at high resolution. Such information will allow to better characterize settlements and identifying processes such as densification that are best represented by continuous LC mapping.

Synthesizing dam-induced land system change.

Dam construction and operation modify land systems. We synthesized 178 observations of dam-induced land system changes from 54 peer-reviewed case studies. Changing extents of forests (23%), agricultural land (21%), and built-up areas (11%) were reported frequently, alongside alterations in land use intensity (23%). Land cover changes were mostly related to hydropower and multi-purpose dams, while irrigation dams dominantly caused land use intensity changes. While a significant share of the changes was caused by reservoir flooding (29%), indirect effects which interact with societal and environmental systems (42%) were of utmost importance. We suggested the distance to the dam and the time since commissioning as potential controls for the direction of land system changes. Our insights provide opportunities for future inductive investigations across large populations of dams at regional to global scales and highlight that multi-disciplinary research perspectives are imperative for the production of generalizable knowledge.

Unravelling the link between global rubber price and tropical deforestation in Cambodia.

Tropical forests continue to undergo a rapid transformation. The expansion of rubber tree (Hevea brasiliensis) plantations has been reported as a major driver of forest loss, linked to a boom in market demand. Distant commodity markets have spurred a surge of large-scale economic land concessions granted throughout tropical Southeast Asia. Using satellite imagery, we show the impact of rubber tree plantations on Cambodian forest cover and analyse how annual forest-to-rubber conversion rates relate to global rubber prices from 2001 to 2015. We found that 23.5 ± 1.8% of national forest cover was cleared in this period, with 23.2 ± 3.6% of cleared forest converted to rubber plantations. Annual forest-to-rubber conversion rates closely correlated with global rubber prices, with a time lag of 8-9 months (Pearson's r = 0.93). Our results reveal a strong link between global commodity markets and tropical forest loss, particularly in countries with land policies geared towards rapid development.

Canopy mortality has doubled in Europe's temperate forests over the last three decades.

Mortality is a key indicator of forest health, and increasing mortality can serve as bellwether for the impacts of global change on forest ecosystems. Here we analyze trends in forest canopy mortality between 1984 and 2016 over more than 30 Mill. ha of temperate forests in Europe, based on a unique dataset of 24,000 visually interpreted spectral trajectories from the Landsat archive. On average, 0.79% of the forest area was affected by natural or human-induced mortality annually. Canopy mortality increased by +2.40% year-1, doubling the forest area affected by mortality since 1984. Areas experiencing low-severity mortality increased more strongly than areas affected by stand-replacing mortality events. Changes in climate and land-use are likely causes of large-scale forest mortality increase. Our findings reveal profound changes in recent forest dynamics with important implications for carbon storage and biodiversity conservation, highlighting the importance of improved monitoring of forest mortality.

Using Landsat time series for characterizing forest disturbance dynamics in the coupled human and natural systems of Central Europe.

Remote sensing is a key information source for improving the spatiotemporal understanding of forest ecosystem dynamics. Yet, the mapping and attribution of forest change remains challenging, particularly in areas where a number of interacting disturbance agents simultaneously affect forest development. The forest ecosystems of Central Europe are coupled human and natural systems, with natural and human disturbances affecting forests both individually and in combination. To better understand the complex forest disturbance dynamics in such systems, we utilize 32-year Landsat time series to map forest disturbances in five sites across Austria, the Czech Republic, Germany, Poland, and Slovakia. All sites consisted of a National Park and the surrounding forests, reflecting three management zones of different levels of human influence (managed, protected, strictly protected). This allowed for a comparison of spectral, temporal, and spatial disturbance patterns across a gradient from natural to coupled human and natural disturbances. Disturbance maps achieved overall accuracies ranging from 81% to 93%. Disturbance patches were generally small, with 95% of the disturbances being smaller than 10 ha. Disturbance rates ranged from 0.29% yr-1 to 0.95% yr-1, and differed substantially among management zones and study sites. Natural disturbances in strictly protected areas were longer in duration (median of 8 years) and slightly less variable in magnitude compared to human-dominated disturbances in managed forests (median duration of 1 year). However, temporal dynamics between natural and human-dominated disturbances showed strong synchrony, suggesting that disturbance peaks are driven by natural events affecting managed and unmanaged areas simultaneously. Our study demonstrates the potential of remote sensing for mapping forest disturbances in coupled human and natural systems, such as the forests of Central Europe. Yet, we also highlight the complexity of such systems in terms of agent attribution, as many natural disturbances are modified by management responding to them outside protected areas.

Remote sensing of forest insect disturbances: Current state and future directions.

Insect disturbance are important agents of change in forest ecosystems around the globe, yet their spatial and temporal distribution and dynamics are not well understood. Remote sensing has gained much attention in mapping and understanding insect outbreak dynamics. Consequently, we here review the current literature on the remote sensing of insect disturbances. We suggest to group studies into three insect types: bark beetles, broadleaved defoliators, and coniferous defoliators. By so doing, we systematically compare the sensors and methods used for mapping insect disturbances within and across insect types. Results suggest that there are substantial differences between methods used for mapping bark beetles and defoliators, and between methods used for mapping broadleaved and coniferous defoliators. Following from this, we highlight approaches that are particularly suited for each insect type. Finally, we conclude by highlighting future research directions for remote sensing of insect disturbances. In particular, we suggest to: 1) Separate insect disturbances from other agents; 2) Extend the spatial and temporal domain of analysis; 3) Make use of dense time series; 4) Operationalize near-real time monitoring of insect disturbances; 5) Identify insect disturbances in the context of coupled human-natural systems; and 6) Improve reference data for assessing insect disturbances. Since the remote sensing of insect disturbances has gained much interest beyond the remote sensing community recently, the future developments identified here will help integrating remote sensing products into operational forest management. Furthermore, an improved spatiotemporal quantification of insect disturbances will support an inclusion of these processes into regional to global ecosystem models.

Carbon emissions from agricultural expansion and intensification in the Chaco.

Carbon emissions from land-use changes in tropical dry forest systems are poorly understood, although they are likely globally significant. The South American Chaco has recently emerged as a hot spot of agricultural expansion and intensification, as cattle ranching and soybean cultivation expand into forests, and as soybean cultivation replaces grazing lands. Still, our knowledge of the rates and spatial patterns of these land-use changes and how they affected carbon emissions remains partial. We used the Landsat satellite image archive to reconstruct land-use change over the past 30 years and applied a carbon bookkeeping model to quantify how these changes affected carbon budgets. Between 1985 and 2013, more than 142 000 km2 of the Chaco's forests, equaling 20% of all forest, was replaced by croplands (38.9%) or grazing lands (61.1%). Of those grazing lands that existed in 1985, about 40% were subsequently converted to cropland. These land-use changes resulted in substantial carbon emissions, totaling 824 Tg C between 1985 and 2013, and 46.2 Tg C for 2013 alone. The majority of these emissions came from forest-to-grazing-land conversions (68%), but post-deforestation land-use change triggered an additional 52.6 Tg C. Although tropical dry forests are less carbon-dense than moist tropical forests, carbon emissions from land-use change in the Chaco were similar in magnitude to those from other major tropical deforestation frontiers. Our study thus highlights the urgent need for an improved monitoring of the often overlooked tropical dry forests and savannas, and more broadly speaking the value of the Landsat image archive for quantifying carbon fluxes from land change.

Using fragmentation to assess degradation of forest edges in Democratic Republic of Congo.

BACKGROUND: Recent studies have shown that fragmentation is an increasing threat to global forests, which has major impacts on biodiversity and the important ecosystem services provided by forested landscapes. Several tools have been developed to evaluate global patterns of fragmentation, which have potential applications for REDD+. We study how canopy height and above ground biomass (AGB) change across several categories of forest edges determined by fragmentation analysis. We use Democratic Republic of Congo (DRC) as an example. RESULTS: An analysis of variance of different edge widths and airborne estimated canopy height found that canopy heights were significantly different in forest edges at a distance of 100 m from the nonforest edge. Biomass was significantly different between fragmentation classes at an edge distance of 300 m. Core forest types were found to have significantly higher canopy height and greater AGB than forest edges and patches, where height and biomass decrease significantly as the level of fragmentation increases. A change analysis shows that deforestation and degradation are increasing over time and biomass loss associated with degradation account for at least one quarter of total loss. We estimate that about 80 % of primary forests are intact, which decreases 3.5 % over the 15 year study period, as primary forest is either deforested or transitioned to forest edge. While the carbon loss per hectare is lower than that of deforestation, degradation potentially affects up to three times more area than deforestation alone. CONCLUSIONS: When defining forest degradation by decreased biomass without any loss in forest area, assessing transitions of core forest to edges over time can contribute an important element to REDD+MRV systems. The estimation of changes between different forest fragmentation types and their associated biomass loss can provide an estimate of degradation carbon emission factors. Forest degradation and emissions due to fragmentation are often underestimated and should comprise an essential component of MRV systems.

Modelling avian biodiversity using raw, unclassified satellite imagery.

Applications of remote sensing for biodiversity conservation typically rely on image classifications that do not capture variability within coarse land cover classes. Here, we compare two measures derived from unclassified remotely sensed data, a measure of habitat heterogeneity and a measure of habitat composition, for explaining bird species richness and the spatial distribution of 10 species in a semi-arid landscape of New Mexico. We surveyed bird abundance from 1996 to 1998 at 42 plots located in the McGregor Range of Fort Bliss Army Reserve. Normalized Difference Vegetation Index values of two May 1997 Landsat scenes were the basis for among-pixel habitat heterogeneity (image texture), and we used the raw imagery to decompose each pixel into different habitat components (spectral mixture analysis). We used model averaging to relate measures of avian biodiversity to measures of image texture and spectral mixture analysis fractions. Measures of habitat heterogeneity, particularly angular second moment and standard deviation, provide higher explanatory power for bird species richness and the abundance of most species than measures of habitat composition. Using image texture, alone or in combination with other classified imagery-based approaches, for monitoring statuses and trends in biological diversity can greatly improve conservation efforts and habitat management.

A high-resolution approach to estimating ecosystem respiration at continental scales using operational satellite data.

A better understanding of the local variability in land-atmosphere carbon fluxes is crucial to improving the accuracy of global carbon budgets. Operational satellite data backed by ground measurements at Fluxnet sites proved valuable in monitoring local variability of gross primary production at highly resolved spatio-temporal resolutions. Yet, we lack similar operational estimates of ecosystem respiration (Re) to calculate net carbon fluxes. If successful, carbon fluxes from such a remote sensing approach would form an independent and sought after measure to complement widely used dynamic global vegetation models (DGVMs). Here, we establish an operational semi-empirical Re model, based only on data from the Moderate Resolution Imaging Spectroradiometer (MODIS) with a resolution of 1 km and 8 days. Fluxnet measurements between 2000 and 2009 from 100 sites across North America and Europe are used for parameterization and validation. Our analysis shows that Re is closely tied to temperature and plant productivity. By separating temporal and intersite variation, we find that MODIS land surface temperature (LST) and enhanced vegetation index (EVI) are sufficient to explain observed Re across most major biomes with a negligible bias [R² = 0.62, RMSE = 1.32 (g C m(-2) d(-1)), MBE = 0.05 (g C m(-2) d(-1))]. A comparison of such satellite-derived Re with those simulated by the DGVM LPJmL reveals similar spatial patterns. However, LPJmL shows higher temperature sensitivities and consistently simulates higher Re values, in high-latitude and subtropical regions. These differences remain difficult to explain and they are likely associated either with LPJmL parameterization or with systematic errors in the Fluxnet sampling technique. While uncertainties remain with Re estimates, the model formulated in this study provides an operational, cross-validated and unbiased approach to scale Fluxnet Re to the continental scale and advances knowledge of spatio-temporal Re variability.