Search NASA⌕ Search

SEARCH · Search NASA

Results for “Variational data assimilation”

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 145 records · Page 8

A VAS-numerical model impact study using the Gal-Chen variational approach

A numerical study based on the use of a variational assimilation technique of Gal-Chen (1983, 1986) was conducted to assess the impact of incorporating temperature data from the VISSR Atmospheric Sounder (VAS) into a regional-scale numerical model. A comparison with the results of a control forecast using only conventional data indicated that the assimilation technique successfully combines actual VAS temperature observations with the dynamically balanced model fields without destabilizing the model during the assimilation cycle. Moreover, increasing the temporal frequency of VAS temperature insertions during the assimilation cycle was shown to enhance the impact on the model forecast through successively longer forecast periods. The incorporation of a nudging technique, whereby the model temperature field is constrained toward the VAS 'updated' values during the assimilation cycle, further enhances the impact of the VAS temperature data.

Aune, Robert M.↗

Interannual Variation in Phytoplankton Primary Production at a Global Scale

We used the NASA Ocean Biogeochemical Model (NOBM) combined with remote sensing data via assimilation to evaluate the contribution of four phytoplankton groups to the total primary production. First, we assessed the contribution of each phytoplankton groups to the total primary production at a global scale for the period 1998-2011. Globally, diatoms contributed the most to the total phytoplankton production ((is)approximately 50%, the equivalent of 20 PgC·y1). Coccolithophores and chlorophytes each contributed approximately 20% ((is) approximately 7 PgC·y1) of the total primary production and cyanobacteria represented about 10% ((is) approximately 4 PgC·y1) of the total primary production. Primary production by diatoms was highest in the high latitudes ((is) greater than 40 deg) and in major upwelling systems (Equatorial Pacific and Benguela system). We then assessed interannual variability of this group-specific primary production over the period 1998-2011. Globally the annual relative contribution of each phytoplankton groups to the total primary production varied by maximum 4% (1-2 PgC·y1). We assessed the effects of climate variability on group-specific primary production using global (i.e., Multivariate El Niño Index, MEI) and "regional" climate indices (e.g., Southern Annular Mode (SAM), Pacific Decadal Oscillation (PDO) and North Atlantic Oscillation (NAO)). Most interannual variability occurred in the Equatorial Pacific and was associated with climate variability as indicated by significant correlation (p (is) less than 0.05) between the MEI and the group-specific primary production from all groups except coccolithophores. In the Atlantic, climate variability as indicated by NAO was significantly correlated to the primary production of 2 out of the 4 groups in the North Central Atlantic (diatoms/cyanobacteria) and in the North Atlantic (chlorophytes and coccolithophores). We found that climate variability as indicated by SAM had only a limited effect on group-specific primary production in the Southern Ocean. These results provide a modeling and data assimilation perspective to phytoplankton partitioning of primary production and contribute to our understanding of the dynamics of the carbon cycle in the oceans at a global scale.

Coccolithophores↗

Estimation of Transpiration and Water Use Efficiency Using Satellite and Field Observations

Structure and function of terrestrial plant communities bring about intimate relations between water, energy, and carbon exchange between land surface and atmosphere. Total evaporation, which is the sum of transpiration, soil evaporation and evaporation of intercepted water, couples water and energy balance equations. The rate of transpiration, which is the major fraction of total evaporation over most of the terrestrial land surface, is linked to the rate of carbon accumulation because functioning of stomata is optimized by both of these processes. Thus, quantifying the spatial and temporal variations of the transpiration efficiency (which is defined as the ratio of the rate of carbon accumulation and transpiration), and water use efficiency (defined as the ratio of the rate of carbon accumulation and total evaporation), and evaluation of modeling results against observations, are of significant importance in developing a better understanding of land surface processes. An approach has been developed for quantifying spatial and temporal variations of transpiration, and water-use efficiency based on biophysical process-based models, satellite and field observations. Calculations have been done using concurrent meteorological data derived from satellite observations and four dimensional data assimilation for four consecutive years (1987-1990) over an agricultural area in the Northern Great Plains of the US, and compared with field observations within and outside the study area. The paper provides substantive new information about interannual variation, particularly the effect of drought, on the efficiency values at a regional scale.

Choudhury, Bhaskar J.↗

Regional Climate and Variability of the Summertime Continental United States in Reanalyses

Understanding climate variability at regional scales is an important for research and societal needs. Atmospheric retrospective-analyses (or reanalyses) integrate multitudes of observing systems with numerical models to produce continuous data that include variables not easily observed, if at all. The breadth of variables as well as observational influence included in reanalyses make them ideal for investigating climate variability. In this paper, we assess NASA s Modern Era Retrospective-analysis for Research and Application (MERRA) regional variability in North America, specifically the United States, in conjunction with current satellite data reanalyses. Emphasis is placed on summertime precipitation because 1) it is a difficult parameter to capture in the most difficult season, 2) significant observational resources exist to benchmark comparisons, and 3) accurate assessment of precipitation variability is crucial to a multitude of sectors and applications. Likewise, we have also begun to evaluate surface air temperature. While precipitation biases are identified, year to year variability of the precipitation variations, in many cases, are quite reasonable. However, some spurious long term trends and sudden shifts in the time series are also identified. In surface air temperature, analysis of station observations provides ERA Interim a clear overall advantage. However, in a number of regions, all the reanalyses are quite comparable in variability and trend. In other regions, significant precipitation biases may occur, which has implications for the ancillary process data in a reanalysis, such as surface fluxes. We also characterize the reanalyses ability to capture variability related to ENSO. In general, the summertime variations of precipitation in the reanalyses are more highly correlated (positively) to ENSO (using ENSO34) than are the observations. The Northwestern US shows the largest positive correlations to ENSO34, and reanalyses agree with that, and correlate highly with the gauge observations there. Additional evaluation of the data assimilation using the MERRA Gridded Innovations and Observations (GIO) data, which consists of the assimilated observations as well as forecast an analysis error fields. This work represents an initial evaluation of using the MERRA reanalysis to study regional climate variations of the United States, and identifying its applied limitations

Bosilovich, M G.↗

Quantifying Errors in TRMM-Based Multi-Sensor QPE Products Over Land in Preparation for GPM

Determining uncertainties in satellite-based multi-sensor quantitative precipitation estimates over land of fundamental importance to both data producers and hydro climatological applications. ,Evaluating TRMM-era products also lays the groundwork and sets the direction for algorithm and applications development for future missions including GPM. QPE uncertainties result mostly from the interplay of systematic errors and random errors. In this work, we will synthesize our recent results quantifying the error characteristics of satellite-based precipitation estimates. Both systematic errors and total uncertainties have been analyzed for six different TRMM-era precipitation products (3B42, 3B42RT, CMORPH, PERSIANN, NRL and GSMap). For systematic errors, we devised an error decomposition scheme to separate errors in precipitation estimates into three independent components, hit biases, missed precipitation and false precipitation. This decomposition scheme reveals hydroclimatologically-relevant error features and provides a better link to the error sources than conventional analysis, because in the latter these error components tend to cancel one another when aggregated or averaged in space or time. For the random errors, we calculated the measurement spread from the ensemble of these six quasi-independent products, and thus produced a global map of measurement uncertainties. The map yields a global view of the error characteristics and their regional and seasonal variations, reveals many undocumented error features over areas with no validation data available, and provides better guidance to global assimilation of satellite-based precipitation data. Insights gained from these results and how they could help with GPM will be highlighted.

Peters-Lidard, Christa D.↗

Insight and Evidence Motivating the Simplification of Dual-Analysis Hybrid Systems into Single-Analysis Hybrid Systems

Many hybrid data assimilation systems currently used for NWP employ some form of dual-analysis system approach. Typically a hybrid variational analysis is responsible for creating initial conditions for high-resolution forecasts, and an ensemble analysis system is responsible for creating sample perturbations used to form the flow-dependent part of the background error covariance required in the hybrid analysis component. In many of these, the two analysis components employ different methodologies, e.g., variational and ensemble Kalman filter. In such cases, it is not uncommon to have observations treated rather differently between the two analyses components; recentering of the ensemble analysis around the hybrid analysis is used to compensated for such differences. Furthermore, in many cases, the hybrid variational high-resolution system implements some type of four-dimensional approach, whereas the underlying ensemble system relies on a three-dimensional approach, which again introduces discrepancies in the overall system. Connected to these is the expectation that one can reliably estimate observation impact on forecasts issued from hybrid analyses by using an ensemble approach based on the underlying ensemble strategy of dual-analysis systems. Just the realization that the ensemble analysis makes substantially different use of observations as compared to their hybrid counterpart should serve as enough evidence of the implausibility of such expectation. This presentation assembles numerous anecdotal evidence to illustrate the fact that hybrid dual-analysis systems must, at the very minimum, strive for consistent use of the observations in both analysis sub-components. Simpler than that, this work suggests that hybrid systems can reliably be constructed without the need to employ a dual-analysis approach. In practice, the idea of relying on a single analysis system is appealing from a cost-maintenance perspective. More generally, single-analysis systems avoid contradictions such as having to choose one sub-component to generate performance diagnostics to another, possibly not fully consistent, component.

Todling, Ricardo↗

Application of an Ensemble Smoother to Precipitation Assimilation

Assimilation of precipitation in a global modeling system poses a special challenge in that the observation operators for precipitation processes are highly nonlinear. In the variational approach, substantial development work and model simplifications are required to include precipitation-related physical processes in the tangent linear model and its adjoint. An ensemble based data assimilation algorithm "Maximum Likelihood Ensemble Smoother (MLES)" has been developed to explore the ensemble representation of the precipitation observation operator with nonlinear convection and large-scale moist physics. An ensemble assimilation system based on the NASA GEOS-5 GCM has been constructed to assimilate satellite precipitation data within the MLES framework. The configuration of the smoother takes the time dimension into account for the relationship between state variables and observable rainfall. The full nonlinear forward model ensembles are used to represent components involving the observation operator and its transpose. Several assimilation experiments using satellite precipitation observations have been carried out to investigate the effectiveness of the ensemble representation of the nonlinear observation operator and the data impact of assimilating rain retrievals from the TMI and SSM/I sensors. Preliminary results show that this ensemble assimilation approach is capable of extracting information from nonlinear observations to improve the analysis and forecast if ensemble size is adequate, and a suitable localization scheme is applied. In addition to a dynamically consistent precipitation analysis, the assimilation system produces a statistical estimate of the analysis uncertainty.

Zhang, Sara↗

Sampling Biases of Space-Based Observations of XCO2 Associated With South American Biomass Burning Events During 2017-2020

- We examine time series of OCO-2 xCO 2 observations coincident with other remotely sensed and in-situ observations to better understand external sampling effects. - The OCO-2 Level 3 assimilation product fills in gaps where OCO-2 does not have observations. - Although the Level 3 clearly shows the increased CO 2 over the Amazon associated with the September 2017 biomass burning event, the total amount of CO 2 may be underestimated because of a sampling bias. - Since TCCON is also a remotely sensed product, it may suffer the same sampling biases associated with biomass burning events as OCO-2. - Aircraft data show that there is significant variation in xCO 2 with altitude. - The OCO-2 fused products created by kriging may be better than the assimilated product at detecting deviations from the average state.

OCO↗

Technical Report Series on Global Modeling and Data Assimilation: A Comparison of GEOS Assimilated Data with FIFE Observations - Volume 14

First ISLSCP Field Experiment (FIFE) observations have been used to validate the near-surface proper- ties of various versions of the Goddard Earth Observing System (GEOS) Data Assimilation System. The site- averaged FIFE data set extends from May 1987 through November 1989, allowing the investigation of several time scales, including the annual cycle, daily means and diurnal cycles. Furthermore, the development of the daytime convective planetary boundary layer is presented for several days. Monthly variations of the surface energy budget during the summer of 1988 demonstrate the affect of the prescribed surface soil wetness boundary conditions. GEOS data comes from the first frozen version of the assimilation system (GEOS-1 DAS) and two experimental versions of GEOS (v. 2.0 and 2.1) with substantially greater vertical resolution and other changes that influence the boundary layer. This report provides a baseline for future versions of the GEOS data assimilation system that will incorporate a state-of-the-art land surface parameterization. Several suggestions are proposed to improve the generality of future comparisons. These include the use of more diverse field experiment observations and an estimate of gridpoint heterogeneity from the new land surface parameterization.

Bosilovich, Michael G.↗

Technical Report Series on Global Modeling and Data Assimilation: Interannual Variability and Potential Predictability in Reanalysis Products - Volume 13

The Data Assimilation Office (DAO) at Goddard Space Flight Center and the National Center for Environmental Prediction and National Center for Atmospheric Research (NCEP/NCAR) have produced multi-year global assimilations of historical data employing fixed analysis systems. These "reanalysis" products are ideally suited for studying short-term climatic variations. The availability of multiple reanalysis products also provides the opportunity to examine the uncertainty in the reanalysis data. The purpose of this document is to provide an updated estimate of seasonal and interannual variability based on the DAO and NCEP/NCAR reanalyses for the 15-year period 1980-1995. Intercomparisons of the seasonal means and their interannual variations are presented for a variety of prognostic and diagnostic fields. In addition, atmospheric potential predictability is re-examined employing selected DAO reanalysis variables.

Min, Wei↗

Assimilation of Cloud- and Land-Affected TOVS/ATOVS Level 1B Radiances in DAO's Next Generation Finite-Volume Data Assimilation System

The Physical-space/Finite-volume Data Assimilation System (fvDAS) is the next generation global atmospheric data assimilation system in development at the Data Assimilation Office (DAO) at NASA's Goddard Space Flight Center. It is based on a new finite-volume general circulation model jointly developed by NASA and NCAR, and on the Physical-Space Statistical Analysis System (PSAS) developed at the DAO. In this talk we will focus on the assimilation of data from the (Advanced) TIROS Operational Vertical Sounder (ATOVS), with emphasis on the impact of cloud- and land-affected level 1B radiances. Recently, it has been shown that the use of observations from satellite-borne microwave and infrared radiometers in data assimilation systems consistently increases forecast skill. Considerable effort has been expended over the past two decades, particularly with the (Advanced) TIROS Operational Vertical Sounder (ATOVS), to achieve this result. The positive impact on forecast skill has resulted from improvements in quality control algorithms, systematic error correction schemes, and more sophisticated data assimilation algorithms. Despite these advances, there are still many issues regarding the use of satellite data in data assimilation systems that remain unresolved. In particular, most operational centers still do not assimilate cloud- and land-affected TOVS data. In this study, we evaluate the impact of assimilating cloud-and land-affected TOVS/ATOVS level 1B data in DAO's next generation fvDAS, using a 1D variational scheme. We will discuss the impact of these data on both tropospheric and stratospheric forecasts, as well as on the general aspects of the earth climate system.

Joiner, J.↗

Apollo 14 mare basalt petrogenesis - Assimilation of KREEP-like components by a fractionating magma

Based on analysis of whole rock and mineral data from 22 samples, it is shown that the trace element characteristics of 14321 high-AL mare basalts exhibit variations that cannot be related to short-range mixing alone. It is suggested that the five groups defined for Apollo 14 high-Al basalts are due to sampling. It is found that the addition of new data results in a continuum of compositions. Detailed trace element modeling demonstrates that a primitive parental magma evolves by assimilation and fractional crystallization (AFC) with KREEP. The ratio of mass assimilated to mass crystallized is found to be 0.22. It is suggested that this AFC process occurs within the lunar crust and that the fractionating phases are effectively removed from the system.

Neal, Clive R.↗

Seasonal and Interannual Variations of Evaporation and their Relations with Precipitation, Net Radiation, and Net Carbon Accumulation for the Gediz Basin Area

A model combining the rate of carbon assimilation with water and energy balance equations has been run using satellite and ancillary data for a period of 60 months (January 1986 to December 1990). Calculations for the Gediz basin area give mean annual evaporation as 395 mm, which is composed of 45% transpiration, 42% soil evaporation and 13% interception. The coefficient of interannual variation of evaporation is found to be 6%, while that for precipitation and net radiation are, respectively, 16% and 2%, illustrating that net radiation has an important effect in modulating interannual variation of evaporation. The mean annual water use efficiency (i.e., the ratio of net carbon accumulation and total evaporation) is ca. 1 g/sq m/mm, and has a coefficient of interannual variation of 5%. A comparison of the mean water use efficiency with field observations suggests that evaporation over the area is utilized well for biomass production. The reference crop evaporation for irrigated areas has annual mean and coefficient of variation as, respectively, 1176 mm and 3%. The total evaporation during three summer months of peak evaporation (June-August) is estimated to be about 575 mm for irrigated crops like maize and cotton. Seasonal variations of the fluxes are presented.

Choudhury, Bhaskar J.↗

The Vertical Error Characteristics of GOES-derived Winds: Description and Impact on Numerical Weather Prediction

Errors in the height assignment of some satellite-derived winds exist because the satellites sense radiation emitted from a finite layer of the atmosphere rather than a specific level. Potential problems in data assimilation may arise because the motion of a measured layer is often represented by a single-level value. In this research, cloud and water vapor motion winds that are derived from the Geostationary Operational Environmental Satellites (GOES winds) are compared to collocated rawinsonde observations (RAOBs). An important aspect of this work is that in addition to comparisons at each assigned height, the GOES winds are compared to the entire profile of the collocated RAOB data to determine the vertical error characteristics of the GOES winds. The impact of these results on numerical weather prediction is then investigated. The comparisons at individual vector height assignments indicate that the error of the GOES winds range from approx. 3 to 10 m/s and generally increase with height. However, if taken as a percentage of the total wind speed, accuracy is better at upper levels. As expected, comparisons with the entire profile of the collocated RAOBs indicate that clear-air water vapor winds represent deeper layers than do either infrared or water vapor cloud-tracked winds. This is because in cloud-free regions the signal from water vapor features may result from emittance over a thicker layer. To further investigate characteristics of the clear-air water vapor winds, they are stratified into two categories that are dependent on the depth of the layer represented by the vector. It is found that if the vertical gradient of moisture is smooth and uniform from near the height assignment upwards, the clear-air water vapor wind tends to represent a relatively deep layer. The information from the comparisons is then used in numerical model simulations of two separate events to determine the forecast impacts. Four simulations are performed for each case: 1) A control simulation that assimilates no satellite wind data, 2) assimilation of all GOES winds according to their assigned single level height, 3) assimilation of all GOES winds spread over multiple levels, and 4) assimilation of all GOES winds spread over multiple levels, but with variations in the vertical influence of clear-air water vapor winds based on the moisture profile in the model. In the first case, a strong mid-latitude cyclone is present and the use of the satellite data results in improved storm tracks during the initial approx. 36 h forecast period. This is because the satellite data improves the analysis of the environment into which the storm progresses. Statistics for mean wind vector and height differences show that, with the exception of the height field at later times in the first case, the use of GOES winds improves the simulation with time. The simulation results suggest that it is beneficial to spread the GOES wind information over multiple levels, particularly when the moisture profile is used to define the vertical influence.

Rao, P. Anil↗

Interannual Variation in Phytoplankton Class-Specific Primary Production at a Global Scale

We used the NASA Ocean Biogeochemical Model (NOBM) combined with remote sensing data via assimilation to evaluate the contribution of 4 phytoplankton groups to the total primary production. First we assessed the contribution of each phytoplankton groups to the total primary production at a global scale for the period 1998-2011. Globally, diatoms were the group that contributed the most to the total phytoplankton production (50, the equivalent of 20 PgC y-1. Coccolithophores and chlorophytes each contributed to 20 (7 PgC y-1 of the total primary production and cyanobacteria represented about 10 (4 PgC y(sub-1) of the total primary production. Primary production by diatoms was highest in high latitude (45) and in major upwelling systems (Equatorial Pacific and Benguela system). We then assessed interannual variability of this group-specific primary production over the period 1998-2011. Globally the annual relative contribution of each phytoplankton groups to the total primary production varied by maximum 4 (1-2 PgC y-1. We assessed the effects of climate variability on the class-specific primary production using global (i.e. Multivariate El Nio Index, MEI) and regional climate indices (e.g. Southern Annular Mode (SAM), Pacific Decadal Oscillation (PDO) and North Atlantic Oscillation (NAO)). Most interannual variability occurred in the Equatorial Pacific and was associated with climate variability as indicated by significant correlation (p 0.05) between the MEI and the class-specific primary production from all groups except coccolithophores. In the Atlantic, climate variability as indicated by NAO was significantly correlated to the primary production of 2 out of the 4 groups in the North Central Atlantic (diatomscyanobacteria) and in the North Atlantic (chlorophytes and coccolithophores). We found that climate variability as indicated by SAM had only a limited effect on the class-specific primary production in the Southern Ocean. These results provide a modeling and data assimilation perspective to phytoplankton partitioning of primary production and contribute to our understanding of the dynamics of the carbon cycle in the oceans at a global scale.

Chl-a↗

The Bi-Modal Pattern of the Summer Circulation Over South America

Submonthly variations in warm-season (January-February) precipitation over South America, in special over the Amazon basin, central southwest Brazil, north Argentina, and Paraguay are shown to be strongly linked to variations in the moisture entering the continent from the Atlantic ocean. Two distinct regimes of lower tropospheric winds (westerlies and easterlies) were observed in Rondonia during the Wet Season Atmospheric Mesoscale Campaign (WETAMC) component of the Large Scale Atmosphere-Biosphere Experiment in Amazonia (LBA) and the Tropical Rainfall Measuring Mission (TRMM) field campaign. The westerly (easterly) winds were associated with the strong (weak) convective activity over the South Atlantic Convergence Zone (SACZ). The whole period of this study (January-February) was divided into SACZ and NSACZ (No SACZ) events. The vertically integrated moisture fluxes over the Amazon and Prata basin from the National Aeronautics and Space Administration/Goddard Data Assimilation Office (NASA/DAO) assimilation show that during SACZ (NSACZ) event strong (weak) convergence occurred over the Amazon basin with divergence (convergence) over the Prata basin. Submonthly variations in the SACZ also can be linked to extreme climate anomalies such as droughts or flooding conditions over the Amazon and Prata basin.

Herdies, Dirceu Luis↗

Effects of the Observed Meridional Flow Variations since 1996 on the Sun's Polar Fields

The cause of the low and extended minimum in solar activity between Sunspot Cycles 23 and 24 was the small size of Sunspot Cycle 24 itself - small cycles start late and leave behind low minima. Cycle 24 is small because the polar fields produced during Cycle 23 were substantially weaker than those produced during the previous cycles and those (weak) polar fields are the seeds for the activity of the following cycle. The polar fields are produced by the latitudinal transport of magnetic flux that emerged in low-latitude active regions. The polar fields thus depend upon the details of both the flux emergence and the flux transport. We have measured the flux transport flows (differential rotation, meridional flow, and supergranules) since 1996 and find systematic and substantial variation in the meridional flow alone. Here we present experiments using a Surface Flux Transport Model in which magnetic field data from SOHO/MDI and SDO/HMI are assimilated into the model only at latitudes between 45-degrees north and south of the equator (this assures that the details of the active region flux emergence are well represented). This flux is then transported in both longitude and latitude by the observed flows. In one experiment the meridional flow is given by the time averaged (and north-south symmetric) meridional flow profile. In the second experiment the time-varying and north-south asymmetric meridional flow is used. Differences between the observed polar fields and those produced in these two experiments allow us to ascertain the effects of these meridional flow variations on the Sun s polar fields.

Hathaway, David↗

Effects of the Observed Meridional Flow Variations since 1996 on the Sun's Polar Fields

The cause of the low and extended minimum in solar activity between Sunspot Cycles 23 and 24 was the small size of Sunspot Cycle 24 itself - small cycles start late and leave behind low minima. Cycle 24 is small because the polar fields produced during Cycle 23 were substantially weaker than those produced during the previous cycles and those (weak) polar fields are the seeds for the activity of the following cycle. The polar fields are produced by the latitudinal transport of magnetic flux that emerged in low-latitude active regions. The polar fields thus depend upon the details of both the flux emergence and the flux transport. We have measured the flux transport flows (differential rotation, meridional flow, and supergranules) since 1996 and find systematic and substantial variation in the meridional flow alone. Here we present experiments using a Surface Flux Transport Model in which magnetic field data from SOHO/MDI and SDO/HMI are assimilated into the model only at latitudes between 45-degrees north and south of the equator (this assures that the details of the active region flux emergence are well represented). This flux is then transported in both longitude and latitude by the observed flows. In one experiment the meridional flow is given by the time averaged (and north-south symmetric) meridional flow profile. In the second experiment the time-varying and north-south asymmetric meridional flow is used. Differences between the observed polar fields and those produced in these two experiments allow us to ascertain the effects of these meridional flow variations on the Sun s polar fields.

Hathaway, David H.↗