SEARCH · Engineering Papers

Results for “land surface temperature”

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.

SourceResearch connectionFrom yearThrough year

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

At least 991 records · Page 55

All sources NASA NTRS DOE OSTI

NASA NTRS2010

Does Aerosol Weaken or Strengthen the South Asian Monsoon?

Aerosols are known to have the ability to block off solar radiation reaching the earth surface, causing it to cool - the so-called solar dimming (SDM) effect. In the Asian monsoon region, the SDM effect by aerosol can produce differential cooling at the surface reducing the meridional thermal contrast between land and ocean, leading to a weakening of the monsoon. On the other hand, absorbing aerosols such as black carbon and dust, when forced up against the steep slopes of the southern Tibetan Plateau can produce upper tropospheric heating, and induce convection-dynamic feedback leading to an advance of the rainy season over northern India and an enhancement of the South Asian monsoon through the "Elevated Heat Pump" (EHP) effect. In this paper, we present modeling results showing that in a coupled ocean-atmosphere-land system in which concentrations of greenhouse gases are kept constant, the response of the South Asian monsoon to dust and black carbon forcing is the net result of the two opposing effects of SDM and EHP. For the South Asian monsoon, if the increasing upper tropospheric thermal contrast between the Tibetan Plateau and region to the south spurred by the EHP overwhelms the reduction in surface temperature contrast due to SDM, the monsoon strengthens. Otherwise, the monsoon weakens. Preliminary observations are consistent with the above findings. We find that the two effects are strongly scale dependent. On interannual and shorter time scales, the EHP effect appears to dominate in the early summer season (May-June). On decadal or longer time scales, the SDM dominates for the mature monsoon (July-August). Better understanding the physical mechanisms underlying the SDM and the EHP effects, the local emission and transport of aerosols from surrounding deserts and arid-regions, and their interaction with monsoon water cycle dynamics are important in providing better prediction and assessment of climate change impacts on precipitation of the Asian monsoon land regions.

Lau, William K.↗

NASA NTRS2015

A Prototype Physical Database for Passive Microwave Retrievals of Precipitation over the US Southern Great Plains

An accurate understanding of the instantaneous, dynamic land surface emissivity is necessary for a physically based, multi-channel passive microwave precipitation retrieval scheme over land. In an effort to assess the feasibility of the physical approach for land surfaces, a semi-empirical emissivity model is applied for calculation of the surface component in a test area of the US Southern Great Plains. A physical emissivity model, using land surface model data as input, is used to calculate emissivity at the 10GHz frequency, combining contributions from the underlying soil and vegetation layers, including the dielectric and roughness effects of each medium. An empirical technique is then applied, based upon a robust set of observed channel covariances, extending the emissivity calculations to all channels. For calculation of the hydrometeor contribution, reflectivity profiles from the Tropical Rainfall Measurement Mission Precipitation Radar (TRMM PR) are utilized along with coincident brightness temperatures (Tbs) from the TRMM Microwave Imager (TMI), and cloud-resolving model profiles. Ice profiles are modified to be consistent with the higher frequency microwave Tbs. Resulting modeled top of the atmosphere Tbs show correlations to observations of 0.9, biases of 1K or less, root-mean-square errors on the order of 5K, and improved agreement over the use of climatological emissivity values. The synthesis of these models and data sets leads to the creation of a simple prototype Tb database that includes both dynamic surface and atmospheric information physically consistent with the land surface model, emissivity model, and atmospheric information.

microwave↗

NASA NTRS2021

Fairfax County Urban Development: Identifying Urban Heat Mitigation Strategies for Climate Adaptation Planning in Fairfax County, Virginia

Extreme high temperatures lead to increased instances of cardiovascular disease, pulmonary disease, and even death, as well as increased energy consumption and infrastructure costs. People in urbanized areas experience higher temperatures than rural areas due to diminished vegetation and increased impervious surfaces which absorb and radiate heat. Fairfax County, Virginia has embarked on Resilient Fairfax, a program aimed at addressing climate adaptation and resilience. The DEVELOP team partnered with the Fairfax County Office of Environmental and Energy Coordination (OEEC) to assess the extent of the urban heat island effect on the county and its most vulnerable populations. The team used data from Landsat 8 Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS), as well as the ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) for the years 2013 to 2021 and found that the hottest spots were in densely urbanized areas, with temperatures as much as 47°F above that of undeveloped reference areas. The team used the Integrated Valuation of Ecosystem Services and Tradeoffs (InVEST) urban cooling model and determined that areas with higher tree canopy cover had greater heat mitigation capacity. Estimates from the InVEST model showed that a 4.5% increase in canopy cover across the county could result in a temperature reduction of up to 2.4°F in some areas. The results will allow partners to assess heat distribution across Fairfax County and implement effective mitigation strategies, including locating prime locations for cooling centers and increasing canopy cover.

W. Pierce Holloway↗

NASA NTRS1967

Second Application of the "Successive Transformations Method" to Predict the Safest Lunar Landing Site for an Astronaut: Report 3

The purpose of the first application of the "Successive Transformations Method" described in the second report was to determine, for a manned spacecraft, the optimum landing site on the moon. The results obtained, made by extrapolating the temperature TO furnished by Surveyor, can be summarized as follows: (1) The method is efficient in allowing selection of the optimum landing area. ( 2 ) The zrea selected is extensive and relatively flat; consequently, it is sufficient to allow trine for the difficult landing operation which results from the high speed of the spacecraft when reaching the moon. ( 3 ) The relative temperature obtained is lower than that where Surveyor landed. In this third report, a second application of the "Successive Transformations" is described to determine the safest conditions of a given landing site for the astronaut, especially in relation to determining the composition of the lunar surface material. As a means of accomplishing this task, it is the intent of this paper to analyze lunar craters since they are present over so much of the lunar surface.

Rojas, Hector R.↗

NASA NTRS1976

Viking 1: early results

A brief outline of the Viking 1 mission to Mars is followed by descriptions of the Martian landing site and the scientific instrumentation aboard Viking 1 orbiter and lander. Measurements of the Martian atmosphere provided data on its molecular composition, temperature and pressure. The detection of nitrogen in the Martian atmosphere indicates the existence of life. Panoramic photographs of the Martian surface were also obtained and are shown. Preliminary chemical and biological investigations on samples of Martian soil indicated the presence of the elements iron, calcium, silicon, titanium and aluminum as major constituents. Observed biochemical reactions were judged conducive of biological activity.

Source record↗

NASA NTRS1994

Two-Layer Variable Infiltration Capacity Land Surface Representation for General Circulation Models

A simple two-layer variable infiltration capacity (VIC-2L) land surface model suitable for incorporation in general circulation models (GCMs) is described. The model consists of a two-layer characterization of the soil within a GCM grid cell, and uses an aerodynamic representation of latent and sensible heat fluxes at the land surface. The effects of GCM spatial subgrid variability of soil moisture and a hydrologically realistic runoff mechanism are represented in the soil layers. The model was tested using long-term hydrologic and climatalogical data for Kings Creek, Kansas to estimate and validate the hydrological parameters. Surface flux data from three First International Satellite Land Surface Climatology Project Field Experiments (FIFE) intensive field compaigns in the summer and fall of 1987 in central Kansas, and from the Anglo-Brazilian Amazonian Climate Observation Study (ABRACOS) in Brazil were used to validate the mode-simulated surface energy fluxes and surface temperature.

GENERAL CIRCULATION MODELS↗

NASA NTRS2000

High Spatial Resolution Airborne Multispectral Thermal Infrared Remote Sensing Data for Analysis of Urban Landscape Characteristics

We have used airborne multispectral thermal infrared (TIR) remote sensing data collected at a high spatial resolution (i.e., 10m) over several cities in the United States to study thermal energy characteristics of the urban landscape. These TIR data provide a unique opportunity to quantify thermal responses from discrete surfaces typical of the urban landscape and to identify both the spatial arrangement and patterns of thermal processes across the city. The information obtained from these data is critical to understanding how urban surfaces drive or force development of the Urban Heat Island (UHI) effect, which exists as a dome of elevated air temperatures that presides over cities in contrast to surrounding non-urbanized areas. The UHI is most pronounced in the summertime where urban surfaces, such as rooftops and pavement, store solar radiation throughout the day, and release this stored energy slowly after sunset creating air temperatures over the city that are in excess of 2-4'C warmer in contrast with non-urban or rural air temperatures. The UHI can also exist as a daytime phenomenon with surface temperatures in downtown areas of cities exceeding 38'C. The implications of the UHI are significant, particularly as an additive source of thermal energy input that exacerbates the overall production of ground level ozone over cities. We have used the Airborne Thermal and Land Applications Sensor (ATLAS), flown onboard a Lear 23 jet aircraft from the NASA Stennis Space Center, to acquire high spatial resolution multispectral TIR data (i.e., 6 bandwidths between 8.2-12.2 (um) over Huntsville, Alabama, Atlanta, Georgia, Baton Rouge, Louisiana, Salt Lake City, Utah, and Sacramento, California. These TIR data have been used to produce maps and other products, showing the spatial distribution of heating and cooling patterns over these cities to better understand how the morphology of the urban landscape affects development of the UHI. In turn, these data have been used by government officials, urban planners, and other decision-makers, to make more informed decisions on how to mitigate the UHI and its subsequent impacts.

Quattrochi, Dale A.↗

NASA NTRS2008

Impact of Urban Growth on Surface Climate: A Case Study in Oran, Algeria

We develop a land use map discriminating urban surfaces from other cover types over a semiarid region in North Africa and use it in a land surface model to assess the impact of urbanized land on surface energy, water and carbon balances. Unlike in temperate climates where urbanization creates a marked heat island effect, this effect is not strongly marked in semiarid regions. During summer, the urban class results in an additional warming of 1.45 C during daytime and 0.81 C at night compared to that simulated for needleleaf trees under similar climate conditions. Seasonal temperatures show urban areas warmer than their surrounding during summer and slightly cooler in winter. The hydrological cycle is practically "shut down" during summer and characterized by relatively large amount of runoff in winter. We estimate the annual amount of carbon uptake to 1.94 million metric tons with only 11.9% assimilated during the rainy season. However, if urbanization expands to reach 50% of the total area excluding forests, the annual total carbon uptake will decline by 35% and the July mean temperature would increase only 0.10 C, compared to current situation. In contrast, if urbanization expands to 50% of the total land excluding forests and croplands but all short vegetation is replaced by native broadleaf deciduous trees, the annual carbon uptake would increase 39% and the July mean temperature would decrease by 0.9 C, compared to current configuration. These results provide guidelines for urban planners and land use managers and indicate possibilities for mitigating the urban heat.

Bounoua, Lahouari↗

NASA NTRS1980

Dryland pasture and crop conditions as seen by HCMM

Ground truth, aircraft, and satellite data were examined in order to: (1) assess the capability for determining wheat and pasture canopy temperatures in a dryland farming region from HCMM data; (2) assess the capability for determining soil moisture from HCMM data in dryland crops (winter wheat) from adjacent range lands; and (3) determine the relationship of HCMM-derived soil moisture and canopy temperature values with the condition of winter wheat and dryland farming areas during the principal growth stages. The IR data were screened to include areas having greater than 60% pasture and surface temperatures were recalculated using the atmospheric correction factor calculated by the modified RADTRA model, and the July 29, 1978 IR data were analyzed. Screening the IR data improved the relationship for July 24/July 13 and October 7/August 31 temperature/API relationship. However the coefficient of determination was not improved in the July 29/July 13 relationship.

Rosenthal, W. D.↗

NASA NTRS1997

Numerical Model Studies of the Martian Mesoscale Circulations

The study objectives were to evaluate by numerical modeling various possible mesoscale circulation on Mars and related atmospheric boundary layer processes. The study was in collaboration with J. Tillman of the University of Washington (who supported the study observationally). Interaction has been made with J. Prusa of Iowa State University in numerical modeling investigation of dynamical effects of topographically-influenced flow. Modeling simulations included evaluations of surface physical characteristics on: (i) the Martian atmospheric boundary layer and (ii) their impact on thermally and dynamically forced mesoscale flows. Special model evaluations were made in support of selection of the Pathfinder landing sites. J. Tillman's finding of VL-2 inter-annual temperature difference was followed by model simulations attempting to point out the forcing for this feature. Publication of the results in the reviewed literature in pending upon completion of the manuscripts in preparation as indicated later.

Segal, Moti↗

NASA NTRS2002

Urban Surface Radiative Energy Budgets Determined Using Aircraft Scanner Data

It is estimated that by the year 2025, 80% of the world's population will live in cities. The extent of these urban areas across the world can be seen in an image of city lights from the Defense Meteorological Satellite Program. In many areas of North America and Europe, it is difficult to separate individual cities because of the dramatic growth and sprawl of urbanized areas. This conversion of the natural landscape vegetation into man-made urban structures such as roads and buildings drastically alter the regional surface energy budgets, hydrology, precipitation patterns, and meteorology. One of the earliest recognized and measured phenomena of urbanization is the urban heat island (UHI) which was reported as early as 1833 for London and 1862 for Paris. The urban heat island results from the energy that is absorbed by man-made materials during the day and is released at night resulting in the heating of the air within the urban area. The magnitude of the air temperature difference between the urban and surrounding countryside is highly dependent on the structure of the urban area, amount of solar immolation received during the day, and atmospheric conditions during the night. These night time air temperature differences can be in the range of 2 to 5 C. or greater. Although day time air temperature differences between urban areas and the countryside exists during the day, atmospheric mixing and stability reduce the magnitude. This phenomena is not limited to large urban areas, but also occurs in smaller metropolitan areas. The UHI has significant impacts on the urban air quality, meteorology, energy use, and human health. The UPI can be mitigated through increasing the amount of vegetation and modification of urban surfaces using high albedo materials for roofs and paved surfaces. To understand why the urban heat island phenomenon exists it is useful to define the surface in terms of the surface energy budget. Surface temperature and albedo is a major component of the surface energy budget. Knowledge of it is important in any attempt to describe the radiative and mass fluxes which occur at the surface. Use of energy terms in modeling surface energy budgets allows the direct comparison of various land surfaces encountered in a urban landscape, from vegetated (forest and herbaceous) to non-vegetated (bare soil, roads, and buildings). These terms are also easily measured using remote sensing from aircraft or satellite platforms allowing one to examine the spacial variability. The partitioning of energy budget terms depends on the surface type. In natural landscapes, the partitioning is dependent on canopy biomass, leaf area index, aerodynamic roughness, and moisture status, all of which are influenced by the development stage of the ecosystem. In urban landscapes, coverage by man-made materials substantially alters the surface face energy budget. The remotely sensed data obtained from aircraft and satellites, when properly calibrated allows the measurement of important terms in the radiative surface energy budget a urban landscape scale.

Luvall, Jeffrey C.↗

NASA NTRS2003

Meridiani Planum Hematite Deposit: Potential for Preservation of Microfossils

Christensen et al., using data from the Mars Global Surveyor Thermal Emission Spectrometer (TES), have identified gray crystalline hematite in a 350 km by 750 km region near Meridiani Planum. The deposit corresponds closely to the low-albedo highlands unit sm, mapped as a wind-eroded, ancient, subaqueous sedimentary deposit. Christensen et al. interpreted the Meridiani Planum deposit to be an in-place, rock-stratigraphic sedimentary unit characterized by smooth, friable layers composed primarily of basaltic sediments with approximately 10 to 15% crystalline gray hematite. The Meridiani Planum hematite deposit has recently been designated as the prime landing site for one of the two Mars Exploration Rover (MER) spacecraft. The MER landings are scheduled for January, 2004. Christensen et al. discussed five possible mechanisms for the formation of this deposit: direct precipitation from standing, oxygenated, Fe-rich water; precipitation from Fe-rich hydrothermal fluids; low-temperature dissolution and precipitation through mobile groundwater leaching; surface weathering and coatings; thermal oxidation of magnetite-rich lavas. Four of these mechanisms involve the interactions of rock with water, and thus have implications in the search for evidence of microbial life.

C C Allen↗

NASA NTRS2008

GPM Pre-Launch Algorithm Development for Physically-Based Falling Snow Retrievals

In this work we compare and correlate the long time series (Nov.-March) n~easurements of precipitation rate from the Parsivels and 2DVD to the passive (89, 150, 183+/-1, +/-3, +/-7 GHz) observations of NOAA's AMSU-B radiometer. There are approximately 5-8 AMSU-B overpass views of the CARE site a day. We separate the comparisons into categories of no precipitation, liquid rain and falling snow precipitation. Scatterplots between the Parsivel snowfall rates and AMSU-B brightness temperatures (TBs) did not show an exploitable relationship for retrievals. We further compared and contrasted brightness temperatures to other surface measurements such as temperature and relative humidity with equally unsatisfying results. We found that there are similar TBs (especially at 89 and 150 GHz) for cases with falling snow and for non-precipitating cases. The comparisons indicate that surface emissivity contributions to the satellite observed TB over land can add uncertainty in detecting and estimating falling snow. The newest results show that the cloud icc scattering signal in the AMSU-B data call be detected by computing clear air TBs based on CARE radiosonde data and a rough estimate of surface emissivity. That is the differences in computed TI3 and AMSU-B TB for precipitating and nonprecipitating cases are unique such that the precipitating versus ~lon-precipitating cases can be identified. These results require that the radiosonde releases are within an hour of the AMSU-B data and allow for three surface types: no snow on the ground, less than 5 cm snow on the ground, and greater than 5 cm on the ground (as given by ground station data). Forest fraction and measured emissivities were combined to calculate the surface emissivities. The above work and future work to incorporate knowledge about falling snow retrievals into the framework of the expected GPM Bayesian retrievals will be described during this presentation.

Gail Skofronick-Jackson↗

NASA NTRS2021

Contribution of Meteorological Downscaling to Skill and Precision of Seasonal Drought Forecasts

Research in meteorological prediction on sub-seasonal to seasonal (S2S) timescales has seen growth in recent years. Concurrent with this, demand for seasonal drought forecasting has risen. While there is obvious synergy between these fields, S2S meteorological forecasting has typically focused on low resolution global models, while the development of drought can be sensitive to the local expression of weather anomalies and their interaction with local surface properties and processes. This suggests that downscaling might play an important role in the application of meteorological S2S forecasts to skillful forecasting of drought. Here, we apply the Generalized Analog Regression Downscaling (GARD) algorithm to downscale meteorological hindcasts from the NASA Goddard Earth Observing System (GEOS) global S2S forecast system. Downscaled meteorological fields are then applied to drive offline simulations with the Catchment Land Surface Model (CLSM) to forecast United States Drought Monitor (USDM) style drought indicators derived from simulated surface hydrology variables. We compare the representation of drought in these downscaled hindcasts to hindcasts that are not downscaled, using the North American Land Data Assimilation System Phase 2 (NLDAS-2) dataset as an observational reference. We find that downscaling using GARD improves hindcasts of temperature and temperature anomalies, but the results for precipitation are mixed and generally small. Overall, GARD downscaling led to improved hindcast skill for total drought across the Contiguous United States (CONUS), and improvements were greatest for extreme (D3) and exceptional (D4) drought categories.

Ryan A Zamora↗

NASA NTRSDate not supplied

Hydrological Data at the NASA GES DISC: Current Capabilities and New Opportunities

The NASA Goddard Earth Sciences Data and Information Services Center (GES DISC) is one of twelve NASA Earth science data centers that document, process, archive and distribute data from Earth observation missions and projects. GES DISC maintains an archive of several hydrology datasets, including the Land Data Assimilation Systems (LDAS) and the Gravity Recovery and Climate Experiment (GRACE) Data Assimilation for Drought Monitoring (GRACE-DA-DM) data products. These datasets include model output of heat fluxes, rain, snow, soil temperature, soil moisture, and runoff; and observational forcing data, including surface pressure, temperature, precipitation, downward shortwave and longwave radiation, humidity, and wind. The temporal resolution of the hydrology data at GES DISC ranges from hourly to monthly, and spatial resolutions range from 0.1° to 1.0°. The GES DISC provides services which enable users to aggregate, temporally and spatially subset, regrid, and visualize archived data including the GES DISC Subsetter, Hydrology Data Rods, and the Geospatial Interactive Online Visualization and Analysis Infrastructure (GIOVANNI). The Hydrology Data Rods service optimally reorganizes large hydrological data sets as extended time series, providing more efficient access for the hydrological community. The time series data (aka “data rods”) were integrated into hydrology community tools, such as the Data Rods Explorer on HydroShare. Furthermore, the GES DISC is in the process of migrating its data and services to the cloud. Hydrological data available at the GES DISC are now available in the Amazon Web Services (AWS) cloud (us-west-2 region) providing users Direct S3 data access and the capability for cloud computing operations. In this presentation, the hydrology data products and services currently available at the GES DISC will be summarized. Also discussed are the migration to the cloud, user support through this transition, and the status of migrating the data rods service to the cloud.

Ashley Heath↗

NASA NTRS2014

Introducing Multisensor Satellite Radiance-Based Evaluation for Regional Earth System Modeling

Earth System modeling has become more complex, and its evaluation using satellite data has also become more difficult due to model and data diversity. Therefore, the fundamental methodology of using satellite direct measurements with instrumental simulators should be addressed especially for modeling community members lacking a solid background of radiative transfer and scattering theory. This manuscript introduces principles of multisatellite, multisensor radiance-based evaluation methods for a fully coupled regional Earth System model: NASA-Unified Weather Research and Forecasting (NU-WRF) model. We use a NU-WRF case study simulation over West Africa as an example of evaluating aerosol-cloud-precipitation-land processes with various satellite observations. NU-WRF-simulated geophysical parameters are converted to the satellite-observable raw radiance and backscatter under nearly consistent physics assumptions via the multisensor satellite simulator, the Goddard Satellite Data Simulator Unit. We present varied examples of simple yet robust methods that characterize forecast errors and model physics biases through the spatial and statistical interpretation of various satellite raw signals: infrared brightness temperature (Tb) for surface skin temperature and cloud top temperature, microwave Tb for precipitation ice and surface flooding, and radar and lidar backscatter for aerosol-cloud profiling simultaneously. Because raw satellite signals integrate many sources of geophysical information, we demonstrate user-defined thresholds and a simple statistical process to facilitate evaluations, including the infrared-microwave-based cloud types and lidar/radar-based profile classifications.

Planetary Boundary Layer↗

NASA NTRS2011

Numerical Model Sensitivity to Heterogeneous Satellite Derived Vegetation Roughness

The sensitivity of a mesoscale weather prediction model to a 1 km satellite-based vegetation roughness initialization is investigated for a domain within the south central United States. Three different roughness databases are employed: i) a control or standard lookup table roughness that is a function only of land cover type, ii) a spatially heterogeneous roughness database, specific to the domain, that was previously derived using a physically based procedure and Moderate Resolution Imaging Spectroradiometer (MODIS) imagery, and iii) a MODIS climatologic roughness database that like (i) is a function only of land cover type, but possesses domain specific mean values from (ii). The model used is the Weather Research and Forecast Model (WRF) coupled to the Community Land Model within the Land Information System (LIS). For each simulation, a statistical comparison is made between modeled results and ground observations within a domain including Oklahoma, Eastern Arkansas, and Northwest Louisiana during a 4-day period within IHOP 2002. Sensitivity analysis compares the impact the three roughness initializations on time-series temperature, precipitation probability of detection (POD), average wind speed, boundary layer height, and turbulent kinetic energy (TKE). Overall, the results indicate that, for the current investigation, replacement of the standard look-up table values with the satellite-derived values statistically improves model performance for most observed variables. Such natural roughness heterogeneity enhances the surface wind speed, PBL height and TKE production up to 10 percent, with a lesser effect over grassland, and greater effect over mixed land cover domains.

Jasinski, Michael↗

NASA NTRS1990

Relationship of multispectral satellite data to land surface evaporation from the Australian continent

Based on NOAA-9 AVHRR and Nimbus-7 SMMR satellite data, satellite indices of vegetation from the Australian continent are calculated for the period of May 1986 to April 1987. Visible (VIS) and near infrared (NIR) reflectances and the normalized difference (ND) vegetation index are calculated from the AVHRR sensor. The microwave polarization difference (PD) is also calculated as the difference between the vertically and horizontally polarized brightness temperatures at 37 GHz. ND, PD, VIS, and NIR indices were plotted against rainfall and water balance estimates of evaporation. It is concluded that direct satellite monitoring of annual evaporation across the Australian continent using PD or VIS satellite indices of vegetation biomass appears possible for areas with evaporation less than 600 mm/y and that use of the ND relationship at continental scale may underpredict monthly evaporation of forests relative to agriculture.

Smith, R. C. G.↗