Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial statistics”

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 235 records · Page 13

Statistics of the geomagnetic secular variation for the past 5 m.y

A new statistical model is proposed for the geomagnetic secular variation over the past 5Ma. Unlike previous models, the model makes use of statistical characteristics of the present day geomagnetic field. The spatial power spectrum of the non-dipole field is consistent with a white source near the core-mantle boundary with Gaussian distribution. After a suitable scaling, the spherical harmonic coefficients may be regarded as statistical samples from a single giant Gaussian process; this is the model of the non-dipole field. The model can be combined with an arbitrary statistical description of the dipole and probability density functions and cumulative distribution functions can be computed for declination and inclination that would be observed at any site on Earth's surface. Global paleomagnetic data spanning the past 5Ma are used to constrain the statistics of the dipole part of the field. A simple model is found to be consistent with the available data. An advantage of specifying the model in terms of the spherical harmonic coefficients is that it is a complete statistical description of the geomagnetic field, enabling us to test specific properties for a general description. Both intensity and directional data distributions may be tested to see if they satisfy the expected model distributions.

Constable, C. G.↗

Optical and Microphysical Retrievals of Marine Stratocumulus Clouds off the Coast of Namibia from Satellite and Aircraft

Though the emphasis of the Southern Africa Regional Science Initiative 2000 (SAFARI-2000) dry season campaign was largely on emission sources and transport, the assemblage of aircraft (including the high altitude NASA ER-2 remote sensing platform and the University of Washington CV-580, UK MRF C-130, and South African Weather Bureau JRA in situ aircrafts) provided a unique opportunity for cloud studies. Therefore, as part of the SAFARI initiative, investigations were undertaken to assess regional aerosol-cloud interactions and cloud remote sensing algorithms. In particular, the latter part of the experiment concentrated on marine boundary layer stratocumulus clouds off the southwest coast of Africa. Associated with cold water upwelling along the Benguela current, the Namibian stratocumulus regime has received limited attention but appears to be unique for several reasons. During the dry season, outflow of continental fires and industrial pollution over this area can be extreme. From below, upwelling provides a rich nutrient source for phytoplankton (a source of atmospheric sulfur through DMS production as well as from decay processes). The impact of these natural and anthropogenic sources on the microphysical and optical properties of the stratocumulus is unknown. Continental and Indian Ocean cloud systems of opportunity were also studied during the campaign. SAFARI 2000 aircraft flights off the coast of Namibia were coordinated with NASA Terra Satellite overpasses for synergy with the Moderate Resolution Imaging Spectroradiometer (MODIS) and other Terra instruments. MODIS was developed by NASA and launched onboard the Terra spacecraft on December 18, 1999 (and Aqua spacecraft on May 4, 2002). Among the remote sensing algorithms developed and applied to this sensor are cloud optical and microphysical properties that include cloud thermodynamic phase, optical thickness, and effective particle radius of both liquid water and ice clouds. The archived products from these algorithms have applications in climate change studies, climate modeling, numerical weather prediction, and fundamental atmospheric research. The archived MODIS Collection 5 cloud products processing stream will be used to analyze low water cloud scenes off the Namibian and Angolan coasts during SAFARI 2000 time period, as well as other years. Pixel-level Terra and Aqua MODIS retrievals (l. km spatial resolution at nadir) and gridded (1' uniform grid) statistics of cloud optical thickness and effective particle radius will be presented, including joint probability distributions between the two quantities. In addition, perspectives from the MODIS Airborne Simulator, which flew on the ER-2 during SAFARI 2000 providing high spatial resolution retrievals (50 m at nadir), will be presented as appropriate. The H-SAF Program requires an experimental operational European-centric Satellite Precipitation Algorithm System (E-SPAS) that produces medium spatial resolution and high temporal resolution surface rainfall and snowfall estimates over the Greater European Region including the Greater Mediterranean Basin. Currently, there are various types of experimental operational algorithm methods of differing spatiotemporal resolutions that generate global precipitation estimates. This address will first assess the current status of these methods and then recommend a methodology for the H-SAF Program that deviates somewhat from the current approach under development but one that takes advantage of existing techniques and existing software developed for the TRMM Project and available through the public domain.

Platnick, Steven E.↗

Far field effects of weak spatially periodic inhomogeneities.

Small variations in the properties of an acoustic medium produce cumulative effects that significantly alter wave forms propagating over long distances. To investigate some aspects of these effects, without statistical models, the one-dimensional wave equation with a weak spatially periodic variation in the sound speed is treated. An approximation is developed valid for such weak inhomogeneities over large distances based upon the approximate treatment of a linear differential equation with a periodic coefficient. It is applied to the propagation of an initial step wave and identifies the dispersive nature of the medium which accounts for the cumulative effects.

Friedman, M. B.↗

Spatially Complete Global Surface Albedos Derived from Terra/MODIS Data

Spectral land surface albedo is an important parameter for describing the radiative properties of the Earth. Accordingly it reflects the consequences of natural and human interactions, such as anthropogenic, meteorological, and phenological effects, on global and local climatological trends. Consequently, albedos are integral parts in a variety of research areas, such as general circulation models (GCMs), energy balance studies, modeling of land use and land use change, and biophysical, oceanographic, and meteorological studies. Recent production of land surface anisotropy, diffuse bihemispherical (white-sky) albedo and direct beam directional hemispherical (black-sky) albedo from observations acquired by the MODIS instruments aboard NASA s Terra and Aqua satellite platforms have provided researchers with unprecedented spatial, spectral, and temporal information on the land surface's radiative characteristics. Cloud cover, which cutails retrievals, and the presence of ephemeral and seasonal snow limit the snow-free data to approximately half the global land surfaces on an annual equal-angle basis. This precludes the MOD43B3 albedo products from being used in some remote sensing and ground-based applications, climate models, and global change research projects. An ecosystem-dependent temporal interpolation technique is described that has been developed to fill missing or seasonally snow-covered data in the official MOD43B3 albedo product. The method imposes pixel-level and local regional ecosystem-dependent phenological behavior onto retrieved pixel temporal data in such a way as to maintain pixel-level spatial and spectral detail and integrity. The phenological curves are derived from statistics based on the MODIS MOD12Q1 IGBP land cover classification product geolocated with the MOD43B3 data. The resulting snow-free value-added products provide the scientific community with spatially and temporally complete global white- and black-sky surface albedo maps and statistics. These products are stored on 1'(approximately 10 km) and coarser resolution equal-angle grids, and are computed for the first seven MODIS wavelengths, ranging from 0.47 through 2.1 microns, and for three broadband wavelengths, 0.3-0.7,0.3-5.0 and 0.7-5.0 microns.

King, Michael D.↗

Multi-resolution processing for fractal analysis of airborne remotely sensed data

Fractal geometry is increasingly becoming a useful tool for modeling natural phenomenon. As an alternative to Euclidean concepts, fractals allow for a more accurate representation of the nature of complexity in natural boundaries and surfaces. Since they are characterized by self-similarity, an ideal fractal surface is scale-independent; i.e. at different scales a fractal surface looks the same. This is not exactly true for natural surfaces. When viewed at different spatial resolutions parts of natural surfaces look alike in a statistical manner and only for a limited range of scales. Images acquired by NASA's Thermal Infrared Multispectral Scanner are used to compute the fractal dimension as a function of spatial resolution. Three methods are used to determine the fractal dimension - Schelberg's line-divider method, the variogram method, and the triangular prism method. A description of these methods and the results of applying these methods to a remotely-sensed image is also presented. Five flights were flown in succession at altitudes of 2 km (low), 6 km (mid), 12 km (high), and then back again at 6 km and 2 km. The area selected was the Ross Barnett reservoir near Jackson, Mississippi. The mission was flown during the predawn hours of 1 Feb. 1992. Radiosonde data was collected for that duration to profile the characteristics of the atmosphere. This corresponds to 3 different pixel sizes - 5m, 15m, and 30m. After, simulating different spatial sampling intervals within the same image for each of the 3 image sets, the results are cross-correlated to compare the extent of detail and complexity that is obtained when data is taken at lower spatial intervals.

Jaggi, S.↗

An X-ray study of five supernova remnants in the Carina spiral arm

The ROSAT Position Sensitive Proportional Counter (PSPC) is used to perform an exploratory study of four fields in the Carina spiral arm containing five radio supernova remnants, only one of which has previously been studied in X-rays. We present upper limits for the detection in X-rays of G298.5 - 0.3, G298.6 - 0.0, and G299.0+0.2, and report a 4 sigma detection of G296.8-0.3. In addition, we present detailed spatial and spectral analysis of the bright X-ray remnant G296.1-0.7, which has previously been studied by both the Einstein IPC and EXOSAT LE/CMA. We detect relatively slight, but statistically significant, variations in the spectrum across the remnant via spatially resolved spectral fits and a study of the spatial variation of hardness ratios. In general, the spectrum is characteristic of a thermal plasma with kT about 0.2 keV and N(sub H) about 1.5 x 10(exp 21/sq. cm). The total X-ray emitting mass is estimated to be about 250 solar mass for an optically estimated distance of 4 kpc to the remnant. At this distance, the linear dimensions of the remnant are roughly 35 - 50 pc, implying an age on the order of 20,000 yr. Assuming that X-ray and radio brightnesses are related by SIGMA(sub R) proportional to SIGMA(exp 0.69)(sub X) and that the four radio remnants have X-ray spectral characteristics similar to G296.1-0.7, we find that the column densities to these sources must be several times 10(exp 22)/sq cm in order to explain their low X-ray count rates. This column density is considerably in excess of the X-ray fitted column density to G296.1-0.7, but is comparable to the total column densities in H I measured via the 21 cm line in the directions to all five remnants. This implies that G296.1 - 0.7 is at a significantly smaller distance than the other remnants.

Hwang, Una↗

Wheat signature modeling and analysis for improved training statistics

The author has identified the following significant results. The spectral, spatial, and temporal characteristics of wheat and other signatures in LANDSAT multispectral scanner data were examined through empirical analysis and simulation. Irrigation patterns varied widely within Kansas; 88 percent of wheat acreage in Finney was irrigated and 24 percent in Morton, as opposed to less than 3 percent for western 2/3's of the State. The irrigation practice was definitely correlated with the observed spectral response; wheat variety differences produced observable spectral differences due to leaf coloration and different dates of maturation. Between-field differences were generally greater than within-field differences, and boundary pixels produced spectral features distinct from those within field centers. Multiclass boundary pixels contributed much of the observed bias in proportion estimates. The variability between signatures obtained by different draws of training data decreased as the sample size became larger; also, the resulting signatures became more robust and the particular decision threshold value became less important.

Nalepka, R. F.↗

Application of Geostatistical Simulation to Enhance Satellite Image Products

With the deployment of Earth Observing System (EOS) satellites that provide daily, global imagery, there is increasing interest in defining the limitations of the data and derived products due to its coarse spatial resolution. Much of the detail, i.e. small fragments and notches in boundaries, is lost with coarse resolution imagery such as the EOS MODerate-Resolution Imaging Spectroradiometer (MODIS) data. Higher spatial resolution data such as the EOS Advanced Spaceborn Thermal Emission and Reflection Radiometer (ASTER), Landsat and airborne sensor imagery provide more detailed information but are less frequently available. There are, however, both theoretical and analytical evidence that burn scars and other fragmented types of land covers form self-similar or self-affine patterns, that is, patterns that look similar when viewed at widely differing spatial scales. Therefore small features of the patterns should be predictable, at least in a statistical sense, with knowledge about the large features. Recent developments in fractal modeling for characterizing the spatial distribution of undiscovered petroleum deposits are thus applicable to generating simulations of finer resolution satellite image products. We will present example EOS products, analysis to investigate self-similarity, and simulation results.

Hlavka, Christine A.↗

The spatial correlation function of rich clusters of galaxies

A series of objective statistical estimators is applied to directly study the three-dimensional distribution of rich clusters, using a recently completed redshift sample of 104 Abell clusters out to distance class four and galactic latitude of 30 deg or more. The sample and the relevant positional and redshift selection effects are discussed, and the two-point angular correlation function and nearest-neighbor distribution of clusters in the sample are determined, as well as in the larger D = 5+6 sample and in several subsamples of different richness classes and latitude zones. The redshift correlation function of clusters for various angular separations on the sky is examined as a test for the reality of the observed angular correlations. The two-point spatial correlation function of clusters in the sample is determined, and the consistency among the three observed correlation functions is tested. A model that describes all the observed correlations is presented.

Bahcall, N. A.↗

Tropospheric Ozone Change from 1980 to 2010 Dominated by Equatorward Redistribution of Emissions

Ozone is an important air pollutant at the surface, and the third most important anthropogenic greenhouse gas in the troposphere. Since 1980, anthropogenic emissions of ozone precursors methane, non-methane volatile organic compounds, carbon monoxide and nitrogen oxides (NOx) have shifted from developed to developing regions. Emissions have thereby been redistributed equatorwards, where they are expected to have a stronger effect on the tropospheric ozone burden due to greater convection, reaction rates and NOx sensitivity. Here we use a global chemical transport model to simulate changes in tropospheric ozone concentrations from 1980 to 2010, and to separate the influences of changes in the spatial distribution of global anthropogenic emissions of short-lived pollutants, the magnitude of these emissions, and the global atmospheric methane concentration. We estimate that the increase in ozone burden due to the spatial distribution change slightly exceeds the combined influences of the increased emission magnitude and global methane. Emission increases in Southeast, East and South Asia may be most important for the ozone change, supported by an analysis of statistically significant increases in observed ozone above these regions. The spatial distribution of emissions dominates global tropospheric ozone, suggesting that the future ozone burden will be determined mainly by emissions from low latitudes.

Zhang, Yuqiang↗

Inversion of data from diffraction-limited multiwavelength remote sensors. I - Linear case

The remote sensing inverse problem is considered in which the sensor does not necessarily view the same area on the earth at each wavelength, and spatial correlations of the geophysical parameters may be present. Under the conditions of linearity and stationary statistics, the minimum mean-square error solution to the problem of inverting such data is a spatial filter of the Wiener-Kolmogorov class. The resulting remote-sensing system can be characterized by an impulse response matrix in ordinary space or by a transfer matrix in frequency space. A signal-to-noise matrix for the geophysical parameters to be sensed is also defined; this matrix depends on the postulated a priori statistics and on the characteristics of the remote-sensing instrument. The system transfer matrix and the signal-to-noise matrix are simultaneously diagonalizable. The optimum transfer matrix filters out of the estimate vector those eigenparameters for which eigenvalues of the signal-to-noise matrix are less than unity.

Rosenkranz, P. W.↗

Substrate Effects on Electronic Properties of Atomic Chains

When the device size is reduced down to 0.07 micrometers, the number of dopant atoms in the channel will no longer be macroscopic, typically less than a hundred. A spatial distribution of these dopant atoms will fluctuate statistically from device to device even in identically designed devices, and this places a serious limitation for integration. It is, however, impractical to control dopant positions within atomic dimension. One fundamental solution to this problem is to create electronics with atomically precise, but preferably simple structures. Atomic chains, precise structures of adatoms created on an atomically regulated surface, are candidates for constituent components in future electronics. All the adatoms will be placed at designated positions on the substrate, and all the device structures will be precise, free from any deviations. It was predicted using the tight-binding calculation with universal parameters that silicon chains were metallic and magnesium chains were semiconducting regardless of the lattice spacing, and a possible doping method was also proposed. In these treatments, the substrate was assumed to serve as a non-interacting template holding the adatoms without a formation of chemical bonding with substrate atoms. However, this scheme may not be easy to implement experimentally. Adatoms will have to be fixed with a van der Waals force on the substrate, but the force is generally weak and an extremely low temperature environment has to be prepared to suppress their unwanted thermal displacement. It may be logical to seek a scheme to allow the adatoms to form chemical bonding with the substrate atoms and secure their positions. The substrate effects are studied in detail.

Yamada, Toshishige↗

The Global Structure of UTLS Ozone in GEOS-5: A Multi-Year Assimilation of EOS Aura Data

Eight years of ozone measurements retrieved from the Ozone Monitoring Instrument (OMI) and the Microwave Limb Sounder, both on the EOS Aura satellite, have been assimilated into the Goddard Earth Observing System version 5 (GEOS-5) data assimilation system. This study thoroughly evaluates this assimilated product, highlighting its potential for science. The impact of observations on the GEOS-5 system is explored by examining the spatial distribution of the observation-minus-forecast statistics. Independent data are used for product validation. The correlation coefficient of the lower-stratospheric ozone column with ozonesondes is 0.99 and the bias is 0.5%, indicating the success of the assimilation in reproducing the ozone variability in that layer. The upper-tropospheric assimilated ozone column is about 10% lower than the ozonesonde column but the correlation is still high (0.87). The assimilation is shown to realistically capture the sharp cross-tropopause gradient in ozone mixing ratio. Occurrence of transport-driven low ozone laminae in the assimilation system is similar to that obtained from the High Resolution Dynamics Limb Sounder (HIRDLS) above the 400 K potential temperature surface but the assimilation produces fewer laminae than seen by HIRDLS below that surface. Although the assimilation produces 5 - 8 fewer occurrences per day (up to approximately 20%) during the three years of HIRDLS data, the interannual variability is captured correctly. This data-driven assimilated product is complementary to ozone fields generated from chemistry and transport models. Applications include study of the radiative forcing by ozone and tracer transport near the tropopause.

GEOS-5↗

Supersonic Disk Gap Band Parachute Performance in the Wake of a Viking-Type Aeroshell from Mach 2 to 2.5

Supersonic wind tunnel testing of 0.813 m diameter Disk-Gap-Band parachutes is being conducted in the NASA Glenn Research Center (GRC) 10' x 10' wind-tunnel. The tests are conducted in support of the Mars Science Laboratory Parachute Decelerator System development and qualification. Four percent of full-scale parachutes were constructed similarly to the flight-article in material and construction techniques. The parachutes are attached to a 4% scale MSL entry-vehicle to simulate the free-flight configuration. The parachutes are tested from Mach 2 to 2.5 over a Reynolds number (Re) range of 1 to 3 x 10(exp 6), representative of the MSL deployment envelope. Constrained and unconstrained test configurations are investigated to quantify the effects of parachute trim, suspension line interaction, and alignment with the capsule wake. The parachute is constrained horizontally through the vent region, to measure canopy breathing and wake interaction for fixed trim angles of 0 and 10 degrees from the velocity vector. In the unconstrained configuration the parachute is permitted to trim and cone, similar to the free-flight varying its alignment relative to the entry-vehicle wake. Test diagnostics were chosen to quantify parachute performance and to provide insight into the flow field structure. An in-line load cell provided measurement of unsteady and mean drag as a function of Mach and Re. High-speed shadowgraph video of the upstream parachute flow field was used to capture bow-shock motion and stand of distance. Particle image velocimetry of the upstream parachute flow field provides spatially and temporally resolved measurement velocity and turbulent statistics. Multiple high speed video views of targets placed in the interior of the canopy enable photo-grammetric measurement of the fabric motion in time and space from reflective. High speed video is also used to document the supersonic inflation and measure trim angle, projected area, and frequency of area oscillations.

EDL↗

Moisture and Temperature Influences on Nonlinear Vegetation Trends in Serengeti National Park

While long-term vegetation greening trends have appeared across large land areas over the late 20th century, uncertainty remains in identifying and attributing finer-scale vegetation changes and trends, particularly across protected areas. Serengeti National Park (SNP) is a critical East African protected area, where seasonal vegetation cycles support vast populations of grazing herbivores and a host of ecosystem dynamics. Previous work has shown how non-climate drivers (e.g. land use) shape the SNP ecosystem, but it is still unclear to what extent changing climate conditions influence SNP vegetation, particularly at finer spatial and temporal scales. We fill this research gap by evaluating long-term (1982–2016) changes in SNP leaf area index (LAI) in relation to both temperature and moisture availability using Ensemble Empirical Mode Decomposition and Principal Component Analysis with regression techniques. We find that SNP LAI trends are nonlinear, display high sub-seasonal variation, and are influenced by lagged changes in both moisture and temperature variables and their interactions. LAI during the long rains (e.g. March) exhibits a greening-to-browning trend reversal starting in the early 2000s, partly due to antecedent precipitation declines. In contrast, LAI during the short rains (e.g. November, December) displays browning-to-greening alongside increasing moisture availability. Rising temperature trends also have important, secondary interactions with moisture variables to shape these SNP vegetation trends. Our findings show complex vegetation-climate interactions occurring at important temporal and spatial scales of the SNP, and our rigorous statistical approaches detect these complex climate-vegetation trends and interactions, while guarding against spurious vegetation signals.

Moisture↗

Lightning Strike Distance Distribution Beyond a Preexisting Lightning Area

The 45th Weather Squadron (45 WS) asked the Applied Meteorology Unit (AMU) to review the 30-year-old, lightning stand-off distances of 5 nautical miles (nmi) for applicability to today's operations. This was based on the realization that previous lightning strike distance studies did not match how 45 WS issues lightning warnings (Roeder, 2008). The previous lightning distance studies were from the point of origin of the lightning or the average starting location that would tend to be in the core of the thunderstorm. However, the 45 WS issues lightning warnings based on the edge of a preexisting lightning area. Before beginning the AMU project, it took several years to develop a method to calculate a distance distribution beyond a preexisting area (Roeder, 2015). The AMU pulled Lightning Detection and Ranging (LDAR) sensor data from 1/1/2013 to 12/31/2013. This dataset consisted of 37 million individual source data points from the LDAR sensors. Only sources within 50 km north, south, east or west of the LDAR grid center were included in the dataset. This limited the use of LDAR data to that with the greatest accuracy of source detection and increased data processing. Points were grouped into flashes based on spatial and temporal criteria. Based on the sensitivity analysis the AMU performed on the flash clustering algorithm, a time value of 0.3 seconds was found to model flashes adequately. Distance parameters were tested from 1,500 to 7,500 meters (m) in 500 m increments. Distance parameters of both 3,000 m and 4,000 m produced results in the plotting tool that were most representative of the physical behavior of lightning. Thus statistics were gathered for the most representative of these spatial and temporal criteria on the flash size and the polygon expansion distance in order to find the correct data distributions. The best fit curves for the LDAR polygon expansion frequency vs. distances for both the 3 kilometer (km) and 4 km distance threshold values were exponential decay functions and had R2 values of > 0.998, indicating good model fits. The equations of the best fit curves were then used to calculate a desired safety radius of 4 nmi for either 3 km or 4 km distance threshold criteria. The AMU analysis concludes the safe reduction of the 5 nmi lightning warning circles to 4 nmi should improve the operational impact by 36% if based on distance from the center of the property area being protected. If based on the edge of the property being protected, then the reduction is 4.5 nmi to 4 nmi and the operational impact is decreased by 21%. For the 6 nmi lightning warning circles, the recommended 4 nmi stand-off distance will result in a safe reduction of operational impact of 31% if based on the center of the area being protected, or 16% if based on the edge of the property.

Flash clustering algorithm↗

Studies in the use of cloud type statistics in mission simulation

A study to further improve NASA's global cloud statistics for mission simulation is reported. Regional homogeneity in cloud types was examined; most of the original region boundaries defined for cloud cover amount in previous studies were supported by the statistics on cloud types and the number of cloud layers. Conditionality in cloud statistics was also examined with special emphasis on temporal and spatial dependencies, and cloud type interdependence. Temporal conditionality was found up to 12 hours, and spatial conditionality up to 200 miles; the diurnal cycle in convective cloudiness was clearly evident. As expected, the joint occurrence of different cloud types reflected the dynamic processes which form the clouds. Other phases of the study improved the cloud type statistics for several region and proposed a mission simulation scheme combining the 4-dimensional atmospheric model, sponsored by MSFC, with the global cloud model.

Fowler, M. G.↗

Ice Mass Change in Greenland and Antarctica Between 1993 and 2013 from Satellite Gravity Measurements

We construct long-term time series of Greenland and Antarctic ice sheet mass change from satellite gravity measurements. A statistical reconstruction approach is developed based on a principal component analysis (PCA) to combine high-resolution spatial modes from the Gravity Recovery and Climate Experiment (GRACE) mission with the gravity information from conventional satellite tracking data. Uncertainties of this reconstruction are rigorously assessed; they include temporal limitations for short GRACE measurements, spatial limitations for the low-resolution conventional tracking data measurements, and limitations of the estimated statistical relationships between low- and high-degree potential coefficients reflected in the PCA modes. Trends of mass variations in Greenland and Antarctica are assessed against a number of previous studies. The resulting time series for Greenland show a higher rate of mass loss than other methods before 2000, while the Antarctic ice sheet appears heavily influenced by interannual variations.

Talpe, Matthieu J.↗