Deep learning estimation of northern hemisphere soil freeze-thaw dynamics using satellite multi-frequency microwave brightness temperature observations

Satellite microwave sensors are well suited for monitoring landscape freeze-thaw (FT) transitions owing to the strong brightness temperature (TB) or backscatter response to changes in liquid water abundance between predominantly frozen and thawed conditions. The FT retrieval is also a sensitive climate indicator with strong biophysical importance. However, retrieval algorithms can have difficulty distinguishing the FT status of soils from that of overlying features such as snow and vegetation, while variable land conditions can also degrade performance. Here, we applied a deep learning model using a multilayer convolutional neural network driven by AMSR2 and SMAP TB records, and trained on surface (~0-5 cm depth) soil temperature FT observations. Soil FT states were classified for the local morning (6 a.m.) and evening (6 p.m.) conditions corresponding to SMAP descending and ascending orbital overpasses, mapped to a 9 km polar grid spanning a five-year (2016-2020) record and Northern Hemisphere domain. Continuous variable estimates of the probability of frozen or thawed conditions were derived using a model cost function optimized against FT observational training data. Model results derived using combined multi-frequency (1.4, 18.7, 36.5 GHz) TBs produced the highest soil FT accuracy over other models derived using only single sensor or single frequency TB inputs. Moreover, SMAP L-band (1.4 GHz) TBs provided enhanced soil FT information and performance gain over model results derived using only AMSR2 TB inputs. The resulting soil FT classification showed favorable and consistent performance against soil FT observations from ERA5 reanalysis (mean percent accuracy, MPA: 92.7%) and in situ weather stations (MPA: 91.0%). The soil FT accuracy was generally consistent between morning and afternoon predictions and across different land covers and seasons. The model also showed better FT accuracy than ERA5 against regional weather station measurements (91.0% vs. 86.1% MPA). However, model confidence was lower in complex terrain where FT spatial heterogeneity was likely beneath the effective model grain size. Our results provide a high level of precision in mapping soil FT dynamics to improve understanding of complex seasonal transitions and their influence on ecological processes and climate feedbacks, with the potential to inform Earth system model predictions.

Carbon uptake in Eurasian boreal forests dominates the high-latitude net ecosystem carbon budget

Arctic-boreal landscapes are experiencing profound warming, along with changes in ecosystem moisture status and disturbance from fire. This region is of global importance in terms of carbon feedbacks to climate, yet the sign (sink or source) and magnitude of the Arctic-boreal carbon budget within recent years remains highly uncertain. Here, we provide new estimates of recent (2003-2015) vegetation gross primary productivity (GPP), ecosystem respiration (R_eco ), net ecosystem CO₂ exchange (NEE; R_eco  - GPP), and terrestrial methane (CH₄ ) emissions for the Arctic-boreal zone using a satellite data-driven process-model for northern ecosystems (TCFM-Arctic), calibrated and evaluated using measurements from >60 tower eddy covariance (EC) sites. We used TCFM-Arctic to obtain daily 1-km² flux estimates and annual carbon budgets for the pan-Arctic-boreal region. Across the domain, the model indicated an overall average NEE sink of -850 Tg CO₂ -C year^-1 . Eurasian boreal zones, especially those in Siberia, contributed to a majority of the net sink. In contrast, the tundra biome was relatively carbon neutral (ranging from small sink to source). Regional CH₄ emissions from tundra and boreal wetlands (not accounting for aquatic CH₄ ) were estimated at 35 Tg CH₄ -C year^-1 . Accounting for additional emissions from open water aquatic bodies and from fire, using available estimates from the literature, reduced the total regional NEE sink by 21% and shifted many far northern tundra landscapes, and some boreal forests, to a net carbon source. This assessment, based on in situ observations and models, improves our understanding of the high-latitude carbon status and also indicates a continued need for integrated site-to-regional assessments to monitor the vulnerability of these ecosystems to climate change.

Pan-Arctic soil moisture control on tundra carbon sequestration and plant productivity

Long-term atmospheric CO₂ concentration records have suggested a reduction in the positive effect of warming on high-latitude carbon uptake since the 1990s. A variety of mechanisms have been proposed to explain the reduced net carbon sink of northern ecosystems with increased air temperature, including water stress on vegetation and increased respiration over recent decades. However, the lack of consistent long-term carbon flux and in situ soil moisture data has severely limited our ability to identify the mechanisms responsible for the recent reduced carbon sink strength. In this study, we used a record of nearly 100 site-years of eddy covariance data from 11 continuous permafrost tundra sites distributed across the circumpolar Arctic to test the temperature (expressed as growing degree days, GDD) responses of gross primary production (GPP), net ecosystem exchange (NEE), and ecosystem respiration (ER) at different periods of the summer (early, peak, and late summer) including dominant tundra vegetation classes (graminoids and mosses, and shrubs). We further tested GPP, NEE, and ER relationships with soil moisture and vapor pressure deficit to identify potential moisture limitations on plant productivity and net carbon exchange. Our results show a decrease in GPP with rising GDD during the peak summer (July) for both vegetation classes, and a significant relationship between the peak summer GPP and soil moisture after statistically controlling for GDD in a partial correlation analysis. These results suggest that tundra ecosystems might not benefit from increased temperature as much as suggested by several terrestrial biosphere models, if decreased soil moisture limits the peak summer plant productivity, reducing the ability of these ecosystems to sequester carbon during the summer.

Respiratory loss during late-growing season determines the net carbon dioxide sink in northern permafrost regions

Warming of northern high latitude regions (NHL, > 50 °N) has increased both photosynthesis and respiration which results in considerable uncertainty regarding the net carbon dioxide (CO₂) balance of NHL ecosystems. Using estimates constrained from atmospheric observations from 1980 to 2017, we find that the increasing trends of net CO₂ uptake in the early-growing season are of similar magnitude across the tree cover gradient in the NHL. However, the trend of respiratory CO₂ loss during late-growing season increases significantly with increasing tree cover, offsetting a larger fraction of photosynthetic CO₂ uptake, and thus resulting in a slower rate of increasing annual net CO₂ uptake in areas with higher tree cover, especially in central and southern boreal forest regions. The magnitude of this seasonal compensation effect explains the difference in net CO₂ uptake trends along the NHL vegetation- permafrost gradient. Such seasonal compensation dynamics are not captured by dynamic global vegetation models, which simulate weaker respiration control on carbon exchange during the late-growing season, and thus calls into question projections of increasing net CO₂ uptake as high latitude ecosystems respond to warming climate conditions.

The northern high latitude permafrost region has been an important contributor to the carbon sink since the 1980s. A new study finds that as tree cover increases, respiratory CO2 loss during late-growing season offsets photosynthetic CO2 uptake and leads to a slower rate of increasing annual net CO2 uptake.

Exceptional heat and atmospheric dryness amplified losses of primary production during the 2020 U.S. Southwest hot drought

Earth's ecosystems are increasingly threatened by "hot drought," which occurs when hot air temperatures coincide with precipitation deficits, intensifying the hydrological, physiological, and ecological effects of drought by enhancing evaporative losses of soil moisture (SM) and increasing plant stress due to higher vapor pressure deficit (VPD). Drought-induced reductions in gross primary production (GPP) exert a major influence on the terrestrial carbon sink, but the extent to which hotter and atmospherically drier conditions will amplify the effects of precipitation deficits on Earth's carbon cycle remains largely unknown. During summer and autumn 2020, the U.S. Southwest experienced one of the most intense hot droughts on record, with record-low precipitation and record-high air temperature and VPD across the region. Here, we use this natural experiment to evaluate the effects of hot drought on GPP and further decompose those negative GPP anomalies into their constituent meteorological and hydrological drivers. We found a 122 Tg C (>25%) reduction in GPP below the 2015-2019 mean, by far the lowest regional GPP over the Soil Moisture Active Passive satellite record. Roughly half of the estimated GPP loss was attributable to low SM (likely a combination of record-low precipitation and warming-enhanced evaporative depletion), but record-breaking VPD amplified the reduction of GPP, contributing roughly 40% of the GPP anomaly. Both air temperature and VPD are very likely to continue increasing over the next century, likely leading to more frequent and intense hot droughts and substantially enhancing drought-induced GPP reductions.

Evaluation of trophic state for inland waters through combining Forel-Ule Index and inherent optical properties

Eutrophication is a severe environmental pollution problem for inland waters and poses significant threats to the water safety. Monitoring trophic state of inland waters using optical remote sensing generally requires the inversion of water quality parameters, such as chlorophyll-a, secchi depth, etc. However, the accurate inversion of these individual indicators remains challenging, while the associated retrieval errors can propagate and degrade the evaluation of trophic state. Hence, we proposed a novel monitoring method by developing a Trophic State Index (TSI) based on optical remote-sensing parameters, i.e., Forel-Ule index (FUI) and non-water absorption coefficient at 674 nm (referred to as a_t-w(674)) retrieved from Sentinel-3 Ocean and Land Color Instrument (OLCI) imagery. The estimated TSI showed favorable correspondence with observed water quality data, including coefficient of determination (r² = 0.91), root mean squared error (RMSE = 5.54), and mean absolute percentage error (MAPE = 10.69%). Using the Sentinel-3 OLCI data, the proposed method also had very good performance in the field spectrum (MAPE = 5.25 % , RMSE = 3.36). The monthly trophic state evaluation also showed congruence (MAPE = 12.51 % , RMSE = 6.41) with surface water quality monthly report (SWQMR) from the Ministry of Environment and Ecology of the People's Republic of China. The monthly TSI showed favorable agreement for 23 ungauged lakes (RMSE = 7.26, MAPE = 12.78%), indicating potential utility for regional lake water quality monitoring. The proposed method was then applied to 47 other large (>50 km²) water bodies in the Middle-and-Lower watershed of Yangtze River and the Huaihe watershed to evaluate the spatial and temporal variation of trophic state from 2016 to 2020. The TSI results revealed several lakes, such as Lake Honghu and Lake Luoma, with rapidly deteriorating water quality during the study period, while other lakes show relative improvement (e.g., Xiashan Reservoir), indicating unbalanced environmental pressure over the region. Overall, this study showed promising performance and potential for satellite-based monitoring of regional aquatic environments.

A framework to integrate innovations in invasion science for proactive management

Invasive alien species (IAS) are a rising threat to biodiversity, national security, and regional economies, with impacts in the hundreds of billions of U.S. dollars annually. Proactive or predictive approaches guided by scientific knowledge are essential to keeping pace with growing impacts of invasions under climate change. Although the rapid development of diverse technologies and approaches has produced tools with the potential to greatly accelerate invasion research and management, innovation has far outpaced implementation and coordination. Technological and methodological syntheses are urgently needed to close the growing implementation gap and facilitate interdisciplinary collaboration and synergy among evolving disciplines. A broad review is necessary to demonstrate the utility and relevance of work in diverse fields to generate actionable science for the ongoing invasion crisis. Here, we review such advances in relevant fields including remote sensing, epidemiology, big data analytics, environmental DNA (eDNA) sampling, genomics, and others, and present a generalized framework for distilling existing and emerging data into products for proactive IAS research and management. This integrated workflow provides a pathway for scientists and practitioners in diverse disciplines to contribute to applied invasion biology in a coordinated, synergistic, and scalable manner.

DroughtCast: A Machine Learning Forecast of the United States Drought Monitor

Drought is one of the most ecologically and economically devastating natural phenomena affecting the United States, causing the U.S. economy billions of dollars in damage, and driving widespread degradation of ecosystem health. Many drought indices are implemented to monitor the current extent and status of drought so stakeholders such as farmers and local governments can appropriately respond. Methods to forecast drought conditions weeks to months in advance are less common but would provide a more effective early warning system to enhance drought response, mitigation, and adaptation planning. To resolve this issue, we introduce DroughtCast, a machine learning framework for forecasting the United States Drought Monitor (USDM). DroughtCast operates on the knowledge that recent anomalies in hydrology and meteorology drive future changes in drought conditions. We use simulated meteorology and satellite observed soil moisture as inputs into a recurrent neural network to accurately forecast the USDM between 1 and 12 weeks into the future. Our analysis shows that precipitation, soil moisture, and temperature are the most important input variables when forecasting future drought conditions. Additionally, a case study of the 2017 Northern Plains Flash Drought shows that DroughtCast was able to forecast a very extreme drought event up to 12 weeks before its onset. Given the favorable forecasting skill of the model, DroughtCast may provide a promising tool for land managers and local governments in preparing for and mitigating the effects of drought.

Testing a Generalizable Machine Learning Workflow for Aquatic Invasive Species on Rainbow Trout (Oncorhynchus mykiss) in Northwest Montana

Biological invasions are accelerating worldwide, causing major ecological and economic impacts in aquatic ecosystems. The urgent decision-making needs of invasive species managers can be better met by the integration of biodiversity big data with large-domain models and data-driven products. Remotely sensed data products can be combined with existing invasive species occurrence data via machine learning models to provide the proactive spatial risk analysis necessary for implementing coordinated and agile management paradigms across large scales. We present a workflow that generates rapid spatial risk assessments on aquatic invasive species using occurrence data, spatially explicit environmental data, and an ensemble approach to species distribution modeling using five machine learning algorithms. For proof of concept and validation, we tested this workflow using extensive spatial and temporal hybridization and occurrence data from a well-studied, ongoing, and climate-driven species invasion in the upper Flathead River system in northwestern Montana, USA. Rainbow Trout (RBT; Oncorhynchus mykiss), an introduced species in the Flathead River basin, compete and readily hybridize with native Westslope Cutthroat Trout (WCT; O. clarkii lewisii), and the spread of RBT individuals and their alleles has been tracked for decades. We used remotely sensed and other geospatial data as key environmental predictors for projecting resultant habitat suitability to geographic space. The ensemble modeling technique yielded high accuracy predictions relative to 30-fold cross-validated datasets (87% 30-fold cross-validated accuracy score). Both top predictors and model performance relative to these predictors matched current understanding of the drivers of RBT invasion and habitat suitability, indicating that temperature is a major factor influencing the spread of invasive RBT and hybridization with native WCT. The congruence between more time-consuming modeling approaches and our rapid machine-learning approach suggest that this workflow could be applied more broadly to provide data-driven management information for early detection of potential invaders.

Satellite Flood Inundation Assessment and Forecast Using SMAP and Landsat

The capability and synergistic use of multisource satellite observations for flood monitoring and forecasts is crucial for improving disaster preparedness and mitigation. Here, surface fractional water cover (FW) retrievals derived from Soil Moisture Active Passive (SMAP) L-band (1.4 GHz) brightness temperatures were used for flood assessment over southeast Africa during the Cyclone Idai event. We then focused on five subcatchments of the Pungwe basin and developed a machine learning based approach with the support of Google Earth Engine for daily (24-h) forecasting of FW and 30-m inundation downscaling and mapping. The Classification and Regression Trees model was selected and trained using retrievals derived from SMAP and Landsat coupled with rainfall forecasts from the NOAA Global Forecast System. Independent validation showed that FW predictions over randomly selected dates are highly correlated (R = 0.87) with the Landsat observations. The forecast results captured the flood temporal dynamics from the Idai event; and the associated 30-m downscaling results showed inundation spatial patterns consistent with independent satellite synthetic aperture radar observations. The data-driven approach provides new capacity for flood monitoring and forecasts leveraging synergistic satellite observations and big data analysis, which is particularly valuable for data sparse regions.

Monitoring Crop Status in the Continental United States Using the SMAP Level-4 Carbon Product

Accurate monitoring of crop condition is critical to detect anomalies that may threaten the economic viability of agriculture and to understand how crops respond to climatic variability. Retrievals of soil moisture and vegetation information from satellite-based remote-sensing products offer an opportunity for continuous and affordable crop condition monitoring. This study compared weekly anomalies in accumulated gross primary production (GPP) from the SMAP Level-4 Carbon (L4C) product to anomalies calculated from a state-scale weekly crop condition index (CCI) and also to crop yield anomalies calculated from county-level yield data reported at the end of the season. We focused on barley, spring wheat, corn, and soybeans cultivated in the continental United States from 2000 to 2018. We found that consistencies between SMAP L4C GPP anomalies and both crop condition and yield anomalies increased as crops developed from the emergence stage (r: 0.4–0.7) and matured (r: 0.6–0.9) and that the agreement was better in drier regions (r: 0.4–0.9) than in wetter regions (r: −0.8–0.4). The L4C provides weekly GPP estimates at a 1-km scale, permitting the evaluation and tracking of anomalies in crop status at higher spatial detail than metrics based on the state-level CCI or county-level crop yields. We demonstrate that the L4C GPP product can be used operationally to monitor crop condition with the potential to become an important tool to inform decision-making and research.

Interannual variability of ecosystem iso/anisohydry is regulated by environmental dryness

●Plants are characterized by the iso/anisohydry continuum depending on how they regulate leaf water potential (Ψ_L ). However, how iso/anisohydry changes over time in response to year-to-year variations in environmental dryness and how such responses vary across different regions remains poorly characterized. ●We investigated how dryness, represented by aridity index, affects the interannual variability of ecosystem iso/anisohydry at the regional scale, estimated using satellite microwave vegetation optical depth (VOD) observations. This ecosystem-level analysis was further complemented with published field observations of species-level Ψ_L . ●We found different behaviors in the directionality and sensitivity of isohydricity (σ) with respect to the interannual variation of dryness in different ecosystems. These behaviors can largely be differentiated by the average dryness of the ecosystem itself: in mesic ecosystems, σ decreases in drier years with a higher sensitivity to dryness; in xeric ecosystems, σ increases in drier years with a lower sensitivity to dryness. These results were supported by the species-level synthesis. ●Our study suggests that how plants adjust their water use across years - as revealed by their interannual variability in isohydricity - depends on the dryness of plants' living environment. This finding advances our understanding of plant responses to drought at regional scales.

Increasing riverine heat influx triggers Arctic sea ice decline and oceanic and atmospheric warming

Arctic river discharge increased over the last several decades, conveying heat and freshwater into the Arctic Ocean and likely affecting regional sea ice and the ocean heat budget. However, until now, there have been only limited assessments of riverine heat impacts. Here, we adopted a synthesis of a pan-Arctic sea ice-ocean model and a land surface model to quantify impacts of river heat on the Arctic sea ice and ocean heat budget. We show that river heat contributed up to 10% of the regional sea ice reduction over the Arctic shelves from 1980 to 2015. Particularly notable, this effect occurs as earlier sea ice breakup in late spring and early summer. The increasing ice-free area in the shelf seas results in a warmer ocean in summer, enhancing ocean-atmosphere energy exchange and atmospheric warming. Our findings suggest that a positive river heat-sea ice feedback nearly doubles the river heat effect.

Increased high-latitude photosynthetic carbon gain offset by respiration carbon loss during an anomalous warm winter to spring transition

Arctic and boreal ecosystems play an important role in the global carbon (C) budget, and whether they act as a future net C sink or source depends on climate and environmental change. Here, we used complementary in situ measurements, model simulations, and satellite observations to investigate the net carbon dioxide (CO₂ ) seasonal cycle and its climatic and environmental controls across Alaska and northwestern Canada during the anomalously warm winter to spring conditions of 2015 and 2016 (relative to 2010-2014). In the warm spring, we found that photosynthesis was enhanced more than respiration, leading to greater CO₂ uptake. However, photosynthetic enhancement from spring warming was partially offset by greater ecosystem respiration during the preceding anomalously warm winter, resulting in nearly neutral effects on the annual net CO₂ balance. Eddy covariance CO₂ flux measurements showed that air temperature has a primary influence on net CO₂ exchange in winter and spring, while soil moisture has a primary control on net CO₂ exchange in the fall. The net CO₂ exchange was generally more moisture limited in the boreal region than in the Arctic tundra. Our analysis indicates complex seasonal interactions of underlying C cycle processes in response to changing climate and hydrology that may not manifest in changes in net annual CO₂ exchange. Therefore, a better understanding of the seasonal response of C cycle processes may provide important insights for predicting future carbon-climate feedbacks and their consequences on atmospheric CO₂ dynamics in the northern high latitudes.

Assessing global surface water inundation dynamics using combined satellite information from SMAP, AMSR2 and Landsat

A method to assess global land surface water (fw) inundation dynamics was developed by exploiting the enhanced fw sensitivity of L-band (1.4 GHz) passive microwave observations from the Soil Moisture Active Passive (SMAP) mission. The L-band fw (fwLBand ) retrievals were derived using SMAP H-polarization brightness temperature (Tb ) observations and predefined L-band reference microwave emissivities for water and land endmembers. Potential soil moisture and vegetation contributions to the microwave signal were represented from overlapping higher frequency Tb observations from AMSR2. The resulting fwLBand global record has high temporal sampling (1-3 days) and 36-km spatial resolution. The fwLBand annual averages corresponded favourably (R=0.85, p-value<0.001) with a 250-m resolution static global water map (MOD44W) aggregated at the same spatial scale, while capturing significant inundation variations worldwide. The monthly fwLBand averages also showed seasonal inundation changes consistent with river discharge records within six major US river basins. An uncertainty analysis indicated generally reliable fwLBand performance for major land cover areas and under low to moderate vegetation cover, but with lower accuracy for detecting water bodies covered by dense vegetation. Finer resolution (30-m) fwLBand results were obtained for three sub-regions in North America using an empirical downscaling approach and ancillary global Water Occurrence Dataset (WOD) derived from the historical Landsat record. The resulting 30-m fwLBand retrievals showed favourable spatial accuracy for water (commission error 31.46%, omission error 30.20%) and land (commission error 0.87%, omission error 0.96%) classifications and seasonal wet and dry periods when compared to independent water maps derived from Landsat-8 imagery. The new fwLBand algorithms and continuing SMAP and AMSR2 operations provide for near real-time, multi-scale monitoring of global surface water inundation dynamics and potential flood risk.

Global Satellite Retrievals of the Near-Surface Atmospheric Vapor Pressure Deficit from AMSR-E and AMSR2

Near-surface atmospheric Vapor Pressure Deficit (VPD) is a key environmental variable affecting vegetation water stress, evapotranspiration, and atmospheric moisture demand. Although VPD is readily derived from in situ standard weather station measurements, more spatially continuous global observations for regional monitoring of VPD are lacking. Here, we document a new method to estimate daily (both a.m. and p.m.) global land surface VPD at a 25-km resolution using a satellite passive microwave remotely sensed Land Parameter Data Record (LPDR) derived from the Advanced Microwave Scanning Radiometer (AMSR) sensors. The AMSR-derived VPD record shows strong correspondence (correlation coefficient ≥ 0.80, p-value < 0.001) and overall good performance (0.48 kPa ≤ Root Mean Square Error ≤ 0.69 kPa) against independent VPD observations from the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) data. The estimated AMSR VPD retrieval uncertainties vary with land cover type, satellite observation time, and underlying LPDR data quality. These results provide new satellite capabilities for global mapping and monitoring of land surface VPD dynamics from ongoing AMSR2 operations. Overall good accuracy and similar observations from both AMSR2 and AMSR-E allow for the development of climate data records documenting recent (from 2002) VPD trends and potential impacts on vegetation, land surface evaporation, and energy budgets.

Future global productivity will be affected by plant trait response to climate

Plant traits are both responsive to local climate and strong predictors of primary productivity. We hypothesized that future climate change might promote a shift in global plant traits resulting in changes in Gross Primary Productivity (GPP). We characterized the relationship between key plant traits, namely Specific Leaf Area (SLA), height, and seed mass, and local climate and primary productivity. We found that by 2070, tropical and arid ecosystems will be more suitable for plants with relatively lower canopy height, SLA and seed mass, while far northern latitudes will favor woody and taller plants than at present. Using a network of tower eddy covariance CO2 flux measurements and the extrapolated plant trait maps, we estimated the global distribution of annual GPP under current and projected future plant community distribution. We predict that annual GPP in northern biomes (≥45 °N) will increase by 31% (+8.1 ± 0.5 Pg C), but this will be offset by a 17.9% GPP decline in the tropics (-11.8 ± 0.84 Pg C). These findings suggest that regional climate changes will affect plant trait distributions, which may in turn affect global productivity patterns.

Characterizing permafrost active layer dynamics and sensitivity to landscape spatial heterogeneity in Alaska

An important feature of the Arctic is large spatial heterogeneity in active layer conditions, which is generally poorly represented by global models and can lead to large uncertainties in predicting regional ecosystem responses and climate feedbacks. In this study, we developed a spatially integrated modelling and analysis framework combining field observations, local scale (~ 50 m resolution) active layer thickness (ALT) and soil moisture maps derived from airborne low frequency (L+P-band) radar measurements, and global satellite environmental observations to investigate the ALT sensitivity to recent climate trends and landscape heterogeneity in Alaska. Modelled ALT results show good correspondence with in situ measurements in higher permafrost probability (PP ≥ 70%) areas (n = 33, R = 0.60, mean bias = 1.58 cm, RMSE = 20.32 cm), but with larger uncertainty in sporadic and discontinuous permafrost areas. The model results also reveal widespread ALT deepening since 2001, with smaller ALT increases in northern Alaska (mean trend = 0.32 ± 1.18 cm yr^-1) and much larger increases (> 3 cm yr^-1) across interior and southern Alaska. The positive ALT trend coincides with regional warming and a longer snow-free season (R = 0.60 ± 0.32). A spatially integrated analysis of the radar retrievals and model sensitivity simulations demonstrated that uncertainty in the spatial and vertical distribution of soil organic carbon (SOC) was the largest factor affecting modeled ALT accuracy, while soil moisture played a secondary role. Potential improvements in characterizing SOC heterogeneity, including better spatial sampling of soil conditions and advances in remote sensing of SOC and soil moisture, will enable more accurate predictions of active layer conditions and refinement of the modelling framework across a larger domain.

Global Assessment of the SMAP Level-4 Surface and Root-Zone Soil Moisture Product Using Assimilation Diagnostics

The Soil Moisture Active Passive (SMAP) mission Level-4 Soil Moisture (L4_SM) product provides 3-hourly, 9-km resolution, global estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture and related land surface variables from 31 March 2015 to present with ~2.5day latency. The ensemble-based L4_SM algorithm assimilates SMAP brightness temperature (Tb) observations into the Catchment land surface model. This study describes the spatially distributed L4_SM analysis and assesses the observation-minus-forecast (O-F) Tb residuals and the soil moisture and temperature analysis increments. Owing to the climatological rescaling of the Tb observations prior to assimilation, the analysis is essentially unbiased, with global mean values of ~0.37 K for the O-F Tb residuals and practically zero for the soil moisture and temperature increments. There are, however, modest regional (absolute) biases in the O-F residuals (under ~3 K), the soil moisture increments (under ~0.01 m³ m^-3), and the surface soil temperature increments (under ~1 K). Typical instantaneous values are ~6 K for O-F residuals, ~0.01 (~0.003) m³ m^-3 for surface (root-zone) soil moisture increments, and ~0.6 K for surface soil temperature increments. The O-F diagnostics indicate that the actual errors in the system are overestimated in deserts and densely vegetated regions and underestimated in agricultural regions and transition zones between dry and wet climates. The O-F auto-correlations suggest that the SMAP observations are used efficiently in western North America, the Sahel, and Australia, but not in many forested regions and the high northern latitudes. A case study in Australia demonstrates that assimilating SMAP observations successfully corrects short-term errors in the L4_SM rainfall forcing.

Global Assessment of the SMAP Level-4 Surface and Root-Zone Soil Moisture Product Using Assimilation Diagnostics

The Soil Moisture Active Passive (SMAP) mission Level-4 Soil Moisture (L4_SM) product provides 3-hourly, 9-km resolution, global estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture and related land surface variables from 31 March 2015 to present with ~2.5day latency. The ensemble-based L4_SM algorithm assimilates SMAP brightness temperature (Tb) observations into the Catchment land surface model. This study describes the spatially distributed L4_SM analysis and assesses the observation-minus-forecast (O-F) Tb residuals and the soil moisture and temperature analysis increments. Owing to the climatological rescaling of the Tb observations prior to assimilation, the analysis is essentially unbiased, with global mean values of ~0.37 K for the O-F Tb residuals and practically zero for the soil moisture and temperature increments. There are, however, modest regional (absolute) biases in the O-F residuals (under ~3 K), the soil moisture increments (under ~0.01 m3 m-3), and the surface soil temperature increments (under ~1 K). Typical instantaneous values are ~6 K for O-F residuals, ~0.01 (~0.003) m3 m-3 for surface (root-zone) soil moisture increments, and ~0.6 K for surface soil temperature increments. The O-F diagnostics indicate that the actual errors in the system are overestimated in deserts and densely vegetated regions and underestimated in agricultural regions and transition zones between dry and wet climates. The O-F auto-correlations suggest that the SMAP observations are used efficiently in western North America, the Sahel, and Australia, but not in many forested regions and the high northern latitudes. A case study in Australia demonstrates that assimilating SMAP observations successfully corrects short-term errors in the L4_SM rainfall forcing.

Climate variables explain neutral and adaptive variation within salmonid metapopulations: the importance of replication in landscape genetics

Understanding how environmental variation influences population genetic structure is important for conservation management because it can reveal how human stressors influence population connectivity, genetic diversity and persistence. We used riverscape genetics modelling to assess whether climatic and habitat variables were related to neutral and adaptive patterns of genetic differentiation (population-specific and pairwise FST ) within five metapopulations (79 populations, 4583 individuals) of steelhead trout (Oncorhynchus mykiss) in the Columbia River Basin, USA. Using 151 putatively neutral and 29 candidate adaptive SNP loci, we found that climate-related variables (winter precipitation, summer maximum temperature, winter highest 5% flow events and summer mean flow) best explained neutral and adaptive patterns of genetic differentiation within metapopulations, suggesting that climatic variation likely influences both demography (neutral variation) and local adaptation (adaptive variation). However, we did not observe consistent relationships between climate variables and FST across all metapopulations, underscoring the need for replication when extrapolating results from one scale to another (e.g. basin-wide to the metapopulation scale). Sensitivity analysis (leave-one-population-out) revealed consistent relationships between climate variables and FST within three metapopulations; however, these patterns were not consistent in two metapopulations likely due to small sample sizes (N = 10). These results provide correlative evidence that climatic variation has shaped the genetic structure of steelhead populations and highlight the need for replication and sensitivity analyses in land and riverscape genetics.

Remote Sensing Derived Fire Frequency, Soil Moisture and Ecosystem Productivity Explain Regional Movements in Emu over Australia

Species distribution modeling has been widely used in studying habitat relationships and for conservation purposes. However, neglecting ecological knowledge about species, e.g. their seasonal movements, and ignoring the proper environmental factors that can explain key elements for species survival (shelter, food and water) increase model uncertainty. This study exemplifies how these ecological gaps in species distribution modeling can be addressed by modeling the distribution of the emu (Dromaius novaehollandiae) in Australia. Emus cover a large area during the austral winter. However, their habitat shrinks during the summer months. We show evidence of emu summer habitat shrinkage due to higher fire frequency, and low water and food availability in northern regions. Our findings indicate that emus prefer areas with higher vegetation productivity and low fire recurrence, while their distribution is linked to an optimal intermediate (~0.12 m3 m(-3)) soil moisture range. We propose that the application of three geospatial data products derived from satellite remote sensing, namely fire frequency, ecosystem productivity, and soil water content, provides an effective representation of emu general habitat requirements, and substantially improves species distribution modeling and representation of the species' ecological habitat niche across Australia.

Vegetation Greening and Climate Change Promote Multidecadal Rises of Global Land Evapotranspiration

Recent studies showed that anomalous dry conditions and limited moisture supply roughly between 1998 and 2008, especially in the Southern Hemisphere, led to reduced vegetation productivity and ceased growth in land evapotranspiration (ET). However, natural variability of Earth’s climate system can degrade capabilities for identifying climate trends. Here we produced a long-term (1982–2013) remote sensing based land ET record and investigated multidecadal changes in global ET and underlying causes. The ET record shows a significant upward global trend of 0.88 mm yr⁻² (P < 0.001) over the 32-year period, mainly driven by vegetation greening (0.018% per year; P < 0.001) and rising atmosphere moisture demand (0.75 mm yr⁻²; P = 0.016). Our results indicate that reduced ET growth between 1998 and 2008 was an episodic phenomenon, with subsequent recovery of the ET growth rate after 2008. Terrestrial precipitation also shows a positive trend of 0.66 mm yr⁻² (P = 0.08) over the same period consistent with expected water cycle intensification, but this trend is lower than coincident increases in evaporative demand and ET, implying a possibility of cumulative water supply constraint to ET. Continuation of these trends will likely exacerbate regional drought-induced disturbances, especially during regional dry climate phases associated with strong El Niño events.

Comparing land surface phenology derived from satellite and GPS network microwave remote sensing

The land surface phenology (LSP) start of season (SOS) metric signals the seasonal onset of vegetation activity, including canopy growth and associated increases in land-atmosphere water, energy and carbon (CO2) exchanges influencing weather and climate variability. The vegetation optical depth (VOD) parameter determined from satellite passive microwave remote sensing provides for global LSP monitoring that is sensitive to changes in vegetation canopy water content and biomass, and insensitive to atmosphere and solar illumination constraints. Direct field measures of canopy water content and biomass changes desired for LSP validation are generally lacking due to the prohibitive costs of maintaining regional monitoring networks. Alternatively, a normalized microwave reflectance index (NMRI) derived from GPS base station measurements is sensitive to daily vegetation water content changes and may provide for effective microwave LSP validation. We compared multiyear (2007-2011) NMRI and satellite VOD records at over 300 GPS sites in North America, and their derived SOS metrics for a subset of 24 homogenous land cover sites to investigate VOD and NMRI correspondence, and potential NMRI utility for LSP validation. Significant correlations (P<0.05) were found at 276 of 305 sites (90.5 %), with generally favorable correspondence in the resulting SOS metrics (r (2)=0.73, P<0.001, RMSE=36.8 days). This study is the first attempt to compare satellite microwave LSP metrics to a GPS network derived reflectance index and highlights both the utility and limitations of the NMRI data for LSP validation, including spatial scale discrepancies between local NMRI measurements and relatively coarse satellite VOD retrievals.

Satellite microwave detection of boreal forest recovery from the extreme 2004 wildfires in Alaska and Canada

The rate of vegetation recovery from boreal wildfire influences terrestrial carbon cycle processes and climate feedbacks by affecting the surface energy budget and land-atmosphere carbon exchange. Previous forest recovery assessments using satellite optical-infrared normalized difference vegetation index (NDVI) and tower CO(2) eddy covariance techniques indicate rapid vegetation recovery within 5-10 years, but these techniques are not directly sensitive to changes in vegetation biomass. Alternatively, the vegetation optical depth (VOD) parameter from satellite passive microwave remote sensing can detect changes in canopy biomass structure and may provide a useful metric of post-fire vegetation response to inform regional recovery assessments. We analyzed a multi-year (2003-2010) satellite VOD record from the NASA AMSR-E (Advanced Microwave Scanning Radiometer for EOS) sensor to estimate forest recovery trajectories for 14 large boreal fires from 2004 in Alaska and Canada. The VOD record indicated initial post-fire canopy biomass recovery within 3-7 years, lagging NDVI recovery by 1-5 years. The VOD lag was attributed to slower non-photosynthetic (woody) and photosynthetic (foliar) canopy biomass recovery, relative to the faster canopy greenness response indicated from the NDVI. The duration of VOD recovery to pre-burn conditions was also directly proportional (P < 0.01) to satellite (moderate resolution imaging spectroradiometer) estimated tree cover loss used as a metric of fire severity. Our results indicate that vegetation biomass recovery from boreal fire disturbance is generally slower than reported from previous assessments based solely on satellite optical-infrared remote sensing, while the VOD parameter enables more comprehensive assessments of boreal forest recovery.

Satellite Microwave Remote Sensing for Environmental Modeling of Mosquito Population Dynamics

Environmental variability has important influences on mosquito life cycles and understanding the spatial and temporal patterns of mosquito populations is critical for mosquito control and vector-borne disease prevention. Meteorological data used for model-based predictions of mosquito abundance and life cycle dynamics are typically acquired from ground-based weather stations; however, data availability and completeness are often limited by sparse networks and resource availability. In contrast, environmental measurements from satellite remote sensing are more spatially continuous and can be retrieved automatically. This study compared environmental measurements from the NASA Advanced Microwave Scanning Radiometer on EOS (AMSR-E) and in situ weather station data to examine their ability to predict the abundance of two important mosquito species (Aedes vexans and Culex tarsalis) in Sioux Falls, South Dakota, USA from 2005 to 2010. The AMSR-E land parameters included daily surface water inundation fraction, surface air temperature, soil moisture, and microwave vegetation opacity. The AMSR-E derived models had better fits and higher forecasting accuracy than models based on weather station data despite the relatively coarse (25-km) spatial resolution of the satellite data. In the AMSR-E models, air temperature and surface water fraction were the best predictors of Aedes vexans, whereas air temperature and vegetation opacity were the best predictors of Cx. tarsalis abundance. The models were used to extrapolate spatial, seasonal, and interannual patterns of climatic suitability for mosquitoes across eastern South Dakota. Our findings demonstrate that environmental metrics derived from satellite passive microwave radiometry are suitable for predicting mosquito population dynamics and can potentially improve the effectiveness of mosquito-borne disease early warning systems.

CLIMATE, HYDROLOGIC DISTURBANCE, AND SUCCESSION: DRIVERS OF FLOODPLAIN PATTERN

Floodplains are among the world's most threatened ecosystems due to the pervasiveness of dams, levee systems, and other modifications to rivers. Few unaltered floodplains remain where we may examine their dynamics over decadal time scales. Our study provides a detailed examination of landscape change over a 60-year period (1945-2004) on the Nyack floodplain of the Middle Fork of the Flathead River, a free-flowing, gravel-bed river in northwest Montana, USA. We used historical aerial photographs and airborne and satellite imagery to delineate habitats (i.e., mature forest, regenerative forest, water, cobble) within the floodplain. We related changes in the distribution and size of these habitats to hydrologic disturbance and regional climate. Results show a relationship between changes in floodplain habitats and annual flood magnitude, as well as between hydrology and the cooling and warming phases of the Pacific Decadal Oscillation (PDO). Large magnitude floods and greater frequency of moderate floods were associated with the cooling phases of the PDO, resulting in a floodplain environment dominated by extensive restructuring and regeneration of floodplain habitats. Conversely, warming phases of the PDO corresponded with decreases in magnitude, duration, and frequency of critical flows, creating a floodplain environment dominated by late successional vegetation and low levels of physical restructuring. Over the 60-year time series, habitat change was widespread throughout the floodplain, though the relative abundances of the habitats did not change greatly. We conclude that the long- and short-term interactions of climate, floods, and plant succession produce a shifting habitat mosaic that is a fundamental attribute of natural floodplain ecosystems.

Simulating effects of fire disturbance and climate change on boreal forest productivity and evapotranspiration

We used a terrestrial ecosystem process model, BIOME-BGC, to investigate historical climate change and fire disturbance effects on regional carbon and water budgets within a 357,500 km(2) portion of the Canadian boreal forest. Historical patterns of increasing atmospheric CO2, climate change, and regional fire activity were used as model drivers to evaluate the relative effects of these impacts to spatial patterns and temporal trends in forest net primary production (NPP) and evapotranspiration (ET). Historical trends of increasing atmospheric CO2 resulted in overall 13% and 5% increases in annual NPP and ET from 1994 to 1996, respectively. NPP was found to be relatively sensitive to changes in air temperature (T(a)), while ET was more sensitive to precipitation (P) change within the ranges of observed climate variability (e.g., +/-2 degrees C for T(a) and +/-20% for P). In addition, the potential effect of climate change related warming on NPP is exacerbated or offset depending on whether these changes are accompanied by respective decreases or increases in precipitation. Historical fire activity generally resulted in reductions of both NPP and ET, which consumed an average of approximately 6% of annual NPP from 1959 to 1996. Areas currently occupied by dry conifer forests were found to be subject to more frequent fire activity, which consumed approximately 8% of annual NPP. The results of this study show that the North American boreal ecosystem is sensitive to historical patterns of increasing atmospheric CO2, climate change and regional fire activity. The relative impacts of these disturbances on NPP and ET interact in complex ways and are spatially variable depending on regional land cover and climate gradients.

Simulating the effects of climate change on the carbon balance of North American high-latitude forests

The large magnitude of predicted warming at high latitudes and the potential feedback of ecosystems to atmospheric CO₂ concentrations make it important to quantify both warming and its effects on high-latitude carbon balance. We analysed long-term, daily surface meteorological records for 13 sites in Alaska and north-western Canada and an 82-y record of river ice breakup date for the Tanana River in interior Alaska. We found increases in winter and spring temperature extrema for all sites, with the greatest increases in spring minimum temperature, average 0.47 °C per 10 y, and a 0.7-day per 10 y advance in ice breakup on the Tanana River. We used the climate records to drive an ecosystem process model, BIOME_BGC, to simulate the effects of climate change on the carbon and water balances of boreal forest ecosystems. The growing season has lengthened by an average of 2.6 days per 10 y with an advance in average leaf onset date of 1.10 days per 10 y. This advance in the start of the active growing season correlates positively with progressively earlier ice breakup on the Tanana River in interior Alaska. The advance in the start of the growing season resulted in a 20% increase in net primary production for both aspen (Populus tremuloides) and white spruce (Picea glauca) stands. Aspen had a greater mean increase in maintenance respiration than spruce, whereas spruce had a greater mean increase in evapotranspiration. Average decomposition rates also increased for both species. Both net primary production and decomposition are enhanced in our simulations, suggesting that productive forest types may not experience a significant shift in net carbon flux as a result of climate warming.

Simulating forest productivity and surface-atmosphere carbon exchange in the BOREAS study region

A process-based, general ecosystem model (BIOME-BGC) was used to simulate daily gross primary production, maintenance and heterotrophic respiration, net primary production and net ecosystem carbon exchange of boreal aspen, jack pine and black spruce stands. Model simulations of daily net carbon exchange of the ecosystem (NEE) explained 51.7% (SE = 1.32 g C m(-2) day(-1)) of the variance in daily NEE derived from stand eddy flux measurements of CO(2) during 1994. Differences between measured and simulated results were attributed to several factors including difficulties associated with measuring nighttime CO(2) fluxes and model assumptions of site homogeneity. However, comparisons between simulations and field data improved markedly at coarser time-scales. Model simulations explained 66.1% (SE = 0.97 g C m(-2) day(-1)) of the variance in measured NEE when 5-day means of daily results were compared. Annual simulations of aboveground net primary production ranged from 0.6-2.4 Mg C ha(-1) year(-1) and were concurrent with results derived from tree increment core measurements and allometric equations. Model simulations showed that all of the sites were net sinks (0.1-4.1 Mg C ha(-1) year(-1)) of atmospheric carbon for 1994. Older conifer stands showed narrow margins between uptake of carbon by net photosynthesis and carbon release through respiration. Younger stands were more productive than older stands, primarily because of lower maintenance respiration costs. However, all sites appeared to be less productive than temperate forests. Productivity simulations were strongly linked to stand morphology and site conditions. Old jack pine and aspen stands showed decreased productivity in response to simulated low soil water contents near the end of the 1994 growing season. Compared with the aspen stand, the jack pine stand appeared better adapted to conserve soil water through lower daily evapotranspiration losses but also exhibited a narrower margin between daily net photosynthesis and respiration. Stands subjected to water stress during the growing season may exist on the edge between being annual sources or sinks for atmospheric carbon.