Search NASA⌕ Search

SEARCH · Search NASA

Results for “circulation modification”

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.

105 records · Page 6

MODIS Aerosol Observations used to Constrain Dust Distributions and Lifecycle in the NASA GEOS-5 Model

Approximately 240 Tg of mineral dust aerosol are transported annually from Saharan Africa to the Atlantic Ocean. Dust affects the Earth radiation budget, and plays direct (through scattering and absorption of radiation) and indirect (through modification of cloud properties and environment) roles in climate. Deposition of dust to the surface provides an important nutrient source to terrestrial and oceanic ecosystems. Dust is additionally a contributor to adverse air quality. Among the tools toward understanding the lifecycle and impacts of mineral dust aerosols are numerical models. Important constraints on these models come from quantitative satellite observations, like those from the space-based Moderate Resolution Imaging Spectroradiometer (MODIS). In particular, Kauhan et al. [2005] used MODIS aerosol observations to infer transport and deposition fluxes of Saharan dust over the Atlantic, Caribbean, and Amazonian basins. Those observations are used here to constrain the transport of dust and its interannual variability simulated in the NASA GEOS-5 general circulation model and data assimilation system. Significant uncertainty exists in the MODIS-derived fluxes, however, due to uncertainty in the wind fields provided by meteorological analyses in this region. That same uncertainty in the wind fields is manifest in our GEOS-5 simulations of dust distributions. Here we use MODIS observations to investigate the seasonality and location of the Saharan dust plume and explore through sensitivity analysis of our model the meteorological controls on the dust distribution, including dust direct radiative effects and sub-gridscale source and sink processes.

Colarco, P.↗

On the structure of the Southern Oscillation

A realistic oceanic general circulation model is forced with winds observed over the tropical Pacific between 1967 and 1979. The structure of the simulated Southern Oscillation is strikingly different in the western and eastern sides of the basin, because the principal interannual zonal-wind fluctuations are confined to the west and are in the form of an equatorial jet. This causes thermocline displacements to have maxima off the equator in the west (where the curl of the wind is large) but on the equator in the east. Zonal phase propagation, both on and off the equator, is at different speeds in the west and east. The phase pattern is complex, and there is, on interannual time scale, no explicit evidence of individual equatorial waves. These results lead to a modification of the 'delayed oscillator' mechanism originally proposed by Schopf and Suarez to explain a continual Southern Oscillation. The results also permit an evaluation of the various coupled ocean-atmosphere models that simulate the Southern Oscillation and indicate which measurements are necessary to determine which models are most relevant to reality.

Chao, YI↗

Method for Impeding Degradation of Porous Silicon Structures

This invention relates to surface modification of porosified silicon (pSi) structures with poly(alkylene) glycols for the purpose of controlled degradation of the silicon matrix and tailored release of encapsulated substances for biomedical applications. The pSi structures are currently used in diverse biomedical applications including bio-molecular screening, optical bio-sensoring, and drug delivery by means of injectable/orally administered carriers and implantable devices. The size of the pores and the surface chemistry of the pSi structure can be controlled during the microfabrication process and thereafter. A fine regulation of the degradation kinetics of mesoporous silicon structures is of fundamental importance. Polyethylene glycols (PEGs) represent the major category of surface modifying agents used in classical drug delivery systems and in pharmaceutical dosage forms. PEGylation enables avoidance of RES uptake, thus prolonging circulation time of intravenously injectable nanovectors. PEG molecules demonstrate little toxicity and immunogenicity, and are cleared from the body through the urine (molecular weight, MW less than 30 kDa) or in the feces (MW greater than 30kDa). The invention focuses on the possibility of finely tuning the degradation kinetics of the pSi nanovectors and other structures through surface conjugation of PEGs with various backbone lengths/MWs. To prove the concept, pSi nanovectors were covalently conjugated to seven PEGs with MW from 245 to 5,000 Da and their degradation kinetics in physiologically relevant media (phosphate buffer saline, PBS pH7.4, and fetal bovine serum) was assessed by the elemental analysis of the Si using inductive coupled plasma atomic emission spectroscopy (ICP-AES). The conjugation of the PEG with lowest MW to the nanovectors surface did not induce any change in the degradation kinetics in serum, but inhibited degradation and consequently the release of orthosilicic acid into buffer. When PEGs with the longer chains were evaluated, Si mass loss from the nanovectors was slowed down, and the PEGylated structures were almost fully degraded within 18 24 hours in serum and within 48 hours in PBS. The most dramatic effect was observed for high MW PEGs 3,400 and 5,000 Da, which prominently inhibited the degradation of the systems, with complete degradation achieved only after four days. For these PEGs, during the early stages of the degradation, there was a lag period of little or no Si mass loss from the nanovector. The obtained profiles were in agreement with the erosion of the nanovector surface as observed by scanning electron microscopy.

Vilentchouk, Biana Godin↗

Evaluation of NESS soundings used in the FGGE special effort

A special effort to produce high quality edited and enhanced global data sets is being conducted for the two special observing periods of the first GARP global experiment. More than 60 cases were selected for enhancement on the basis of meteorological interest. These cases include situations of blocking, cut off low development, cyclogenesis, and tropical circulations. The sounding data enhancement process consists of supplementing the operational satellite sounding data set with higher resolution sounding in meteorologically active regions and with new soundings where data voids or soundings of questionable quality exist. After the enhancement process is completed, a final editing and quality assessment of the enhanced data are performed. The data is archived for further analyses. Subjective evaluation of the effects of sounding data enhancement indicate that: (1) enhanced and operational IR retrievals are similar in most situations; (2) large improvements in 1000 to 300 mb thickness and mandatory level temperature and intensification of atmospheric thermal gradients occasionally occur; (3) interactive processing can fill in major gaps in the TIROS-N coverage; and (4) assimilation of special effort soundings can result in major modifications to large scale analyses and prognoses.

Atlas, R.↗

Detailed studies of aviation fuel flowability

Six Jet A fuels, with varying compositions, were tested for low temperature flowability in a 190-liter simulator tank that modeled a section of a wing tank of a wide-body commercial airplane. The insulated tank was chilled by circulating coolant through the upper and lower surfaces. Flow-ability was determined as a function of fuel temperature by holdup, the fraction of unflowable fuel remaining in the tank after otherwise complete withdrawal. In static tests with subfreezing tank conditions, hold up varied with temperature and fuel composition. However, a general correlation of two or three classes of fuel type was obtained by plotting holdup as a function of the difference between freezing point and boundary-layer temperature, measured 0.6 cm above the bottom tank surface. Dynamic conditions of vibrations and slosh or rate of fuel withdrawal had very minor effects on holdup. Tests with cooling schedules to represent extreme, cold-day flights showed, at most, slight holdup for any combination of fuel type or dynamic conditions. Tests that superimposed external fuel heating and recirculation during the cooldown period indicates reduced hold up by modification of the low-temperature boundary layer. Fuel heating was just as effective when initiated during the later times of the tests as when applied continuously.

Mehta, H. K.↗

Stars and linear dunes on Mars

A field containing 11 star and incipient star dunes occurs on Mars at 8.8 deg S, 270.9 deg W. Examples of linear dunes are found in a crater at 59.4 deg S, 343 deg W. While rare, dune varieties that form in bi- and multidirectional wind regimes are not absent from the surface of Mars. The occurence of both of these dune fields offers new insight into the nature of martian wind conditions and sand supply. The linear dunes appears to have formed through modification of a formerly transverse aeolian deposit, suggesting a relatively recent change in local wind direction. The 11 dunes in the star dune locality show a progressive change from barchan to star form as each successive dune has traveled up into a valley, into a more complex wind regime. The star dunes corroborate the model of N. Lancaster (1989), for the formation of star dunes by projection of transverse dunes into a complex, topographically influenced wind regime. The star dunes have dark streaks emanating from them, providing evidence that the dunes were active at or near the time the relevant image was obtained by the Viking 1 orbiter in 1978. The star and linear dunes described here are located in different regions on the martian surface. Unlike most star and linear dunes on Earth, both martian examples are isolated occurrences; neither is part of a major sand sea. Previously published Mars general circulation model results suggest that the region in which the linear dune field occurs should be a bimodal wind regime, while the region in which the star dunes occur should be unimodal. The star dunes are probably the result of localized complication of the wind regime owing to topographic confinement of the dunes. Local topographic influence on wind regime is also evident in the linear dune field, as there are transverse dunes in close proximity to the linear dunes, and their occurrence is best explained by funneling of wind through a topographic gap in the upwind crater wall.

Edgett, Kenneth S.↗

Interpretation of lidar and satellite data sets using a global photochemical model

The status in the beginning of the report period was that the existing General Circulation Model (GCM) was running with a chemistry module compiled for stratospheric simulation studies. The chemistry simulation was not working sufficiently in the troposphere and any tropospheric trace gas sources or dry deposition sinks were not yet incorporated. The current status concerning the chemistry module is that the chemistry simulation has been modified to also simulate the chemistry in the troposphere with resulting mixing ratios close to other model simulations as described in Olson et al. (1996). The mechanism to incorporate trace gas source and dry deposition sinks, testing for H202, CH300H, 03, HCHO, HN03, and N02, are incorporated and is currently being tested. Existing model and development versions include: the full GCM model, currently still running with the original stratospheric chemistry module; an off-line version of the GCM, i.e. wind and photolysis rates are pre-calculated and prescribed in read-in arrays; a box model version of the modified chemistry module for developments and first tests of new modifications; and a box model with the same chemistry simulated but flexible partitioning and integration methods for test purposes of those. The grantee's work focused mainly on three areas: trace gas sources and dry deposition; method to introduce a NO source instead of a NO(x) source; and investigating integration methods.

Chyba, Thomas↗

Large Eddy Simulations of Upper-Ocean Turbulence Across the Equatorial Pacific Cold Tongue

Turbulence in the equatorial Pacific cold tongue plays multiple roles in the sea-surface temperature budget. However, observations of turbulence are limited to a handful of field experiments lasting from a few days to one month at 0 N , 140 W, and the key turbulent fluxes can only be inferred. Although explicit simulations of turbulence resolve these fluxes, the required large eddy simulations (LES) are costly and thus also limited to a few weeks. A two-fold hypothesis motivates this contribution to the PUMP pre-field modelling experiments in support of TPOS: 1) that there exist useful relationships between the surface wind, wave, and buoyancy forcing, oceanic shear and stratification, and turbulent mixing, but 2) these relationships differ between regions of equatorial dynamics (f~0) and higher latitudes, such that the state space is insufficiently sampled to understand these relationships and represent them in ocean general circulation models (OGCM). To address the lack of observations, we simulated 35 continuous days of upper ocean turbulence using LES during October-November 1985 at sites in two regions (70 days in total) of the cold tongue: 0 N, 140 W and 3 N, 140 W. The model reproduces many features of the observed turbulence at 0,140W, including the “deep-cycle” observed by Gregg et al. (1985) Nature and the turbulent pulses every few hours during the night and in the early morning. Furthermore, the results at 3 N support the hypothesis advanced by Cherian et al. (2021) JPO that strong, but intermittent deep-cycle turbulence is observed to the north of the equator in the presence of strong shear associated with Tropical Instability Waves. Comparisons with the K-profile parameterization (KPP) scheme show that additional tuning of the KPP parameters can reduce differences between the turbulent fluxes produced by KPP and LES at both sites. The results will suggest modifications to OGCMs and already emphasize the need for a more extensive measurements of vertical profiles of shear, stratification and turbulent mixing to reveal the spatio-temporal structure of turbulence and how it relates to forcing, shear and stratification across the equatorial Pacific Cold tongue.

Large Eddy Simulations↗

Dependence of tropospheric temperature on the parameterization of cumulus convection in the GLAS model of the general circulation

Analysis of the simulation of seasonal change by the GLAS model of the general circulation reveals deficiencies in the simulation of tropospheric temperature and of convective cloud cover. These interrelated deficiencies are due to a spurious doubling from January to July in the convective cloud cover of the Northern Hemisphere. The spurious doubling, in turn, is due to the oversensitivity of cumulus convection, in the GLAS model, to the specific humidity of the lower atmosphere. The oversensitivity is enhanced by a feedback mechanism which perpetuates the existence of deep, penetrative convective clouds at certain preferred locations. The cumulus parameterization scheme has been modified to more realistically relate the onset of cumulus convection to the relative humidity of the lower atmosphere. The modified parameterization has improved the simulation of tropospheric temperature, planetary albedo and convective cloud cover as well as their seasonal variations. Comparison of this experiment with its control has shown a high degree of interrelation among these fields in the GLAS model and has demonstrated the sensitivity of the atmospheric heat budget to the design of the cumulus parameterization scheme. Also, the modification to the cumulus scheme has demonstrated a plausible mechanism to explain the correlation between convective cloud cover and relative humidity in the real atmosphere.

Helfand, H. M.↗

A System of Conservative Regridding for Ice-Atmosphere Coupling in a General Circulation Model (GCM)

The method of elevation classes, in which the ice surface model is run at multiple elevations within each grid cell, has proven to be a useful way for a low-resolution atmosphere inside a general circulation model (GCM) to produce high-resolution downscaled surface mass balance fields for use in one-way studies coupling atmospheres and ice flow models. Past uses of elevation classes have failed to conserve mass and energy because the transformation used to regrid to the atmosphere was inconsistent with the transformation used to downscale to the ice model. This would cause problems for two-way coupling. A strategy that resolves this conservation issue has been designed and is presented here. The approach identifies three grids between which data must be regridded and five transformations between those grids required by a typical coupled atmosphere-ice flow model. This paper develops a theoretical framework for the problem and shows how each of these transformations may be achieved in a consistent, conservative manner. These transformations are implemented in Glint2, a library used to couple atmosphere models with ice models. Source code and documentation are available for download. Confounding real-world issues are discussed, including the use of projections for ice modeling, how to handle dynamically changing ice geometry, and modifications required for finite element ice models.

Ice-atmosphere↗

Calculating Soil Wetness, Evapotranspiration and Carbon Cycle Processes Over Large Grid Areas Using a New Scaling Technique

Soil wetness typically shows great spatial variability over the length scales of general circulation model (GCM) grid areas (approx 100 km ), and the functions relating evapotranspiration and photosynthetic rate to local-scale (approx 1 m) soil wetness are highly non-linear. Soil respiration is also highly dependent on very small-scale variations in soil wetness. We therefore expect significant inaccuracies whenever we insert a single grid area-average soil wetness value into a function to calculate any of these rates for the grid area. For the particular case of evapotranspiration., this method - use of a grid-averaged soil wetness value - can also provoke severe oscillations in the evapotranspiration rate and soil wetness under some conditions. A method is presented whereby the probability distribution timction(pdf) for soil wetness within a grid area is represented by binning. and numerical integration of the binned pdf is performed to provide a spatially-integrated wetness stress term for the whole grid area, which then permits calculation of grid area fluxes in a single operation. The method is very accurate when 10 or more bins are used, can deal realistically with spatially variable precipitation, conserves moisture exactly and allows for precise modification of the soil wetness pdf after every time step. The method could also be applied to other ecological problems where small-scale processes must be area-integrated, or upscaled, to estimate fluxes over large areas, for example in treatments of the terrestrial carbon budget or trace gas generation.

Sellers, Piers↗

A meteorological interpretation of the Arctic Boundary Layer Expedition (ABLE) 3B flight series

The Arctic Boundary Layer Expedition (ABLE) 3B was conducted to determine the summertime tropospheric distribution, sources, and sinks of important trace gas and aerosol species over the wetlands and boreal forests of central and eastern Canada. Isentropic trajectories and analyzed midtropospheric circulation patterns were used to group flights according to the transport histories of polar, midlatitude, or tropical air masses which were sampled. These data were then divided into bands of potential temperature levels representing the low, middle, and maximum aircraft altitudes to assess the effects of both local and long distance transport and natural and man-made pollutants to the measured chemical species. Detailed case studies are provided to depict the complex three-dimensional airflow regimes that transported air with differing chemical signatures to the study area. Mission 6 details the large-scale movement of smoke in the generally prevailing west to northwesterly airflow that was observed on the majority of flights. Mission 1 analyzes the horizontal and vertical motions of maritime Pacific air in the upper troposphere that was routinely mixed downward to the aircraft altitude. Finally, mission 14 tracks the far northward excursion of tropical air that had been associated with a Pacific typhoon. The following three factors all had important influences on the collected chemical data sets: (1) local and distant stratospheric in puts into the upper and middle troposphere; (2) biomass-burning plumes from active fires in Alaska and Canada; (3) a band of 'low ozone' upper tropospheric air that was observed by airborne differential absorption lidar (DIAL) above the aircraft maximum altitude. Other modification factors observed on some flights included urban pollution from U.S. and Canadian cities, tropical air that had been associated with a Pacific typhoon, and precipitation scavenging by clouds and rain. Many flights were affected by several of the above factors which led to complex chemical signatures that will be discussed in other companion papers.

Shipham, Mark C.↗

Description, Usage, and Validation of the MVL-15 Modified Vortex Lattice Analysis Capability

MVL-15 is the most recent version of the Modified Vortex-Lattice (MVL) code developed within the Aerodynamics Systems Analysis Branch (ASAB) at NASA LaRC. The term "modified" refers to the primary modification of the core vortex-lattice methodology: inclusion of viscous aerodynamics tables that are linked to the linear solution via iterative processes. The inclusion of the viscous aerodynamics inherently converts the MVL-15 from a purely analytic linearized method to a semi-empirical blend which retains the rapid execution speed of the linearized method while empirically characterizing the section aerodynamics at all spanwise lattice points. The modification provides a means to assess non-linear effects on lift that occur at angles of attack near stall, and provides a means to determine the drag associated with the application of design strategies for lift augmentation such as the use of flaps or blowing. The MVL-15 code is applicable to the analyses of aircraft aerodynamics during cruise, but it is most advantageously applied to the analysis of aircraft operating in various high-lift configurations. The MVL methodology has been previously conceived and implemented; the initial concept version was delivered to the ASAB in 2001 (van Dam, C.), subsequently revised (Gelhausen, P. and Ozoroski, T. 2002 / AVID Inc., Gelhausen, P., and Roberts, M. 2004), and then overhauled (Ozoroski, T., Hahn, A. 2008). The latest version, MVL-15 has been refined to provide analysis transparency and enhanced to meet the analysis requirements of the Environmentally Responsible Aviation (ERA) Project. Each revision has been implemented with reasonable success. Separate applications of the methodology are in use, including a similar in-house capability, developed by Olson, E. that is tailored for structural and acoustics analyses. A central premise of the methodology is that viscous aerodynamic data can be associated with analytic inviscid aerodynamic results at each spanwise wing section, thereby providing a pathway to map viscous data to the inviscid results. However, a number of factors can sidetrack the analysis consistency during various stages of this process. For example, it should be expected that the final airplane lift curve and drag polar results depend strongly on the geometry and aerodynamics of the airfoil section; however, flap deflections and flap chord extensions change the local reference geometry of the input airfoil, the airplane wing, the tabulated non-dimensional viscous aerodynamics, and the spanwise links between the linear and the viscous aerodynamics. These changes also affect the bound circulation and therefore, calculation and integration of the induced angle of attack and induced drag. MVL-15 is configured to ensure these types of challenges are properly addressed. This report is a comprehensive manual describing the theory, use, and validation of the MVL-15 analysis tool. Section 3 summarizes theoretical, procedural, and characteristic features of MVL-15, and includes a list of the files required to setup, execute, and summarize an analysis. Section 4, Section 5, Section 6, and Section 7 combine to comprise the User's Guide portions of this report. The MVL-15 input and output files are described in Section 4 and Section 5, respectively; the descriptions are supplemented with example files and information about the file formats, parameter definitions, and typical parameter values. Section 6 describes the Wing Geometry Setup Utility and the 2d-Variants Utility files that simplify and assist setting up a consistent set of MVL-15 geometry and aerodynamics input parameters and input files. Section 7 describes the use of the 3d-Results Presentation Utility file that can be used to automatically create summary tables and charts from the MVL-15 output files. Section 8 documents the Validation Results of an extensive and varied validation test matrix, including results of an airplane analysis representative of the ERA Program. A start-to-finish example of the airplane analysis procedure is described in Section 7.

Ozoroski, Thomas A.↗

Geology and landscape evolution of the Hellas region of Mars

Hellas basin on Mars has been the site of volcanism, tectonism, and modification by fluvial, mass-wasting, and eolian processes over its more than 4-b.y. existence. Our detailed geologic mapping and related studies have resulted in the following new interpretations. The asymmetric distribution of highland massifs and other structures that define the uplifted basin rim suggest a formation of the basin by the impact of a low-angle bolide having a trajectory heading S60E. During the Late Noachian, the basin was infilled, perhaps by lava flows, that were sufficiently thick (greater than 1 km) to produce wrinkle ridges on the fill material and extensional faulting along the west rim of the basin. At about the same time, deposits buried northern Malea Planum, which are interpreted to be pyroclastic flows from Amphitrites and Peneus Paterae on the basis of their degraded morphology, topology, and the application of a previous model for pyroclastic volcanism on Mars. Peneus forms a distinctive caldera structure that indicates eruption of massive volumes of magma, whereas Amphitrites is a less distinct circular feature surrounded by a broad, low, dissected shield that suggests generally smaller volume eruptions. During the Early Hesperian, an approximately 1-to 2km-thick sequence of primarily fined-grained, eolian material was deposited on the floor of Hellas basin. Subsequently, the deposit was deeply eroded, except where armored by crater ejecta, and it retreated as much as 200-300 km along its western margin, leaving behind pedestal craters and knobby outliers of the deposit. Local debris flows within the deposit attest to concentrations of groundwater, perhaps in part brought in by outflow floods along the east rim of the basin. These floods may have deposited approximately 100-200m of sediment, subduing wrinkle ridges in the eastern part of the basin floor. During the Late Hesperian and Amazonian, eolian mantles were emplaced on the basin rim and floor and surrounding highlands. Their subsequent erosion resulted in pitted and etched plains and crater fill, irregular mesas, and pedestal craters. Local evidence occurs for the possible former presence of ground ice or ice sheets approximately 100 km across; however, we disagree with a hypothesis that suggest that the entire south rim and much of the floor of Hellas have been glaciated. Orientations of dune fields and yardangs in lower parts of Hellas basin follow directions of the strongest winds predicted by a recently published general circulation model (GCM). Transient frost and dust splotches in the region are, by contrast, related to the GCM prediction for the season in which the images they appear in were taken.

Tanaka, Kenneth L.↗

Geology and landscape evolution of the Hellas region of Mars

Hellas basin on Mars has been the site of volcanism, tectonism, and modification by fluvial, mass-wasting, and eolian processes over its more than 4-b.y. existence. Our detailed geologic mapping and related studies have resulted in the following new interpretations. The asymmetric distribution of highland massifs and other structures that define the uplifted basin rim suggest a formation of the basin by the impact of a low-angle bolide having a trajectory heading S 60 deg E. During the Late Noachian, the basin was infilled, perhaps by lava flows, that were sufficiently thick (>1 km) to produce wrinkle ridges on the fill material and extensional faulting along the west rim of the basin. At about the same time, deposits buried northern Malea Planum, which are interpreted to be pyroclastic flows from Amphitrites and Peneus Paterae on the basis of their degraded morphology, topography, and the application of a previous model for pyroclastic volcanism on Mars. Peneus forms a distinctive caldera structure that indicates eruption of massive volumes of magma, whereas Amphitrites is a less distinct circular feature surrounded by a broad, low, dissected shield that suggests generally smaller volume eruptions. During the Early Hesperian, an approx. 1- to 2-km-thick sequence of primarily fined-grained, eolian material was deposited on the floor of Hellas basin. Subsequently, the deposit was deeply eroded, except where armored by crater ejecta, and it retreated as much as 200-300 km along its western margin, leaving behind pedestal craters and knobby outliers of the deposit. Local debris flows within the deposit attest to concentrations of groundwater, perhaps in part brought in by outflow floods along the east rim of the basin. These floods may have deposited approx. 100-200 m of sediment, subduing wrinkle ridges in the eastern part of the basin floor. During the Late Hesperian and Amazonian, eolian mantles were emplaced on the basin rim and floor and surrounding highlands. Their subsequent erosion resulted in pitted and etched plains and crater fill, irregular mesas, and pedestal craters. Local evidence occurs for the possible former presence of ground ice or ice sheets approx. 100 km across; however, we disagree with a hypothesis that suggests that the entire south rim and much of the floor of Hellas have been glaciated. Orientations of dune fields and yardangs in lower parts of Hellas basin follow directions of the strongest winds predicted by a recently published general circulation model (GCM). Transient frost and dust splotches in the region are, by contrast, related to the GCM prediction for the season in which the images they appear in were taken.

Tanaka, Kenneth L.↗