Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bias correction”

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 361 records · Page 20

Improved Satellite Retrievals of NO2 and SO2 over the Canadian Oil Sands and Comparisons with Surface Measurements

Satellite remote sensing is increasingly being used to monitor air quality over localized sources such as the Canadian oil sands. Following an initial study, significantly low biases have been identified in current NO2 and SO2 retrieval products from the Ozone Monitoring Instrument (OMI) satellite sensor over this location resulting from a combination of its rapid development and small spatial scale. Air mass factors (AMFs) used to convert line-of-sight "slant" columns to vertical columns were re-calculated for this region based on updated and higher resolution input information including absorber profiles from a regional-scale (15 km × 15 km resolution) air quality model, higher spatial and temporal resolution surface reflectivity, and an improved treatment of snow. The overall impact of these new Environment Canada (EC) AMFs led to substantial increases in the peak NO2 and SO2 average vertical column density (VCD), occurring over an area of intensive surface mining, by factors of 2 and 1.4, respectively, relative to estimates made with previous AMFs. Comparisons are made with long-term averages of NO2 and SO2 (2005-2011) from in situ surface monitors by using the air quality model to map the OMI VCDs to surface concentrations. This new OMI-EC product is able to capture the spatial distribution of the in situ instruments (slopes of 0.65 to 1.0, correlation coefficients of greater than 0.9). The concentration absolute values from surface network observations were in reasonable agreement, with OMI-EC NO2 and SO2 biased low by roughly 30%. Several complications were addressed including correction for the interference effect in the surface NO2 instruments and smoothing and clear-sky biases in the OMI measurements. Overall these results highlight the importance of using input information that accounts for the spatial and temporal variability of the location of interest when performing retrievals.

surface concentration↗

Status of Middle Atmosphere-Climate Models: Results SPARC-GRIPS

The middle atmosphere is an important component of the climate system, primarily because of the radiative forcing of ozone. Middle atmospheric ozone can change, over long times, because of changes in the abundance of anthropogenic pollutants which catalytically destroy it, and because of the temperature sensitivity of kinetic reaction rates. There is thus a complex interaction between ozone, involving chemical and climatic mechanisms. One question of interest is how ozone will change over the next decades , as the "greenhouse-gas cooling" of the middle atmosphere increases but the concentrations of chlorine species decreases (because of policy changes). concerns the climate biases in current middle atmosphere-climate models, especially their ability to simulate the correct seasonal cycle at high latitudes, and the existence of temperature biases in the global mean. A major obstacle when addressing this question This paper will present a summary of recent results from the "GCM-Reality Intercomparison Project for SPARC" (GRIPS) initiative. A set of middle atmosphere-climate models has been compared, identifying common biases. Mechanisms for these biases are being studied in some detail, including off-line assessments of the radiation transfer codes and coordinated studies of the impacts of gravity wave drag due to sub-grid-scale processes. ensemble of models will be presented, along with numerical experiments undertaken with one or more models, designed to investigate the mechanisms at work in the atmosphere. The discussion will focus on dynamical and radiative mechanisms in the current climate, but implications for coupled ozone chemistry and the future climate will be assessed.

Pawson, Steven↗

ICESAT GLAS Altimetry Measurements: Received Signal Dynamic Range and Saturation Correction

NASAs Ice, Cloud, and land Elevation Satellite (ICESat), which operated between 2003 and 2009, made the first satellite-based global lidar measurement of earths ice sheet elevations, sea-ice thickness, and vegetation canopy structure. The primary instrument on ICESat was the Geoscience Laser Altimeter System (GLAS), which measured the distance from the spacecraft to the earth's surface via the roundtrip travel time of individual laser pulses. GLAS utilized pulsed lasers and a direct detection receiver consisting of a silicon avalanche photodiode and a waveform digitizer. Early in the mission, the peak power of the received signal from snow and ice surfaces was found to span a wider dynamic range than anticipated, often exceeding the linear dynamic range of the GLAS 1064-nm detector assembly. The resulting saturation of the receiver distorted the recorded signal and resulted in range biases as large as approximately 50 cm for ice- and snow-covered surfaces. We developed a correction for this saturation range bias based on laboratory tests using a spare flight detector, and refined the correction by comparing GLAS elevation estimates with those derived from Global Positioning System surveys over the calibration site at the salar de Uyuni, Bolivia. Applying the saturation correction largely eliminated the range bias due to receiver saturation for affected ICESat measurements over Uyuni and significantly reduced the discrepancies at orbit crossovers located on flat regions of the Antarctic ice sheet.

Sun, Xiaoli↗

An Analysis of Degradation in Low-Cost Particulate Matter Sensors

Low-cost sensors (LCSs) are increasingly being used to measure fine particulate matter (PM 2.5 ) concentrations in cities around the world. One of the most commonly deployed LCSs is PurpleAir with ∼15 000 sensors deployed in the United States, alone. PurpleAir measurements are widely used by the public to evaluate PM 2.5 levels in their neighborhoods. PurpleAir measurements are also increasingly being integrated into models by researchers to develop large-scale estimates of PM 2.5 . However, the change in sensor performance over time has not been well studied. It is important to understand the lifespan of these sensors to determine when they should be serviced or replaced, and when measurements from these devices should or should not be used for various applications. This paper fills this gap by leveraging the fact that: (1) each PurpleAir sensor is composed of two identical sensors and the divergence between their measurements can be observed, and (2) there are numerous PurpleAir sensors within 50 meters of regulatory monitors allowing for the comparison of measurements between these instruments. We propose empirically-derived degradation outcomes for the PurpleAir sensors and evaluate how these outcomes change over time. On average, we find that the number of ‘flagged’ measurements, where the two sensors within each PurpleAir sensor disagree, increases with time to ∼4% after 4 years of operation. Approximately, 2 percent of all PurpleAir sensors were permanently degraded. The largest fraction of permanently degraded PurpleAir sensors appeared to be in the hot and humid climate zone, suggesting that sensors in these locations may need to be replaced more frequently. We also find that the bias of PurpleAir sensors, or the difference between corrected PM2.5 levels and the corresponding reference measurements, changed over time by −0.12 μg m −3 (95% CI: −0.13 μg m −3 , −0.10 μg m −3 ) per year. The average bias increases dramatically after 3.5 years. Further, climate zone is a significant modifier of the association between degradation outcomes and time.

Priyanka deSouza↗

Soil Moisture Estimation in South Asia via Assimilation of SMAP Retrievals

A soil moisture retrieval assimilation framework is implemented across South Asia in an attempt to improve regional soil moisture estimation as well as to provide a consistent regional soil moisture dataset. This study aims to improve the spatiotemporal variability of soil moisture estimates by assimilating Soil Moisture Active Passive (SMAP) near-surface soil moisture retrievals into a land surface model. The Noah-MP (v4.0.1) land surface model is run within the NASA Land Information System software framework to model regional land surface processes. NASA Modern-Era Retrospective Analysis for Research and Applications (MERRA2) and Global Precipitation Measurement (GPM) Integrated Multi-satellitE Retrievals (IMERG) provide the meteorological boundary conditions to the land surface model. Assimilation is carried out using both cumulative distribution function (CDF)-corrected (DA-CDF) and uncorrected SMAP retrievals (DA-NoCDF). CDF matching is applied to correct the statistical moments of the SMAP soil moisture retrieval relative to the land surface model. Comparison of assimilated and model-only soil moisture estimates with publicly available in situ measurements highlights the relative improvement in soil moisture estimates by assimilating SMAP retrievals. Across the Tibetan Plateau, DA-NoCDF reduced the mean bias and RMSE by 8.4 % and 9.4 %, even though assimilation only occurred during less than 10 % of the study period due to frozen (or partially frozen) soil conditions. The best goodness-of-fit statistics were achieved for the IMERG DA-NoCDF soil moisture experiment. The general lack of publicly available in situ measurements across irrigated areas limited a domain-wide direct model validation. However, comparison with regional irrigation patterns suggested correction of biases associated with an unmodeled hydrologic phenomenon (i.e., anthropogenic influence via irrigation) as a result of SMAP soil moisture retrieval assimilation. The greatest sensitivity to assimilation was observed in cropland areas. Improvements in soil moisture potentially translate into improved spatiotemporal patterns of modeled evapotranspiration, although limited influence from soil moisture assimilation was observed on modeled processes within the carbon cycle such as gross primary production. Improvement in fine-scale modeled estimates by assimilating coarse-scale retrievals highlights the potential of this approach for soil moisture estimation over data-scarce regions.

Jawairia Ahmad↗

Drag bias feedback for the analytic drag control entry guidance system

The Analytic Drag Control (ADC) entry guidance has been developed and baselined for the space shuttle orbiter entry. A method is presented which corrects the orbiter entry guidance commanded bank angle for biases between navigated drag and guidance computed reference drag. This is accomplished by an integral feedback technique, which uses the drag bias information to adjust the difference between navigated and reference altitude rate used by the ADC guidance. The method improves the capability of the ADC guidance system by compensating for any error source which causes a bias between the navigated drag and reference drag profile. These errors include navigated altitude rate errors, atmosphere dispersions, and roll attitude deadband effects. A discussion of the method and results of digital computer entry simulations is presented.

Kyle, H. C.↗

GMI High Frequency Antenna Pattern Correction Update Based on GPM Inertial Hold and Comparison with ATMS

In an inertial hold, the spacecraft does not attempt to maintain geodetic pointing, but rather maintains the same inertial position throughout the orbit. The result is that the spacecraft appears to pitch from 0 to 360 degrees around the orbit. Two inertial holds were performed with the GPM spacecraft: 1) May 20, 2014 16:48:31 UTC-18:21:04 UTC, spacecraft flying forward +X (0yaw), pitch from 55 degrees (FCS) to 415 degrees (FCS) over the orbit2) Dec 9, 2014 01:30:00 UTC-03:02:32 UTC, spacecraft flying backward X (180yaw), pitch from 0 degrees (FCS) to 360 degrees (FCS) over the orbitThe inertial hold affords a view of the earth through the antenna backlobe. The antenna spillover correction may be evaluated based on the inertial hold data.The current antenna pattern correction does not correct for spillover in the 166 and 183 GHz channels. The two inertial holds both demonstrate that there is significant spillover from the 166 and 183 GHz channels. By not correcting the spillover, the 166 and 183 GHz channels are biased low by about 1.8 to 3K. We propose to update the GMI calibration algorithm with the spill-over correction presented in this document for 166 GHz and 183 GHz.

calibration↗

A New Method to Correct the Electrochemical Concentration Cell (ECC) Ozonesonde Time Response and its Implications for "Background Current" and Pump Efficiency

The electrochemical concentration cell (ECC) ozonesonde has been the main instrument for in situ profiling of ozone worldwide; yet, some details of its operation, which contribute to the ozone uncertainty budget, are not well understood. Here, we investigate the time response of the chemical reactions inside the ECC and how corrections can be used to remove some systematic biases. The analysis is based on the understanding that two reaction pathways involving ozone occur inside the ECC that generate electrical currents on two very different timescales. The main fast-reaction pathway with a time constant of about 20 s is due the conversion of iodide to molecular iodine and the generation of two free electrons per ozone molecule. A secondary slow-reaction pathway involving the buffer generates an excess current of about 2 %–10 % with a time constant of about 25 min. This excess current can be interpreted as what has conventionally been considered the “background current”. This contribution can be calculated and removed from the measured current instead of the background current. Here we provide an algorithm to calculate and remove the contribution of the slow-reaction pathway and to correct for the time lag of the fast-reaction pathway. This processing algorithm has been applied to ozonesonde profiles at Costa Rica and during the Central Equatorial Pacific Experiment (CEPEX) as well as to laboratory experiments evaluating the performance of ECC ozonesondes. At Costa Rica, where a 1 % KI, 1/10th buffer solution is used, there is no change in the derived total ozone column; however, in the upper troposphere and lower stratosphere, average reported ozone concentrations increase by up to 7 % and above 30 km decrease by up to 7 %. During CEPEX, where a 1 % KI, full-buffer solution was used, ozone concentrations are increased mostly in the upper troposphere, with no change near the top of the profile. In the laboratory measurements, the processing algorithms have been applied to measurements using the majority of current sensing solutions and using only the stronger pump efficiency correction reported by Johnson et al. (2002). This improves the accuracy of the ECC sonde ozone profiles, especially for low ozone concentrations or large ozone gradients and removes systematic biases relative to the reference instruments. In the surface layer, operational procedures prior to launch, in particular the use of filters, influence how typical gradients above the surface are detected. The correction algorithm may report gradients that are steeper than originally reported, but their uncertainty is strongly influenced by the prelaunch procedures.

electrochemical concentration cell (ECC) ozonesond↗

Shortcomings with Tree-Structured Edge Encodings for Neural Networks

In evolutionary algorithms a common method for encoding neural networks is to use a tree structured assembly procedure for constructing them. Since node operators have difficulties in specifying edge weights and these operators are execution-order dependent, an alternative is to use edge operators. Here we identify three problems with edge operators: in the initialization phase most randomly created genotypes produce an incorrect number of inputs and outputs; variation operators can easily change the number of input/output (I/O) units; and units have a connectivity bias based on their order of creation. Instead of creating I/O nodes as part of the construction process we propose using parameterized operators to connect to preexisting I/O units. Results from experiments show that these parameterized operators greatly improve the probability of creating and maintaining networks with the correct number of I/O units, remove the connectivity bias with I/O units and produce better controllers for a goal-scoring task.

Hornby, Gregory S.↗

Statistics of Delta v magnitude for a trajectory correction maneuver containing deterministic and random components

A number of interplanetary missions now being planned involve placing deterministic maneuvers along the flight path to alter the trajectory. Lee and Boain (1973) examined the statistics of trajectory correction maneuver (TCM) magnitude with no deterministic ('bias') component. The Delta v vector magnitude statistics were generated for several values of random Delta v standard deviations using expansions in terms of infinite hypergeometric series. The present investigation uses a different technique (Monte Carlo simulation) to generate Delta v magnitude statistics for a wider selection of random Delta v standard deviations and also extends the analysis to the case of nonzero deterministic Delta v's. These Delta v magnitude statistics are plotted parametrically. The plots are useful in assisting the analyst in quickly answering questions about the statistics of Delta v magnitude for single TCM's consisting of both a deterministic and a random component. The plots provide quick insight into the nature of the Delta v magnitude distribution for the TCM.

Bollman, W. E.↗

Comparison and testing of extended Kalman filters for attitude estimation of the Earth Radiation Budget Satellite

This paper presents the testing and comparison of two extended Kalman filters (EKFs) developed for the Earth Radiation Budget Satellite (ERBS). One EKF updates the attitude quaternion using a four-component additive error quaternion. This technique is compared to that of a second EKF, which uses a multiplicative error quaternion. Brief developments of the algorithms are included. The convergence properties, singularity problems, and normalization techniques of the two filters are addressed. Both filters are tested with simulated ERBS sensor data in addition to real ERBS sensor data. The results of the two filters are also compared to those from the ERBS operational ground support software, which uses a batch differential correction algorithm to estimate attitude and gyro biases.

Bar-Itzhack, Itzhack Y.↗

Comparison and testing of extended Kalman filters for attitude estimation of the Earth radiation budget satellite

The testing and comparison of two Extended Kalman Filters (EKFs) developed for the Earth Radiation Budget Satellite (ERBS) is described. One EKF updates the attitude quaternion using a four component additive error quaternion. This technique is compared to that of a second EKF, which uses a multiplicative error quaternion. A brief development of the multiplicative algorithm is included. The mathematical development of the additive EKF was presented in the 1989 Flight Mechanics/Estimation Theory Symposium along with some preliminary testing results using real spacecraft data. A summary of the additive EKF algorithm is included. The convergence properties, singularity problems, and normalization techniques of the two filters are addressed. Both filters are also compared to those from the ERBS operational ground support software, which uses a batch differential correction algorithm to estimate attitude and gyro biases. Sensitivity studies are performed on the estimation of sensor calibration states. The potential application of the EKF for real time and non-real time ground attitude determination and sensor calibration for future missions such as the Gamma Ray Observatory (GRO) and the Small Explorer Mission (SMEX) is also presented.

Deutschmann, Julie↗

A new look at Wolf sunspot numbers in the late 1700's

A procedure is described for the derivation of the Wolf sunspot number from the number of sunspot groups, and the approach is used to obtain a 'Group Wolf number'. It is shown that this technique is superior to the classical method of determining the Wolf number, because corrections for differences due to the observers' bias are reduced and self-consistent long-term time series can be developed; the level of solar activity can be calculated with an accuracy of +/- 5 percent. The procedure was used to determine Wolf sunspot numbers for the solar cycles 1, 2, and 3 (1761-1777). It is found that the standard Wolf numbers are nearly homogeneous with sunspot numbers measured from 1875 to 1976, but the peak of solar cycle 2 is too low by 30 percent.

Hoyt, Douglas V.↗

A deep x ray survey in the Lockman hole and the soft x ray log N - log S

The longest pointed observation (152 ksec) with the Rosat position sensitive proportional counter in the direction of the absolutely lowest neutral hydrogen column density is discussed. 26 shallower fields from the Rosat medium sensivity survey are analyzed. 1176 X-ray sources were detected in at least one Rosat energy band in these fields covering a total solid angle of 9.3 deg; 661 of these sources constitute a statistically complete sample detected in the hard band with 0.5 to 2 keV fluxes greater than 2.5 times 10 to the power of minus 15 erg/sq cm s. The faintest limiting flux of the survey is analyzed. Detailed simulations show that confusion effects and other selection biases are relatively small and can be corrected for in the sample. From an analysis in the Lockman field, a best fit slope of approximately 1.8 is found for the extrapolation of the differential X-ray counts below 2.5 times 10 to the power of minus 15 erg/sq cm s. On the basis of this analysis an upper limit of approximately 25% can be found for a truly diffuse background component in the Rosat hard band.

Hasinger, G.↗

Estimates of Ground Temperature and Atmospheric Moisture from CERES Observations

A method is developed to retrieve surface ground temperature (Tg) and atmospheric moisture using clear sky fluxes (CSF) from CERES-TRMM observations. In general, the clear sky outgoing long-wave radiation (CLR) is sensitive to upper level moisture (q(sub h)) over wet regions and Tg over dry regions The clear sky window flux from 800 to 1200 /cm (RadWn) is sensitive to low level moisture (q(sub j)) and Tg. Combining these two measurements (CLR and RadWn), Tg and q(sub h) can be estimated over land, while q(sub h) and q(sub t) can be estimated over the oceans. The approach capitalizes on the availability of satellite estimates of CLR and RadWn and other auxiliary satellite data. The basic methodology employs off-line forward radiative transfer calculations to generate synthetic CSF data from two different global 4-dimensional data assimilation products. Simple linear regression is used to relate discrepancies in CSF to discrepancies in Tg, q(sub h) and q(sub t). The slopes of the regression lines define sensitivity parameters that can be exploited to help interpret mismatches between satellite observations and model-based estimates of CSF. For illustration, we analyze the discrepancies in the CSF between an early implementation of the Goddard Earth Observing System Data Assimilation System (GEOS-DAS) and a recent operational version of the European Center for Medium-Range Weather Prediction data assimilation system. In particular, our analysis of synthetic total and window region SCF differences (computed from two different assimilated data sets) shows that simple linear regression employing (Delta)Tg and broad layer (Delta)q(sub l) from 500 hPa to surface and (Delta)q(sub h) from 200 to 500 hPa provides a good approximation to the full radiative transfer calculations, typically explaining more than 90% of the 6-hourly variance in the flux differences. These simple regression relations can be inverted to "retrieve" the errors in the geophysical parameters. Uncertainties (normalized by standard deviation) in the monthly mean retrieved parameters range from 7% for (Delta)T to about 20% for (Delta)q(sub t). Our initial application of the methodology employed an early CERES-TRMM data set (CLR and Radwn) to assess the quality of the GEOS2 data. The results showed that over the tropical and subtropical oceans GEOS2 is, in general, too wet in the upper troposphere (mean bias of 0.99 mm) and too dry in the lower troposphere (mean bias of -4.7 mm). We note that these errors, as well as a cold bias in the Tg, have largely been corrected in the current version of GEOS-2 with the introduction of a land surface model, a moist turbulence scheme and the assimilation of SSTM/I total precipitable water.

Wu, Man Li C.↗

High Frequency, Long Time Series Measurements from the Bermuda Testbed Mooring in Support of SIMBIOS

It has been recognized that optical moorings are important platforms for the validation of Sea-Viewing Wide Field-of-view Sensor (SeaWiFS). It was recommended that optical moorings be maintained in order to: (1) provide long-term time series comparisons between in situ and SeaWIFS measurements of normalized water-leaving radiance; (2) develop and test algorithms for pigment biomass and phytoplankton primary productivity; and (3) provide long-term, virtually continuous in situ observations which can be used to determine and optimize the accuracy of derived satellite products. These applications require the use of in situ radiometers for long periods of time to evaluate and correct for inherent satellite undersampling (aliasing and biasing) and degradation of satellite color sensors (e.g., drifts as experienced by the Coastal Zone Color Scanner). The Bermuda Testbed Mooring (BTM) program was initiated in 1994 at a site located about 80km southeast of Bermuda in waters of about 4530 m depth. In August 1997, with NASA's support, we started to provide the Sensor Intercomparison and Merger for Biological and Interdisciplinary Oceanic Studies (SIMBIOS) program with large volumes of high frequency, long-term time-series bio-optical data from the BTM for SeaWiFS satellite ocean color groundtruthing and algorithm development. This NASA supported portion of the BTM activity spanned three years and covered five BTM deployments. During these three years, the quality of radiometric data has improved dramatically. Excellent agreement between BTM moored data and both SeaWiFS and nearby ship profile radiometric data demonstrate that technical advances in the moored optical observations have reduced the major difficulties that moored platforms face: biofouling and less frequent calibration.

Dickey, Tommy↗

View angle dependence of cloud optical thicknesses retrieved by MODIS

This study examines whether cloud inhomogeneity influences the view angle dependence of MODIS cloud optical thickness (tau) retrieval results. The degree of cloud inhomogeneity is characterized through the local gradient in 11 microns brightness temperature. The analysis of liquid phase clouds in a one year long global dataset of Collection 4 MODIS data reveals that while optical thickness retrievals give remarkably consistent results for all view directions if clouds are homogeneous, they give much higher tau-values for oblique views than for overhead views if clouds are inhomogeneous and the sun is fairly oblique. For solar zenith angles larger than 55deg, the mean optical thickness retrieved for the most inhomogeneous third of cloudy pixels is more than 30% higher for oblique views than for overhead views. After considering a variety of possible scenarios, the paper concludes that the most likely reason for the increase lies in three-dimensional radiative interactions that are not considered in current, one-dimensional retrieval algorithms. Namely, the radiative effect of cloud sides viewed at oblique angles seems to contribute most to the enhanced tau-values. The results presented here will help understand cloud retrieval uncertainties related to cloud inhomogeneity. They complement the uncertainty estimates that will start accompanying MODIS cloud products in Collection 5 and may eventually help correct for the observed view angle dependent biases.

Marshak, Alexander↗

Validation of SO2 Retrievals from the Ozone Monitoring Instrument over NE China

The Dutch-Finnish Ozone Monitoring Instrument (OMI) launched on the NASA Aura satellite in July 2004 offers unprecedented spatial resolution, coupled with contiguous daily global coverage, for space-based UV measurements of sulfur dioxide (SO2). We present a first validation of the OMI SO2 data with in situ aircraft measurements in NE China in April 2005. The study demonstrates that OMI can distinguish between background SO2 conditions and heavy pollution on a daily basis. The noise (expressed as the standard deviation,sigma) is approximately 1.5 DU (Dobson units; 1 DU = 2.69 10 (exp 16) molecules/cm (exp 2)) for instantaneous field of view boundary layer (PBL) SO2 data. Temporal and spatial averaging can reduce the noise to sigma approximetly 0.3 DU over a remote region of the South Pacific; the long-term average over this remote location was within 0.1 DU of zero. Under polluted conditions collection 2 OMI data are higher than aircraft measurements by a factor of two. Improved calibrations of the radiance and irradiance data (collection 3) result in better agreement with aircraft measurements on polluted days. The air mass corrected collection 3 data still show positive bias and sensitivity to UV absorbing aerosols. The difference between the in situ data and the OMI SO2 measurements within 30 km of the aircraft profiles was about 1 DU, equivalent to approximately 5 ppb from 0 to 3000 m altitude. Quantifying the SO2 and aerosol profiles and spectral dependence of aerosol absorption between 310 and 330 nm are critical for an accurate estimate of SO2 from satellite UV measurements.

Krotkov, Nickolay A.↗