Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatial patterns”

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

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

At least 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Pre-industrial to End 21st Century Projections of Tropospheric Ozone from the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP)

Present day tropospheric ozone and its changes between 1850 and 2100 are considered, analysing 15 global models that participated in the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP). The ensemble mean compares well against present day observations. The seasonal cycle correlates well, except for some locations in the tropical upper troposphere. Most (75 %) of the models are encompassed with a range of global mean tropospheric ozone column estimates from satellite data, but there is a suggestion of a high bias in the Northern Hemisphere and a low bias in the Southern Hemisphere, which could indicate deficiencies with the ozone precursor emissions. Compared to the present day ensemble mean tropospheric ozone burden of 337+/-23 Tg, the ensemble mean burden for 1850 time slice is approx. 30% lower. Future changes were modelled using emissions and climate projections from four Representative Concentration Pathways (RCPs). Compared to 2000, the relative changes in the ensemble mean tropospheric ozone burden in 2030 (2100) for the different RCPs are: −4% (−16 %) for RCP2.6, 2% (−7%) for RCP4.5, 1% (−9%) for RCP6.0, and 7% (18 %) for RCP8.5. Model agreement on the magnitude of the change is greatest for larger changes. Reductions in most precursor emissions are common across the RCPs and drive ozone decreases in all but RCP8.5, where doubled methane and a 40-150% greater stratospheric influx (estimated from a subset of models) increase ozone. While models with a high ozone burden for the present day also have high ozone burdens for the other time slices, no model consistently predicts large or small ozone changes; i.e. the magnitudes of the burdens and burden changes do not appear to be related simply, and the models are sensitive to emissions and climate changes in different ways. Spatial patterns of ozone changes are well correlated across most models, but are notably different for models without time evolving stratospheric ozone concentrations. A unified approach to ozone budget specifications and a rigorous investigation of the factors that drive tropospheric ozone is recommended to help future studies attribute ozone changes and inter-model differences more clearly.

Troposphere↗
NASA NTRS2013

Tropical Atlantic Dust and Smoke Aerosol Variations Related to the Madden-Julian Oscillation in MODIS and MISR Observations

In this study, Moderate Resolution Imaging Spectroradiometer (MODIS) fine mode fraction and Multi-angle Imaging SpectroRadiometer (MISR) nonspherical fraction data are used to derive dust and smoke aerosol optical thickness (T(sub dust) and T(sub smoke)) over the tropical Atlantic in a complementary way: due to its wider swath, MODIS has 3-4 times greater sampling than MISR, but MISR dust discrimination is based on particle shape retrievals, whereas an empirical scheme is used for MODIS. MODIS and MISR show very similar dust and smoke winter climatologies. T(sub dust) is the dominant aerosol component over the tropical Atlantic, accounting for 40-70 percent of the total aerosol optical thickness (AOT), whereas T(sub smoke) is significantly smaller than T(sub dust). The consistency and high correlation between these climatologies and their daily variations lends confidence to their use for investigating the relative dust and smoke contributions to the total AOT variation associated with the Madden-Julian Oscillation (MJO). The temporal evolution and spatial patterns of the τdus anomalies associated with the MJO are consistent between MODIS and MISR: the magnitude of MJO-realted T(sub dust) anomalies is comparable to or even larger than that of the total T, while the T(sub smoke) anomaly represents about 15 percent compared to the total, which is quite different from their relative magnitudes to the total T on the climatological time scale. This suggests that dust and smoke are not influenced by the MJO in the same way. Based on correlation analysis, dust is strongly influenced by the MJO-modulated trade wind and precipitation anomalies, and can last as long as one MJO phase, whereas smoke is less affected.

Madden-Julian Oscillation↗
NASA NTRS2013

Estimating Root Mean Square Errors in Remotely Sensed Soil Moisture over Continental Scale Domains

Root Mean Square Errors (RMSE) in the soil moisture anomaly time series obtained from the Advanced Scatterometer (ASCAT) and the Advanced Microwave Scanning Radiometer (AMSR-E; using the Land Parameter Retrieval Model) are estimated over a continental scale domain centered on North America, using two methods: triple colocation (RMSETC ) and error propagation through the soil moisture retrieval models (RMSEEP ). In the absence of an established consensus for the climatology of soil moisture over large domains, presenting a RMSE in soil moisture units requires that it be specified relative to a selected reference data set. To avoid the complications that arise from the use of a reference, the RMSE is presented as a fraction of the time series standard deviation (fRMSE). For both sensors, the fRMSETC and fRMSEEP show similar spatial patterns of relatively highlow errors, and the mean fRMSE for each land cover class is consistent with expectations. Triple colocation is also shown to be surprisingly robust to representativity differences between the soil moisture data sets used, and it is believed to accurately estimate the fRMSE in the remotely sensed soil moisture anomaly time series. Comparing the ASCAT and AMSR-E fRMSETC shows that both data sets have very similar accuracy across a range of land cover classes, although the AMSR-E accuracy is more directly related to vegetation cover. In general, both data sets have good skill up to moderate vegetation conditions.

Remotely sensed soil moisture validation↗
NASA NTRS2014

Revisiting AVHRR Tropospheric Aerosol Trends Using Principal Component Analysis

The advanced very high resolution radiometer (AVHRR) satellite instruments provide a nearly 25 year continuous record of global aerosol properties over the ocean. It offers valuable insights into the long-term change in global aerosol loading. However, the AVHRR data record is heavily influenced by two volcanic eruptions, El Chichon on March 1982 and Mount Pinatubo on June 1991. The gradual decay of volcanic aerosols may last years after the eruption, which potentially masks the estimation of aerosol trends in the lower troposphere, especially those of anthropogenic origin. In this study, we show that a principal component analysis approach effectively captures the bulk of the spatial and temporal variability of volcanic aerosols into a single mode. The spatial pattern and time series of this mode provide a good match to the global distribution and decay of volcanic aerosols. We further reconstruct the data set by removing the volcanic aerosol component and reestimate the global and regional aerosol trends. Globally, the reconstructed data set reveals an increase of aerosol optical depth from 1985 to 1990 and decreasing trend from 1994 to 2006. Regionally, in the 1980s, positive trends are observed over the North Atlantic and North Arabian Sea, while negative tendencies are present off the West African coast and North Pacific. During the 1994 to 2006 period, the Gulf of Mexico, North Atlantic close to Europe, and North Africa exhibit negative trends, while the coastal regions of East and South Asia, the Sahel region, and South America show positive trends.

volcanoes↗
NASA NTRS2013

Siberian Pine Decline and Mortality in Southern Siberian Mountains

The causes and resulting spatial patterns of Siberian pine mortality in eastern Kuznetzky Alatau Mountains, Siberia were analyzed based on satellite (Landsat, MODIS) and dendrochronology data. Climate variables studied included temperature, precipitation and Standardized Precipitation-Evapotranspiration Index (SPEI) drought index. Landsat data analysis showed that stand mortality was first detected in the year 2006 at an elevation of 650 m, and extended up to 900 m by the year 2012. Mortality was accompanied by a decrease in MODIS derived vegetation index (EVI).. The area of dead stands and the upper mortality line were correlated with increased drought. The uphill margin of mortality was limited by elevational precipitation gradients. Dead stands (i.e., >75% tree mortality) were located mainly on southern slopes. With respect to slope, mortality was observed within a 7 deg - 20 deg range with greatest mortality occurring on convex terrain. Tree radial incrementmeasurements correlate and were synchronous with SPEI (r sq = 0.37, r(sub s) = 80). Increasing synchrony between tree ring growth and SPEI indicates that drought has reduced the ecological niche of Siberian pine. The results also showed the primary role of drought stress on Siberian pine mortality. A secondary role may be played by bark beetles and root fungi attacks. The observed Siberian pine mortality is part of a broader phenomenon of "dark needle conifers" (DNC, i.e., Siberian pine, fir and spruce) decline and mortality in European Russia, Siberia, and the Russian Far East. All locations of DNC decline coincided with areas of observed drought increase. The results obtained are one of the first observations of drought-induced decline and mortality of DNC at the southern border of boreal forests. Meanwhile if model projections of increased aridity are correct DNC, within the southern part of its range may be replaced by drought-resistant Pinus silvestris and Larix sibirica.

tree die-off↗
NASA NTRS2014

Global Free-tropospheric NO2 Abundances Derived Using a Cloud Slicing Technique from AURA OMI

We derive free-tropospheric NO2 volume mixing ratios (VMRs) by applying a cloud-slicing technique to data from the Ozone Monitoring Instrument (OMI) on the Aura satellite. In the cloud-slicing approach, the slope of the above-cloud NO2 column versus the cloud scene pressure is proportional to the NO2 VMR. In this work, we use a sample of nearby OMI pixel data from a single orbit for the linear fit. The OMI data include cloud scene pressures from the rotational-Raman algorithm and above-cloud NO2 vertical column density (VCD) (defined as the NO2 column from the cloud scene pressure to the top of the atmosphere) from a differential optical absorption spectroscopy (DOAS) algorithm. We compare OMI-derived NO2 VMRs with in situ aircraft profiles measured during the NASA Intercontinental Chemical Transport Experiment Phase B (INTEX-B) campaign in 2006. The agreement is generally within the estimated uncertainties when appropriate data screening is applied. We then derive a global seasonal climatology of free-tropospheric NO2 VMR in cloudy conditions. Enhanced NO2 in the free troposphere commonly appears near polluted urban locations where NO2 produced in the boundary layer may be transported vertically out of the boundary layer and then horizontally away from the source. Signatures of lightning NO2 are also shown throughout low and middle latitude regions in summer months. A profile analysis of our cloud-slicing data indicates signatures of lightning-generated NO2 in the upper troposphere. Comparison of the climatology with simulations from the global modeling initiative (GMI) for cloudy conditions (cloud optical depth less than10) shows similarities in the spatial patterns of continental pollution outflow. However, there are also some differences in the seasonal variation of free-tropospheric NO2 VMRs near highly populated regions and in areas affected by lightning-generated NOx.

NO2↗
NASA NTRS2015

Global Soil Moisture from the Aquarius/SAC-D Satellite: Description and Initial Assessment

Aquarius satellite observations over land offer a new resource for measuring soil moisture from space. Although Aquarius was designed for ocean salinity mapping, our objective in this investigation is to exploit the large amount of land observations that Aquarius acquires and extend the mission scope to include the retrieval of surface soil moisture. The soil moisture retrieval algorithm development focused on using only the radiometer data because of the extensive heritage of passive microwave retrieval of soil moisture. The single channel algorithm (SCA) was implemented using the Aquarius observations to estimate surface soil moisture. Aquarius radiometer observations from three beams (after bias/gain modification) along with the National Centers for Environmental Prediction model forecast surface temperatures were then used to retrieve soil moisture. Ancillary data inputs required for using the SCA are vegetation water content, land surface temperature, and several soil and vegetation parameters based on land cover classes. The resulting global spatial patterns of soil moisture were consistent with the precipitation climatology and with soil moisture from other satellite missions (Advanced Microwave Scanning Radiometer for the Earth Observing System and Soil Moisture Ocean Salinity). Initial assessments were performed using in situ observations from the U.S. Department of Agriculture Little Washita and Little River watershed soil moisture networks. Results showed good performance by the algorithm for these land surface conditions for the period of August 2011-June 2013 (rmse = 0.031 m(exp 3)/m(exp 3), Bias = −0.007 m(exp 3)/m(exp 3), and R = 0.855). This radiometer-only soil moisture product will serve as a baseline for continuing research on both active and combined passive-active soil moisture algorithms. The products are routinely available through the National Aeronautics and Space Administration data archive at the National Snow and Ice Data Center.

Soil MOisture↗
NASA NTRS2014

Diagnosis of Middle Atmosphere Climate Sensitivity by the Climate Feedback Response Analysis Method

We present a new method to diagnose the middle atmosphere climate sensitivity by extending the Climate Feedback-Response Analysis Method (CFRAM) for the coupled atmosphere-surface system to the middle atmosphere. The Middle atmosphere CFRAM (MCFRAM) is built on the atmospheric energy equation per unit mass with radiative heating and cooling rates as its major thermal energy sources. MCFRAM preserves the CFRAM unique feature of an additive property for which the sum of all partial temperature changes due to variations in external forcing and feedback processes equals the observed temperature change. In addition, MCFRAM establishes a physical relationship of radiative damping between the energy perturbations associated with various feedback processes and temperature perturbations associated with thermal responses. MCFRAM is applied to both measurements and model output fields to diagnose the middle atmosphere climate sensitivity. It is found that the largest component of the middle atmosphere temperature response to the 11-year solar cycle (solar maximum vs. solar minimum) is directly from the partial temperature change due to the variation of the input solar flux. Increasing CO2 always cools the middle atmosphere with time whereas partial temperature change due to O3 variation could be either positive or negative. The partial temperature changes due to different feedbacks show distinctly different spatial patterns. The thermally driven globally averaged partial temperature change due to all radiative processes is approximately equal to the observed temperature change, ranging from 0.5 K near 70 km from the near solar maximum to the solar minimum.

MCFRAM↗
NASA NTRS2013

Natural Air-Sea Flux of CO2 in Simulations of the NASA-GISS Climate Model: Sensitivity to the Physical Ocean Model Formulation

Results from twin control simulations of the preindustrial CO2 gas exchange (natural flux of CO2) between the ocean and the atmosphere are presented here using the NASA-GISS climate model, in which the same atmospheric component (modelE2) is coupled to two different ocean models, the Russell ocean model and HYCOM. Both incarnations of the GISS climate model are also coupled to the same ocean biogeochemistry module (NOBM) which estimates prognostic distributions for biotic and abiotic fields that influence the air-sea flux of CO2. Model intercomparison is carried out at equilibrium conditions and model differences are contrasted with biases from present day climatologies. Although the models agree on the spatial patterns of the air-sea flux of CO2, they disagree on the strength of the North Atlantic and Southern Ocean sinks mainly because of kinematic (winds) and chemistry (pCO2) differences rather than thermodynamic (SST) ones. Biology/chemistry dissimilarities in the models stem from the different parameterizations of advective and diffusive processes, such as overturning, mixing and horizontal tracer advection and to a lesser degree from parameterizations of biogeochemical processes such as gravitational settling and sinking. The global meridional overturning circulation illustrates much of the different behavior of the biological pump in the two models, together with differences in mixed layer depth which are responsible for different SST, DIC and nutrient distributions in the two models and consequently different atmospheric feedbacks (in the wind, net heat and freshwater fluxes into the ocean).

Gas exchange↗
NASA NTRS2013

Decadal Prediction Skill in the GEOS-5 Forecast System

A suite of decadal predictions has been conducted with the NASA Global Modeling and Assimilation Office's (GMAO's) GEOS-5 Atmosphere-Ocean general circulation model. The hind casts are initialized every December 1st from 1959 to 2010, following the CMIP5 experimental protocol for decadal predictions. The initial conditions are from a multivariate ensemble optimal interpolation ocean and sea-ice reanalysis, and from GMAO's atmospheric reanalysis, the modern-era retrospective analysis for research and applications. The mean forecast skill of a three-member-ensemble is compared to that of an experiment without initialization but also forced with observed greenhouse gases. The results show that initialization increases the forecast skill of North Atlantic sea surface temperature compared to the uninitialized runs, with the increase in skill maintained for almost a decade over the subtropical and mid-latitude Atlantic. On the other hand, the initialization reduces the skill in predicting the warming trend over some regions outside the Atlantic. The annual-mean Atlantic meridional overturning circulation index, which is defined here as the maximum of the zonally-integrated overturning stream function at mid-latitude, is predictable up to a 4-year lead time, consistent with the predictable signal in upper ocean heat content over the North Atlantic. While the 6- to 9-year forecast skill measured by mean squared skill score shows 50 percent improvement in the upper ocean heat content over the subtropical and mid-latitude Atlantic, prediction skill is relatively low in the sub-polar gyre. This low skill is due in part to features in the spatial pattern of the dominant simulated decadal mode in upper ocean heat content over this region that differ from observations. An analysis of the large-scale temperature budget shows that this is the result of a model bias, implying that realistic simulation of the climatological fields is crucial for skillful decadal forecasts.

Decadal Prediction↗
NASA NTRS2015

An Assessment of Southern Ocean Water Masses and Sea Ice During 1988-2007 in a Suite of Interannual CORE-II Simulations

We characterise the representation of the Southern Ocean water mass structure and sea ice within a suite of 15 global ocean-ice models run with the Coordinated Ocean-ice Reference Experiment Phase II (CORE-II) protocol. The main focus is the representation of the present (1988-2007) mode and intermediate waters, thus framing an analysis of winter and summer mixed layer depths; temperature, salinity, and potential vorticity structure; and temporal variability of sea ice distributions. We also consider the interannual variability over the same 20 year period. Comparisons are made between models as well as to observation-based analyses where available. The CORE-II models exhibit several biases relative to Southern Ocean observations, including an underestimation of the model mean mixed layer depths of mode and intermediate water masses in March (associated with greater ocean surface heat gain), and an overestimation in September (associated with greater high latitude ocean heat loss and a more northward winter sea-ice extent). In addition, the models have cold and fresh/warm and salty water column biases centred near 50 deg S. Over the 1988-2007 period, the CORE-II models consistently simulate spatially variable trends in sea-ice concentration, surface freshwater fluxes, mixed layer depths, and 200-700 m ocean heat content. In particular, sea-ice coverage around most of the Antarctic continental shelf is reduced, leading to a cooling and freshening of the near surface waters. The shoaling of the mixed layer is associated with increased surface buoyancy gain, except in the Pacific where sea ice is also influential. The models are in disagreement, despite the common CORE-II atmospheric state, in their spatial pattern of the 20-year trends in the mixed layer depth and sea-ice.

Sea ice↗
NASA NTRS2015

Late-Glacial to Late-holocene Shifts in Global Precipitation Delta(sup 18)O

Reconstructions of Quaternary climate are often based on the isotopic content of paleo-precipitation preserved in proxy records. While many paleo-precipitation isotope records are available, few studies have synthesized these dispersed records to explore spatial patterns of late-glacial precipitation delta(sup 18)O. Here we present a synthesis of 86 globally distributed groundwater (n 59), cave calcite (n 15) and ice core (n 12) isotope records spanning the late-glacial (defined as 50,000 to 20,000 years ago) to the late-Holocene (within the past 5000 years). We show that precipitation delta(sup 18)O changes from the late-glacial to the late-Holocene range from -7.1% (delta(sup 18)O(late-Holocene) > delta(sup 18)O(late-glacial) to +1.7% (delta(sup 18)O(late-glacial) > delta(sup 18)O(late-Holocene), with the majority (77) of records having lower late-glacial delta(sup 18)O than late-Holocene delta(sup 18)O values. High-magnitude, negative precipitation delta(sup 18)O shifts are common at high latitudes, high altitudes and continental interiors.

isotopes↗
NASA NTRS2015

Extension and Statistical Analysis of the GACP Aerosol Optical Thickness Record.

The primary product of the Global Aerosol Climatology Project (GACP) is a continuous record of the aerosol optical thickness (AOT) over the oceans. It is based on channel-1 and -2 radiance data from the Advanced Very High Resolution Radiometer (AVHRR) instruments flown on successive National Oceanic and Atmospheric Administration (NOAA) platforms. We extend the previous GACP dataset by four years through the end of 2009 using NOAA-17 and -18 AVHRR radiances recalibrated against MODerate resolution Imaging Spectroradiometer (MODIS) radiance data, thereby making the GACP record almost three decades long. The temporal overlap of over three years of the new NOAA-17 and the previous NOAA-16 record reveals an excellent agreement of the corresponding global monthly mean AOT values, thereby confirming the robustness of the vicarious radiance calibration used in the original GACP product. The temporal overlap of the NOAA-17 and -18 instruments is used to introduce a small additive adjustment to the channel-2 calibration of the latter resulting in a consistent record with increased data density. The Principal Component Analysis (PCA) of the newly extended GACP record shows that most of the volcanic AOT variability can be isolated into one mode responsible for ~12% of the total variance. This conclusion is confirmed by a combined PCA analysis of the GACP, MODIS, andMulti-angle Imaging SpectroRadiometer (MISR) AOTs during the volcano-free period fromFebruary 2000 to December 2009.We show that the modes responsible for the tropospheric AOT variability in the three datasets agree well in terms of correlation and spatial patterns. A previously identified negative AOT trend which started in the late 1980s and continued into the early 2000s is confirmed. Its magnitude and duration indicate that it was caused by changes in tropospheric aerosols. The latest multi-satellite segment of the GACP record shows that this trend tapered off, with no noticeable AOT change after 2002. This result is consistent with the MODIS andMISR AOT records as well as with the recent gradual reversal frombrightening to dimming revealed by surface flux measurements in many aerosol producing regions. Thus the robustness of the GACP record is confirmed, increasing our confidence in the validity of the negative trend. Although the nominal negative GACP AOT trend could partially be an artifact of increasing aerosol absorption, we argue that the time dependence of the GACP record, including the latest flat period, is more consistent with the actual decrease in the tropospheric AOT.

aerosols↗
NASA NTRS2015

Comparison of GFED3, QFED2 and FEER1 Biomass Burning Emissions Datasets in a Global Model

Biomass burning contributes about 40% of the global loading of carbonaceous aerosols, significantly affecting air quality and the climate system by modulating solar radiation and cloud properties. However, fire emissions are poorly constrained in models on global and regional levels. In this study, we investigate 3 global biomass burning emission datasets in NASA GEOS5, namely: (1) GFEDv3.1 (Global Fire Emissions Database version 3.1); (2) QFEDv2.4 (Quick Fire Emissions Dataset version 2.4); (3) FEERv1 (Fire Energetics and Emissions Research version 1.0). The simulated aerosol optical depth (AOD), absorption AOD (AAOD), angstrom exponent and surface concentrations of aerosol plumes dominated by fire emissions are evaluated and compared to MODIS, OMI, AERONET, and IMPROVE data over different regions. In general, the spatial patterns of biomass burning emissions from these inventories are similar, although the strength of the emissions can be noticeably different. The emissions estimates from QFED are generally larger than those of FEER, which are in turn larger than those of GFED. AOD simulated with all these 3 databases are lower than the corresponding observations in Southern Africa and South America, two of the major biomass burning regions in the world.

biomass burning↗
NASA NTRS2016

Comparison of AIRS Version-6 OLR Climatologies and Anomaly Time Series with Those of CERES and MERRA-2

RCs of AIRS and MERRA-2 500 mb specific humidity agree very well in terms of spatial patterns, but MERRA-2 ARCs are larger in magnitude and show a spurious moistening globally and over Central Africa. AIRS and MERRA-2 fractional cloud cover ARCs agree less well with each other. MERRA-2 shows a spurious global mean increase in cloud cover that is not found in AIRS, including a large spurious cloud increase in Central Africa. AIRS and MERRA-2 ARCs of surface skin and surface air temperatures are all similar to each other in patterns. AIRS shows a small global warming over the 13 year period, while MERRA-2 shows a small global cooling. This difference results primarily from spurious MERRA-2 temperature trends at high latitudes and over Central Africa. These differences all contribute to the spurious negative global MERRA-2 OLR trend. AIRS Version-6 confirms that 2015 is the warmest year on record and that the Earth's surface is continuing to warm.

AIRS↗
NASA NTRS2016

Have GRACE Satellites Overestimated Groundwater Depletion in the Northwest India Aquifer?

The Northwest India Aquifer (NWIA) has been shown to have the highest groundwater depletion (GWD) rate globally, threatening crop production and sustainability of groundwater resources. Gravity Recovery and Climate Experiment (GRACE) satellites have been emerging as a powerful tool to evaluate GWD with ancillary data. Accurate GWD estimation is, however, challenging because of uncertainties in GRACE data processing. We evaluated GWD rates over the NWIA using a variety of approaches, including newly developed constrained forward modeling resulting in a GWD rate of 3.1 plus or minus 0.1 centimeters per acre (or 14 plus or minus 0.4 cubic kilometers per acre) for Jan 2005-Dec 2010, consistent with the GWD rate (2.8 centimeters per acre or 12.3 cubic kilometers per acre) from groundwater-level monitoring data. Published studies (e.g., 4 plus or minus 1 centimeter per acre or 18 plus or minus 4.4 cubic kilometers per acre) may overestimate GWD over this region. This study highlights uncertainties in GWD estimates and the importance of incorporating a priori information to refine spatial patterns of GRACE signals that could be more useful in groundwater resource management and need to be paid more attention in future studies.

Long, Di↗
NASA NTRS2015

Climate-Induced Mortality of Spruce Stands in Belarus

The aim of this work is an analysis of the causes of spruce (Picea abies L.) decline and mortality in Belarus. The analysis was based on forest inventory and Landsat satellite (land cover classification, climate variables (air temperature, precipitation, evaporation, vapor pressure deficit, SPEI drought index)), and GRACE-derived soil moisture estimation (equivalent of water thickness anomalies, EWTA). We found a difference in spatial patterns between dead stands and all stands (i.e., before mortality). Dead stands were located preferentially on relief features with higher water stress risk (i.e., higher elevations, steeper slopes, south and southwestern exposure). Spruce mortality followed a series of repeated droughts between 1990 and 2010. Mortality was negatively correlated with air humidity (r = -0.52), and precipitation (r = -0.57), and positively correlated with the prior year vapor pressure deficit (r = 0.47), and drought increase (r = 0.57). Mortality increased with the increase in occurrence of spring frosts (r = 0.5), and decreased with an increase in winter cloud cover (r = -0.37). Spruce mortality was negatively correlated with snow water accumulation (r = -0.81) and previous year anomalies in water soil content (r = -0.8). Weakened by water stress, spruce stands were attacked by pests and phytopathogens. Overall, spruce mortality in Belarussian forests was caused by drought episodes and drought increase in synergy with pest and phytopathogen attacks. Vast Picea abies mortality in Belarus and adjacent areas of Russia and Eastern Europe is a result of low adaptation of that species to increased drought. This indicates the necessity of spruce replacement by drought-tolerant indigenous (e.g., Pinus sylvestris, Querqus robur) or introduced (e.g., Larix sp. or Pseudotsuga menzieslii) species to obtain sustainable forest growth management.

spruce decline↗
NASA NTRS2016

Combining Livestock Production Information in a Process-Based Vegetation Model to Reconstruct the History of Grassland Management

Grassland management type (grazed or mown) and intensity (intensive or extensive) play a crucial role in the greenhouse gas balance and surface energy budget of this biome, both at field scale and at large spatial scale. However, global gridded historical information on grassland management intensity is not available. Combining modelled grass-biomass productivity with statistics of the grass-biomass demand by livestock, we reconstruct gridded maps of grassland management intensity from 1901 to 2012. These maps include the minimum area of managed vs. maximum area of unmanaged grasslands and the fraction of mown vs. grazed area at a resolution of 0.5deg by 0.5deg. The grass-biomass demand is derived from a livestock dataset for 2000, extended to cover the period 19012012. The grass-biomass supply (i.e. forage grass from mown grassland and biomass grazed) is simulated by the process-based model ORCHIDEE-GM driven by historical climate change, risingCO2 concentration, and changes in nitrogen fertilization. The global area of managed grassland obtained in this study increases from 6.1 x 10(exp 6) km(exp 2) in 1901 to 12.3 x 10(exp 6) kmI(exp 2) in 2000, although the expansion pathway varies between different regions. ORCHIDEE-GM also simulated augmentation in global mean productivity and herbage-use efficiency over managed grassland during the 20th century, indicating a general intensification of grassland management at global scale but with regional differences. The gridded grassland management intensity maps are model dependent because they depend on modelled productivity. Thus specific attention was given to the evaluation of modelled productivity against a series of observations from site-level net primary productivity (NPP) measurements to two global satellite products of gross primary productivity (GPP) (MODIS-GPP and SIF data). Generally, ORCHIDEE-GM captures the spatial pattern, seasonal cycle, and inter-annual variability of grassland productivity at global scale well and thus is appropriate for global applications presented here.

Livestock dataset↗