Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “net present value”

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 163 records · Page 9

Multimodel Estimates of Atmospheric Lifetimes of Long-lived Ozone-Depleting Substances: Present and Future

We have diagnosed the lifetimes of long-lived source gases emitted at the surface and removed in the stratosphere using six three-dimensional chemistry-climate models and a two-dimensional model. The models all used the same standard photochemical data. We investigate the effect of different definitions of lifetimes, including running the models with both mixing ratio (MBC) and flux (FBC) boundary conditions. Within the same model, the lifetimes diagnosed by different methods agree very well. Using FBCs versus MBCs leads to a different tracer burden as the implied lifetime contained in the MBC value does not necessarily match a model's own calculated lifetime. In general, there are much larger differences in the lifetimes calculated by different models, the main causes of which are variations in the modeled rates of ascent and horizontal mixing in the tropical midlower stratosphere. The model runs have been used to compute instantaneous and steady state lifetimes. For chlorofluorocarbons (CFCs) their atmospheric distribution was far from steady state in their growth phase through to the 1980s, and the diagnosed instantaneous lifetime is accordingly much longer. Following the cessation of emissions, the resulting decay of CFCs is much closer to steady state. For 2100 conditions the model circulation speeds generally increase, but a thicker ozone layer due to recovery and climate change reduces photolysis rates. These effects compensate so the net impact on modeled lifetimes is small. For future assessments of stratospheric ozone, use of FBCs would allow a consistent balance between rate of CFC removal and model circulation rate

CFC↗

Structure and Dynamics of Coronal Plasma

Brief summaries of the four published papers produced within the present performance period of NASA Grant NAGW-4081 are presented. The full text of the papers are appended to the report. The first paper titled "Coronal Structures Observed in X-rays and H-alpa Structures" was published in the Kofu Symposium proceedings. The study analyzes cool and hot behavior of two x-ray events, a small flare and a surge. It was found that a large H-alpha surge appears in x-rays as a very weak event, while a weak H-alpha feature corresponds to the brightest x-ray emission on the disk at the time of the observation. Calculations of the heating necessary to produce these signatures, and implications for the driving and heating mechanisms of flares vs. surges are presented. The second paper "Differential Magnetic Field Shear in an Active Region" has been published in The Astrophysical Journal. The study compared the three dimensional extrapolation of magnetic fields with the observed coronal structure in an active region. Based on the fit between observed coronal structure throughout the volume of the region and the calculated magnetic field configurations, the authors propose a differential magnetic field shear model for this active region. The decreasing field shear in the outer portions of the AR may indicate a continual relaxation of the magnetic field with time, corresponding to a net transport of helicity outward. The third paper "Difficulties in Observing Coronal Structure" has been published in the journal Solar Physics. This paper discusses the evidence that the temperature and density structure of the corona are far more complicated than had previously been thought. The discussion is based on five studies carried out by the group on coronal plasma properties, showing that any one x-ray instrument does see all of the plasma present in the corona, that hot and cool material may appear to be co-spatial at a given location in the corona, and that simple magnetic field extrapolations provide only a poor fit to the observed structure. The fourth paper "Analysis and Comparison of Loop Structures Imaged with NIXT and Yohkoh/SXT" has been published in Astronomy and Astrophysics. This paper analyzes and compares a variety of coronal loops, deriving loop pressure and emission measure from loop models. They are able to determine the volume filling factor in the corona, which is found to be in the range 0.001 to 0.01 for compact loops, and of order 1 for large structures. The small values suggest highly filamented structures, especially at lower temperatures.

Golub, Leon↗

Thermal balance of the atmospheres of Jupiter and Uranus

Two-dimensional, radiative-convective-dynamical models of the visible atmospheres of Jupiter and Uranus are presented. Zonally-averaged temperatures and heat fluxes are calculated numerically as functions of pressure and latitude. In addition to radiative heat fluxes, the dynamical heat flux due to large-scale baroclinic eddies is included and is parametrized using a mixing length theory which gives heat fluxes similar to those of Stone. The results for Jupiter indicate that the internal heat flow is non-uniform in latitude and nearly balances the net radiative flux leaving the atmosphere. The thermal emission is found to be uniform in latitude in agreement with Pioneer and Voyager observations. Baroclinic eddies are calculated to transport only a small amount of the meridional heat flow necessary to account for the uniformity of thermal emission with latitude. The bulk of the meridional heat transfer is found to occur very deep in the stable interior of Jupiter as originally proposed by Ingersoll and Porco. The relative importance of baroclinic eddies vs. internal heat flow in the thermal balance of Uranus depends on the ratio of emitted thermal power to absorbed solar power. The thermal balance of Uranus is compared to that of Jupiter for different values of this ratio.

Friedson, A. J.↗

Directly Measured Heating Rates of a Tropical Subvisible Cirrus Cloud

We present the first direct measurements of the infrared and solar heating rates of a tropical subvisible cirrus (SVC) cloud sampled off the east coast of Nicaragua on 25 July 2007 by the NASA ER-2 aircraft during the Tropical Composition, Cloud and Climate Coupling Experiment (TC4). On this day a persistent thin cirrus layer, with mostly clear skies underneath, was detected in real time by the cloud lidar on the ER-2, and the aircraft was directed to profile down through the SVC. Measurements of the net broadband infrared irradiance and spectrally integrated solar irradiance above, below, and through the SVC are used to determine the infrared and solar heating rates of the cloud. The lidar measurements show that the variable SVC layer was located between approximately 13 and 15 km. Its midvisible optical depth varied from 0.01 to 0.10 with a mean of 0.034 +/- 0.033. Its depolarization ratio was approximately 0.4, indicative of ice clouds. From the divergence of the measured net irradiances the infrared heating rate of the SVC was determined to be approximately 2.50 - 3.24 K/d and the solar heating rate was found to be negligible. These values are consistent with previous indirect observations of other SVC and with model-generated heating rates of SVC with similar optical depths. This study illustrates the utility and potential of the profiling sampling strategy employed here. A more fully instrumented high-altitude aircraft that also included in situ cloud and aerosol probes would provide a comprehensive data set for characterizing both the radiative and microphysical properties of these ubiquitous tropical clouds

Bucholtz, Anthongy↗

Effects of CO2 concentration and light intensity on photosynthesis of a rootless submerged plant, Ceratophyllum demersum L., used for aquatic food production in bioregenerative life support systems

In addition to green microalgae, aquatic higher plants are likely to play an important role in aquatic food production modules in bioregenerative systems for producing feed for fish, converting CO2 to O2 and remedying water quality. In the present study, the effects of culture conditions on the net photosynthetic rate of a rootless submerged plant, Ceratophyllum demersum L., was investigated to determine the optimum culture conditions for maximal function of plants in food production modules including both aquatic plant culture and fish culture systems. The net photosynthetic rate in plants was determined by the increase in dissolved O2 concentrations in a closed vessel containing a plantlet and water. The water in the vessel was aerated sufficiently with a gas containing a known concentration of CO2 gas mixed with N2 gas before closing the vessel. The CO2 concentrations in the aerating gas ranged from 0.3 to 10 mmol mol-1. Photosynthetic photon flux density (PPFD) in the vessel ranged from 0 (dark) to 1.0 mmol m-2 s-1, which was controlled with a metal halide lamp. Temperature was kept at 28 degrees C. The net photosynthetic rate increased with increasing PPFD levels and was saturated at 0.2 and 0.5 mmol m-2 s-1 PPFD under CO2 levels of 1.0 and 3.0 mmol mol-1, respectively. The net photosynthetic rate increased with increasing CO2 levels from 0.3 to 3.0 mmol mol-1 showing the maximum value, 75 nmol O2 gDW-1 s-1, at 2-3 mmol mol-1 CO2 and gradually decreased with increasing CO2 levels from 3.0 to 10 mmol mol-1. The results demonstrate that C. demersum could be an efficient CO2 to O2 converter under a 2.0 mmol mol-1 CO2 level and relatively low PPFD levels in aquatic food production modules. c2003 COSPAR. Published by Elsevier Science Ltd. All rights reserved.

Ecological Systems, Closed↗

On the Global Oxygen Anomaly and Air-Sea Flux

A new climatology of monthly air-sea oxygen fluxes throughout the ice-free surface global ocean is presented. The climatology is based on weighted linear least squares regressions using heat flux monthly anomalies for spatial and temporal interpolation of historical O2 data. The seasonal oceanic variations show that the tropical belt (20 S - 20 N) is characterized by relatively small air-sea fluxes when compared to the middle to high latitudes (40 deg - 70 deg). The largest and lowest seasonal fluxes occur during summer and winter in both hemispheres. By means of an atmospheric transport model we show that our climatology is in better agreement with the observed amplitude and phasing of the variations in atmospheric O2/N2 ratios because of seasonal air-sea exchanges at baseline stations in the Pacific Ocean than with previous air-sea O2 climatologies. Our study indicates that the component of the air-sea O2 flux that correlates with heat flux dominates the large-scale air-sea O2 exchange on seasonal timescales. The contribution of each major oceanic basin to the atmospheric observations is described. The seasonal net thermal (SNO(sub T)) and biological (SNO(sub B)) outgassing components of the flux are examined in relation to latitudinal bands, basin-wide, and hemispheric contributions. The Southern Hemisphere's SNO(sub B) (approximately 0.26 Pmol) and SNO(sub T) (approximately 0.29 Pmol) values are larger than the Northern Hemisphere's SNO(sub B) (approximately 0.15 Pmol) and SNO(sub T) (approximately 0.16 Pmol) values (1 Pmol = 10(exp 15) mol). We estimate a global extratropical carbon new production during the outgassing season of 3.7 Pg C (1 Pg = 10(exp 15) g), lower than previous estimates with air-sea O2 climatologies.

Garcia, Hernan E.↗

Free convection in the Matian atmosphere

The 'free convective' regime for the Martian atmospheric boundary layer (ABL) was investigated. This state occurs when the mean windspeed at the top of the ABL drops below some critical value U(sub c) and positive buoyant forces are present. Such forces can arise either from vertical temperature or water vapor gradients across the atmospheric surface layer. During free convection, buoyant forces drive narrow plumes that ascend to the inversion height with a return circulation consisting of broad slower-moving downdraughts. Horizontal pressure, temperature, windspeed, and water vapor fluctuations resulting form this circulation pattern can be quite large adjacent to the ground (within the surface layer). The local turbulent fluctuations cause non-zero mean surface stresses, sensible heat fluxes, and latent heat fluxes, even when the mean regional windspeed is zero. Although motions above the surface layer are insensitive to the nature of the surface, the sensible and latent heat fluxes are primarily controlled by processes within the interfacial sublayer immediately adjacent to the ground during free convection. Thus the distinction between aerodynamically smooth and rough airflow within the interfacial sublayer is more important than for the more typical situation where the mean regional windspeed is greater than U(sub c). Buoyant forces associated with water vapor gradients are particularly large on Mars at low pressures and high temperatures when the surface relative humidity is 100 percent, enhancing the likelihood of free convection under these conditions. On this basis, Ingersol postulated the evaporative heat losses from an icy surface on Mars at 237 K and current pressures would exceed the available net radiative flux at the surface, thus prohibiting ice from melting at low atmospheric pressures. Schumann has developed equations describing the horizontal fluctuations and mean vertical gradients occurring during free convection. Schumann's model was generalized to include convection driven by water vapor gradients and to include the effects of circulation above both aerodynamically smooth and rough surfaces.

Clow, G. D.↗

The Atmospheric Impact of the 1991 Mount Pinatubo Eruption

The 1991 eruption of Pinatubo produced about 5 cubic kilometers of dacitic magma and may be the second largest volcanic eruption of the century. Eruption columns reached 40 kilometers in altitude and emplaced a giant umbrella cloud in the middle to lower stratosphere that injected about 17 megatons of SO2, slightly more than twice the largest yielded by the 1982 eruption of El Chichon, Mexico. The SO2 formed sulfate aerosols that produced the largest perturbation to the stratospheric aerosol layer since the eruption of Krakatau in 1883. The aerosol cloud spread rapidly around the Earth in about 3 weeks and attained global coverage by about 1 year after the eruption. Peak local midvisible optical depths of up to 0.4 were measured in late 1992, and globally averaged values were about 0.1 to 0.15 for 2 years. The large aerosol cloud caused dramatic decreases in the amount of net radiation reaching the Earth's surface, producing a climate forcing that was two times stronger than the aerosols of El Chichon. Effects on climate were an observed surface cooling in the northern hemisphere of up to 0.5 to 0.6 C, equivalent to a hemispheric-wide reduction in net radiation of 4 watts per square meter and a cooling of perhaps as large as -0.4 C over large parts of the earth in 1992-93. Climate models seem to have predicted the cooling with a reasonable degree of accuracy. The Pinatubo climate forcing was stronger than the opposite warming of either the El Nino event or anthropogenic greenhouse gases in the period 1991-93. As a result of the presence of the aerosol particles, midlatitude ozone concentrations reached their lowest levels on record during 1992-93, the southern hemisphere 'ozone hole' increased in 1992 to an unprecedented size and ozone depletion rates were observed to be faster than ever before recorded. The atmospheric impact of the Pinatubo eruption has been profound, and it has sparked a lively interest in the role that volcanic aerosols play in climate change. This event has shown that a powerful eruption providing a 15 to 20 megaton release of SO2 into the stratosphere can produce sufficient aerosols to offset the present global warming trends and severely impact the ozone budget.

Self, Stephen↗

Modeling Variability of Plasma Conditions in the Io Torus

Telescopic observations an in situ measurements of the Io plasma torus show the density, temperature and composition to vary over time, sometimes up to a factor of 2. While previous models of the physical and chemical processes in the Io plasma torus have reasonably modeled the conditions of the Voyager 1 era, their authors have not addressed the observed variability nor explored the sensitivity of torus conditions to input parameters. In this paper we present a homogeneous torus model parameterized by five variables (transport timescale, neutral source strength, ratio of oxygen sulfur to atoms in the source, fraction of superthermal electrons, temperature of these hot electrons). The model incorporates the most recent data for ionization, recombination, charge exchange and radiative energy losses for the major torus species (S, S(sup +), S(sup ++), S(sup +++), O, O(sup +), O(sup ++). We solve equations of conservation of mass and energy to find equilibrium conditions for a set of input parameters. We compare model plasma conditions with those observed by Voyager 1 Voyager 2, and Cassini. Furthermore, we explore the sensitivity of torus conditions to each parameter. We find that (1) torus conditions are distinctly different for the Voyager 1, Voyager 2 and Cassini eras, (2) unique torus input parameters for any given era are poorly constrained given the wide range of solution space that is consistent with the range of observed torus conditions, (3) ion composition is highly sensitive to the specification of a non-thermal electron distribution, (4) neutral O/S source ratio is highly variable with model values ranging between 1.7 for Cassini to 4.0 for Voyager 1 conditions, (5) transport times range between 23 days for Voyager 2 to 50 days for Voyager 1 and Cassini, (6) neutral source strengths range between 7 to 30 x 10(sup -4) cm (sup -3) s(sup -1) which corresponds to a net production of 0.4 to 1.3 tons/s for a torus volume of 1.4 x 10(sup 31) cm(sup 3), or 38 R(sub j)(sup 3).

Delamere, P. A.↗

Measuring and modeling near-surface reflected and emitted radiation fluxes at the FIFE site

Information is presented pertaining to the measurement and estimation of reflected and emitted components of the radiation balance. Information is included about reflectance and transmittance of solar radiation from and through the leaves of some grass and forb prairie species, bidirectional reflectance from a prairie canopy is discussed and measured and estimated fluxes are described of incoming and outgoing longwave and shortwave radiation. Results of the study showed only very small differences in reflectances and transmittances for the adaxial and abaxial surfaces of grass species in the visible and infrared wavebands, but some differences in the infrared wavebands were noted for the forbs. Reflectance from the prairie canopy changed as a function of solar and view zenith angles in the solar principal plane with definite asymmetry about nadir. The surface temperature of prairie canopies was found to vary by as much as 5 C depending on view zenith and azimuth position and on the solar azimuth. Aerodynamic temperature calculated from measured sensible heat fluxes ranged from 0 to 3 C higher than nadir-viewed temperatures. Models were developed to estimate incoming and reflected shortwave radiation from data collected with a Barnes Modular Multiband Radiometer. Several algorithms for estimating incoming longwave radiation were evaluated and compared to actual measures of that parameter. Net radiation was calculated using the estimated components of the shortwave radiation streams, determined from the algorithms developed, and from the longwave radiation streams provided by the Brunt, modified Deacon, and the Stefan-Boltzmann models. Estimates of net radiation were compared to measured values and found to be within the measurement error of the net radiometers used in the study.

Blad, Blaine L.↗

Towards a Marine Stratus Climatology on Drizzle Occurrence from CALIPSO

Marine stratus are a predominant feature of our planet with the annual mean coverage exceeding 20%. They strongly reflect sunlight, yet exert only a modest effect on outgoing infrared radiation, providing a significant net cooling to the Earth’s radiative balance. Their formation is coupled to boundary layer circulations that are driven, in part, by cloud top radiative cooling and evaporative cooling from precipitation in downdrafts. Understanding how these cloud systems evolve as the climate changes is a key question that requires additional information on their lifecycle and microphysical properties to accurately represent their behavior in global circulation models. From a large-scale perspective, insight into the microphysical properties of marine stratus at cloud top can be realized through estimates of the effective radius (Re) of the droplet size distributions derived from MODIS observations. Estimates on the occurrence of rain/drizzle are available from CloudSat. Together these observations indicate that precipitation frequently occurs in clouds with higher cloud top Re. This relationship is consistent with the well documented shift in cloud top droplet size distributions towards fewer, yet larger droplets prior the onset of precipitation. Here we report on a new and complementary set observations from the CALIPSO mission. The approach derives an extinction-to-backscatter ratio (Sc, also known as the cloud lidar ratio) using an established relationship that depends on observations of the lidar attenuated backscatter and volume depolarization ratio within the cloud. Because Sc is strongly and inversely related to Re, a change in the derived Sc from higher to lower values corresponds to a change in the droplet size distribution as seen by MODIS. This change in the lidar signals at cloud top clearly identifies clouds that are capable of precipitation. The presentation provides a brief overview of the approach for deriving Sc and compares CALIOP-derived Sc with observations from other techniques. CALIOP classifications of drizzling clouds, based on the retrieved values Sc, are compared to independent, collocated assessments of drizzle occurrence reported in the standard CloudSat data products. Regional and seasonal comparisons highlight the strengths and weaknesses of the two sensors. A machine learning approach that combines information from both CALIOP and CloudSat showcases possible improvements in the global identification of scenes likely to contain rain-bearing clouds.

CALIPSO↗

Airborne Sunphotometer Measurements of Aerosol Optical Depth and Columnar Water Vapor During the Puerto Rico Dust Experiment, and Comparison with Land, Aircraft, and Satellite Measurements

Analyses of aerosol optical depth (AOD) and columnar water vapor (CWV) measurements obtained with the six-channel NASA Ames Airborne Tracking Sunphotometer (AATS-6) mounted on a twin-engine aircraft during the summer 2000 Puerto Rico Dust Experiment are presented. In general, aerosol extinction values calculated from AATS-6 AOD measurements acquired during aircraft profiles up to 5 km ASL reproduce the vertical structure measured by coincident aircraft in-situ measurements of total aerosol number and surface area concentration. Calculations show that the spectral dependence of AOD was small (mean Angstrom wavelength exponents of approximately 0.20) within three atmospheric layers defined as the total column beneath the top of each aircraft profile, the region beneath the trade wind inversion, and the region within the Saharan Air Layer (SAL) above the trade inversion. This spectral behavior is consistent with attenuation of incoming solar radiation by large dust particles or by dust plus sea salt. Values of CWV calculated from profile measurements by AATS-6 at 941.9 nm and from aircraft in-situ measurements by a chilled mirror dewpoint hygrometer agree to within approximately 4% (0.13 g/sq cm). AATS-6 AOD values measured on the ground at Roosevelt Roads Naval Air Station and during low altitude aircraft runs over the adjacent Cabras Island aerosol/radiation ground site agree to within 0.004 to 0.030 with coincident data obtained with an AERONET Sun/sky Cimel radiometer located at Cabras Island. For the same observation times, AERONET retrievals of CWV exceed AATS-6 values by a mean of 0.74 g/sq cm (approximately 21 %) for the 2.9-3.9 g/sq cm measured by AATS-6. Comparison of AATS-6 aerosol extinction values obtained during four aircraft ascents over Cabras Island with corresponding values calculated from coincident aerosol backscatter measurements by a ground-based micro-pulse lidar (MPL-Net) located at Cabras yields a similar vertical structure above the trade inversion. Finally, AATS-6 AOD values measured during low altitude aircraft traverses over the ocean are compared with corresponding AOD values retrieved over water from upwelling radiance measurements by the MODIS, TOMS, and GOES-8 Imager satellite sensors, with mixed results. These exercises highlight the need for continued satellite sensor comparison/validation studies to improve satellite AOD retrieval algorithms, and the usefulness of airborne sunphotometer measurements in the validation process.

Livingston, John M.↗

Validation Assessment for the Soil Moisture Active Passive (SMAP) Level 4 Carbon (L4_C) Data Product Version 5

The post-launch Cal/Val phase of the SMAP mission is guided by two primary objectives for each science product team: 1) to calibrate, verify, and improve the performance of the science algorithms, and 2) validate accuracies of the science data products as specified in the SMAP Level-1 mission science requirements. Algorithm science and product maintenance activities during the SMAP extended mission phase have also involved periodic algorithm calibration and product refinements to maintain or enhance product consistency and performance as well as science utility. This report provides an assessment of the latest (Version 5) SMAP Level 4 Carbon (L4_C) product. The L4_C Version 5 (v5) global record now spans more than six years (March 2015 – present) of SMAP operations and has benefited from five major reprocessing updates to the operational product. These reprocessing events and L4_C product release updates have incorporated various algorithm refinements and calibration adjustments to account for similar refinements to the upstream GEOS land model assimilation system, SMAP brightness temperatures, and MODIS vegetation inputs used for L4_C processing. The SMAP L4_C algorithms utilize a terrestrial carbon flux model informed by daily surface and root zone soil moisture information contributed from the SMAP Level 4 Soil Moisture (L4_SM) product along with optical remote sensing-based (e.g. MODIS-based) land cover and canopy fractional photosynthetic active radiation (fPAR), and other ancillary biophysical data. The carbon flux model estimates global daily net ecosystem CO2 exchange (NEE) and the component carbon fluxes, namely, vegetation gross primary production (GPP) and soil heterotrophic respiration (Rh). Other L4_C product elements include surface (~0-5 cm depth) soil organic carbon (SOC) stocks and associated environmental constraints to these processes, including soil moisture-related controls on GPP and ecosystem respiration (Kimball et al. 2014, Jones et al. 2017). The L4_C product addresses SMAP carbon cycle science objectives by: 1) providing a direct link between terrestrial carbon fluxes and underlying freeze/thaw and soil moisture-related constraints to these processes, 2) documenting primary connections between terrestrial water, energy and carbon cycles, and 3) improving understanding of terrestrial carbon sink activity. The SMAP L4_C algorithms and operational product are mature and at a CEOS Validation Stage 4 level (Jackson et al. 2012) based on extensive validation of the multi-year record against a diverse array of independent benchmarks, well characterized global performance, and systematic refinements gained from five major reprocessing events. There are no Level-1 mission science requirements for the L4_C product; however, self-imposed requirements have been established focusing on NEE as the primary product field for validation, and on demonstrating L4_C accuracy and success in meeting product science requirements (Jackson et al. 2012). The other L4_C product fields also have strong utility for carbon science applications (e.g., Liu et al. 2019, Endsley et al. 2020); however, analysis of these other fields is considered secondary relative to primary validation activities focusing on NEE. The L4_C targeted accuracy requirements are to meet or exceed a mean unbiased root-mean-square error (ubRMSE, or standard deviation of the error) for NEE of 1.6 g C m-2 d-1 and 30 g C m-2 yr-1, emphasizing northern (≥45°N) boreal and arctic ecosystems; this accuracy is similar to that of tower eddy covariance measurement-based observations (Baldocchi 2008). Methods used for the latest v5 L4_C product performance and validation assessment have been established from the SMAP Cal/Val plan and previous studies (Jackson et al. 2012, Jones et al. 2017) and include: 1) consistency evaluations of the product fields against earlier product releases (version 4 or earlier); 2) comparisons of daily carbon flux estimates with independent tower eddy 4covariance measurement-based daily carbon (CO2) flux observations from core tower validation sites (CVS); and 3) consistency checks against other global carbon products, including soil carbon inventory records, global GPP records derived from tower observation upscaling methods, and satellite-based observations of canopy solar induced chlorophyll fluorescence (SIF) as a surrogate for GPP. Metrics used to evaluate relative agreement between L4_C product fields and observational benchmarks include correlation (r-value), RMSE differences, bias and model sensitivity diagnostics. Following these validation criteria, the present report provides a validation assessment of the latest L4_C product release (v5). Detailed descriptions of the L4_C algorithm and additional global product accuracy and performance results are given elsewhere (Jones et al. 2017, Endsley et al. 2020). The v5 L4_C product replaces earlier product versions and continues to show: (i) accuracy and performance levels meeting or exceeding SMAP L4_C science requirements; (ii) improvement over the previous product version (version 4); and (iii) suitability for a diversity of science applications. Example L4_C applications from the recent literature include clarifying environmental trends and controls on the northern terrestrial carbon sink (Liu et al. 2019), diagnosing drought-related impacts on ecosystem productivity (Li et al. 2020), and regional monitoring of cropland conditions for projecting annual yields (Wurster et al. 2020). 2EXPECTED L4_C ALGORITHM AND PRODUCTPERFORMANCE The L4_C algorithm performance, including variance and uncertainty estimates of model outputs, was determined during the mission pre-launch phase through spatially explicit model sensitivity studies using available model inputs similar to those currently being used for operational production and evaluating the resulting model simulations over the observed range of northern (≥45 °N) and global conditions (Kimball et al. 2012, Entekhabi et al. 2014). The L4_C algorithm options were also evaluated during the mission prelaunch phase, including deriving canopy fPAR from lower order NDVI (Normalized Difference Vegetation Index) inputs in lieu of using MODIS (MOD15) fPAR; and including an explicit model representation of boreal fire disturbance recovery impacts. These results indicated that the L4_C accuracy requirements (i.e., NEE ubRMSE ≤ 30 g C m-2 yr-1or ≤ 1.6 g C m-2 d-1) could be met from the baseline algorithms over more than 82% and 89% of global and northern vegetated land areas, respectively (Yi et al. 2013, Kimball et al. 2014). The global L4_C algorithm error budget for NEE derived during the mission prelaunch phase indicated that the estimated NEE ubRMSE uncertainty is proportional to GPP and is therefore larger in higher biomass productivity areas, including forests and croplands (Kimball et al. 2014). Likewise, NEE ubRMSE uncertainty is expected to be lower in less-productive areas, including grasslands and shrublands. Expected model NEE ubRMSE levels were also generally within targeted accuracy levels for characteristically less-productive boreal and Arctic biomes, even though relative model error as a proportion of total productivity (NEE RMSE / GPP) may be large in these areas. The estimated NEE uncertainty was lower than expected in some warmer tropical high biomass productivity areas (e.g. Amazon rainforest) because of reduced low temperature and moisture constraints to the L4_C respiration calculations so that the bulk of model uncertainty is contributed by GPP in these areas. Model NEE uncertainty in the African Congo was estimated to

SMAP↗

GC13I-0857: Designing a Frost Forecasting Service for Small Scale Tea Farmers in East Africa

Kenya is the third largest tea exporter in the world, producing 10% of the world's black tea. Sixty percent of this production occurs largely by small scale tea holders, with an average farm size of 1.04 acres, and an annual net income of $1,075. According to a recent evaluation, a typical frost event in the tea growing region causes about $200 dollars in losses which can be catastrophic for a small holder farm. A 72-hour frost forecast would provide these small-scale tea farmers with enough notice to reduce losses by approximately 80 USD annually. With this knowledge, SERVIR, a joint NASA-USAID initiative that brings Earth observations for improved decision making in developing countries, sought to design a frost monitoring and forecasting service that would provide farmers with enough lead time to react to and protect against a forecasted frost occurrence on their farm. SERVIR Eastern and Southern Africa, through its implementing partner, the Regional Centre for Mapping of Resources for Development (RCMRD), designed a service that included multiple stakeholder engagement events whereby stakeholders from the tea industry value chain were invited to share their experiences so that the exact needs and flow of information could be identified. This unique event allowed enabled the design of a service that fit the specifications of the stakeholders. The monitoring service component uses the MODIS Land Surface Temperature product to identify frost occurrences in near-real time. The prediction component, currently under testing, uses the 2-m air temperature, relative humidity, and 10-m wind speed from a series of high-resolution Weather Research and Forecasting (WRF) numerical weather prediction model runs over eastern Kenya as inputs into a frost prediction algorithm. Accuracy and sensitivity of the algorithm is being assessed with observations collected from the farmers using a smart phone app developed specifically to report frost occurrences, and from data shared through our partner network developed at the stakeholder engagement meeting. This presentation will illustrate the efficacy of our frost forecasting algorithm, and a way forward for incorporating these forecasts in a meaningful way to the key decision makers - the small-scale farmers of East Africa.

frost↗

Cloud Scattering Impact on Thermal Radiative Transfer and Global Longwave Radiation

The potential importance of longwave (LW) cloud scattering has been recognized but the actual estimate of this effect on thermal radiation varies greatly among different studies. General circulation models (GCMs) generally neglect or simplify the multiple scattering in the LW. In this study, we use a rigorous radiative transfer algorithm to explicitly consider LW multiple-scattering and apply the GCM to quantify the impact of cloud LW scattering on thermal radiation fluxes. Our study shows that the cloud scattering effect on downward thermal radiation at the surface is concentrated in the infrared atmospheric window spectrum (800–1250 cm9exp −1)). The scattering effect on the outgoing longwave radiation (OLR) is also present in the window region over low clouds but it is mainly in the far-infrared spectrum (300–600 cm(exp −1)) over high clouds. For clouds with small to moderate optical depth (τ < 10), the scattering effect on thermal fluxes shows large variation with the cloud τ and has a maximum at an optical depth of ∼3. For opaque clouds, the scattering effect approaches an asymptote and is smaller and less important. The 2-stream radiative transfer scheme could have an error over 10% with an RMS error around 3.5%–4.0% in the calculated LW flux. This algorithm error of the 2-stream approximation could readily exceed the no-scattering error in the LW, and thus it is worthless to include the time-consuming computation of multiple scattering in a 2-stream radiative transfer scheme. However, the calculation error rapidly decreases as stream number increases and the RMS error in LW flux using the 4-stream scheme is under 0.3%, an accuracy sufficient for most climate studies. We implement the 4-stream discrete-ordinate algorithm in the GISS GCM and run the GCM for 20 years with and without the LW scattering effect, respectively. When cloud LW scattering is included, we find that the global annual mean OLR is reduced by 2.7 W/m(exp 2), and the downward surface flux and the net atmospheric absorption are increased by 1.6 W/m2 and 1.8 W/m(exp 2), respectively. Using one year of ISCCP clouds and running the standalone radiative transfer offline, the global annual mean non-scattering errors in OLR, surface LW downward flux and net atmospheric absorption are 3.6 W/m(exp 2), −1.1 W/m(exp 2), and −2.5 W/m(exp 2), respectively. The global scattering impact of 2.7 W/m(exp 2) on the OLR is small when compared to the typical global OLR value of 240 W/m2, but it is significant when compared to cloud LW radiative forcing (30 W/m2) and net cloud forcing (−14 W/m(exp 2)). Overall, the effect of neglecting scattering on the thermal fluxes is comparable to the reported clear sky radiative effect of doubling CO2.

longwave cloud scattering↗

Motion of Drops on Surfaces with Wettability Gradients

A liquid drop present on a solid surface can move because of a gradient in wettability along the surface, as manifested by a gradient in the contact angle. The contact angle at a given point on the contact line between a solid and a liquid in a gaseous medium is the angle between the tangent planes to the liquid and the solid surfaces at that point and is measured within the liquid side, by convention. The motion of the drop occurs in the direction of increasing wettability. The cause of the motion is the net force exerted on the drop by the solid surface because of the variation of the contact angle around the periphery. This force causes acceleration of an initially stationary drop, and leads to its motion in the direction of decreasing contact angle. The nature of the motion is determined by the balance between the motivating force and the resisting hydrodynamic force from the solid surface and the surrounding gaseous medium. A wettability gradient can be chemically induced as shown by Chaudhury and Whitesides who provided unambiguous experimental evidence that drops can move in such gradients. The phenomenon can be important in heat transfer applications in low gravity, such as when condensation occurs on a surface. Daniel et al have demonstrated that the velocity of a drop on a surface due to a wettability gradient in the presence of condensation can be more than two orders of magnitude larger than that observed in the absence of condensation. In the present research program, we have begun to study the motion of a drop in a wettability gradient systematically using a model system. Our initial efforts will be restricted to a system in which no condensation occurs. The experiments are performed as follows. First, a rectangular strip of approximate dimensions 10 x 20 mm is cut out of a silicon wafer. The strip is cleaned thoroughly and its surface is exposed to the vapor from an alkylchlorosilane for a period lasting between one and two minutes inside a desiccator. This is done using an approximate line source of the vapor in the form of a string soaked in the alkylchlorosilane. Ordinarily, many fluids, including water, wet the surface of silicon quite well. This means that the contact angle is small. But the silanized surface resists wetting, with contact angles that are as large as 100 degs. Therefore, a gradient of wettability is formed on the silicon surface. The region near the string is highly hydrophobic, and the contact angle decreases gradually toward a small value at the hydrophilic end away from this region. The change in wettability occurs over a distance of several mm. The strip is placed on a platform within a Plexiglas cell. Drops of a suitable liquid are introduced on top of the strip near the hydrophobic end. An optical system attached to a video camera is trained on the drop so that images of the moving drop can be captured on videotape for subsequent analysis. We have performed preliminary experiments with water as well as ethylene glycol drops. Results from these experiments will be presented in the poster. Future plans include the refinement of the experimental system so as to permit images to be recorded from the side as well as the top, and the conduct of a systematic study in which the drop size is varied over a good range. Experiments will be conducted with different fluids so as to obtain the largest possible range of suitably defined Reynolds and Capillary numbers. Also, an effort will be initiated on theoretical modeling of this motion. The challenges in the development of the theoretical description lie in the proper analysis of the region in the vicinity of the contact line, as well as in the free boundary nature of the problem. It is known that continuum models assuming the no slip condition all the way to the contact line fail by predicting that the stress on the solid surface becomes singular as the contact line is approached. One approach for dealing with this issue has been to relax the no-slip boundary condition using the Navier model. Molecular dynamics simulations of the contact line region show that for a non-polar liquid on a solid surface, the no-slip boundary condition is in fact incorrect near the contact line. Furthermore, the same simulations also show that the usual relationship between stress and the rate of deformation breaks down in the vicinity of the contact line. In developing continuum theoretical models of the system, we shall accommodate this knowledge to the extent possible.

Subramanian, R. Shankar↗

Model of Image Artifacts from Dust Particles

A mathematical model of image artifacts produced by dust particles on lenses has been derived. Machine-vision systems often have to work with camera lenses that become dusty during use. Dust particles on the front surface of a lens produce image artifacts that can potentially affect the performance of a machine-vision algorithm. The present model satisfies a need for a means of synthesizing dust image artifacts for testing machine-vision algorithms for robustness (or the lack thereof) in the presence of dust on lenses. A dust particle can absorb light or scatter light out of some pixels, thereby giving rise to a dark dust artifact. It can also scatter light into other pixels, thereby giving rise to a bright dust artifact. For the sake of simplicity, this model deals only with dark dust artifacts. The model effectively represents dark dust artifacts as an attenuation image consisting of an array of diffuse darkened spots centered at image locations corresponding to the locations of dust particles. The dust artifacts are computationally incorporated into a given test image by simply multiplying the brightness value of each pixel by a transmission factor that incorporates the factor of attenuation, by dust particles, of the light incident on that pixel. With respect to computation of the attenuation and transmission factors, the model is based on a first-order geometric (ray)-optics treatment of the shadows cast by dust particles on the image detector. In this model, the light collected by a pixel is deemed to be confined to a pair of cones defined by the location of the pixel s image in object space, the entrance pupil of the lens, and the location of the pixel in the image plane (see Figure 1). For simplicity, it is assumed that the size of a dust particle is somewhat less than the diameter, at the front surface of the lens, of any collection cone containing all or part of that dust particle. Under this assumption, the shape of any individual dust particle artifact is the shape (typically, circular) of the aperture, and the contribution of the particle to the attenuation factor for a given pixel is the fraction of the cross-sectional area of the collection cone occupied by the particle. Assuming that dust particles do not overlap, the net transmission factor for a given pixel is calculated as one minus the sum of attenuation factors contributed by all dust particles affecting that pixel. In a test, the model was used to synthesize attenuation images for random distributions of dust particles on the front surface of a lens at various relative aperture (F-number) settings. As shown in Figure 2, the attenuation images resembled dust artifacts in real test images recorded while the lens was aimed at a white target.

Willson, Reg↗

Using ISCCP Weather States to Decompose Cloud Radiative Effects

The presentation will examine the shortwave (SW) and longwave (LW) cloud radiative effect CRE (aka "cloud radiative forcing") at the top-of-the-atmosphere and surface of ISCCP weather states (aka "cloud regimes") in three distinct geographical zones, one tropical and two mid-latitude. Our goal is to understand and quantify the contribution of the different cloud regimes to the planetary radiation budget. In the tropics we find that the three most convectively active states are the ones with largest SW, LW and net TOA CRE contributions to the overall daytime tropical CRE budget. They account for 59%, 71% and 55% of the total CRE, respectively. The boundary layer-dominated weather states account for only 34% of the total SW CRE and 41% of the total net CRE, so to focus only on them in cloud feedback studies may be imprudent. We also find that in both the northern and southern midlatitude zones only two weather states, the first and third most convectively active with large amounts of nimbostratus-type clouds, contribute ",40% to both the SW and net TOA CRE budgets, highlighting the fact that cloud regimes associated with frontal systems are not only important for weather (precipitation) but also for climate (radiation budget). While all cloud regimes in all geographical zones have a slightly larger SFC than TOA SW CRE, implying cooling of the surface and slight warming of the atmosphere, their LW radiative effects are more subtle: in the tropics the weather states with plentiful high clouds warm the atmosphere while those with copious amounts of low clouds cool the atmosphere. In both midlatitude zones only the weather states with peak cloud fractions at levels above 440 mbar warm the atmosphere while all the rest cool it. These results make the connection of the contrasting CRE effects to the atmospheric dynamics more explicit - "storms" tend to warm the atmosphere whereas fair weather clouds cool it, suggesting a positive feedback of clouds on weather systems. The breakdown of CRE by cloud regime are however not entirely similar between the two midlatitude zones. Despite the existence of an additional state in the nort!lern midlatitudes, only four weather states have net daytime CREs with absolute values above 100 Watts per square meter compared to six in the south. This reminds us that the environment where clouds occur also has a crucial role in determining their radiative effects. All the above make evident that reproducing grand averages of current CRE by climate models in only part of the challenge. If existing cloud regimes and shifts in their distributions and frequency of occurrence in a changed climate are not properly simulated, the radiative role of clouds will not be adequately predicted.

Oreopoulos, L.↗