Search NASA⌕ Search

SEARCH · Search NASA

Results for “spatially explicit modeling”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

The Impact of ARM on Climate Modeling

Climate models are among humanity's most ambitious and elaborate creations. They are designed to simulate the interactions of the atmosphere, ocean, land surface, and cryosphere on time scales far beyond the limits of deterministic predictability, and including the effects of time-dependent external forcings. The processes involved include radiative transfer, fluid dynamics, microphysics, and some aspects of geochemistry, biology, and ecology. The models explicitly simulate processes on spatial scales ranging from the circumference of the Earth down to one hundred kilometers or smaller, and implicitly include the effects of processes on even smaller scales down to a micron or so. The atmospheric component of a climate model can be called an atmospheric global circulation model (AGCM). In an AGCM, calculations are done on a three-dimensional grid, which in some of today's climate models consists of several million grid cells. For each grid cell, about a dozen variables are time-stepped as the model integrates forward from its initial conditions. These so-called prognostic variables have special importance because they are the only things that a model remembers from one time step to the next; everything else is recreated on each time step by starting from the prognostic variables and the boundary conditions. The prognostic variables typically include information about the mass of dry air, the temperature, the wind components, water vapor, various condensed-water species, and at least a few chemical species such as ozone. A good way to understand how climate models work is to consider the lengthy and complex process used to develop one. Lets imagine that a new AGCM is to be created, starting from a blank piece of paper. The model may be intended for a particular class of applications, e.g., high-resolution simulations on time scales of a few decades. Before a single line of code is written, the conceptual foundation of the model must be designed through a creative envisioning that starts from the intended application and is based on current understanding of how the atmosphere works and the inventory of mathematical methods available.

atmospheric circulation↗
NASA NTRS2016

Regional Disparities in the Beneficial Effects of Rising CO2 Emissions on Crop Water Productivity

Rising atmospheric carbon dioxide concentrations are expected to enhance photosynthesis and reduce crop water use. However, there is high uncertainty about the global implications of these effects for future crop production and agricultural water requirements under climate change. Here we combine results from networks of field experiments and global crop models to present a spatially explicit global perspective on crop water productivity (CWP, the ratio of crop yield to evapotranspiration) for wheat, maize, rice and soybean under elevated carbon dioxide and associated climate change projected for a high-end greenhouse gas emissions scenario. We find carbon dioxide effects increase global CWP by 10[0;47]%-27[7;37]% (median[interquartile range] across the model ensemble) by the 2080s depending on crop types, with particularly large increases in arid regions (by up to 48[25;56]% for rain fed wheat). If realized in the fields, the effects of elevated carbon dioxide could considerably mitigate global yield losses whilst reducing agricultural consumptive water use (4-17%). We identify regional disparities driven by differences in growing conditions across agro-ecosystems that could have implications for increasing food production without compromising water security. Finally, our results demonstrate the need to expand field experiments and encourage greater consistency in modeling the effects of rising carbon dioxide across crop and hydrological modeling communities.

Deryng, Delphine↗
NASA NTRS2017

Daily Ambient Air Pollution Metrics for Five Cities: Evaluation of Data Fusion-Based Estimates and Uncertainties

Spatiotemporal characterization of ambient air pollutant concentrations is increasingly relying on the combination of observations and air quality models to provide well-constrained, spatially and temporally complete pollutant concentration fields. Air quality models, in particular, are attractive, as they characterize the emissions, meteorological, and physiochemical process linkages explicitly while providing continuous spatial structure. However, such modeling is computationally intensive and has biases. The limitations of spatially sparse and temporally incomplete observations can be overcome by blending the data with estimates from a physically and chemically coherent model, driven by emissions and meteorological inputs. We recently developed a data fusion method that blends ambient ground observations and chemical transport-modeled (CTM) data to estimate daily, spatially resolved pollutant concentrations and associatedcorrelations. In this study, we assess the ability of the data fusion method to produce daily metrics (i.e., 1-hr max, 8-hr max, and 24-hr average) of ambient air pollution that capture spatiotemporal air pollution trends for 12 pollutants (CO, NO2, NOx, O3, SO2, PM (sub10), PM (sub 2.5), and five PM (sub 2.5) components) across five metropolitan areas (Atlanta, Birmingham, Dallas, Pittsburgh, and St. Louis), from 2002 to 2008. Three sets of comparisons are performed: (1) the CTM concentrations are evaluated for each pollutant and metropolitan domain, (2) the data fusion concentrations are compared with the monitor data, (3) a comprehensive cross-validation analysis against observed data evaluates the quality of the data fusion model simulations across multiple metropolitan domains. The resulting daily spatial field estimates of air pollutant concentrations and uncertainties are not only consistent with observations, emissions, andmeteorology, but substantially improve CTM-derived results for nearly all pollutants and all cities, with the exception of NO2 for Birmingham. The greatest improvements occur for O3 and PM (sub 2.5). Squared spatiotemporal correlation coefficients range between simulations and observations determined using cross-validation across all cities for air pollutants of secondary and mixed origins are R-squared equal to 0.88-0.93 (O3), 0.81-0.89 (SO4), 0.67-0.83 (PM (sub 2.5)), 0.52-0.72 (NO3), 0.43-0.80 (NH4), 0.32-0.51 (OC), and 0.14-0.71 (PM (sub 10)). Results for relatively homogeneous pollutants of secondary origin, tend to be better than those for more spatially heterogeneous (larger spatial gradients) pollutants of primary origin (NOx, CO, SO2 and EC). Generally, background concentrations and spatial concentration gradients reflect interurban airshed complexity and the effects of regional transport, whereas daily spatial pattern variability shows intra urban consistency in the fused data. With sufficiently high CTM spatial resolution, traffic-related pollutants exhibit gradual concentration gradients that peak toward the urban centers. Ambient pollutant concentration uncertainty estimates for the fused data are both more accurate and smaller than those for either the observations or the model simulations alone.

spatiotemporal fusion↗
NASA NTRS1996

The barberplaid illusion: plaid motion is biased by elongated apertures

The perceived direction of motion of plaids windowed by elongated spatial Gaussians is biased toward the window's long axis. The bias increases as the relative angle between the plaid motion and the long axis of the window increases, peaks at a relative angle of approximately 45 deg, and then decreases. The bias increases as the window is made narrower (at fixed height) and decreases as the component spatial frequency increases (at fixed aperture size). We examine several models of human motion processing (cross-correlation, motion-energy, intersection-of-constraints, and vector-sum), and show that none of these standard models can predict our data. We conclude that spatial integration of motion signals plays a crucial role in plaid motion perception and that current models must be explicitly expanded to include such spatial interactions.

NASA Center ARC↗
NASA NTRS2007

Alaska at the Crossroads of Migration: Space Based Ornithology

Understanding bird migration on a global scale is one of the most compelling and challenging problems of modern biology with major implications for human health and conservation biology. Revolutionary advances in remote sensing now provide us with near real-time measurements of atmospheric and land surface conditions at high spatial resolution over entire continents. We use spatially-explicit, individual based bird migration models driven by numerical weather prediction models of atmospheric conditions, dynamic habitat suitability maps derived from remotely sensed land surface conditions, biophysiological models, and biological field data to simulate migration routes, timing, energy budgets, and survival of individual birds and populations. Long-distance migratory birds travel annually between breeding grounds in Alaska and wintering grounds in Latin Amierica. Approximately 25% of these species are potential vectors of Avian Influenza. Alaska is at the crossroads of Asian and New World migratory flyways and is likely to be a point of introduction of Asian H5N1 AI into the western hemisphere. If/when an infected bird is detected, a pressing question will be where was this bird several days ago, and where is it likely to go after it was released from the survey site? Answers to such questions will increase effectiveness of AI surveillance and mitigation measures. From a conservation perspective, Alaska's diverse landscape provides breeding sites for many migrants, and climatic and land surface changes along migratory flyways in the western hemisphere may reduce bird survival and physical condition upon arrival at Alaskan breeding territories, success and migrant populations.

Deppe, Jill↗
NASA NTRS2019

Parameterization-Induced Uncertainties and Impacts of Crop Management Harmonization in a Global Gridded Crop Model Ensemble

Global gridded crop models (GGCMs) combine agronomic or plant growth models with gridded spatial input data to estimate spatially explicit crop yields and agricultural externalities at the global scale. Differences in GGCM outputs arise from the use of different biophysical models, setups, and input data. GGCM ensembles are frequently employed to bracket uncertainties in impact studies without investigating the causes of divergence in outputs. This study explores differences in maize yield estimates from five GGCMs based on the public domain field-scale model Environmental Policy Integrated Climate (EPIC) that participate in the AgMIP Global Gridded Crop Model Intercomparison initiative. Albeit using the same crop model, the GGCMs differ in model version, input data, management assumptions, parameterization, and selection of subroutines affecting crop yield estimates via cultivar distributions, soil attributes, and hydrology among others. The analyses reveal inter-annual yield variability and absolute yield levels in the EPIC-based GGCMs to be highly sensitive to soil parameterization and crop management. All GGCMs show an intermediate performance in reproducing reported yields with a higher skill if a static soil profile is assumed or sufficient plant nutrients are supplied. An in-depth comparison of setup domains for two EPIC-based GGCMs shows that GGCM performance and plant stress responses depend substantially on soil parameters and soil process parameterization, i.e. hydrology and nutrient turnover, indicating that these often neglected domains deserve more scrutiny. For agricultural impact assessments, employing a GGCM ensemble with its widely varying assumptions in setups appears the best solution for coping with uncertainties from lack of comprehensive global data on crop management, cultivar distributions and coefficients for agro-environmental processes. However, the underlying assumptions require systematic specifications to cover representative agricultural systems and environmental conditions. Furthermore, the interlinkage of parameter sensitivity from various domains such as soil parameters, nutrient turnover coefficients, and cultivar specifications highlights that global sensitivity analyses and calibration need to be performed in an integrated manner to avoid bias resulting from disregarded core model domains. Finally, relating evaluations of the EPIC-based GGCMs to a wider ensemble based on individual core models shows that structural differences outweigh in general differences in configurations of GGCMs based on the same model, and that the ensemble mean gains higher skill from the inclusion of structurally different GGCMs. Although the members of the wider ensemble herein do not consider crop-soil-management interactions, their sensitivity to nutrient supply indicates that findings for the EPIC-based sub-ensemble will likely become relevant for other GGCMs with the progressing inclusion of such processes.

Folberth, Christian↗
NASA NTRS1996

Spectra of High-Ionization Seyfert 1 Galaxies: Implications for the Narrow-Line Region

We present line profiles and profile parameters for the Narrow-Line Regions (NLRs) of six Seyfert 1 galaxies with high-ionization lines: MCG 8-11-11, Mrk 79, Mrk 704, Mrk 841, NGC 4151, and NGC 5548. The sample was chosen primarily with the goal of obtaining high-quality [Fe VII] lambda6087 and, when possible, [Fe X] lambda6374 profiles to determine if these lines are more likely formed in a physically distinct 'coronal line region' or are formed throughout the NLR along with lines of lower critical density (n(sub cr)) and/or Ionization Potential (IP). We discuss correlations of velocity shift and width with n(sub cr) and IP. In some objects, lines of high IP and/or n(sub cr) are systematically broader than those of low IP/n(sub cr). Of particular interest, however, are objects that show no correlations of line width with either IP or n(sub cr). In these objects, lines of high and low IP/n(sub cr), are remarkably similar, which is difficult to reconcile with the classical picture of the NLR, in which lines of high and low IP/n(sub cr) are formed in physically distinct regions. We argue for similar spatial extents for the flux in lines with similar profiles. Here, as well as in a modeling-oriented companion paper, we develop further an idea suggested by Moore & Cohen that objects that do and do not show line width correlations with IP/n(sub cr) can both be explained in terms of a single NLR model with only a small difference in the cloud column density distinguishing the two types of object. Overall, our objects do not show correlations between the Full Width at Half-Maximum (FWHM) and IP and/or n(sub cr). The width must be defined by a parameter that is sensitive to extended profile wings in order for the correlations to result. We present models in which FWHM correlations with IP and/or n(sub cr) result only after simulating the lower spectral resolution used in previous observational studies. The models that simulate the higher spectral resolution of our observational study produce line width correlations only if the width is defined by a parameter that is more sensitive to extended profile wings than is the FWHM. Our sample of six objects is in effect augmented by incorporating the larger sample (16 objects) of Veilleux into some of our discussion. This paper focuses on new interpretations of NLR emission-line spectra and line profiles that stem directly from the observations. Paper 2 focuses on modeling and complements this paper by illustrating explicitly the effects that spatial variations in electron density, ionization parameter, and column density have on model profiles. By comparing model profiles with the observed profiles presented here, as well as with those presented by Veilleux, Paper 2 yields insight into how the electron density, ionization parameter, and column density likely vary throughout the NLR.

Moore, David↗
NASA NTRS1984

The numerical simulation of the evolution of a Saharan dust outbreak by using a general prognostic aerosol model

A three-dimensional model is developed that explicitly predicts the evolution of the spatial distributions of aerosols utilizing a time-dependent, initial value prognostic aerosol model based on aerosol physicochemical processes, a knowledge of aerosol sinks and sources, and the predicted motions of the atmosphere. These models serve a useful purpose in determining standardized aerosol spectra for utilization in radiative transfer calculations. The first numerical prediction of a Saharan dust outbreak evolution is described. It is shown that the results are in qualitative agreement with measurements despite the fact that in this preliminary simulation the generation of the dust depended upon several empirically based assumptions. IR satellite photographs taken during GATE are being examined in an attempt to determine the sources of the dust.

Westphal, Douglas L.↗
NASA NTRS2016

Modeling Forest Biomass and Growth: Coupling Long-Term Inventory and Lidar Data

Combining spatially-explicit long-term forest inventory and remotely sensed information from Light Detection and Ranging (LiDAR) datasets through statistical models can be a powerful tool for predicting and mapping above-ground biomass (AGB) at a range of geographic scales. We present and examine a novel modeling approach to improve prediction of AGB and estimate AGB growth using LiDAR data. The proposed model accommodates temporal misalignment between field measurements and remotely sensed data-a problem pervasive in such settings-by including multiple time-indexed measurements at plot locations to estimate AGB growth. We pursue a Bayesian modeling framework that allows for appropriately complex parameter associations and uncertainty propagation through to prediction. Specifically, we identify a space-varying coefficients model to predict and map AGB and its associated growth simultaneously. The proposed model is assessed using LiDAR data acquired from NASA Goddard's LiDAR, Hyper-spectral & Thermal imager and field inventory data from the Penobscot Experimental Forest in Bradley, Maine. The proposed model outperformed the time-invariant counterpart models in predictive performance as indicated by a substantial reduction in root mean squared error. The proposed model adequately accounts for temporal misalignment through the estimation of forest AGB growth and accommodates residual spatial dependence. Results from this analysis suggest that future AGB models informed using remotely sensed data, such as LiDAR, may be improved by adapting traditional modeling frameworks to account for temporal misalignment and spatial dependence using random effects.

Babcock, Chad↗
NASA NTRS1993

Effects of spatial variability and scale on areal -average evapotranspiration

This paper explores the effect of spatial variability and scale on areally-averaged evapotranspiration. A spatially-distributed water and energy balance model is employed to determine the effect of explicit patterns of model parameters and atmospheric forcing on modeled areally-averaged evapotranspiration over a range of increasing spatial scales. The analysis is performed from the local scale to the catchment scale. The study area is King's Creek catchment, an 11.7 sq km watershed located on the native tallgrass prairie of Kansas. The dominant controls on the scaling behavior of catchment-average evapotranspiration are investigated by simulation, as is the existence of a threshold scale for evapotranspiration modeling, with implications for explicit versus statistical representation of important process controls. It appears that some of our findings are fairly general, and will therefore provide a framework for understanding the scaling behavior of areally-averaged evapotranspiration at the catchment and larger scales.

Famiglietti, J. S.↗
NASA NTRS2013

Representation of Nucleation Mode Microphysics in a Global Aerosol Model with Sectional Microphysics

In models, nucleation mode (1 nm<Dp <10 nm) particle microphysics can be represented explicitly with aerosol microphysical processes or can be parameterized to obtain the growth and survival of nuclei to the model's lower size boundary. This study investigates how the representation of nucleation mode microphysics impacts aerosol number predictions in the TwO-Moment Aerosol Sectional (TOMAS) aerosol microphysics model running with the GISS GCM II-prime by varying its lowest diameter boundary: 1 nm, 3 nm, and 10 nm. The model with the 1 nm boundary simulates the nucleation mode particles with fully resolved microphysical processes, while the model with the 10 nm and 3 nm boundaries uses a nucleation mode dynamics parameterization to account for the growth of nucleated particles to 10 nm and 3 nm, respectively.We also investigate the impact of the time step for aerosol microphysical processes (a 10 min versus a 1 h time step) to aerosol number predictions in the TOMAS models with explicit dynamics for the nucleation mode particles (i.e., 3 nm and 1 nm boundary). The model with the explicit microphysics (i.e., 1 nm boundary) with the 10 min time step is used as a numerical benchmark simulation to estimate biases caused by varying the lower size cutoff and the time step. Different representations of the nucleation mode have a significant effect on the formation rate of particles larger than 10 nm from nucleated particles (J10) and the burdens and lifetimes of ultrafinemode (10 nm≤Dp ≤70 nm) particles but have less impact on the burdens and lifetimes of CCN-sized particles. The models using parameterized microphysics (i.e., 10 nm and 3 nm boundaries) result in higher J10 and shorter coagulation lifetimes of ultrafine-mode particles than the model with explicit dynamics (i.e., 1 nm boundary). The spatial distributions of CN10 (Dp ≥10 nm) and CCN(0.2 %) (i.e., CCN concentrations at 0.2%supersaturation) are moderately affected, especially CN10 predictions above ∼700 hPa where nucleation contributes most strongly to CN10 concentrations. The lowermost-layer CN10 is substantially improved with the 3 nm boundary (compared to 10 nm) in most areas. The overprediction in CN10 with the 3 nm and 10 nm boundaries can be explained by the overprediction of J10 or J3 with the parameterized microphysics, possibly due to the instantaneous growth rate assumption in the survival and growth parameterization. The errors in CN10 predictions are sensitive to the choice of the lower size boundary but not to the choice of the time step applied to the microphysical processes. The spatial distribution of CCN(0.2 %) with the 3 nm boundary is almost identical to that with the 1 nm boundary, but that with the 10 nm boundary can differ more than 10-40% in some areas. We found that the deviation in the 10 nm simulations is partly due to the longer time step (i.e., 1 h time step used in the 10 nm simulations compared to 10 min time step used in the benchmark simulations), but, even with the same time step, the 10 nm cutoff showed noticeably higher errors than the 3 nm cutoff. In conclusion, we generally recommend using a lower diameter boundary of 3 nm for studies focused on aerosol indirect effects but down to 1 nm boundary for studies focused on CN10 predictions or nucleation.

microphysics↗
NASA NTRS2020

A 3D Model to Predict Explicit Morphologies and Volume Fraction of Lack-of-Fusion Pores Generated in Selective Laser Melting Processes

The performance of an additively manufactured (AM) component is dependent on the distribution of process-induced defects in addition to the complex microstructure of the material, surface roughness of the component and the process-induced residual stresses. For instance, it has been well demonstrated that lack-of-fusion (LoF) pores produced in the selective laser melting (SLM) AM process can significantly limit the fatigue performance of the material. Although two-dimensional (2D) models exist to predict the 2D profiles of LoF pores, the 2D pore profiles cannot be directly inserted into a three-dimensional (3D) microstructure domain that is output from several prevailing process simulation packages. A few commercial packages that simulate the SLM process can predict LoF pores in a 3D domain, but the morphologies of LoF pores are voxelated and hence do not capture sharp corners of the pores, thereby obviating their use in fatigue crack initiation studies. In order to address the aforementioned gaps, a high-fidelity model that predicts not only the volume fraction, but also the explicit 3D morphologies and spatial distributions of LoF pores has been developed using a computer aided design-based environment. The model has been partially validated for Ti-6Al-4V alloy by comparing the predictions of the volume fraction of LoF pores predicted by the model with experimental data obtained from the literature. Absolute error in predicted volume fraction of LoF pores varied between 5.16% and 1.87% for energy density values between 13 J/mm3 and 45 J/mm3 where a significant amount (over 3% volume fraction) of LoF porosity was measured. The absolute error was within 1.87% for energy density values greater than 45 J/mm3.

Saikumar R. Yeratapally↗
NASA NTRS2018

The Use of Gridded Fossil Fuel CO2 Emissions (FFCO2) Inventory for Climate Mitigation Applications: Errors, Uncertainties, and Current and Future Challenges

Emission Inventory (EI) is a fundamental tool to monitor global compliance of greenhouse gases (GHGs) emissions reduction actions. Inventory guidelines provide a best practice to help EI compilers to make comparable national emission estimates, in spite of the differences in data availability across countries and regions. There are a variety of sources of errors and uncertainties, however, that originate beyond what the inventory guidelines define. For example, spatially-explicit EIs, which are a key product for atmospheric modeling applications, are often developed for research purposes, and there are no specific guidelines to disaggregate emission estimates from country scale. On top of that, EIs are fundamentally prone to systematic biases due to the simple calculation methodology and thus an objective evaluation (e.g. atmospheric top-down estimates) is needed to assure the accuracy of the estimates. ODIAC is a global high-resolution (1x1 km) fossil fuel carbon dioxide (CO2) gridded EI that is now often used in atmospheric CO2 modeling. ODIAC is based on disaggregation of national emission estimates made by CDIAC, which is the well accepted standard in the community. The ODIAC emission data product is updated on an annual basis using best available statistical data. Subnational spatial emission patterns are estimated using power plant profiles and satellite-observations of nighttime lights. In addition to the conventional CDIAC gridded data product, ODIAC carries international bunker emissions (shipping and aviation), which allows flux inversion modelers to accurately impose the global total fossil fuel emissions and their horizontal and vertical distribution. We have extensively evaluated ODIAC emissions using fine-grained EIs as well as a high-resolution atmospheric model simulation across different scales (national, subnational/regional, and urban policy relevant) with a focus on the uncertainties associated with the emission disaggregation. We have examined the use of NASA's Black Marble Suomi-NPP/VIIRS nightlight data.

Oda, Tomohiro↗
NASA NTRS1991

Nonlinear influence of mesoscale land use on weather and climate

It is shown that the influence of mesoscale landscape spatial variability on the atmosphere must be parameterized or explicitly modeled in larger-scale atmospheric model simulations including general circulation models. The mesoscale fluxes of heat that result from this variability are shown to be of the same order of magnitude but with a different vertical structure than found for the turbulent fluxes. These conclusions are based on experiments in which no phase changes of water were permitted. To parameterize surface thermal inhomogeneities, the influence of landscape must be evaluated using spectral analysis or an equivalent procedure. To include the nonlinear contribution of each scale, numerical model simulations for the range of observed surface and overlying atmospheric conditions must be performed.

Pielke, R. A.↗
NASA NTRS2024

Quantifying Food Security and Mitigation Risks Consequential to Climate Change Impacts on Crop Yields

Climate change is expected to impact crop yields globally, with some regions benefiting from favorable conditions and CO2 fertilization, while others face adverse effects from altered precipitation and higher temperatures. Changes in crop yields can destabilize the global food system and pose challenges to food security. Moreover, crop production is crucial, as biofuels are becoming increasingly important contributors to climate change mitigation measures aimed at limiting global warming. This study uses the Integrated Model to Assess the Global Environment integrated assessment model framework to analyze different indicators related to food security and climate change mitigation under varying climate change impacts on crop yields. Twelve spatially explicit crop productivity projections were taken from the full archive of the Global Gridded Crop Model Intercomparison of 120 climate-crop model combinations, forced by CMIP6-based climate scenarios. The selection includes two average-performing climate-crop model combinations, two pessimistic combinations that perform one standard deviation below the mean, and two optimistic model combinations that perform one standard deviation above the mean. To single out the effect of climate change on productivity changes, we drew samples from two representative concentration pathways (RCP2.6 and RCP8.5). These productivity projections were applied within an otherwise uniform scenario (SSP2) and analyzed for their effect on total calorie demand, crop prices, and number of people at risk of undernourishment to quantify food security. Risks to climate change mitigation targets were explored by modeling the total bioenergy supply, emissions, and global mean temperature. The results revealed significant differences in the risk of food security and mitigation potential between different regions and climate change scenarios. Across scenarios, the crop area extent can vary up to 2 million km2 due to changing crop yields. The projected change in global hunger ranges from 60 to 160 million undernourished people, indicating uncertainty between climate and crop model combinations. Low-income regions are especially impacted because of their high sensitivity to changes in food prices. Global climate change mitigation ambitions can also deviate by the latter part of the 21st century, as changes in yields will impact biofuel production as well as agriculture, forestry and other land use emissions. The quantitative insights generated by this study highlight the need for global policy efforts to make the agricultural system more adaptive to climate change to handle potential negative impacts.

crop yields↗
NASA NTRS2005

Spatial Predictive Modeling and Remote Sensing of Land Use Change in the Chesapeake Bay Watershed

This project was focused on modeling the processes by which increasing demand for developed land uses, brought about by changes in the regional economy and the socio-demographics of the region, are translated into a changing spatial pattern of land use. Our study focused on a portion of the Chesapeake Bay Watershed where the spatial patterns of sprawl represent a set of conditions generally prevalent in much of the U.S. Working in the region permitted us access to (i) a time-series of multi-scale and multi-temporal (including historical) satellite imagery and (ii) an established network of collaborating partners and agencies willing to share resources and to utilize developed techniques and model results. In addition, a unique parcel-level tax assessment database and linked parcel boundary maps exists for two counties in the Maryland portion of this region that made it possible to establish a historical cross-section time-series database of parcel level development decisions. Scenario analyses of future land use dynamics provided critical quantitative insight into the impact of alternative land management and policy decisions. These also have been specifically aimed at addressing growth control policies aimed at curbing exurban (sprawl) development. Our initial technical approach included three components: (i) spatial econometric modeling of the development decision, (ii) remote sensing of suburban change and residential land use density, including comparisons of past change from Landsat analyses and more traditional sources, and (iii) linkages between the two through variable initialization and supplementation of parcel level data. To these we added a fourth component, (iv) cellular automata modeling of urbanization, which proved to be a valuable addition to the project. This project has generated both remote sensing and spatially explicit socio-economic data to estimate and calibrate the parameters for two different types of land use change models and has undertaken analyses of these models. One (the CA model) is driven largely by observations on past patterns of land use change, while the other (the EC model) is driven by mechanisms of the land use change decision at the parcel level. Our project may be the first serious attempt at developing both types of models for the same area, using as much common data as possible. We have identified the strengths and weaknesses of the two approaches and plan to continue to revise each model in the light of new data and new lessons learned through continued collaboration. Questions, approaches, findings, publication and presentation lists concerning the research are also presented.

Goetz, Scott J.↗
NASA NTRS2018

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

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

Drought↗
NASA NTRS1982

ADI on staggered mesh - A method for the calculation of compressible convection

An alternating direction implicit (ADI) method has been applied to a staggered grid for the computation of convection in a highly stratified fluid. Since artificial viscosity is not needed, subtle effects like the onset of convection can be studied. These computations compare well with the 2-D results by Graham (1975) and also agree with standard Boussinesq results when taken to that limit. Good efficiency has been achieved with a time step hundreds of times larger than the stability limit imposed by the explicit treatment of diffusion and the Courant number is not restricted to be below 1. The Navier-Stokes equation contains cross spatial derivatives which are treated explicitly in most ADI schemes. The destabilizing effect of such a practice on a 2-D model system with second-order spatial derivative terms only was analyzed and found to be not excessive.

Chan, K. L.↗