Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “depth profiling”

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 1,153 records · Page 64

Airborne lidar measurements of El Chichon stratospheric aerosols, January 1984

A lidar-equipped NASA Electra aircraft was flown in January 1984 between the latitude of 38 and 90 deg N. One of the primary purposes of this mission was to determine the spatial distribution and aerosol characteristics of El Chichon produced stratospheric material. Lidar data from that portion of the flight mission between 38 deg N and 77 deg N is presented. Representative profiles of lidar backscatter ratio, a plot of the integral backscattering function versus latitude, and contours of backscatter mixing ratio versus altitude and latitude are given. In addition, tables containing numerical values of the backscatter ratio and backscattering function versus altitude are applied for each profile. These data clearly show that material produced by the El Chichon eruptions of late March-early April 1982 had spread throughout the latitudes covered by this mission, and that the most massive portion of the material resided north of 55 deg N and was concentrated below 17 km in a layer that peaked at 13 to 15 km. In this latitude region, peak backscatter ratios at a wavelength of 0.6943 microns were approximately 3 and the peak integrated backscattering function was about 15 X 10 to the -4/sr corresponding to a peak optical depth of approximately 0.07. This report presents the results of this mission in a ready-to-use format for atmospheric and climatic studies.

Mccormick, M. Patrick↗

Sources, Sinks, and Transatlantic Transport of North African Dust Aerosol: A Multimodel Analysis and Comparison With Remote Sensing Data

This study evaluates model-simulated dust aerosols over North Africa and the North Atlantic from five global models that participated in the Aerosol Comparison between Observations and Models phase II model experiments. The model results are compared with satellite aerosol optical depth (AOD) data from Moderate Resolution Imaging Spectroradiometer (MODIS), Multiangle Imaging Spectroradiometer (MISR), and Sea-viewing Wide Field-of-view Sensor, dust optical depth (DOD) derived from MODIS and MISR, AOD and coarse-mode AOD (as a proxy of DOD) from ground-based Aerosol Robotic Network Sun photometer measurements, and dust vertical distributions/centroid height from Cloud Aerosol Lidar with Orthogonal Polarization and Atmospheric Infrared Sounder satellite AOD retrievals. We examine the following quantities of AOD and DOD: (1) the magnitudes over land and over ocean in our study domain, (2) the longitudinal gradient from the dust source region over North Africa to the western North Atlantic, (3) seasonal variations at different locations, and (4) the dust vertical profile shape and the AOD centroid height (altitude above or below which half of the AOD is located). The different satellite data show consistent features in most of these aspects; however, the models display large diversity in all of them, with significant differences among the models and between models and observations. By examining dust emission, removal, and mass extinction efficiency in the five models, we also find remarkable differences among the models that all contribute to the discrepancies of model-simulated dust amount and distribution. This study highlights the challenges in simulating the dust physical and optical processes, even in the best known dust environment, and stresses the need for observable quantities to constrain the model processes.

Africa↗

Analysis of CO Transport from Australia Fires and Its Impact to Ozone, Water Vapor and Clouds Using CrIS Single-Field-of-View Sounder Products and Other Remote Sensing and MERRA-2 Reanalysis Data

Australia’s unprecedented fire disasters at the end of 2019 to early 2020 emitted huge amounts of carbon monoxide (CO) and fire aerosol particles, and the long-range transport of the smoke plumes poses great impact to climate change, such as altering the Antarctic ozone (O3) hole. Spaceborne observations from satellites have been proven successful in monitoring CO from wildfires in the past decade and can be used to examine the associated CO transport. The primary remote sensing data used in this study is the Single Field of View (SFOV) Sounder Atmospheric Product (SiFSAP) derived from the Cross-track Infrared Sounder (CrIS) on SNPP, a new product that will be generated at NASA and released to the public in a few months. SiFSAP has a high spatial resolution of ~14 km and includes atmospheric profiles of temperature, water vapor, CO, O3, as well as cloud products, and as the SiFSAP algorithm performs both day- and night-time retrievals, it provides an additional observation for night-time in comparison to the day-time-only solar polar satellite instruments such as MOPPIT, OMPS, and TROPOMI. Most recent significant improvements in the quality of this product plus its good spatial resolution (better than most global weather and climate models) enable us to conduct process-oriented analysis of CO transport, and the variation of water vapor, O3 and clouds (optical depth, effective radius and cloud height) within and around the fire plumes. The wind fields from MERRA-2 reanalysis will be used to help characterize the transport. Other remote sensing data, such as CO from AIRS, MOPPIT and TROPOMI, O3 from OMPS and TROPOMI, as well as the UV aerosol index from OMI and OMPS, will be used and compared with SiFSAP sounder products. Cross-validation of CO and O3 from these different sensors with the corresponding MERRA-2 products will also be carried out. Along the plume transport path, some enhancement of O3 and correlation between fire aerosol, cloud optical depth and water vapor have been found but need further investigation. Using the transport of fire plumes over the Southern Pacific Ocean as a unique testbed, this study will demonstrate the trace species monitoring capability of current satellite observations, in particular, for CO and O3, and hopefully, it will provide some insight on the impacts on fire emissions on air quality and climate.

Xiaozhen (Shawn) Xiong↗

ICESat-2 Atmospheric Channel Description, Data Processing and First Results

The Advanced Topographic Laser Altimeter System (ATLAS) was launched aboard the Ice Cloud and land-Elevation Satellite-2 (ICESat-2) satellite in September 2018. ATLAS is a single wavelength (532 nm) lidar system designed to acquire high resolution measurements of the earth's surface while also obtaining atmospheric backscatter from molecules, clouds, and aerosols. Because ATLAS is optimized for altimetry, the atmospheric data acquired is unique in many respects and requires non-standard analysis techniques. For example, the high repetition rate laser limits the vertical extent of the profiles to just 14 km and causes atmospheric scattering from above 15 km to be added to the scattering in the lower 0–14 km profile. In addition, the limited vertical range of the acquired profiles renders it difficult to compute the magnitude of the solar background and hinders the application of standard calibration techniques. Despite these limitations, methods have been developed to successfully produce data products that have value to the atmospheric community for cloud and aerosol research and are currently available at the National Snow and Ice Data Center (NSIDC). In this paper we describe the ICESat-2 atmospheric channel and the methods used to process the ATLAS raw photon count data to obtain calibrated backscatter and higher level products such as layer heights and type, blowing snow, column optical depth and apparent surface reflectance.

Stephen P Palm↗

Assessing Uncertainties of a Geophysical Approach to Estimate Surface Fine Particulate Matter Distributions from Satellite-Observed Aerosol Optical Depth

Health impact analyses are increasingly tapping the broad spatial coverage of satellite aerosol optical depth (AOD) products to estimate human exposure to fine particulate matter (PM2.5). We use a forward geophysical approach to derive ground-level PM2.5 distributions from satellite AOD at 1 km2(exp) resolution for 2011 over the northeastern US by applying relationships between surface PM2.5 and column AOD (calculated offline from speciated mass distributions) from a regional air quality model (CMAQ; 12×12 km2(exp) horizontal resolution). Seasonal average satellite-derived PM2.5 reveals more spatial detail and best captures observed surface PM2.5 levels during summer. At the daily scale, however, satellite-derived PM2.5 is not only subject to measurement uncertainties from satellite instruments, but more importantly to uncertainties in the relationship between surface PM2.5 and column AOD. Using 11 ground-based AOD measurements within 10 km of surface PM2.5 monitors, we show that uncertainties in modeled PM2.5∕AOD can explain more than 70 % of the spatial and temporal variance in the total uncertainty in daily satellite-derived PM2.5 evaluated at PM2.5 monitors. This finding implies that a successful geophysical approach to deriving daily PM2.5 from satellite AOD requires model skill at capturing day-to-day variations in PM2.5∕AOD relationships. Overall, we estimate that uncertainties in the modeled PM2.5∕AOD lead to an error of 11 µg m−3(exp) in daily satellite-derived PM2.5, and uncertainties in satellite AOD lead to an error of 8 µg m−3(exp). Using multi-platform ground, airborne, and radiosonde measurements, we show that uncertainties of modeled PM2.5∕AOD are mainly driven by model uncertainties in aerosol column mass and speciation, while model representation of relative humidity and aerosol vertical profile shape contributes some systematic biases. The parameterization of aerosol optical properties, which determines the mass extinction efficiency, also contributes to random uncertainty, with the size distribution being the largest source of uncertainty and hygroscopicity of inorganic salt the second largest. Future efforts to reduce uncertainty in geophysical approaches to derive surface PM2.5 from satellite AOD would thus benefit from improving model representation of aerosol vertical distribution and aerosol optical properties, to narrow uncertainty in satellite-derived PM2.5.

aerosol optical depth↗

Effect of Smoke on Cloud Formation during the Biomass Burning Season over the Amazon Basin

Aerosol absorption of sunlight reduces surface irradiation and heats the aerosol layer. The consequent changes in the temperature and humidity profiles can affect cloud formation extent and life time, which is called the semi-direct effect. We evaluate this aerosol semi-direct effect using data collected during the 2002 biomass burning season over the Amazon basin from the MODIS instrument on the Aqua satellite. MODIS measures the cloud coverage and the aerosol optical thickness among the clouds. We found that the radiative heating of the atmosphere and cooling of the surface due to the presence of the smoke decreases the cloud coverage. A very clear negative correlation emerges between the cloud fraction and the smoke optical depth. The results are compared to calculations using 1-D radiation model (M.D. Chou), and used to calculate this regional semi direct effect on climate forcing.

Koren, I.↗

Connecting the Interstellar Gas and Dust Properties in Distant Galaxies Using Quasar Absorption Systems

Gas and dust grains are fundamental components of the interstellar medium and significantly impact many of the physical processes driving galaxy evolution, such as star-formation, and the heating, cooling, and ionization of the interstellar material. Quasar absorption systems (QASs), which trace intervening galaxies along the sightlines to luminous quasars, provide a valuable tool to directly study the properties of the interstellar gas and dust in distant, normal galaxies. We have established the presence of silicate dust grains in at least some gas-rich QASs, and find that they exist at higher optical depths than expected for diffuse gas in the Milky Way. Differences in the absorption feature shapes additionally suggest variations in the silicate dust grain properties, such as in the level of grain crystallinity, from system-to-system. We present results from a study of the gas and dust properties of QASs with adequate archival IR data to probe the silicate dust grain properties. We discuss our measurements of the strengths of the 10 and 18 micron silicate dust absorption features in the QASs, and constraints on the grain properties (e.g., composition, shape, crystallinity) based on fitted silicate profile templates. We investigate correlations between silicate dust abundance, reddening, and gas metallicity, which will yield valuable insights into the history of star formation and chemical enrichment in galaxies.

Aller, Monique C.↗

The chromospheres of late-type stars. I - Eridani as a test case of multiline modelling

A new model of the lower chromosphere of the dwarf K2 star Epsilon Eridani is derived by matching flux profiles of the Ca IR triplet lines 8498 and 8542 A H-alpha and H-beta lines and the Na D lines (all observed simultaneously at the AAT), and the Ca II K line. The coupled non-LTE equations of statistical equilibrium and radiative transfer are solved under the constraint of hydrostatic equilibrium using the Carlsson (1986) code. Within the framework of the model, the Na D lines are an important photospheric diagnostic, and the Ca IR triplet lines can be used to locate the temperature minimum. The computed H-alpha and H-beta depths are highly sensitive constraints on the transition zone gradients and base pressures allowing us to derive a pressure at the base of the transition zone of 0.9 dyn/cm.

Thatcher, John D.↗

NiH2 Cycle Life Study

Cycle life studies have been performed at Eagle Picher Technologies (EPT), on HST Mantech design cells with various pedigrees of slurry and dry sinter processed electrodes, to evaluate peak load voltage performance during generic load profile testing. These tests provide information for determining voltage and capacity fade (degradation) mechanisms, and their impact on nickel hydrogen cell cycle life. Comparison of peak load voltage fade, as a function of State of Charge and cycle life, with capacity data from HST indicates that the cycle life limiting mechanism is due to impedance growth, and formation of a second discharge plateau. With a second plateau on discharge, capacity from the cell is still available, but at an unacceptable low voltage of 0.8 V per cell (17.6 V battery). Data shows that cell impedance increases with cycle number and depth of discharge, as expected.

Hollandsworth, Roger P.↗

Applications of ArcticDEM for measuring volcanic dynamics, landslides, retrogressive thaw slumps, snowdrifts, and vegetation heights

Topographical changes are of fundamental interest to a wide range of Arctic science disciplines faced with the need to anticipate, monitor, and respond to the effects of climate change, including geohazard management, glaciology, hydrology, permafrost, and ecology. This study demonstrates several geomorphological, cryospheric, and biophysical applications of ArcticDEM – a large collection of publicly available, time-dependent digital elevation models (DEMs) of the Arctic. Our study illustrates ArcticDEM's applicability across different disciplines and five orders of magnitude of elevation derivatives, including measuring volcanic lava flows, ice cauldrons, post-failure landslides, retrogressive thaw slumps, snowdrifts, and tundra vegetation heights. We quantified surface elevation changes in different geological settings and conditions using the time series of ArcticDEM. Following the 2014–2015 Bárðarbunga eruption in Iceland, ArcticDEM analysis mapped the lava flow field, and revealed the post-eruptive ice flows and ice cauldron dynamics. The total dense-rock equivalent (DRE) volume of lava flows is estimated to be (1431 ± 2) million m 3 . Then, we present the aftermath of a landslide in Kinnikinnick, Alaska, yielding a total landslide volume of (400 ± 8) × 103 m 3 and a total area of 0.025 km 2 . ArcticDEM is further proven useful for studying retrogressive thaw slumps (RTS). The ArcticDEM-mapped RTS profile is validated by ICESat-2 and drone photogrammetry resulting in a standard deviation of 0.5 m. Volume estimates for lake-side and hillslope RTSs range between 40,000 ± 9000 m 3 and 1,160,000 ± 85,000 m 3 , highlighting applicability across a range of RTS magnitudes. A case study for mapping tundra snow demonstrates ArcticDEM's potential for identifying high-accumulation, late-lying snow areas. The approach proves effective in quantifying relative snow accumulation rather than absolute values (standard deviation of 0.25 m, bias of −0.41 m, and a correlation coefficient of 0.69 with snow depth estimated by unmanned aerial systems photogrammetry). Furthermore, ArcticDEM data show its feasibility for estimating tundra vegetation heights with a standard deviation of 0.3 m (no bias) and a correlation up to 0.8 compared to the light detection and ranging (LiDAR). The demonstrated capabilities of ArcticDEM will pave the way for the broad and pan-Arctic use of this new data source for many disciplines, especially when combined with other imagery products. The wide range of signals embedded in ArcticDEM underscores the potential challenges in deciphering signals in regions affected by various geological processes and environmental influences.

Chunli Dai↗

The effects of artificially-produced defects on the film thickness distribution in sliding EHD point contacts

Effects of artificially produced dents and grooves on the elastohydrodynamic film thickness profile in a sliding point contact are investigated using optical interferometry. Defects on the surface of a highly polished ball are held stationary at various locations within the vicinity of the contact region while the disk is rotating. It is shown that defects can dramatically change the film thickness which exists when no defects are present in or near the contact, and it is shown that the change in thickness is mainly a function of the position of the defects in the inlet region, the geometry of the defects, the orientation of the defects in the case of grooves, and the depth of the defect relative to the central film thickness.

Cusano, C.↗

HST-STIS Spectra of Saturn's Rings and Implications for Their Reddening Agent

We obtained HST-STIS spectra of Saturn's main rings in May 2011, using the G230L (and G430L) gratings, with final averaged radial resolution of 160 (and 330) km/pixel. The dataset filled a previous 200-330nm "spectral gap" between Cassini and ground-based spectra. The data provide radial profiles as a function of wavelength, but our most basic product at this point is a set of very low-noise spectra, radially averaged over broad regions of the rings (A, B, C, and Cassini Division). The raw spectra required special processing to remove artifacts due to extended-source grating scatter. We have modeled the spectra using a new particle surface model, which corrects for on-surface shadowing due to the likely very rough ring particle surfaces, and avoids overestimation of intra-mixed "neutral absorber". We correct for non-classical layer effects and finite ring optical depth, and relate our observed reflectivities to the spherical albedos of individual smooth particles. We model these smooth particle albedos using standard Hapke theory for regolith grain mixtures that are either homogeneous and "intramixed" (nonicy absorbers dispersed in water ice regolith grains) or heterogeneous "intimate" mixtures. As candidates for the nonicy contaminants we have considered amorphous carbon, aromatic-rich and aliphatic-rich organic tholins, silicates, hematite and iron metal. For the A and B rings, we find that iron metal (including a new theoretical estimate of the refractive indices of nanometer-sized grains of iron) is not spectrally steep enough in the 200-300nm range, and that aliphatic-rich tholins are either too steep at short wavelengths or too flat at long wavelengths. However, less than 1% by mass of aromatic-rich tholins provides a very good fit across the entire spectral range with no gratuitous "neutral absorber" needed, and a minimum of additional free parameters. The best fits require forward-scattering regolith grains. For the C Ring and Cassini Division, additional absorbers are needed (updated results will be given).

Saturn's Rings↗

Remote sensing of snow properties by passive microwave radiometry: GSFC truck experiment

Recent results indicate that microwave radiometry has the potential for inferring the snow depth and water equivalent information from snowpacks. In order to assess this potential for determining the water equivalent of a snowpack, it is necessary to understand the microwave emission and scattering behavior of the snow at various wavelengths under carefully controlled conditions. Truck-mounted microwave instrumentation was used to study the microwave characteristics of the snowpack in the Colorado Rocky Mountain region during the winters of 1977 to 78 and 7978 to 79. The spectral signatures of C, X, K sub u, and K sub a band radiometers with dual polarization were used, together with measurements of snowpack density, temperature an ram profiles, liquid water content, and rough characterization of the crystal sizes. These data compared favorably with calculated results based on recent microscopic scattering models.

Chang, A. T. C.↗

Airborne Sunphotometer, Airborne in-situ, Space-borne, and Ground-Based Measurements of Troposoheric Aerosol in Ace-2

We report on clear-sky column closure experiments performed in the Canary Islands during the second Aerosol Characterization Experiment (ACE-2) in June/July 1997. We present results obtained by combining airborne sunphotometer and in-situ aerosol measurements taken aboard the Pelican aircraft, space-borne NOAA/AVHRR data and ground-based lidars A wide range of aerosol types was encountered throughout the ACE-2 area, including background Atlantic marine, European pollution-derived, and African mineral dust. During !he two days discussed here, vertical profiles flown in cloud free air masses revealed three distinctly different layers: a marine boundary layer (MBL) with varying pollution levels, an elevated dust layer, and a very clean layer between the MBL and the dust layer. We found that the presence of the elevated dust layer removes the good agreement between satellite and sunphotometer AOD usually found in the absence of the dust layer. Using size-resolved composition information we have computed optical properties of the ambient aerosol from the in-situ measurements and subsequently compared those to the sunphotometer results. In the dust, the agreement in layer aerosol optical depth (380-1060 nm) is 3-8%. In the MBL there is tendency for the in-situ results to be slightly lower than the sunphotometer measurements (10-17% at 525 nm), but these differences are within the combined error bars of the measurements and computations.

Schmid, Beat↗

Lunar seismic profiling experiment natural activity study

The Lunar Seismic Experiment Natural Activity Study has provided a unique opportunity to study the high frequency (4-20 Hz) portion to the seismic spectrum on the moon. The data obtained from the LSPE was studied to evaluate the origin and importance of the process that generates thermal moonquakes and the characteristics of the seismic scattering zone at the lunar surface. The detection of thermal moonquakes by the LSPE array made it possible to locate the sources of many events and determine that they are definitely not generated by astronaut activities but are the result of a natural process on the moon. The propagation of seismic waves in the near-surface layers was studied in a qualitative manner. In the absence of an adequate theoretical model for the propagation of seismic waves in the moon, it is not possible to assign a depth for the scattering layer. The LSPE data does define several parameters which must be satisfied by any model developed in the future.

Duennebier, F. K.↗

Numerical simulations of galactic wakes

We have performed a series of numerical simulations which examine the simple but realistic case of a galaxy with a smooth potential corresponding to a King profile, and with no interstellar medium. This corresponds to the case of an early-type cluster galaxy that has previously been stripped of its gas. The simulations provide a numerical basis for future work involving galaxies with interstellar media. We use a time-explicit finite difference code (Lufkin & Hawley 1992, ApJ submitted) to obtain solutions to the nonlinear equations for adiabatic flow in axisymmetry. The computational grid covers a region extending 500 kpc from a stationary galaxy radially and along the symmetry axis. We use a graded mesh, allowing for full resolution near the galaxy. Because the assumed potential has a finite depth, no artificial inner boundary condition is necessary; reflection symmetry at the axis is assumed. With a total grid size of 64 radial x 128 vertical, each run requires approximately 5 minutes of cpu time on the NCSA Cray-2.

Lufkin, Eric A.↗

Infrared spectroscopy of Seyfert 2 galaxies: A look through the obscuring Torus?

We present both high-resolution (R = 1260) and low-resolution (R = 345 and 425) J-band spectra of a sample of 15 Seyfert 2 galaxies. Our goal is to look for broad Pa beta lines, indicating broad-line regions which are hidden by dust from our view at optical wavelengths. Of the 15 objects studied here, three have broad Pa beta lines: MCG-05.23.16, Mrk 463E, and NGC 2992. Mrk 176 and NGC 5728 may also have weak broad lines. In NGC 5506, previously reported to have broad Pa beta and hydrogen alpha lines, we find that the Pa beta line profile is continuous and has the same shape as the nearby line (Fe II) lambda 1.2567, which should not have a broad component. We interpret these observations as gas from the narrow-line region (NLR) with no broad component. In NGC 5506, however, the NLR profiles become broader with increasing wavelength, indicating that highly reddened wings are becoming more readily visible at the longer wavelengths. We confirm the correlation of (O I) lambda 6300/hydrogen alpha and (Fe II) lambda 1.644/Br gamma (the latter transformed to (Fe II) lambda 1.2567/Pa beta to compare with our data) found by previous authors when comparing active galactic nuclei (AGNs), supernova remnants, starbursts, and H II regions. The correlation confirms that in all of these objects both (O I) lambda 6300 and the (Fe II) lines come from partially ionized regions in which hydrogen is mostly neutral. Comparison of the infrared optical depths with column depths determined from X-ray data show a general tendency for the objects with detected broad Pa beta to have lower X-ray columns.

Goodrich, Robert W.↗

Evaluating Mars Science Laboratory Landing Sites with the Mars Global Reference Atmospheric Model (Mars-GRAM 2005)

The Mars Global Reference Atmospheric Model (Mars-GRAM) is an engineering-level atmospheric model widely used for diverse mission applications. Mars-GRAM s perturbation modeling capability is commonly used, in a Monte-Carlo mode, to perform high fidelity engineering end-to-end simulations for entry, descent, and landing (EDL) [1]. From the surface to 80 km altitude, Mars-GRAM is based on the NASA Ames Mars General Circulation Model (MGCM). Mars-GRAM and MGCM use surface topography from Mars Global Surveyor Mars Orbiter Laser Altimeter (MOLA), with altitudes referenced to the MOLA areoid, or constant potential surface. Traditional Mars-GRAM options for representing the mean atmosphere along entry corridors include: (1) Thermal Emission Spectrometer (TES) mapping years 1 and 2, with Mars-GRAM data coming from NASA Ames Mars General Circulation Model (MGCM) results driven by observed TES dust optical depth or (2) TES mapping year 0, with user-controlled dust optical depth and Mars-GRAM data interpolated from MGCM model results driven by selected values of globally-uniform dust optical depth. Mars-GRAM 2005 has been validated [2] against Radio Science data, and both nadir and limb data from TES [3]. There are several new features included in Mars-GRAM 2005. The first is the option to use input data sets from MGCM model runs that were designed to closely simulate conditions observed during the first two years of TES observations at Mars. The TES Year 1 option includes values from April 1999 through January 2001. The TES Year 2 option includes values from February 2001 through December 2002. The second new feature is the option to read and use any auxiliary profile of temperature and density versus altitude. In exercising the auxiliary profile Mars-GRAM option, values from the auxiliary profile replace data from the original MGCM databases. Some examples of auxiliary profiles include data from TES nadir or limb observations and Mars mesoscale model output at a particular location and time. The final new feature is the addition of two Mars-GRAM parameters that allow standard deviations of Mars-GRAM perturbations to be adjusted. The parameter rpscale can be used to scale density perturbations up or down while rwscale can be used to scale wind perturbations.

Justh, H. L.↗