Search NASA⌕ Search

SEARCH · Search NASA

Results for “dust emission scheme”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Assessment of State-of-the-Art Dust Emission Scheme in GEOS

The GEOS modeling system has been extended with state of the art parameterization of dust emissions based on the vertical flux formulation described in Kok et al 2014. The new dust scheme was coupled with the GOCART and MAM aerosol models. In the present study we compare dust emissions, aerosol optical depth (AOD) and radiative fluxes from GEOS experiments with the standard and new dust emissions. AOD from the model experiments are also compared with AERONET and satellite based data. Based on this comparative analysis we concluded that the new parameterization improves the GEOS capability to model dust aerosols originating from African sources, however it lead to overestimation of dust emissions from Asian and Arabian sources. Further regional tuning of key parameters controlling the threshold friction velocity may be required in order to achieve more definitive and uniform improvement in the dust modeling skill.

Darmenov, Anton↗
NASA NTRS2022

Importance of Different Parameterization Changes for the Updated Dust Cycle Modeling in the Community Atmosphere Model (Version 6.1)

The Community Atmosphere Model (CAM6.1), the atmospheric component of the Community Earth System Model (CESM; version 2.1), simulates the life cycle (emission, transport, and deposition) of mineral dust and its interactions with physio-chemical components to quantify the impacts of dust on climate and the Earth system. The accuracy of such quantifications relies on how well dust-related processes are represented in the model. Here we update the parameterizations for the dust module, including those on the dust emission scheme, the aerosol dry deposition scheme, the size distribution of transported dust, and the treatment of dust particle shape. Multiple simulations were undertaken to evaluate the model performance against diverse observations, and to understand how each update alters the modeled dust cycle and the simulated dust direct radiative effect. The model–observation comparisons suggest that substantially improved model representations of the dust cycle are achieved primarily through the new more physically-based dust emission scheme. In comparison, the other modifications induced small changes to the modeled dust cycle and model–observation comparisons, except the size distribution of dust in the coarse mode, which can be even more influential than that of replacing the dust emission scheme. We highlight which changes introduced here are important for which regions, shedding light on further dust model developments required for more accurately estimating interactions between dust and climate.

mineral dust parameterizations↗
NASA NTRS2024

Evaluation of 10‐m Wind Speed From ISD Meteorological Stations and the MERRA‐2 Reanalysis: Impacts on Dust Emission in the Arabian Peninsula

Mineral dust is one of the most important aerosols when studying the radiative balance and climate of the planet. There are different dust emission schemes utilized by the atmospheric modeling communities, many of which disagree on basic output quantities such as mass of dust emitted and distribution of mass among size bins. In this work, we examined mineral dust emission from a leading model scheme, the Goddard Chemistry Aerosol Radiation and Transport (GOCART), as utilized in the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) Reanalysis and compared it to dust emissions calculated using wind measurements from ground based weather stations located in the Arabian Peninsula that are included in the National Oceanic and Atmospheric Administration’s (NOAA) integrated surface database (ISD). An intercomparison of 10-m wind speed is shown for the Arabian Peninsula region, differences of the observed and modeled wind field are quantified, and impacts of differences on dust emissions are calculated. This analysis shows 10-m winds in the ISD were generally lower than MERRA-2 winds, which propagated to dust emissions errors. Our estimate of one of the most significant mass impacts in dust emission is 0.178 Tg/year/grid box with a percent change of over 200% to the recalculated dust emissions from MERRA-2. These differences in wind speed propagated to a difference in dust mass emitted by the use of a static source function which aids in scaling the mass emitted by the availability of dust in each grid. Additionally, the magnitude of these differences varies seasonally.

winds↗
NASA NTRS2017

Assimilation of MODIS Dark Target and Deep Blue Observations in the Dust Aerosol Component of NMMB-MONARCH version 1.0

A data assimilation capability has been built for the NMMB-MONARCH chemical weather prediction system, with a focus on mineral dust, a prominent type of aerosol. An ensemble-based Kalman filter technique (namely the local ensemble transform Kalman filter - LETKF) has been utilized to optimally combine model background and satellite retrievals. Our implementation of the ensemble is based on known uncertainties in the physical parametrizations of the dust emission scheme. Experiments showed that MODIS AOD retrievals using the Dark Target algorithm can help NMMB-MONARCH to better characterize atmospheric dust. This is particularly true for the analysis of the dust outflow in the Sahel region and over the African Atlantic coast. The assimilation of MODIS AOD retrievals based on the Deep Blue algorithm has a further positive impact in the analysis downwind from the strongest dust sources of the Sahara and in the Arabian Peninsula. An analysis-initialized forecast performs better (lower forecast error and higher correlation with observations) than a standard forecast, with the exception of underestimating dust in the long-range Atlantic transport and degradation of the temporal evolution of dust in some regions after day 1. Particularly relevant is the improved forecast over the Sahara throughout the forecast range thanks to the assimilation of Deep Blue retrievals over areas not easily covered by other observational datasets.The present study on mineral dust is a first step towards data assimilation with a complete aerosol prediction system that includes multiple aerosol species.

atmospheric dust↗
NASA NTRS2013

A 4-D Climatology (1979-2009) of the Monthly Tropospheric Aerosol Optical Depth Distribution over the Mediterranean Region from a Comparative Evaluation and Blending of Remote Sensing and Model Products

Since the 1980s several spaceborne sensors have been used to retrieve the aerosol optical depth (AOD) over the Mediterranean region. In parallel, AOD climatologies coming from different numerical model simulations are now also available, permitting to distinguish the contribution of several aerosol types to the total AOD. In this work, we perform a comparative analysis of this unique multiyear database in terms of total AOD and of its apportionment by the five main aerosol types (soil dust, seasalt, sulfate, black and organic carbon). We use 9 different satellite-derived monthly AOD products: NOAA/AVHRR, SeaWiFS (2 products), TERRA/MISR, TERRA/MODIS, AQUA/MODIS, ENVISAT/MERIS, PARASOL/POLDER and MSG/SEVIRI, as well as 3 more historical datasets: NIMBUS7/CZCS, TOMS (onboard NIMBUS7 and Earth- Probe) and METEOSAT/MVIRI. Monthly model datasets include the aerosol climatology from Tegen et al. (1997), the climate-chemistry models LMDz-OR-INCA and RegCM-4, the multi-model mean coming from the ACCMIP exercise, and the reanalyses GEMS and MACC. Ground-based Level- 2 AERONET AOD observations from 47 stations around the basin are used here to evaluate the model and satellite data. The sensor MODIS (on AQUA and TERRA) has the best average AOD scores over this region, showing a relevant spatio-temporal variability and highlighting high dust loads over Northern Africa and the sea (spring and summer), and sulfate aerosols over continental Europe (summer). The comparison also shows limitations of certain datasets (especially MERIS and SeaWiFS standard products). Models reproduce the main patterns of the AOD variability over the basin. The MACC reanalysis is the closest to AERONET data, but appears to underestimate dust over Northern Africa, where RegCM-4 is found closer to MODIS thanks to its interactive scheme for dust emissions. The vertical dimension is also investigated using the CALIOP instrument. This study confirms differences of vertical distribution between dust aerosols showing a large vertical spread, and other continental and marine aerosols which are confined in the boundary layer. From this compilation, we propose a 4-D blended product from model and satellite data, consisting in monthly time series of 3-D aerosol distribution at a 50 km horizontal resolution over the Euro-Mediterranean marine and continental region for the 2003-2009 period. The product is based on the total AOD from AQUA/MODIS, apportioned into sulfates, black and organic carbon from the MACC reanalysis, and into dust and sea-salt aerosols from RegCM-4 simulations, which are distributed vertically based on CALIOP climatology.We extend the 2003-2009 reconstruction to the past up to 1979 using the 2003-2009 average and applying the decreasing trend in sulfate aerosols from LMDz-OR-INCA, whose AOD trends over Europe and the Mediterranean are median among the ACCMIP models. Finally optical properties of the different aerosol types in this region are proposed from Mie calculations so that this reconstruction can be included in regional climate models for aerosol radiative forcing and aerosol-climate studies.

troposphere↗
NASA NTRS1986

A cloud collision model for water maser excitation

High-velocity collisions between small, dense, neutral clouds or between a dense cloud and a dense shell can provide the energy source required to excite H2O maser emission. The radiative precursor from the surface of the collisional shock front rapidly diffuses through the cloud, heating the dust grains but leaving the H2 molecules cool. Transient maser emission occurs as the conditions for the Goldreich and Kwan (1974) 'hot-dust cold-gas' maser pump scheme are realized locally within the cloud. In time, the local maser action quenches due to the heating of the H2 molecules by collisions against the grains. Although this model cannot explain the very long-lived steady maser features, it is quite successful in explaining a number of the observed properties of the high-velocity features in such sources as Orion, W51, and W49. In particular, it provides a natural explanation for the rapid time variations, the narrow line widths, juxtaposition of high- and low-velocity features, and the short lifetimes which are frequently observed for the so-called high-velocity maser 'bullets' thought to be accelerated by strong stellar winds.

Tarter, T. C.↗
NASA NTRS2003

The ISO View of Palomar-Green Quasars

Mining the ISO data archive we provide the complete ISO view of PG quasars containing 64 infrared spectral energy distributions between 5 and 200 mu m. About half of the sample was supplemented by MAMBO and SCUBA (sub-)millimeter data. Since the PG quasars were selected optically, the high infrared detection rate of more than 80% suggests that every quasar possesses luminous to hyper-luminous dust emission with dust masses comparable to Seyferts and ultra-luminous IR galaxies (ULIRGs). The gas to-dust mass ratio (of those sources where CO measurements are available in the literature) is consistent with the galactic value providing further evidence for the thermal nature of the IR emission of radio quiet quasars. The SEDs represent templates of unprecedented detail and sensitivity. We suggest that the diversity of the SEDs reflects largely the evolution of the dust distribution, and we propose a classification of the SED shapes as well as an evolutionary scheme in which this variety can be understood. During the evolution the surrounding dust redistributes, settling more and more into a torus/disk like configuration, while the SEDs show an initial FIR bump, then an increasing MIR emission and a steeper near- to mid-infrared slope, both of which finally also decrease. Regarding cosmic evolution, our hyper-luminous quasars in the "local" universe at z=l do not show the hyper-luminous (LFIR >? 10(exp 13) L(sub sun)) starburst activity inferred for z=4 quasars detected in several (sub-)millimeter surveys. In view of several caveats this difference should be established further, but it already suggests that in the early dense universe stronger merger events led to more powerful starbursts accompanying the quasar phenomenon, while at later cosmic epochs any coeval starbursts obviously do not reach that high power and are outshone by the AGN. Additional information is included in the original extended abstract.

Haas, M.↗
NASA NTRS1995

The 11 Micron Emissions of Carbon Stars

A new classification scheme of the IRAS LRS carbon stars is presented. It comprises the separation of 718 probable carbon stars into 12 distinct self-similar spectral groupings. Continuum temperatures are assigned and range from 470 to 5000 K. Three distinct dust species are identifiable: SiC, alpha:C-H, and MgS. In addition to the narrow 11 + micron emission feature that is commonly attributed to SiC, a broad 11 + micron emission feature, that is correlated with the 8.5 and 7.7 micron features, is found and attributed to alpha:C-H. SiC and alpha:C-H band strengths are found to correlate with the temperature progression among the Classes. We find a spectral sequence of Classes that reflects the carbon star evolutionary sequence of spectral types, or alternatively developmental sequences of grain condensation in carbon-rich circumstellar shells. If decreasing temperature corresponds to increasing evolution, then decreasing temperature corresponds to increasing C/O resulting in increasing amounts of carbon rich dust, namely alpha:C-H. If decreasing the temperature corresponds to a grain condensation sequence, then heterogeneous, or induced nucleation scenarios are supported. SiC grains precede alpha:C-H and form the nuclei for the condensation of the latter material. At still lower temperatures, MgS appears to be quite prevalent. No 11.3 micron PAH features are identified in any of the 718 carbon stars. However, one of the coldest objects, IRAS 15048-5702, and a few others, displays an 11.9 micron emission feature characteristic of laboratory samples of coronene. That feature corresponds to the C-H out of plane deformation mode of aromatic hydrocarbon. This band indicates the presence of unsaturated, sp(sup 3), hydrocarbon bonds that may subsequently evolve into saturated bonds, sp(sup 2), if, and when, the star enters the planetary nebulae phase of stellar evolution. The effusion of hydrogen from the hydrocarbon grain results in the evolution in wavelength of this 11.9 micron emission feature to the 11.3 micron feature.

Goebel, J. H.↗
NASA NTRS1995

The 11 Micron Emissions of Cabon Stars

A new classification scheme of the IRAS LRS carbon stars is presented. It comprises the separation of 718 probable carbon stars into 12 distinct self-similar spectral groupings. Continuum temperatures are assigned and range from 470 to 5000 K. Three distinct dust species are identifiable: SiC, alpha:C-H, and MgS. In addition to the narrow 11 + micron emission feature that is commonly attributed to SiC, a broad 11 + micron emission feature, that is correlated with the 8.5 and 7.7 micron features, is found and attributed to alpha:C-H. SiC and alpha:C-H band strengths are found to correlate with the temperature progression among the Classes. We find a spectral sequence of Classes that reflects the carbon star evolutionary sequence of spectral types, or alternatively developmental sequences of grain condensation in carbon-rich circumstellar shells. If decreasing temperature corresponds to increasing evolution, then decreasing temperature corresponds to increasing CIO resulting in increasing amounts of carbon rich dust, namely alpha:C-H. If decreasing the temperature corresponds to a grain condensation sequence, then heterogeneous, or induced nucleation scenarios are supported. SiC grains precede alpha:C-H and form the nuclei for the condensation of the latter material. At still lower temperatures, MgS appears to be quite prevalent. No 11.3 micron PAH features are identified in any of the 718 carbon stars. However, one of the coldest objects, IRAS 15048-5702, and a few others, displays an 11.9 micron emission feature characteristic of laboratory samples of coronene. That feature corresponds to the C-H out of plane deformation mode of aromatic hydrocarbon. This band indicates the presence of unsaturated, sp(sup 3), hydrocarbon bonds that may subsequently evolve into saturated bonds, sp(sup 2), if, and when, the star enters the planetary nebulae phase of stellar evolution. The effusion of hydrogen from the hydrocarbon grain results in the evolution in wavelength of this 11.9 micron emission feature to the 11.3 micron feature.

Goebel, J. H.↗
NASA NTRS2013

Impact of Surface Roughness and Soil Texture on Mineral Dust Emission Fluxes Modeling

Dust production models (DPM) used to estimate vertical fluxes of mineral dust aerosols over arid regions need accurate data on soil and surface properties. The Laboratoire Inter-Universitaire des Systemes Atmospheriques (LISA) data set was developed for Northern Africa, the Middle East, and East Asia. This regional data set was built through dedicated field campaigns and include, among others, the aerodynamic roughness length, the smooth roughness length of the erodible fraction of the surface, and the dry (undisturbed) soil size distribution. Recently, satellite-derived roughness length and high-resolution soil texture data sets at the global scale have emerged and provide the opportunity for the use of advanced schemes in global models. This paper analyzes the behavior of the ERS satellite-derived global roughness length and the State Soil Geographic data base-Food and Agriculture Organization of the United Nations (STATSGO-FAO) soil texture data set (based on wet techniques) using an advanced DPM in comparison to the LISA data set over Northern Africa and the Middle East. We explore the sensitivity of the drag partition scheme (a critical component of the DPM) and of the dust vertical fluxes (intensity and spatial patterns) to the roughness length and soil texture data sets. We also compare the use of the drag partition scheme to a widely used preferential source approach in global models. Idealized experiments with prescribed wind speeds show that the ERS and STATSGO-FAO data sets provide realistic spatial patterns of dust emission and friction velocity thresholds in the region. Finally, we evaluate a dust transport model for the period of March to July 2011 with observed aerosol optical depths from Aerosol Robotic Network sites. Results show that ERS and STATSGO-FAO provide realistic simulations in the region.

textures↗
NASA NTRS1993

The NASA/GISS Mars general circulation model: Preliminary experiments

The NASA/GISS Mars General Circulation Model (GCM) is an adapted version of the GISS Global Climate/Middle Atmosphere Model, specifically developed for the diagnostic validation and objective analysis of measured atmospheric temperatures from the Mars Observer Pressure Modulator Infrared Radiometer (PMIRR) experiment. The GISS Mars GCM has 23 vertical layers extending from the surface to approximately 80 km altitude, representing a vertical resolution of about 0.3 scale heights. The primitive (vertically hydrostatic) equations are solved in finite difference form on the Krakawa B grid, with a horizontal resolution of 8 deg x 10 deg (latitude-longitude). The model includes a diurnal solar cycle, heat transport within a two-layer ground, and a high-order 'slopes-scheme' for the advection of heat in the upper atmosphere. The radiative transfer scheme is based on the correlated k distribution method for the treatment of nongray gaseous absorption thermal emission, and multiple scattering, including options for suspended dust. A special feature of the model of particular importance for Mars is a parameterization of gravity-wave-induced drag incorporating orographic forcing, wind shear, convection, and radiative damping. The implementation of the GISS Mars model includes global maps of topography, roughness, and albedo.

Allison, Michael↗
NASA NTRSDate not supplied

FUTURE LEVELS OF AEROSOLS & RADIATIVE FORCING OVER THE ARCTIC

The Arctic Monitoring and Assessment Programme (AMAP) is currently assessing the impacts of Short-Lived Climate Forcers (SLCF) on Arctic climate and air quality. In support of the assessment, we used the NASA Goddard Institute of Space Sciences (GISS) Earth System Model (modelE2.1), with prescribed sea surface temperature and sea-ice fraction, to simulate SLCF concentrations globally between 1995 and 2015. Two simulations were conducted, using the One-Moment Aerosol (OMA) and the Multiconfiguration Aerosol TRacker of mIXing state (MATRIX) aerosol modules. OMA is a mass-based scheme in which aerosols are assumed to remain externally mixed and have a prescribed and constant size distribution, while MATRIX is an aerosol microphysics scheme based on the quadrature method of moments, which is able to explicitly simulate the mixing state of aerosols. Anthropogenic emissions from the ECLIPSE v6b emissions database were used, along with emissions from aircrafts and open biomass burning from the Coupled Model Intercomparison Project Phase 6 (CMIP6), while the natural emissions of sea salt, DMS, isoprene and dust are calculated interactively. The simulated monthly surface concentrations of sulfate (SO4), black carbon (BC), organic carbon (OA), and ozone (O3) are compared with observations from a set of Arctic stations, extracted from the EBAS and IMPROVE databases, as well as a few additional locations. Simulated aerosol optical depths (AOD) are also compared with Advanced Very-High Resolution Radiometer (AVHRR). The study will present the evaluation of the modelE2.1 in simulating SLCF levels over the Arctic using different aerosol schemes, along with observed and simulated trends of SLCFs over the Arctic between 1995 and 2015.

Aersosols↗
NASA NTRS1998

Simulation of Tropical Biomass Burning

The work proposed was carried out as planned. The work described in this final report formed the basis for a follow-on research grant research grant from NASA Ames Research Center. The research objectives that were achieved during the course of our studies include the following: (1) the evaluation of several components of MM5 (Meteorological Model 5 version 2) and the Global/Regional Atmospheric Chemistry Event Simulator (GRACES) combined modeling system; (2) improved calculations of the transport of tracers for both NASA airborne missions, Study of Ozone and Nitrogen oxides experiment (SONEX) and Pacific Exploratory MIssion in the Tropics (PEM-Tropics); (3) improved source strength estimates for isoprene, dust and similar emissions from the Earth's surface. This required the use of newly available databases on the Earth's surface and vegetation; (4) completed atmospheric chemistry simulations of radicals and nitrogen oxide species; (5)improved the handling of cumulonimbus convection by modifying the existing scheme; (6) identified the role of the African Intertropical Front, using MM5's nesting capability to refine model resolution in crucial areas; modified the MM5 trajectory program to allow it to work much better for a parcel crossing the west/east boundaries.

Hamill, Patrick↗
NASA NTRS1996

Simulation of Tropical Biomass Burning

The research objectives that were achieved during the course of our studies include the following: (1) Over the last few years, a model has been developed in the Atmospheric Chemistry and Dynamics Branch at Ames Research Center in collaboration with the Physics Department at San Jose State University. It is referred to as the Global/Regional Atmospheric Chemistry Event Simulator (GRACES). Currently, the GRACES model system combines an atmospheric chemistry and transport model, and a regional mesoscale meteorological model. Therefore this system is suitable for simulating the conditions observed by the tropical observation missions, such as the Pacific Exploratory Mission in the 'Tropics (PEM-Tropics), Study of Ozone and Nitrogen oxides Experiment (SONEX), and other periods. Specifically, the research carried out included the evaluation of the behavior of several components of the MM5 (I.e., Meteorological Model 5, version 2) and the GRACES combined modeling system. We initiated research on (a) the ability of the MM5 model to assimilate downward vertical velocities at least as high as the analyses, (b) the ability of the Graces model to incorporate the vertical velocities from MM5, and (c) other factors related to transport patterns required to transport CO in the observed manner. We carried out improved calculations of the transport of tracers for both NASA airborne missions, SONEX and PEM-Tropics. We also made improved source strength estimates fopr isoprene dust, and similar emissions from the Earths surface. This required the use of newly available databases on the Earth's surface and vegetation. We completed atmospheric chemistry simulations of radicals and nitrogen oxide species. We have greatly improved the handling of cumulonimnbus convection by modifing an existing scheme.

Hamill, Patrick↗
NASA NTRS2003

Modeling the Absorbing Aerosol Index

We propose a scheme to model the absorbing aerosol index and improve the biomass carbon inventories by optimizing the difference between TOMS aerosol index (AI) and modeled AI with an inverse model. Two absorbing aerosol types are considered, including biomass carbon and mineral dust. A priori biomass carbon source was generated by Liousse et al [1996]. Mineral dust emission is parameterized according to surface wind and soil moisture using the method developed by Ginoux [2000]. In this initial study, the coupled CCM1 and GRANTOUR model was used to determine the aerosol spatial and temporal distribution. With modeled aerosol concentrations and optical properties, we calculate the radiance at the top of the atmosphere at 340 nm and 380 nm with a radiative transfer model. The contrast of radiance at these two wavelengths will be used to calculate AI. Then we compare the modeled AI with TOMS AI. This paper reports our initial modeling for AI and its comparison with TOMS Nimbus 7 AI. For our follow-on project we will model the global AI with aerosol spatial and temporal distribution recomputed from the IMPACT model and DAO GEOS-1 meteorology fields. Then we will build an inverse model, which applies a Bayesian inverse technique to optimize the agreement of between model and observational data. The inverse model will tune the biomass burning source strength to reduce the difference between modelled AI and TOMS AI. Further simulations with a posteriori biomass carbon sources from the inverse model will be carried out. Results will be compared to available observations such as surface concentration and aerosol optical depth.

Penner, Joyce↗
NASA NTRS2005

Trace Gas/Aerosol Interactions and GMI Modeling Support

Current global aerosol models use different physical and chemical schemes and parameters, different meteorological fields, and often different emission sources. Since the physical and chemical parameterization schemes are often tuned to obtain results that are consistent with observations, it is difficult to assess the true uncertainty due to meteorology alone. Under the framework of the NASA global modeling initiative (GMI), the differences and uncertainties in aerosol simulations (for sulfate, organic carbon, black carbon, dust and sea salt) solely due to different meteorological fields are analyzed and quantified. Three meteorological datasets available from the NASA DAO GCM, the GISS-II' GCM, and the NASA finite volume GCM (FVGCM) are used to drive the same aerosol model. The global sulfate and mineral dust burdens with FVGCM fields are 40% and 20% less than those with DAO and GISS fields, respectively due to its heavier rainfall. Meanwhile, the sea salt burden predicted with FVGCM fields is 56% and 43% higher than those with DAO and GISS, respectively, due to its stronger convection especially over the Southern Hemispheric Ocean. Sulfate concentrations at the surface in the Northern Hemisphere extratropics and in the middle to upper troposphere differ by more than a factor of 3 between the three meteorological datasets. The agreement between model calculated and observed aerosol concentrations in the industrial regions (e.g., North America and Europe) is quite similar for all three meteorological datasets. Away from the source regions, however, the comparisons with observations differ greatly for DAO, FVGCM and GISS, and the performance of the model using different datasets varies largely depending on sites and species. Global annual average aerosol optical depth at 550 nm is 0.120-0.131 for the three meteorological datasets.

Penner, Joyce E.↗