Search NASA⌕ Search

SEARCH · Search NASA

Results for “uncertainty distribution”

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 163 records · Page 9

Optimal Estimation Framework for Ocean Color Atmospheric Correction and Pixel-level Uncertainty Quantification

Ocean color remote sensing requires compensation for atmospheric scattering and absorption (aerosol, Rayleigh, and trace gases), referred to as atmospheric correction (AC). AC allows inference of parameters such as spectrally resolved remote sensing reflectance ( R rs )(λ) ; sr 1 ) at the ocean surface from the top-of-atmosphere reflectance. Often, the uncertainty of this process is not fully explored. Bayesian inference techniques provide a simultaneous AC and uncertainty assessment via a full posterior distribution of the relevant variables, given the prior distribution of those variables and the radiative transfer (RT) likelihood function. Given uncertainties in the algorithm inputs, the Bayesian framework enables better constraints on the AC process by using the complete spectral information compared to traditional approaches that use only a subset of bands for AC. This paper investigates a Bayesian inference research method (Optimal Estimation, OE) for ocean color AC by simultaneously retrieving atmospheric and ocean properties using all visible and near-infrared spectral bands. The OE algorithm analytically approximates the posterior distribution of parameters based on normality assumptions and provides a potentially viable operational algorithm with a reduced computational expense. We developed a Neural Network (NN) RT forward model look-up-table-based emulator to increase algorithm efficiency further and thus speed up the likelihood computations. We then applied the OE algorithm to synthetic data and observations from the MODerate resolution Imaging Spectroradiometer (MODIS) on NASA’s Aqua spacecraft. We compared the R rs )(λ) retrieval and its uncertainty estimates from the OE method with in-situ validation data from the SeaWiFS Bio-optical Archive and Storage System (SeaBASS) and Aerosol Robotic Network Ocean Color (AERONET-OC) datasets. The OE algorithm improved R rs )(λ) estimates relative to the NASA standard operational algorithm by improving all statistical metrics at 443, 555, and 667 nm. Unphysical negative R rs )(λ) , which often appear in complex water conditions, was reduced by a factor of 3. The OE-derived pixel-level R rs )(λ) uncertainty estimates were also assessed relative to in-situ data and were shown to have skill.

Atmospheric correction↗

Atmospheric processing and aerosol aging responsible for observed increase in absorptivity of long-range-transported smoke over the southeast Atlantic

Biomass burning aerosol (BBA) from agricultural fires in southern Africa contributes about one-third of the global carbonaceous aerosol load. These particles have strong radiative effects in the southeast Atlantic (SEA), which depend in part on the radiative contrast between the aerosol layer in the free troposphere (FT) and the underlying cloud layer. However, there is large disagreement in model estimates of aerosol-driven climate forcing due to uncertainties in the vertical distribution, optical properties, and life cycle of these particles. This study applies a novel method combining remote sensing observations with regional model outputs to investigate the aging of the BBA and its impact on the optical properties during transatlantic transport from emission sources in Africa to the SEA. Results show distinct variations in extinction Ångström exponent (EAE) and single-scattering albedo (SSA) as aerosols age. Near the source, fresh aerosols are characterized by low mean SSA (0.84) and high EAE (1.85), indicating smaller, highly absorbing particles. By isolating marine contributions from the total column during BBA transport across the SEA, our analysis reveals an initial decrease in BBA absorptivity, with mean FT SSA of 0.87 after 6–7 d, followed by increased absorptivity with mean FT SSA of 0.84 after 10 d, suggesting enhanced absorption due to chemical aging. These findings indicate that BBA becomes more absorbing during extended transport across the SEA, with implications for reducing model uncertainties. Our remote-sensing-based results agree well with previous in situ studies and offer new insights into aerosol–radiation interactions and the energy balance over the SEA.

54 ENVIRONMENTAL SCIENCES↗

Finding Uncertainties that Cause the Age Dependence of Dose Limits to Be Immature

Space radiation permissible exposure limits (PEL) are intended to set acceptable levels of cancer risks, and avoid any clinical significant non-cancer effects. The 1989 recommendation of the National Council of Radiation Protection and Measurements (NCRP) recommended a strong age dependence of dose limits that departed drastically from the then mature 1970 dose limits recommendations from the National Academy of Science, which were independent of age. In 2000, the NCRP recommended revised limits that showed a similar trend of risk with age to the 1989 report. In this model, the cancer risk per Sv varies by more than 2-fold for ages between 30- and 50-yr. Therefore for galactic cosmic rays exposure, astronaut age has a larger influence on risk then radiation shielding mass or material composition, vehicle propulsion method, or position in the solar cycle. For considering the control of mission costs and resources, the possibility of using astronaut age as a trade variable in mission design could be considered. However, the uncertainties in describing the age dependence on risk have not been fully explored. We discuss biological factors that influence the age dependence of radiation risks, including susceptibility, expression and latency, and radiation quality. These factors depend not only on the individual s age, but also their genetic sensitivity and interaction with other environmental factors. Epidemiological data is limited in describing the age dependence on risk. The 2005, BEIR VII report recommends an age dependence for cancer risk attributable solely to the life-table disagreeing strongly with the NCRP model. However, BEIR VII also noted the limited power of human data for concomitantly describing both age and age after exposure dependences of cancer risks. Many experimental studies have shown that high LET radiation (e.g., high charge and energy (HZE) nuclei and neutrons) display reduced latency compared to low LET radiation, suggesting distinct biological factors are important. We discuss potential molecular mechanisms that would influence the age dependence of radiation risks. A probability distribution function for the uncertainties in age-dependence of risk models is described and predictions for Mars missions discussed. Our report suggests that theoretical considerations based on new experimental studies are needed to ensure the correct age dependence in space radiation risk models and the resulting Astronaut PEL.

Cucinotta, Francis A.↗

On the appropriateness of applying chi-square distribution based confidence intervals to spectral estimates of helicopter flyover data

The validity of applying chi-square based confidence intervals to far-field acoustic flyover spectral estimates was investigated. Simulated data, using a Kendall series and experimental acoustic data from the NASA/McDonnell Douglas 500E acoustics test, were analyzed. Statistical significance tests to determine the equality of distributions of the simulated and experimental data relative to theoretical chi-square distributions were performed. Bias and uncertainty errors associated with the spectral estimates were easily identified from the data sets. A model relating the uncertainty and bias errors to the estimates resulted, which aided in determining the appropriateness of the chi-square distribution based confidence intervals. Such confidence intervals were appropriate for nontonally associated frequencies of the experimental data but were inappropriate for tonally associated estimate distributions. The appropriateness at the tonally associated frequencies was indicated by the presence of bias error and noncomformity of the distributions to the theoretical chi-square distribution. A technique for determining appropriate confidence intervals at the tonally associated frequencies was suggested.

Rutledge, Charles K.↗

Accelerating Hamiltonian Monte Carlo for Bayesian inference in neural networks and neural operators

Hamiltonian Monte Carlo (HMC) is a powerful and accurate method to sample from the posterior distribution in Bayesian inference. However, HMC techniques are computationally demanding for Bayesian neural networks due to the high dimensionality of the network’s parameter space and the non-convexity of their posterior distributions. Therefore, various approximation techniques, such as variational inference (VI) or stochastic gradient MCMC, are often employed to infer the posterior distribution of the network parameters. Such approximations introduce inaccuracies in the inferred distributions, resulting in unreliable uncertainty estimates. In this work, we propose a hybrid approach that combines inexpensive VI and accurate HMC methods to efficiently and accurately quantify uncertainties in neural networks and neural operators. The proposed approach leverages an initial VI training on the full network. We examine the influence of individual parameters on the prediction uncertainty, which shows that a large proportion of the parameters do not contribute substantially to uncertainty in the network predictions. This information is then used to significantly reduce the dimension of the parameter space, and HMC is performed only for the subset of network parameters that strongly influence prediction uncertainties. This yields a framework for accelerating the full batch HMC for posterior inference in neural networks. We demonstrate the efficiency and accuracy of the proposed framework on deep neural networks and operator networks, showing that inference can be performed for large networks with tens to hundreds of thousands of parameters. Finally, we show that this method can effectively learn surrogates for complex physical systems by modeling the operator that maps from upstream conditions to wall-pressure data on a cone in hypersonic flow.

Bayesian inference↗

Constrained variational optimization of counting-time allocation in sequential scattering measurements: Application to Bonse–Hart USANS

Sequential scattering measurements are often performed under a fixed experimental-time budget, even though the expected count rate varies strongly across the measured coordinate. When the dwell time at each measurement position can be controlled independently, this variation creates a general resource-allocation problem: how should the available time be distributed to minimize the uncertainty of the reconstructed profile? We formulate this problem as a constrained variational optimization for measurements governed by Poisson counting statistics. When each measurement is treated independently, minimizing the averaged squared relative uncertainty yields an inverse-square-root intensity allocation. The formulation is then generalized to include correlations between neighboring measurements and an instrumental resolution operator, leading to an allocation criterion that equalizes the marginal reduction in posterior uncertainty per unit measurement time. Bonse–Hart ultra-small-angle neutron scattering (USANS), in which reciprocal space is sampled sequentially through analyzer-angle stepping, provides an experimentally grounded application. Computational benchmarking shows that the optimized allocation outperforms uniform-time and constant-relative-error strategies, while application to an experimentally measured graphite USANS profile from the Spallation Neutron Source, using Poisson resampling under alternative schedules, demonstrates how counting time should be redistributed toward weak-intensity regions under an identical total duration. The resulting framework applies to sequential scattering and related scanning measurements whenever local dwell times are adjustable and directly determine the measurement uncertainties, and when the relevant correlation and instrumental-response models are available.

Tung, Chi-Huan [ORNL] (ORCID:0000000221972074)↗

Electromagnetic energy calibration of the SoLid detector with horizontal muons

SoLid is a neutrino experiment at very-short baselinesearching for active-to-sterile oscillations of reactorantineutrinos. The detection principle is based on the pairing oftwo types of solid scintillators: polyvinyl toluene and $^{6}$Li:ZnS(Ag), which is a new technology used in this field ofPhysics. In addition to good neutron-gamma discrimination, thissetup allows the detector to be highly segmented; the basicdetection unit is a 5 cm cube. High segmentation provides numerousadvantages including precise localisation of the Inverse Beta Decay(IBD) products, the derivation of an antineutrino energy estimatorbased on the isolated positron energy, and a powerful backgroundreduction tool that relies on the topological signature of thesignal. Finally, the system is read out by a network ofwavelength-shifting fibres coupled to photosensors. A relativeelectromagnetic calibration is performed with horizontal cosmicmuons. This source poses the simplest calibration problem in which asingle detection unit is involved. In addition, large muon energydeposits allow us to perform a calibration at the most detailedlevel (i.e. per fibre) and to accurately define the fraction ofenergy escaping to neighbouring detection cells. A statisticalprecision at the sub-percent level is reached. The paper alsodiscusses two methods to calibrate the absolute energy scale andpresents their implementation and results. The first method relieson horizontal muons, though the precision is limited to around 10%because of the uncertainty in the energy distribution of suchmuons. A novel, alternative method based on the radioactiveamericium-beryllium source is proposed. It takes advantage of theelectron-positron pair-production process and provides a calibrationpoint at 3.4 MeV (i.e. in the core of the IBD positronspectrum). The paper is concluded with various cross-check includinga determination of the energy spectrum of the standard cosmogenicbackground candle: $^{12}$B.

Abreu, Y. [Antwerp U.]↗

Demand Response Under Stochastic, Price-Dependent User Behavior

This letter focuses on price-based demand response (DR) implemented through dynamic adjustments of electricity prices. It extends existing DR models to a stochastic framework in which customer response is represented by price-dependent random variables, leveraging models and tools from the theory of stochastic optimization with decision-dependent distributions. The inherent epistemic uncertainty in the customers' responses renders open-loop, model-based DR strategies impractical. We propose a stochastic, feedback-based pricing strategy to compensate for estimation errors and uncertainty in customer response, establish theoretical results demonstrating the stability and near-optimality of the proposed approach, and validate its effectiveness through numerical simulations.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Consistent detection of global predicates

A fundamental problem in debugging and monitoring is detecting whether the state of a system satisfies some predicate. If the system is distributed, then the resulting uncertainty in the state of the system makes such detection, in general, ill-defined. Three algorithms are presented for detecting global predicates in a well-defined way. These algorithms do so by interpreting predicates with respect to the communication that has occurred in the system.

Cooper, Robert↗

Spectroscopic measurement of HO2, H2O2, and OH in the stratosphere

Stratospheric concentrations of HO2, H2O2, and OH have been retrieved simultaneously from the far-infrared emission spectra obtained with a balloon-borne Fourier transform spectrometer in June 1983 at 32 deg N latitude. Retrieved concentrations of HO2 and H2O2 are reported, along with vertical distributions of OH which were reported in an earlier paper for the afternoon, sunset, and nighttime periods for altitudes from 26 to 38 km. HO2 distributions are obtained with uncertainties that are about the same as OH for the same vertical range and for the afternoon and sunset periods. H2O2 concentration is obtained at an altitude of 30 km for the period that covers afternoon and sunset hours. The retrieved concentrations of these HO(x) species agree well with other individually measured results and the steady state photochemical predictions. The ratio HO2/OH at around 32 km seems to increase from the afternoon period to the sunset period.

Park, J. H.↗

A comparison between general circulation model simulations using two sea surface temperature datasets for January 1979

Simulations with the UCLA atmospheric general circulation model (AGCM) using two different global sea surface temperature (SST) datasets for January 1979 are compared. One of these datasets is based on Comprehensive Ocean-Atmosphere Data Set (COADS) (SSTs) at locations where there are ship reports, and climatology elsewhere; the other is derived from measurements by instruments onboard NOAA satellites. In the former dataset (COADS SST), data are concentrated along shipping routes in the Northern Hemisphere; in the latter dataset High Resolution Infrared Sounder (HIRS SST), data cover the global domain. Ensembles of five 30-day mean fields are obtained from integrations performed in the perpetual-January mode. The results are presented as anomalies, that is, departures of each ensemble mean from that produced in a control simulation with climatological SSTs. Large differences are found between the anomalies obtained using COADS and HIRS SSTs, even in the Northern Hemisphere where the datasets are most similar to each other. The internal variability of the circulation in the control simulation and the simulated atmospheric response to anomalous forcings appear to be linked in that the pattern of geopotential height anomalies obtained using COADS SSTs resembles the first empirical orthogonal function (EOF 1) in the control simulation. The corresponding pattern obtained using HIRS SSTs is substantially different and somewhat resembles EOF 2 in the sector from central North America to central Asia. To gain insight into the reasons for these results, three additional simulations are carried out with SST anomalies confined to regions where COADS SSTs are substantially warmer than HIRS SSTs. The regions correspond to warm pools in the northwest and northeast Pacific, and the northwest Atlantic. These warm pools tend to produce positive geopotential height anomalies in the northeastern part of the corresponding oceans. Both warm pools in the Pacific produce large-scale circulation anomalies with a pattern that resembles that obtained using COADS SSTs as well as EOF 1 of the control simulation; the warm pool in the Atlantic does not. These results suggest that the differences obtained with COADS SSTs and HIRS SSTs are mostly due to the differences in the datasets over the northern Pacific. There was a blocking episode near Greenland in late January 1979. Both simulations with warm SST anomalies over the northwest and northeast Pacific show a tendency toward increased incidence of North Atlantic blocking; the simulation with warm SST anomalies over the northwest Atlantic shows a tendency toward decreased incidence. These results suggest that features in both SST datasets that do not have a counterpart in the other dataset contribute signficantly to the differences between the simulated and observed fields. The results of this study imply that uncertainties in current SST distributions for the world oceans can be as important as the SST anomalies themselves in terms of their impact on the atmospheric circulation. Caution should be exercised, therefore, when linking anomalous circulation and SST patterns, especially in long-range prediction.

Ose, Tomoaki↗

Submillimeter-Wave Cloud Ice Radiometry

Submillimeter-wave cloud ice radiometry is a new and innovative technique for characterizing cirrus ice clouds. Cirrus clouds affect Earth's climate and hydrological cycle by reflecting incoming solar energy, trapping outgoing IR radiation, sublimating into vapor, and influencing atmospheric circulation. Since uncertainties in the global distribution of cloud ice restrict the accuracy of both climate and weather models, successful development of this technique could provide a valuable tool for investigating how clouds affect climate and weather. Cloud ice radiometry could fill an important gap in the observational capabilities of existing and planned Earth-observing systems. Using submillimeter-wave radiometry to retrieve properties of ice clouds can be understood with a simple model. There are a number of submillimeter-wavelength spectral regions where the upper troposphere is transparent. At lower tropospheric altitudes water vapor emits a relatively uniform flux of thermal radiation. When cirrus clouds are present, they scatter a portion of the upwelling flux of submillimeter-wavelength radiation back towards the Earth as shown in the diagram, thus reducing the upward flux o f energy. Hence, the power received by a down-looking radiometer decreases when a cirrus cloud passes through the field of view causing the cirrus cloud to appear radiatively cool against the warm lower atmospheric thermal emissions. The reduction in upwelling thermal flux is a function of both the total cloud ice content and mean crystal size. Radiometric measurements made at multiple widely spaced frequencies permit flux variations caused by changes in crystal size to be distinguished from changes in ice content, and polarized measurements can be used to constrain mean crystal shape. The goal of the cloud ice radiometry program is to further develop and validate this technique of characterizing cirrus. A multi-frequency radiometer is being designed to support airborne science and spacecraft validation missions. This program has already extended the initial millimeter-wave modeling studies to submillimeter-wavelengths and has improved the realism of the cloud scattering models. Additionally a proof-of-concept airborne submillimeter-wave radiometer was constructed and fielded. It measured a radiometric signal from cirrus confirming the basic technical feasibility of this technique. This program is a cooperative effort of the University of Colorado, Colorado State University, Swales Aerospace, and Jet Propulsion Laboratory. Additional information is contained in the original.

Walter, Steven J.↗

Modeling and Visualizing Uncertainty in Continuous Variables Predicted using Remotely Sensed Data

The use of remotely sensed images to map continuous biophysical variables, such as those related to terrestrial vegetation amount, sea surface temperature, and many other targets of NASA s Earth Observing System (EOS), includes variable, parametric, positional, spatial support and structural sources of uncertainty. A complete description of uncertainty will lead to a probability distribution at each location, allowing the exploration of the spatial dimension of uncertainty, that is, where the field is not well quantified. To achieve this purpose, convenient visualization tools are required. We have produced such a tool, called PDFVis, that facilitates the display of probability density functions (pdfs) on a per-grid-cell basis. The density estimate from Monte-Carlo generated realizations is interactively displayed as well as parametric and non-parametric summaries of the pdf field (such as mean, median, quartiles, standard deviation, number of modes, and locations of modes). Shaded surface renderings of pdfs along a transect can also be projected onto a plane. This tool will become more useful as richer descriptions of spatial uncertainty become available.

Dungan, Jennifer L.↗

Development of an Expert Judgement Elicitation and Calibration Methodology for Risk Analysis in Conceptual Vehicle Design

A comprehensive expert-judgment elicitation methodology to quantify input parameter uncertainty and analysis tool uncertainty in a conceptual launch vehicle design analysis has been developed. The ten-phase methodology seeks to obtain expert judgment opinion for quantifying uncertainties as a probability distribution so that multidisciplinary risk analysis studies can be performed. The calibration and aggregation techniques presented as part of the methodology are aimed at improving individual expert estimates, and provide an approach to aggregate multiple expert judgments into a single probability distribution. The purpose of this report is to document the methodology development and its validation through application to a reference aerospace vehicle. A detailed summary of the application exercise, including calibration and aggregation results is presented. A discussion of possible future steps in this research area is given.

Unal, Resit↗

Cratering Rates in the Jovian System

We use several independent constraints on the number of ecliptic comets to estimate impact cratering rates on the Jupiter moons. The impact rate on Jupiter by 1.5-km diameter ecliptic comets is currently NY(d > 1.5km) = 0.005(+0.006)(-0.003) per annum. Asteroids and long period comets are currently unimportant. The size-number distribution of ecliptic comets smaller than 20 km is inferred from size-number distributions of impact craters on Europa, Ganymede, and Triton. For comets bigger than 50 km we use the size-number distribution of Kuiper Belt Objects. The overview of the impact rate at Jupiter in general and at Europa in particular are given. These impact rates imply cratering rates on Europa of 0.5 per Ma per 10(exp 6) sq km for impact craters bigger than 1 km, and of 0.015 per Ma per 10(exp 6) sq km for impact craters bigger than 20 km. The latter corresponds to an average recurrence time of 2.2 Ma for 20 km craters. The best current estimates for the number of 20 km craters on Europa appear to range between about twelve to thirty. This implies that the average age of Europa's surface is between 30 and 70 Ma. The average density of craters with diameter greater than 1 km on well-mapped swaths on Europa is 30 per 10(exp 6) sq km. The corresponding nominal surface age would be 60 Ma. These two estimates are not truly independent because we have used size-number distribution of the Europan craters to help generate the size-number distribution of comets. The uncertainty of the best estimate - call it 42 Ma for specificity - is at least a factor of 3.

Zahnle, K.↗

Evaluating Shielding Effectiveness for Reducing Space Radiation Cancer Risks

We discuss calculations of probability distribution functions (PDF) representing uncertainties in projecting fatal cancer risk from galactic cosmic rays (GCR) and solar particle events (SPE). The PDF s are used in significance tests of the effectiveness of potential radiation shielding approaches. Uncertainties in risk coefficients determined from epidemiology data, dose and dose-rate reduction factors, quality factors, and physics models of radiation environments are considered in models of cancer risk PDF s. Competing mortality risks and functional correlations in radiation quality factor uncertainties are treated in the calculations. We show that the cancer risk uncertainty, defined as the ratio of the 95% confidence level (CL) to the point estimate is about 4-fold for lunar and Mars mission risk projections. For short-stay lunar missions (<180 d), SPE s present the most significant risk, however one that is mitigated effectively by shielding, especially for carbon composites structures with high hydrogen content. In contrast, for long duration lunar (>180 d) or Mars missions, GCR risks may exceed radiation risk limits, with 95% CL s exceeding 10% fatal risk for males and females on a Mars mission. For reducing GCR cancer risks, shielding materials are marginally effective because of the penetrating nature of GCR and secondary radiation produced in tissue by relativistic particles. At the present time, polyethylene or carbon composite shielding can not be shown to significantly reduce risk compared to aluminum shielding based on a significance test that accounts for radiobiology uncertainties in GCR risk projection.

Cucinotta, Francis A.↗

Aeroheating Testing and Predictions for Project Orion CEV at Turbulent Conditions

An investigation of the aeroheating environment of the Project Orion Crew Exploration Vehicle was performed in the Arnold Engineering Development Center Hypervelocity Wind Tunnel No. 9 Mach 8 and Mach 10 nozzles and in the NASA Langley Research Center 20 - Inch Mach 6 Air Tunnel. Heating data were obtained using a thermocouple-instrumented approx.0.035-scale model (0.1778-m/7-inch diameter) of the flight vehicle. Runs were performed in the Tunnel 9 Mach 10 nozzle at free stream unit Reynolds numbers of 1x10(exp 6)/ft to 20x10(exp 6)/ft, in the Tunnel 9 Mach 8 nozzle at free stream unit Reynolds numbers of 8 x 10(exp 6)/ft to 48x10(exp 6)/ft, and in the 20-Inch Mach 6 Air Tunnel at free stream unit Reynolds numbers of 1x10(exp 6)/ft to 7x10(exp 6)/ft. In both facilities, enthalpy levels were low and the test gas (N2 in Tunnel 9 and air in the 20-Inch Mach 6) behaved as a perfect-gas. These test conditions produced laminar, transitional and turbulent data in the Tunnel 9 Mach 10 nozzle, transitional and turbulent data in the Tunnel 9 Mach 8 nozzle, and laminar and transitional data in the 20- Inch Mach 6 Air Tunnel. Laminar and turbulent predictions were generated for all wind tunnel test conditions and comparisons were performed with the experimental data to help define the accuracy of computational method. In general, it was found that both laminar data and predictions, and turbulent data and predictions, agreed to within less than the estimated 12% experimental uncertainty estimate. Laminar heating distributions from all three data sets were shown to correlate well and demonstrated Reynolds numbers independence when expressed in terms of the Stanton number based on adiabatic wall-recovery enthalpy. Transition onset locations on the leeside centerline were determined from the data and correlated in terms of boundary-layer parameters. Finally turbulent heating augmentation ratios were determined for several body-point locations and correlated in terms of the boundary-layer momentum Reynolds number.

Hollis, Brian R.↗