Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ordered 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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

SAR backscatter from coniferous forest gaps

A study is in progress comparing Airborne Synthetic Aperture Radar (AIRSAR) backscatter from coniferous forest plots containing gaps to backscatter from adjacent gap-free plots. Issues discussed are how do gaps in the range of 400 to 1600 sq m (approximately 4-14 pixels at intermediate incidence angles) affect forest backscatter statistics and what incidence angles, wavelengths, and polarizations are most sensitive to forest gaps. In order to visualize the slant-range imaging of forest and gaps, a simple conceptual model is used. This strictly qualitative model has led us to hypothesize that forest radar returns at short wavelengths (eg., C-band) and large incidence angles (e.g., 50 deg) should be most affected by the presence of gaps, whereas returns at long wavelengths and small angles should be least affected. Preliminary analysis of 1989 AIRSAR data from forest near Mt. Shasta supports the hypothesis. Current forest backscatter models such as MIMICS and Santa Barbara Discontinuous Canopy Backscatter Model have in several cases correctly predicted backscatter from forest stands based on inputs of measured or estimated forest parameters. These models do not, however, predict within-stand SAR scene texture, or 'intrinsic scene variability' as Ulaby et al. has referred to it. For instance, the Santa Barbara model, which may be the most spatially coupled of the existing models, is not truly spatial. Tree locations within a simulated pixel are distributed according to a Poisson process, as they are in many natural forests, but tree size is unrelated to location, which is not the case in nature. Furthermore, since pixels of a simulated stand are generated independently in the Santa Barbara model, spatial processes larger than one pixel are not modeled. Using a different approach, Oliver modeled scene texture based on an hypothetical forest geometry. His simulated scenes do not agree well with SAR data, perhaps due to the simple geometric model used. Insofar as texture is the expression of biological forest processes, such as succession and disease, and physical ones, such as fire and wind-throw, it contains useful information about the forest, and has value in image interpretation and classification. Forest gaps are undoubtedly important contributors to scene variance. By studying the localized effects of gaps on forest backscatter, guided by our qualitative model, we hope to understand more clearly the manner in which spatial heterogeneities in forests produce variations in backscatter, which collectively give rise to scene texture.

Day, John L.↗
NASA NTRS1995

Satellite diversity and its implications on the RAKE receiver architecture for CDMA-based S-PCN's

In this paper we examine the applicability of RAKE receivers in a mobile LEO satellite channel and identify the potential problem areas. We then proceed to investigate the possibility of a coherent combining architecture (downlink) in the presence of satellite diversity. We closely examine the path delay difference statistics of a diversity channel and propose a delay compensation scheme for the downlink in order to reduce the complexity of the user terminal. Finally, the required modifications to the conventional RAKE receiver are proposed and discussed.

Taaghol, P.↗
NASA NTRS2001

Mg II Spectral Atlas and Flux Catalog for Late-Type Stars in the Hyades Cluster

In the course of a long-running IUE Guest Observer program, UV spectral images were obtained for more than 60 late-type members of the Hyades Cluster in order to investigate their chromospheric emissions. The emission line fluxes extracted from those observations were used to study the dependence of stellar dynamo activity upon age and rotation (IUE Observations of Rapidly Rotating Low-Mass Stars in Young Clusters: The Relation between Chromospheric Activity and Rotation). However, the details of those measurements, including a tabulation of the line fluxes, were never published. The purpose of the investigation summarized here was to extract all of the existing Hyades long-wavelength Mg II spectra in the IUE public archives in order to survey UV chromospheric emission in the cluster, thereby providing a consistent dataset for statistical and correlative studies of the relationship between stellar dynamo activity, rotation, and age over a broad range in mass.

Simon, Theodore↗
NASA NTRS2003

Linearly Forced Isotropic Turbulence

Stationary isotropic turbulence is often studied numerically by adding a forcing term to the Navier-Stokes equation. This is usually done for the purpose of achieving higher Reynolds number and longer statistics than is possible for isotropic decaying turbulence. It is generally accepted that forcing the Navier-Stokes equation at low wave number does not influence the small scale statistics of the flow provided that there is wide separation between the largest and smallest scales. It will be shown, however, that the spectral width of the forcing has a noticeable effect on inertial range statistics. A case will be made here for using a broader form of forcing in order to compare computed isotropic stationary turbulence with (decaying) grid turbulence. It is shown that using a forcing function which is directly proportional to the velocity has physical meaning and gives results which are closer to both homogeneous and non-homogeneous turbulence. Section 1 presents a four part series of motivations for linear forcing. Section 2 puts linear forcing to a numerical test with a pseudospectral computation.

Lundgren, T. S.↗
NASA NTRS2013

Detection of Extremes with AIRS and CrIS

Climate change is expected to be detected first as changes in extreme values rather than in mean values. The availability of data of from two instruments in the same orbit, AIRS data for the past eleven years and AIRS and CrIS data from the past year, provides an opportunity to evaluate this using examples of climate relevance: Desertification, seen as changes in hot extremes, severe storm, seen as a change in extremely cold clouds and the warming of the polar zone. We use AIRS to establish trends for the 1%tile, the mean and 99%tile brightness temperatures measured with the 900 cm(exp -1) channel from AIRS for the past 11 years. This channel is in the clearest part of the 11 micron atmospheric window. Substantial trends are seen for land and ocean, which in the case of the 1%tile (cold) extremes are related to the current shift of deep convection from ocean to land. Changes are also seen in the 99%tile for day tropical land, but their interpretation is at present unclear. We also see dramatic changes for the mean and 99%tile of the North Polar area. The trends are an order of magnitude larger than the instrument trend of about 3 mK/year. We use the statistical distribution from the past year derived from AIRS to evaluate the accuracy of continuing the trends established with AIRS with CrIS data. We minimize the concern about differences in the spectral response functions by limiting the analysis to the channel at 900 cm(exp -1).While the two instruments agree within 100 mK for the global day/night land/ocean mean values, there are significant differences when evaluating the1% and 99%tiles. We see a consistent warm bias in the CrIS data relative to AIRS for the 1%tile (extremely cold, cloudy) data in the tropical zone, particularly for tropical land, but the bias is not day/night land/ocean consistent. At this point the difference appears to be due to differences in the radiometric response of AIRS and CrIS to differences in the day/night land/ocean cloud types. Unless the effect can be mitigated by a future reprocessing the CrIS data, it will significantly complicate the concatenation of the AIRS and CrIS data records for the continuation of trends in extreme values.

sampling bias↗
NASA NTRSDate not supplied

Evaluating The Relationships Between Supine Propriorception Assessments with Upright Functional Mobility and Balance Tests

INTRODUCTION Upon return to Earth, spaceflight sensorimotor adaptations can result in impaired posture and locomotion [1], [2], and current exercise countermeasures in the International Space Station are not sufficient to maintain sensorimotor function. In-flight countermeasures and assessment tools for sensorimotor function are needed to mitigate risks associated with mission-critical task performance upon return to Earth or arrival to 0.38 G on Mars. One proposed countermeasure for proprioceptive deconditioning includes training on a tilt board device [3], which will be tested in an upcoming bed rest study. In preparation for the said study, the objective of this pilot study is to compare performance in supine proprioceptive assessments to performance in functional upright activities. The results are intended to guide the selection of supine assessment tools for the upcoming bedrest study and the development of future inflight proprioceptive training and assessment tools for exploration missions. METHODS Seventeen healthy participants (8 males and 9 females, 27.9 ± 8.5 years) provided informed consent as approved by the Institutional Review Board (IRB) at NASA. A horizontal air-bearing sled was used to provide a proprioceptive challenge in a supine body orientation, allowing for mediolateral motion with minimal friction during supine stance [4]. Participants were loaded axially (30 to 60% body weight) with their feet on a vertically oriented tilt board and instructed to perform supine assessment activities using custom software that displays a cursor controlled by tilting the tilt board. In a single-leg static activity on the tilt board (TB-St), performance was measured as the percent time spent on a center target during a 30-second trial. During a two-feet dynamic activity on the tilt board (TB-Dy), performance was measured as the number of targets captured in 30 seconds. These supine tasks were compared to performance during two upright activities used to represent upright performance: completion time for a functional mobility task (FMT) [2],[5] and a performance score from force plate sway data during 30-second single-leg upright standing (USL). RESULTS AND DISCUSSION A Spearman's rank-order correlation was run to assess the relationship between upright and supine assessments. There were statistically significant, strong correlations between scores for the FMT and each TB activity, including TB-St, rs(17) = 0.50, p<0.05, TB-Dy3, rs(17) = 0.73, p<0.001, and TB-Dy9, rs(15) = 0.53, p<0.05. There were no statistically significant correlations between USL and each TB activity. Despite a lack of correlation with upright single-leg balance, the ability of the TB to indicate upright performance in the FMT is promising, as this task is a standard post-flight measure designed to characterize locomotor dysfunction [2], [5]. CONCLUSION These efforts will inform the selection of an appropriate sensorimotor assessment method for the upcoming bed rest study and similar future studies. While current work uses a bed rest analog to develop these technologies, future work aims to prepare capabilities for future in-flight sensorimotor training and assessment to mitigate the risks of proprioception impairment after spaceflight. ACKNOWLEDGEMENTS This work was supported by the University Space Research Association Internship Program. REFERENCES [1] Mulavara A.P. et al (2018) Med Sci Sports Exerc, 50(9), 1961-1980. [2] Mulavara A.P. et al (2010) Exp Brain Res, 202(3), 649-659. [3] Macaulay T.R. et al (2021) Front. Syst. Neurosci., 15, 658-985. [4] Goel R. et al (2017) Front. Syst. Neurosci, 11. [5] Koppelmans V. et al (2013) BMC Neurology, 12.

R Bellisle↗
DOE OSTI2023

Developing and applying quantifiable metrics for diagnostic and experiment design on Z

This project applies methods in Bayesian inference and modern statistical methods to quantify the value of new experimental data, in the form of new or modified diagnostic configurations and/or experiment designs. We demonstrate experiment design methods that can be used to identify the highest priority diagnostic improvements or experimental data to obtain in order to reduce uncertainties on critical inferred experimental quantities and select the best course of action to distinguish between competing physical models. Bayesian statistics and information theory provide the foundation for developing the necessary metrics, using two high impact experimental platforms on Z as exemplars to develop and illustrate the technique. We emphasize that the general methodology is extensible to new diagnostics (provided synthetic models are available), as well as additional platforms. We also discuss initial scoping of additional applications that began development in the last year of this LDRD.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Propagation of partially spatially coherent laser beams in instantaneous Kerr media

The propagation of intense, partially spatially coherent laser beams in a medium with instantaneous third-order susceptibility is studied analytically and numerically. For sufficiently high power relative to that required for nonlinear self-focusing, the propagation initially proceeds in two stages. In the first stage, spatial coherence builds up, and in the second stage, the number of speckles reduces. Once the degree of coherence is sufficiently high, whole-beam self-focusing occurs. The beam power is mostly confined within the initial spot radius. Two analytical approaches for describing the evolution of the beam are presented. The method of moments leads to an analytical solution for the rms spot radius that is in excellent agreement with simulations. This method does not require any knowledge of the field statistics beyond the initial conditions and provides no information about the evolution of the individual speckles. The other approach employs a self-similar solution for the second-order coherence function of the field and assumes that the fourth-order coherence function is factorizable and obeys complex circular Gaussian random statistics. The latter method also leads to an analytical expression for the spot radius, but its predictions for the qualitative evolution of the speckles disagree with wave-optics simulations.

lasers↗
DOE OSTI2023

Combining cosmic shear data with correlated photo-z uncertainties: constraints from DESY1 and HSC-DR1

An accurate calibration of the source redshift distribution p(z) is a key aspect in the analysis of cosmic shear data. This, one way or another, requires the use of spectroscopic or high-quality photometric samples. However, the difficulty to obtain colour-complete spectroscopic samples matching the depth of weak lensing catalogs means that the analyses of different cosmic shear datasets often use the same samples for redshift calibration. This introduces a source of statistical and systematic uncertainty that is highly correlated across different weak lensing datasets, and which must be accurately characterised and propagated in order to obtain robust cosmological constraints from their combination. In this paper we introduce a method to quantify and propagate the uncertainties on the source redshift distribution in two different surveys sharing the same calibrating sample. The method is based on an approximate analytical marginalisation of the p(z) statistical uncertainties and the correlated marginalisation of residual systematics. We apply this method to the combined analysis of cosmic shear data from the DESY1 data release and the HSC-DR1 data, using the COSMOS 30-band catalog as a common redshift calibration sample. We find that, although there is significant correlation in the uncertainties on the redshift distributions of both samples, this does not change the final constraints on cosmological parameters significantly. The same is true also for the impact of residual systematic uncertainties from the errors in the COSMOS 30-band photometric redshifts. Additionally, we show that these effects will still be negligible in Stage-IV datasets. Finally, the combination of DESY1 and HSC-DR1 allows us to constrain the "clumpiness" parameter to S 8 = ${0.768}_{-0.017}^{+0.021}$. This corresponds to a ~√(2) improvement in uncertainties with respect to either DES or HSC alone.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Direct cross section measurement of 102 Pd ⁢(𝛾,𝑝) and 102 Pd ⁢(𝛾,𝛼) for the astrophysical 𝑝 process

Background: A handful of neutron-deficient stable nuclei, known as the “p nuclei,” cannot be produced through astrophysical neutron capture processes. Instead, some of these nuclei are proposed to be produced by 𝛾-induced reactions on existing r- and s-process seeds. The specific astrophysical site or sites are not yet identified, however, with uncertainties in the cross sections of these 𝛾-induced reactions playing a role. Databases of reaction rates for astrophysical simulations often rely on theoretical statistical model calculations, such as Hauser-Feshbach, for rates where no experimental information is known. However, reasonable variations in the choice of parametrizations of various nuclear properties can create order-of-magnitude variations in the final predicted cross sections and reaction rates, which are then propagated through the models to the predicted final abundances. Purpose: To better constrain these statistical model calculations and ultimately reduce the uncertainties from the nuclear physics on our understanding of the p nuclei, a measurement of the cross sections of 𝛾-induced reactions on the p-nucleus 102 Pd was undertaken. This work represents the first measurement of its kind, using segmented silicon detectors to measure prompt charged particle emission from 𝛾-induced reactions. Methods: Quasimonoenergetic gamma beams from the High Intensity 𝛾 Source facility bombarded an enriched 102 Pd target. A segmented silicon array was arranged to detect the particles emitted from (𝛾,𝑝) and (𝛾,𝛼) reactions. Results: Reaction cross sections were deduced at multiple 𝛾-beam energies between 10 and 19 MeV, and compared to statistical model calculations using talys-1.96. The 102 Pd ⁢(𝛾,𝑝)⁢ 101 Rh reaction cross section was reasonably well reproduced by a subset of photon strength functions and level densities, though the strength to the ground state of 101 Rh was underestimated at higher incident gamma energies. The 102 Pd ⁢(𝛾,𝛼)⁢ 98 Ru was in general overpredicted by the various alpha-nucleus optical model potentials. Conclusions: While the theoretical cross sections used to model the (𝛾,𝑝) reactions for the p process may be reasonable, a more careful approach is needed in the case of (𝛾,𝛼). Further work to probe gamma-induced reaction cross sections at and near the p nuclei is warranted.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1994

Voids and constraints on nonlinear clustering of galaxies

Void statistics of the galaxy distribution in the Center for Astrophysics Redshift Survey provide strong constraints on galaxy clustering in the nonlinear regime, i.e., on scales R equal to or less than 10/h Mpc. Computation of high-order moments of the galaxy distribution requires a sample that (1) densely traces the large-scale structure and (2) covers sufficient volume to obtain good statistics. The CfA redshift survey densely samples structure on scales equal to or less than 10/h Mpc and has sufficient depth and angular coverage to approach a fair sample on these scales. In the nonlinear regime, the void probability function (VPF) for CfA samples exhibits apparent agreement with hierarchical scaling (such scaling implies that the N-point correlation functions for N greater than 2 depend only on pairwise products of the two-point function xi(r)) However, simulations of cosmological models show that this scaling in redshift space does not necessarily imply such scaling in real space, even in the nonlinear regime; peculiar velocities cause distortions which can yield erroneous agreement with hierarchical scaling. The underdensity probability measures the frequency of 'voids' with density rho less than 0.2 -/rho. This statistic reveals a paucity of very bright galaxies (L greater than L asterisk) in the 'voids.' Underdensities are equal to or greater than 2 sigma more frequent in bright galaxy samples than in samples that include fainter galaxies. Comparison of void statistics of CfA samples with simulations of a range of cosmological models favors models with Gaussian primordial fluctuations and Cold Dark Matter (CDM)-like initial power spectra. Biased models tend to produce voids that are too empty. We also compare these data with three specific models of the Cold Dark Matter cosmogony: an unbiased, open universe CDM model (omega = 0.4, h = 0.5) provides a good match to the VPF of the CfA samples. Biasing of the galaxy distribution in the 'standard' CDM model (omega = 1, b = 1.5; see below for definitions) and nonzero cosmological constant CDM model (omega = 0.4, h = 0.6 lambda(sub 0) = 0.6, b = 1.3) produce voids that are too empty. All three simulations match the observed VPF and underdensity probability for samples of very bright (M less than M asterisk = -19.2) galaxies, but produce voids that are too empty when compared with samples that include fainter galaxies.

Vogeley, Michael S.↗
NASA NTRS1981

Temporal correlation measurements of pulsed dual CO2 lidar returns

A pulsed dual-laser direct-detection differential-absorption lidar DIAL system, operating near 10.6 microns, is used to measure the temporal correlation and statistical properties of backscattered returns from specular and diffuse topographic targets. Results show that atmospheric-turbulence fluctuations can effectively be frozen for pulse separation times on the order of 1-3 msec or less. The diffuse target returns, however, yielded a much lower correlation than that obtained with the specular targets; this being due to uncorrelated system noise effects and different statistics for the two types of target returns.

Menyuk, N.↗
NASA NTRS1984

GASP cloud encounter statistics - Implications for laminar flow control flight

The cloud observation archive from the NASA Global Atmospheric Sampling Program (GASP) is analyzed in order to derive the probability of cloud encounter at altitudes normally flown by commercial airliners, for application to a determination of the feasability of Laminar Flow Control (LFC) on long-range routes. The probability of cloud encounter is found to vary significantly with season. Several meteorological circulation features are apparent in the latitudinal distribution of cloud cover. The cloud encounter data are shown to be consistent with the classical midlatitude cyclone model with more clouds encountered in highs than in lows. Aircraft measurements of route-averaged time-in-clouds fit a gamma probability distribution model which is applied to estimate the probability of extended cloud encounter, and the associated loss of LFC effectiveness along seven high-density routes. The probability is demonstrated to be low.

Jasperson, W. H.↗
NASA NTRS2018

More Efficient Precipitation Rate Data Acquisition with the STORM Swath-Based Analysis Tool

A multitude of locations exist for acquiring satellite-derived precipitation rate data. Because the volume of data exceeds what can reasonably be analyzed manually, online tools have been developed to help researchers make more informed data acquisition decisions. Some interfaces provide sample images generated from prioritized variables to enable efficient visual inspection, while others include the ability to subset geographically to reduce the volume of the data to the area of interest. Building upon such value-added services is the STORM (the NASA Precipitation Processing System (PPS) data product ordering system) Swath-Based Analysis Tool (https://storm.pps.eosdis.nasa.gov/storm/Analysis.jsp). With it, users can compare overflights of multiple precipitation satellites simultaneously and from this make determinations about which swaths will best serve their research. They can preview these overflights with two distinct visualization tools. Finally, they can then download geographically subset swaths that have precipitation rate characteristics matching their research needs. A user starts by making up to six selections. They then select their date range of interest going back to 1997. Finally, they select a geographical region either by inputting latitude-longitude corners of a rectangle or drawing that rectangle on an interactive map. The system processes each file, determining aggregate statistics about the precipitation rates within the selected region as observed by each instrument. These statistical values are added to a time series graph as a "time remaining" counter is continuously updated. The chart is fully interactive, with detailed values displayed upon mouseover, and numerous options available upon clicking each point. Since the user is likely to need the relevant files for their own research, they can generate an order from their request. All files are geographically subset to the region of interest, and the user can use the statistical values to limit which files are added to the order. This flexibility helps to ensure researchers are only downloading the files they need for their work, making for a more efficient process both from a time-consumption and a storage perspective. This paper will introduce the Swath-Based Analysis Tool and explore its numerous features. It will discuss several use cases for precipitation scientists.

remote sensing↗
DOE OSTI2021

Identification of Life Models for Li-Ion Batteries Using Penalized Regression and Bilevel Optimization

Reduced-order physics-based life models are extremely useful for rapidly predicting battery state-of-health and for simulating battery lifetime in arbitrary aging conditions. However, identification of well-parameterized models is difficult. This is because, for maximum usefulness in predicting lifetime under a variety of conditions, aging test data exhibits many degradation mechanisms, which all need to be accurately modeled. However, because aging tests are time-consuming and expensive, especially for large-format batteries, a minimum of tests are conducted while probing many stress factors. Building a well-parameterized model is then very challenging: an under-parameterized model will neglect critical degradation modes, and an over-parameterized model will extrapolate poorly to new testing conditions. To complicate this matter, the functional form of the model for any individual degradation rate can be very difficult to identify. In this work, the statistical tools of penalized regression and bilevel optimization are used to help identify both the functional forms of and optimize the parameters of reduced-order life models, accelerating identification of robust models. Model robustness is demonstrated through traditional statistics methods of cross-validation and sensitivity analysis, uncertainty quantification through bootstrap resampling and Monte-Carlo simulation, and simulation of real-world use cases.

47 OTHER INSTRUMENTATION↗
NASA NTRS1978

A statistical technique for determining rainfall over land employing Nimbus-6 ESMR measurements

At 37 GHz, the frequency at which the Nimbus 6 Electrically Scanning Microwave Radiometer (ESMR 6) measures upwelling radiance, it was shown theoretically that the atmospheric scattering and the relative independence on electromagnetic polarization of the radiances emerging from hydrometers make it possible to monitor remotely active rainfall over land. In order to verify experimentally these theoretical findings and to develop an algorithm to monitor rainfall over land, the digitized ESMR 6 measurements were examined statistically. Horizontally and vertically polarized brightness temperature pairs (TH, TV) from ESMR 6 were sampled for areas of rainfall over land as determined from the rain recording stations and the WSR 57 radar, and areas of wet and dry ground (whose thermodynamic temperatures were greater than 5 C) over the Southeastern United States. These three categories of brightness temperatures were found to be significantly different in the sense that the chances that the mean vectors of any two populations coincided were less than 1 in 100.

Rodgers, E.↗
DOE OSTI2022

Priors on red galaxy stochasticity from hybrid effective field theory

ABSTRACT We investigate the stochastic properties of typical red galaxy samples in a controlled numerical environment. We use halo occupation distribution (HOD) modelling to create mock realizations of three separate bright red galaxy samples consistent with data sets used for clustering and lensing analyses in modern galaxy surveys. Second-order Hybrid Effective Field Theory (HEFT) is used as a field-level forward model to describe the full statistical distribution of these tracer samples, and their stochastic power spectra are directly measured and compared to the Poisson shot-noise prediction. While all of the galaxy samples we consider are hosted within haloes with sub-Poisson stochasticity, we observe that the galaxy samples themselves possess stochasticities that range from sub-Poisson to super-Poisson, in agreement with predictions from the halo model. As an application of our methodology, we place priors on the expected degree of non-Poisson stochasticity in cosmological analyses using such samples. We expect these priors will be useful in reducing the complexity of the full parameter space for future analyses using second-order Lagrangian bias models. More generally, the techniques outlined here present the first application of HEFT methods to characterize models of the galaxy–halo connection at the field level, revealing new connections between once-disparate modelling frameworks.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

Empirical convection models for northward IMF

It is clear that polar cap convection during times of northward Interplanetary Magnetic Field (IMF) is more structured and of lower mean speed than at times of southward IMF. This, coupled with the fact that the polar cap is smaller, means that empirical models are more difficult to construct with certainty. It is also clear that sunward flow deep in the polar cap is often observed, but its connection with the rest of the flow pattern is controversial. At present, empirical models are of three types: 'statistical' models wherein data from different days but with similar IMF conditions are averaged together; 'pattern recognition' models, which are built up by examining individually hundreds of passes to derive a 'typical' pattern which embodies features frequently observed; and 'assimilative' models, which use data of different types and from as many locations as possible, but all taken at the same time, in order to derive a snapshot (or series of snapshots) of the entire pattern. Each type of model has its own difficulties. Statistical models, by their very nature, smooth out flow features (e.g. the convection reversal, and the locus of sunward flow deep in the polar cap) which are not found at precisely the same invariant latitudes and magnetic local times on different days. Pattern recognition models are better at reproducing small-scale features, but the large-scale pattern can be a matter of interpretation. Assimilative models (such as AMIE) hold out the best hope for creating instantaneous, global convection patterns; however, the analysis technique tends to be most irregular (and least reliable) in the regions which are not well covered by in situ data. It appears that, at least at times, a four cell model with sunward flow at the highest and lowest latitudes, and antisunward flow in between, is consistent with the observations. At other times, the observations may be consistent with a two-cell convection pattern, but which includes significant meanders within the polar cap.

Moses, Julie J.↗