Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear Regression”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Re-Greening Sahel: 30 Years of Remote Sensing Data and Field Observations (Mali, Niger)

Desertification of the Sahel region has been debated for decades, while the concept of a "re-greening" Sahel appeared with satellite remote sensing data that allowed vegetation monitoring across wide regions and over increasingly long series of years (nowadays 30 years with the GIMMS-3g dataset). However, the scarcity of long-term field observations of vegetation in the Sahel prevents ground validation and deeper analysis of such trends. After assessing the consistency of the new GIMMS-3g NDVI product by comparison to three other AVHRR-NDVI datasets and MODIS NDVI, regional GIMMS-3g NDVI trends over 1981-2011 are analyzed. Trends are found positive and statistically significant almost everywhere in Sahel over the 1981-2011 period. Long-term field observations of the aboveground herbaceous layer mass have been collected within the Gourma region in Mali (1984-2011) and within the Fakara region in western Niger (1994-2011). These observations sample ecosystem and soil diversity, thus enabling estimation of averaged values representative of the Gourma and Fakara. NDVI measurements are found in good agreement with field observations, both over the Gourma and Fakara regions where re-greening and negative trends are observed respectively. A linear regression analysis performed between spatially averaged seasonal NDVI and a weighted average of field measurements explains 59% of the variability for the Gourma region over 1984-2011, and 38% for the Fakara region over 1994-2011. In the Gourma, which is a pastoral region, the re-greening trend is mainly observed over sandy soils, and attests for the ecosystem's resilience to the 1980s' drought, able to react to the more favorable rainfall of the 1990s and 2000s. However, contrasted changes in the landscape's functioning have occurred locally. An increase in erosion and run-off processes in association with decreasing or stable vegetation cover was observed over shallow soils, which occupy 30% of the area. In the agro-pastoral Fakara, the decreasing trends observed both from satellite NDVI and field assessments of herbaceous mass are hardly explained by rainfall. These results give confidence in the dominant positive trends in Sahelian greenness, but indicate that degradation trends can also be observed, both in situ and from satellite time series.

Sahel↗
NASA NTRS2015

Solar Signals in CMIP-5 Simulations: The Stratospheric Pathway

The 11 year solar-cycle component of climate variability is assessed in historical simulations of models taken from the Coupled Model Intercomparison Project, phase 5 (CMIP-5). Multiple linear regression is applied to estimate the zonal temperature, wind and annular mode responses to a typical solar cycle, with a focus on both the stratosphere and the stratospheric influence on the surface over the period approximately 1850-2005. The analysis is performed on all CMIP-5 models but focuses on the 13 CMIP-5 models that resolve the stratosphere (high-top models) and compares the simulated solar cycle signature with reanalysis data. The 11 year solar cycle component of climate variability is found to be weaker in terms of magnitude and latitudinal gradient around the stratopause in the models than in the reanalysis. The peak in temperature in the lower equatorial stratosphere (approximately 70 hPa) reported in some studies is found in the models to depend on the length of the analysis period, with the last 30 years yielding the strongest response. A modification of the Polar Jet Oscillation (PJO) in response to the 11 year solar cycle is not robust across all models, but is more apparent in models with high spectral resolution in the short-wave region. The PJO evolution is slower in these models, leading to a stronger response during February, whereas observations indicate it to be weaker. In early winter, the magnitude of the modeled response is more consistent with observations when only data from 1979-2005 are considered. The observed North Pacific high-pressure surface response during the solar maximum is only simulated in some models, for which there are no distinguishing model characteristics. The lagged North Atlantic surface response is reproduced in both high- and low-top models, but is more prevalent in the former. In both cases, the magnitude of the response is generally lower than in observations.

stratosphere↗
NASA NTRS2015

A Combined Solar and Geomagnetic Index for Thermospheric Climate

Infrared radiation from nitric oxide (NO) at 5.3 Â is a primary mechanism by which the thermosphere cools to space. The SABER instrument on the NASA TIMED satellite has been measuring thermospheric cooling by NO for over 13 years. Physically, changes in NO emission are due to changes in temperature, atomic oxygen, and the NO density. These physical changes however are driven by changes in solar irradiance and changes in geomagnetic conditions. We show that the SABER time series of globally integrated infrared power (Watts) radiated by NO can be replicated accurately by a multiple linear regression fit using the F10.7, Ap, and Dst indices. This fit enables several fundamental properties of NO cooling to be determined as well as their variability with time, permitting reconstruction of the NO power time series back nearly 70 years with extant databases of these indices. The relative roles of solar ultraviolet and geomagnetic processes in determining the NO cooling are derived and shown to be solar cycle dependent. This reconstruction provides a long-term time series of an integral radiative constraint on thermospheric climate that can be used to test climate models.

Hunt, Linda↗
NASA NTRS2017

Carotid Intima Media Thickness in the Astronaut Corps: Association to Spacecraft

Background: Carotid Intima Media Thickness (CIMT) has been demonstrated to be predictive of future cardiovascular events. Within various populations, radiation exposure, stress, and physical confinement have all been linked to an increased CIMT. Recent research discovered CIMT was significantly increased in ten long duration astronauts from pre-flight to four days post flight. The relationship between spaceflight and CIMT is not understood and trends in CIMT within the larger astronaut population are unknown. Methods: In 2010, CIMT was offered as part of the astronaut annual exam at the JSC Flight Medicine Clinic using a standardized CIMT screening protocol and professional sonographers. Between 2010 and 2016, CIMT measurements were collected on 213 NASA astronauts and payload specialists. The values used in this retrospective chart review are the mean of the CIMT from the right and left. Spaceflight exposure was categorized based on the total number of days spent in space at the time of the ground-based ultrasound (0, 1-29, 30-100, 101-200, ≥200). Linear regression with generalized estimating equations were used to estimate the association between spaceflight exposures and CIMT. Results: 530 studies were completed among 213 astronauts with a mean of 2.5 studies (range 1-6) per astronaut over the six year period. As in other populations, CIMT was significantly associated with age; however, gender was not. While there was no significant direct correlation between total spaceflight exposure and CIMT found, astronauts with 30-100 spaceflight days and astronauts with greater than 100 spaceflight days had significantly increased CIMT over astronauts who had never flown (p=0.002 and p=<0.0001 respectively) after adjustment for age. Conclusion: Further work is needed to fully understand CIMT and its association to spaceflight. Current occupational surveillance activities are under way to study CIMT values in conjunction with other cardiovascular risk factors among astronauts as compared to the general population.

Suffredini, John↗
NASA NTRS2016

Quantifying PM2.5-Meteorology Sensitivities in a Global Climate Model

Climate change can influence fine particulate matter concentrations (PM2.5) through changes in air pollution meteorology. Knowledge of the extent to which climate change can exacerbate or alleviate air pollution in the future is needed for robust climate and air pollution policy decision-making. To examine the influence of climate on PM2.5, we use the Geophysical Fluid Dynamics Laboratory Coupled Model version 3 (GFDL CM3), a fully-coupled chemistry-climate model, combined with future emissions and concentrations provided by the four Representative Concentration Pathways (RCPs). For each of the RCPs, we conduct future simulations in which emissions of aerosols and their precursors are held at 2005 levels while other climate forcing agents evolve in time, such that only climate (and thus meteorology) can influence PM2.5 surface concentrations. We find a small increase in global, annual mean PM2.5 of about 0.21 micro-g/cu m3 (5%) for RCP8.5, a scenario with maximum warming. Changes in global mean PM2.5 are at a maximum in the fall and are mainly controlled by sulfate followed by organic aerosol with minimal influence of black carbon. RCP2.6 is the only scenario that projects a decrease in global PM2.5 with future climate changes, albeit only by -0.06 micro-g/cu m (1.5%) by the end of the 21st century. Regional and local changes in PM2.5 are larger, reaching upwards of 2 micro-g/cu m for polluted (eastern China) and dusty (western Africa) locations on an annually averaged basis in RCP8.5. Using multiple linear regression, we find that future PM2.5 concentrations are most sensitive to local temperature, followed by surface wind and precipitation. PM2.5 concentrations are robustly positively associated with temperature, while negatively related with precipitation and wind speed. Present-day (2006-2015) modeled sensitivities of PM2.5 to meteorological variables are evaluated against observations and found to agree reasonably well with observed sensitivities (within 10e50% over the eastern United States for several variables), although the modeled PM2.5 is less sensitive to precipitation than in the observations due to weaker convective scavenging. We conclude that the hypothesized "climate penalty" of future increases in PM2.5 is relatively minor on a global scale compared to the influence of emissions on PM2.5 concentrations.

PM2.5↗
NASA NTRS2016

Modeling Relationships Between Flight Crew Demographics and Perceptions of Interval Management

The Interval Management Alternative Clearances (IMAC) human-in-the-loop simulation experiment was conducted to assess interval management system performance and participants' acceptability and workload while performing three interval management clearance types. Twenty-four subject pilots and eight subject controllers flew ten high-density arrival scenarios into Denver International Airport during two weeks of data collection. This analysis examined the possible relationships between subject pilot demographics on reported perceptions of interval management in IMAC. Multiple linear regression models were created with a new software tool to predict subject pilot questionnaire item responses from demographic information. General patterns were noted across models that may indicate flight crew demographics influence perceptions of interval management.

Remy, Benjamin↗
NASA NTRS2016

GOCI Yonsei Aerosol Retrieval (YAER) Algorithm and Validation During the DRAGON-NE Asia 2012 Campaign

The Geostationary Ocean Color Imager (GOCI) onboard the Communication, Ocean, and Meteorological Satellite (COMS) is the first multi-channel ocean color imager in geostationary orbit. Hourly GOCI top-of-atmosphere radiance has been available for the retrieval of aerosol optical properties over East Asia since March 2011. This study presents improvements made to the GOCI Yonsei Aerosol Retrieval (YAER) algorithm together with validation results during the Distributed Regional Aerosol Gridded Observation Networks - Northeast Asia 2012 campaign (DRAGONNE Asia 2012 campaign). The evaluation during the spring season over East Asia is important because of high aerosol concentrations and diverse types of Asian dust and haze. Optical properties of aerosol are retrieved from the GOCI YAER algorithm including aerosol optical depth (AOD) at 550 nm, fine-mode fraction (FMF) at 550 nm, single-scattering albedo (SSA) at 440 nm, Angstrom exponent (AE) between 440 and 860 nm, and aerosol type. The aerosol models are created based on a global analysis of the Aerosol Robotic Networks (AERONET) inversion data, and covers a broad range of size distribution and absorptivity, including nonspherical dust properties. The Cox-Munk ocean bidirectional reflectance distribution function (BRDF) model is used over ocean, and an improved minimum reflectance technique is used over land. Because turbid water is persistent over the Yellow Sea, the land algorithm is used for such cases. The aerosol products are evaluated against AERONET observations and MODIS Collection 6 aerosol products retrieved from Dark Target (DT) and Deep Blue (DB) algorithms during the DRAGON-NE Asia 2012 campaign conducted from March to May 2012. Comparison of AOD from GOCI and AERONET resulted in a Pearson correlation coefficient of 0.881 and a linear regression equation with GOCI AOD = 1.083 x AERONET AOD - 0.042. The correlation between GOCI and MODIS AODs is higher over ocean than land. GOCI AOD shows better agreement with MODIS DB than MODIS DT. The other GOCI YAER products (AE, FMF, and SSA) show lower correlation with AERONET than AOD, but still show some skills for qualitative use.

Choi, Myungje↗
NASA NTRS2017

Changes in the Structure and Propagation of the MJO with Increasing CO2

Changes in the Madden-Julian Oscillation (MJO) with increasing CO2 concentrations are examined using the Goddard Institute for Space Studies Global Climate Model (GCM). Four simulations performed with fixed CO2 concentrations of 0.5, 1, 2 and 4 times pre-industrial levels using the GCM coupled with a mixed layer ocean model are analyzed in terms of the basic state, rainfall and moisture variability, and the structure and propagation of the MJO.The GCM simulates basic state changes associated with increasing CO2 that are consistent with results from earlier studies: column water vapor increases at approximately 7.1% K(exp -1), precipitation also increases but at a lower rate (approximately 3% K(exp -1)), and column relative humidity shows little change. Moisture and rainfall variability intensify with warming. Total moisture and rainfall variability increases at a rate that is similar to that of the mean state change. The intensification is faster in the MJO-related anomalies than in the total anomalies, though the ratio of the MJO band variability to its westward counterpart increases at a much slower rate. On the basis of linear regression analysis and space-time spectral analysis, it is found that the MJO exhibits faster eastward propagation, faster westward energy dispersion, a larger zonal scale and deeper vertical structure in warmer climates.

Adames, Angel F.↗
NASA NTRS2016

An Evaluation of MODIS-Retrieved Aerosol Optical Depth over a Mountainous AERONET Site in the Southeastern US

The literature shows that aerosol optical depth (AOD) derived from the MODIS Collection 5 (C5) dark target algorithm has been extensively validated by spatiotemporal collocation with AERONET sites on both global and regional scales.Although generally comparing well over the eastern US region, poor performance over mountains in other regions indicate the need to evaluate the MODIS product over a mountain site. This study compares MODIS C5 AOD at 550nm to AOD measured at the Appalachian State University AERONET site in Boone, NC over 30 months between August 2010 and September 2013. For the combined Aqua and Terra datasets, although more than 70% of the 500 MODIS AOD measurements agree with collocated AERONET AOD to within error envelope of +/- (0.05 + 15%), MODIS tends to have a low bias (0.02-0.03). The agreement between MODIS and AERONET AOD does not depend on MODIS quality assurance confidence (QAC) value. However, when stratified by satellite, MODIS-Terra data does not perform as well as Aqua, with especially poor correlation (r = 0.39) for low aerosol loading conditions (AERONET AOD less than 0.15).Linear regressions between Terra and AERONET possess statistically-different slopes for AOD < 0.15 and AOD > or = 0.15. AERONET AOD measured only during MODIS overpass hours is highly correlated with daily-averaged AERONET AOD. MODIS monthly-averaged AOD also tracks that of AERONET over the study period. These results indicate that MODIS is sensitive to the day-to-day variability, as well as the annual cycle of AOD over the Appalachian State AERONET site. The complex topography and high seasonality in AOD and vegetation indices allow us to specifically evaluate MODIS dark target algorithm surface albedo and aerosol model assumptions at a regionally-representative SE US mountain site.

air quality↗
NASA NTRS2015

NASA Software Cost Estimation Model: An Analogy Based Estimation Model

The cost estimation of software development activities is increasingly critical for large scale integrated projects such as those at DOD and NASA especially as the software systems become larger and more complex. As an example MSL (Mars Scientific Laboratory) developed at the Jet Propulsion Laboratory launched with over 2 million lines of code making it the largest robotic spacecraft ever flown (Based on the size of the software). Software development activities are also notorious for their cost growth, with NASA flight software averaging over 50% cost growth. All across the agency, estimators and analysts are increasingly being tasked to develop reliable cost estimates in support of program planning and execution. While there has been extensive work on improving parametric methods there is very little focus on the use of models based on analogy and clustering algorithms. In this paper we summarize our findings on effort/cost model estimation and model development based on ten years of software effort estimation research using data mining and machine learning methods to develop estimation models based on analogy and clustering. The NASA Software Cost Model performance is evaluated by comparing it to COCOMO II, linear regression, and K-­ nearest neighbor prediction model performance on the same data set.

Hihn, Jairus↗
NASA NTRS2018

Calibration of the DSCOVR EPIC Visible and NIR Channels using MODIS Terra and Aqua Data and EPIC Lunar Observations

The unique position of the Deep Space Climate Observatory (DSCOVR) Earth Polychromatic Imaging Camera (EPIC) at the Lagrange 1 point makes an important addition to the data from currently operating low Earth orbit observing instruments. EPIC instrument does not have an onboard calibration facility. One approach to its calibration is to compare EPIC observations to the measurements from polar-orbiting radiometers. Moderate Resolution Imaging Spectroradiometer (MODIS) is a natural choice for such comparison due to its well-established calibration record and wide use in remote sensing. We use MODIS Aqua and Terra L1B 1km reflectances to infer calibration coefficients for four EPIC visible and NIR channels: 443, 551, 680 and 780 nm. MODIS and EPIC measurements made between June 2015 and 2016 are employed for comparison. We first identify favorable MODIS pixels with scattering angle matching temporarily collocated EPIC observations. Each EPIC pixel is then spatially collocated to a subset of the favorable MODIS pixels within 25 km radius. Standard deviation of the selected MODIS pixels as well as of the adjacent EPIC pixels is used to find the most homogeneous scenes. These scenes are then used to determine calibration coefficients using a linear regression between EPIC counts/sec and reflectances in the close MODIS spectral channels. We present thus inferred EPIC calibration coefficients and discuss sources of uncertainties. The lunar EPIC observations are used to calibrate EPIC O2 absorbing channels (688 and 764 nm), assuming that there is a small difference between moon reflectances separated by approx.10 nm in wavelength provided the calibration factors of the red (680 nm) and near-IR (780 nm) are known from comparison between EPIC and MODIS.

Geogdzhayev, Igor V.↗
NASA NTRS2017

A Consistent EPIC Visible Channel Calibration using VIIRS and MODIS as a Reference

The Earth Polychromatic Imaging Camera (EPIC) aboard the Deep Space Climate Observatory (DSCOVR) satellite constantly images the sunlit disk of Earth from the Lagrange-1 (L1) point in 10 spectral channels spanning the UV, VIS, and NIR spectrums. Recently, the DSCOVR EPIC team has publicly released version 2 dataset, which has implemented improved navigation, stray-light correction, and flat-fielding of the CCD array. The EPIC 2-year data record must be well-calibrated for consistent cloud, aerosol, trace gas, land use and other retrievals. Because EPIC lacks onboard calibrators, the observations made by EPIC channels must be calibrated vicariously using the coincident measurements from radiometrically stable instruments that have onboard calibration systems. MODIS and VIIRS are best-suited instruments for this task as they contain similar spectral bands that are well-calibrated onboard using solar diffusers and lunar tracking. We have previously calibrated the EPIC version 1 dataset by using EPIC and VIIRS angularly matched radiance pairs over both all-sky ocean and deep convective clouds (DCC). We noted that the EPIC image required navigations adjustments, and that the EPIC stray-light correction provided an offset term closer to zero based on the linear regression of the EPIC and VIIRS ray-matched radiance pairs. We will evaluate the EPIC version 2 navigation and stray-light improvements using the same techniques. In addition, we will monitor the EPIC channel calibration over the two years for any temporal degradation or anomalous behavior. These two calibration methods will be further validated using desert and DCC invariant Earth targets. The radiometric characterization of the selected invariant targets is performed using multiple years of MODIS and VIIRS measurements. Results of these studies will be shown at the conference.

Haney, C. O.↗
NASA NTRS2018

Monitoring Springs in the Mojave Desert Using Landsat Time Series Analysis

The purpose of this study, based on Landsat satellite data was to characterize variations and trends over 30 consecutive years (1985-2016) in perennial vegetation green cover at over 400 confirmed Mojave Desert spring locations. These springs were surveyed between in 2015 and 2016 on lands managed in California by the U.S. Bureau of Land Management (BLM) and on several land trusts within the Barstow, Needles, and Ridgecrest BLM Field Offices. The normalized difference vegetation index (NDVI) from July Landsat images was computed at each spring location and a trend model was first fit to the multi-year NDVI time series using least squares linear regression.Â

Time Series↗
NASA NTRS2018

Total Ozone Trends from 1979 to 2016 Derived from Five Merged Observational Datasets - The Emergence into Ozone Recovery

We report on updated trends using different merged datasets from satellite and ground-based observations for the period from 1979 to 2016. Trends were determined by applying a multiple linear regression (MLR) to annual mean zonal mean data. Merged datasets used here include NASA MOD v8.6 and National Oceanic and Atmospheric Administration (NOAA) merge v8.6, both based on data from the series of Solar Backscatter UltraViolet (SBUV) and SBUV-2 satellite instruments (1978–present) as well as the Global Ozone Monitoring Experiment (GOME)-type Total Ozone (GTO) and GOME-SCIAMACHY-GOME-2 (GSG) merged datasets (1995-present), mainly comprising satellite data from GOME, the Scanning Imaging Absorption Spectrometer for Atmospheric Chartography (SCIAMACHY), and GOME-2A. The fifth dataset consists of the monthly mean zonal mean data from ground-based measurements collected at World Ozone and UV Data Center (WOUDC). The addition of four more years of data since the last World Meteorological Organization (WMO) ozone assessment (2013-2016) shows that for most datasets and regions the trends since the stratospheric halogen reached its maximum (approximately 1996 globally and approximately 2000 in polar regions) are mostly not significantly different from zero. However, for some latitudes, in particular the Southern Hemisphere extratropics and Northern Hemisphere subtropics, several datasets show small positive trends of slightly below +1 percent decade(exp. -1) that are barely statistically significant at the 2 Sigma uncertainty level. In the tropics, only two datasets show significant trends of +0.5 to +0.8 percent(exp.-1), while the others show near-zero trends. Positive trends since 2000 have been observed over Antarctica in September, but near-zero trends are found in October as well as in March over the Arctic. Uncertainties due to possible drifts between the datasets, from the merging procedure used to combine satellite datasets and related to the low sampling of ground-based data, are not accounted for in the trend analysis. Consequently, the retrieved trends can be only considered to be at the brink of becoming significant, but there are indications that we are about to emerge into the expected recovery phase. However, the recent trends are still considerably masked by the observed large year-to-year dynamical variability in total ozone.

column ozone trends↗
NASA NTRS2018

Preliminary Error Characterization and Parametric Error Model for the Automatic Dependent Surveillance - Contract Extended Projected Profile Message

A critical component of Trajectory-Based Operations (TBO) is the ability for a consistent and accurate 4-dimensional trajectory (4DT) to be shared and synchronized between airborne and ground systems as well as amongst various ground automation systems. The Aeronautical Telecommunication Network—Baseline 2 (ATN-B2) standard defines the Extended Projected Profile (EPP) trajectory that can be sent via Automatic Dependent Surveillance-Contract (ADS-C) from an aircraft to ground automation. The EPP trajectory message contains a representation of the reference trajectory from an aircraft’s Flight Management System (FMS). In this work, a set of scenarios were run in a high-fidelity aircraft and FMS simulation to perform an initial characterization of EPP trajectory errors under a given set of conditions. The parameters investigated were the route length, route type, wind magnitude error, wind direction error, and with and without a required time-of-arrival (RTA) constraint. In addition, linear regression was used to identify a set of EPP error models for the cross-track, vertical, and time errors.

Guerreiro, Nelson M.↗
NASA NTRS2018

Agricultural Burning and Air Quality over Northern India: A Synergistic Analysis Using NASA's A-Train Satellite Data and Ground Measurements

In the recent years, New Delhi, the capital city of India, has ranked among the most polluted cities in the world regarding its air quality related to the submicron Particulate Matter (PM(sub 2.5)). Using NASA's A-train satellite data (MODIS, OMI, and CALIOP), ground-level PM2.5 measured in New Delhi (2013-2016), and back-trajectory calculations, we show that the PM(sub 2.5). over New Delhi is strongly affected by the agricultural fires in the northwestern Indian states of Punjab and Haryana during the post-monsoon season (October and November). The mass concentration of PM2.5 escalates from ~50 μg m(exp -3) measured prior to the onset of residue burning in early October to as high as 300 μg m(exp -3) (24-hour averaged, 7- day running mean) during the peak burning period in early November. A linear regression analysis reveals that the variations in PM(sub 2.5) over New Delhi can be attributed to the concurrent changes in the satellite retrievals of fire counts and aerosols over the crop burning area. The back-trajectory analysis shows that most clusters (> 80%) of the northwesterly flow near the ground intercepted the crop burning region before arriving at the receptor location in New Delhi; this further corroborates the transport patterns inferred from the satellite data. A 15-year long satellite record (2002-2016) reveals an increasing trend in agricultural fires (~617 per year) and aerosol loading (0.031 and 0.04 per year in aerosol optical depth and UV aerosol index) in November. Increasing levels of crop residue burning and resulting particulate matter pollution at an alarming rate over northern India is a pressing concern demanding corrective measures to substantially reduce or completely diminish the crop burning through an effective residue management system.

Crop Residue Fires↗
NASA NTRS2016

The Viability of Trajectory Analysis for Diagnosing Dynamical and Chemical Influences on Ozone Concentrations in the UTLS

The viability of trajectory analysis for diagnosing the interplay between chemistry and dynamics is investigated by comparing ozone mixing ratios modelled using air-parcel pathways to values observed along flight tracks during ATTREX (Airborne Tropical TRopopause EXperiment). Trajectories are initiated at the locations of ozone observations and tracked backward in time to their sources at termini of backward trajectories. The modelled values of ozone utilize 3-dimensional analysis fields from WACCM (Whole Atmosphere Community Climate Model) (a chemical-climate model with dynamical fields nudged towards MERRA (Modern-Era Retrospective Analysis and Research Applications) reanalysis) and ERA-interim (product of ECMWF - the European Centre for Medium-Range Weather Forecasts) to determine source mixing ratios with chemical production and loss terms derived from the ozone chemistry used in WACCM. A statistical base of modelled ozone is constructed with 6 trajectory platforms (adiabatic, diabatic, and kinematic forced by ERA-interim and MERRA), two chemical models (WACCM chemistry and no chemistry), and 4 trajectory lengths (5, 10, 20, and 30 days). Linear regression is employed to separate systematic errors from random errors and to characterize the impact of source mixing ratios, path length, vertical motion, and chemistry on modelled ozone errors. Errors in the analysis ozone fields are large, if not dominant, contributors to model error. Random errors are particularly large for point-by-point comparisons, however averaging over 800 km (75 minutes) flight segments substantially reduces random error and exposes systematic errors. Of the two analysis ozone data sets, WACCM, which incorporates detailed chemistry, provides the smaller systematic errors while ERA-interim, which has crude chemistry but assimilates observational data, has the smaller random errors. Of the different trajectory platforms, adiabatic calculations produce the smaller random errors (irrespective of the use of chemistry) but both vertical motion and chemistry are required to optimally reduce systematic errors. These results suggest that meaningful analysis of dynamical and chemical interactions that control ozone mixing ratios are viable on spatial scales larger than a few reanalysis grid spaces, that errors in the analyzed ozone data sets are large but not prohibitively so, and that vertical velocities and heating rates from reanalysis data, while problematic, contain useful information [on the ozone concentrations in the UTLS (Upper Troposphere/Lower Stratosphere)].

ozone↗
NASA NTRS2018

Understanding Climate-Vegetation Interactions in Global Rainforests Through a GP-Tree Analysis

The tropical rainforests are the largest reserves of terrestrial carbon and, therefore, the future of these rainforests is a question that is of immense importance in the geoscience research community. With the recent severe Amazonian droughts in 2005 and 2010 and on-going drought in the Congo region for more than two decades, there is growing concern that these forests could succumb to precipitation reduction, causing extensive carbon release and feedback to the carbon cycle. However, there is no single ecosystem model that quantifies the relationship between vegetation health in these rainforests and climatic factors. Small scale studies have used statistical correlation measure and simple linear regression to model climate-vegetation interactions, but suffer from the lack of comprehensive data representation as well as simplistic assumptions about dependency of the target on the covariates. In this paper we use genetic programming (GP) based symbolic regression for discovering equations that govern the vegetation climate dynamics in the rainforests. Expecting micro-regions within the rainforests to have unique characteristics compared to the overall general characteristics, we use a modified regression-tree based hierarchical partitioning of the space to build individual models for each partition. The discovery of these equations reveal very interesting characteristics about the Amazon and the Congo rainforests. Our method GP-tree shows that the rainforests exhibit tremendous resiliency in the face of extreme climatic events by adapting to changing conditions.

genetic programming↗