Search NASA⌕ Search

SEARCH · Search NASA

Results for “extended least squares”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Tracking and Data Relay Satellite (TDRS) Orbit Estimation Using an Extended Kalman Filter

Alternatives to the Tracking and Data Relay Satellite (TDRS) orbit estimation procedure were studied to develop a technique that both produces more reliable results and is more amenable to automation than the prior procedure. The Earth Observing System (EOS) Terra mission has TDRS ephemeris prediction 3(sigma) requirements of 75 meters in position and 5.5 millimeters per second in velocity over a 1.5-day prediction span. Meeting these requirements sometimes required reruns of the prior orbit determination (OD) process, with manual editing of tracking data to get an acceptable solution. After a study of the available alternatives, the Flight Dynamics Facility (FDF) began using the Real-Time Orbit Determination (RTOD(Registered TradeMark)) Kalman filter program for operational support of TDRSs in February 2007. This extended Kalman filter (EKF) is used for daily support, including within hours after most thrusting, to estimate the spacecraft position, velocity, and solar radiation coefficient of reflectivity (C(sub R)). The tracking data used are from the Bilateration Ranging Transponder System (BRTS), selected TDRS System (TDRSS) User satellite tracking data, and Telemetry, Tracking, and Command (TT&C) data. Degraded filter results right after maneuvers and some momentum unloads provided incentive for a hybrid OD technique. The results of combining EKF strengths with the Goddard Trajectory Determination System (GTDS) Differential Correction (DC) program batch-least-squares solutions, as recommended in a 2005 paper on the chain-bias technique, are also presented.

Ward, Douglas T.↗
NASA NTRS1994

Separation of lithospheric, external, and core components of the south polar geomagnetic field at satellite altitudes

We present a new approach to producing scalar Magsat magnetic anomaly maps based on correlation coefficient filtering and the use of almost all of the available orbits. Our method differs from earlier techniques with respect to the following: (1) Passes are selected based on their variance properties rather than planetary indices such as Kp. (2) The core field model is least squares fit to individual passes and subsequently removed instead of substracting the model directly. This technique replaces band pass filtering or polynomial trend removal methods. (3) Each selected pass is sorted geographically and by local time, placed into one of four different altitude bands, and correlation coefficient filtered with the two adjacent passes. The filtering is the second step toward isolating the static lithospheric signal from the more dynamic external field signals. (4) Least squares collocation is used to grid the correlated passes; subsequently, the dawn and dusk maps are also correlation filtered providing another step toward removal of external fields. (5) The four resultant total field maps are continued to a common altitude and again correlation filtered for the static lithospheric anomalies. (6) The filtered results are then averaged together to provide a new total field map of the lithosphere south of 40 deg S latitude. Our total field map differs from previous efforts over the crustal blocks of West Antarctica. We obtained a positive anomaly over Edward VII Peninsula, extending into the Byrd subglacial basin and obtained a negative anomaly over the Ellsworth Mountains and parts of the Byrd subglacial basin. Also, a positive anomaly extending from the Ross Sea to offshore Wilkes Land is present in our map; however, this feature is absent in other maps. Positive anomalies marking the Weddell Sea in previous efforts are not present in our map. Prominent external field anomalies in the quadrant offshore of Wilkes Land are present in all previous efforts; however, these signatures are reduced in our total field map. This map significantly advances efforts for the identification and interpretation of magnetic anomalies of the south polar lithosphere.

Alsdorf, Douglas E.↗
NASA NTRS2014

LightForce Photon-Pressure Collision Avoidance: Updated Efficiency Analysis Utilizing a Highly Parallel Simulation Approach

This paper provides an updated efficiency analysis of the LightForce space debris collision avoidance scheme. LightForce aims to prevent collisions on warning by utilizing photon pressure from ground based, commercial off the shelf lasers. Past research has shown that a few ground-based systems consisting of 10 kilowatt class lasers directed by 1.5 meter telescopes with adaptive optics could lower the expected number of collisions in Low Earth Orbit (LEO) by an order of magnitude. Our simulation approach utilizes the entire Two Line Element (TLE) catalogue in LEO for a given day as initial input. Least-squares fitting of a TLE time series is used for an improved orbit estimate. We then calculate the probability of collision for all LEO objects in the catalogue for a time step of the simulation. The conjunctions that exceed a threshold probability of collision are then engaged by a simulated network of laser ground stations. After those engagements, the perturbed orbits are used to re-assess the probability of collision and evaluate the efficiency of the system. This paper describes new simulations with three updated aspects: 1) By utilizing a highly parallel simulation approach employing hundreds of processors, we have extended our analysis to a much broader dataset. The simulation time is extended to one year. 2) We analyze not only the efficiency of LightForce on conjunctions that naturally occur, but also take into account conjunctions caused by orbit perturbations due to LightForce engagements. 3) We use a new simulation approach that is regularly updating the LightForce engagement strategy, as it would be during actual operations. In this paper we present our simulation approach to parallelize the efficiency analysis, its computational performance and the resulting expected efficiency of the LightForce collision avoidance system. Results indicate that utilizing a network of four LightForce stations with 20 kilowatt lasers, 85% of all conjunctions with a probability of collision Pc > 10 (sup -6) can be mitigated.

conjunction analysis↗
NASA NTRS2003

Termolecular ion-molecule reactions in Titan's atmosphere. IV. A search made at up to 1 micron in pure hydrocarbons

The results of a study of ion-molecule reactions occurring in pure methane, acetylene, ethylene, ethane, propyne, propene, propane, and diacetylene at pressures up to 40 microns of pressure are reported. A variety of experimental methods are used: The standard double resonance in an ICR, for determination of the precursor ions and the modulated double resonance ejection in an ICR, for the determination of the daughter ions. The FA-SIFT technique was used for validation and examination of termolecular reactions with rate coefficients that are less than 10(-26) cm(6) s(-1). An extensive database of reaction kinetics already exists for many of these reactions. The main point of this study was the determination of the accuracy of this database and to search for any missing reactions and reaction channels that may have been omitted from earlier investigations. A specific objective of this work was to extend the study to the highest pressures possible to find out if there were any important termolecular reaction channels occurring. A new approach was used here. In the pure hydrocarbon gases the mass spectra were followed as a function of the pressure changes of the gas. An initial guess was first made using the current literature as a source of the reaction kinetics that were expected. A model of the ion abundances was produced from the solution of the partial differential equations in terms of reaction rate coefficients and initial abundances. The experimental data was fitted to the model for all of the pressures by a least squares minimization to the reaction rate coefficients and initial abundances. The reaction rate coefficients obtained from the model were then compared to the literature values. Several new channels and reactions were discovered when the modeled fits were compared to the actual data. This is all explained in the text and the implications of these results are discussed for the Titan atmosphere.

Hydrocarbons, Acyclic/analysis↗
NASA NTRS2000

Preliminary development and evaluation of an algae-based air regeneration system

The potential of air regeneration system based on the growth of microalgae on the surface of porous ceramic tubes is evaluated. The algae have been maintained in the system for extended periods, up to 360 days. Preliminary measurements of the photosynthetic capacity have been made for Chlorella vulgaris (UTEX 259), Neospongiococcum punctatum (UTEX 786), Stichococcus sp., and Gloeocapsa sp. Under standard test conditions (photosynthetic photon flux approximately 66 micromoles m-2 s-1, initial CO2 concentration approximately 450 micromoles mol-1), mature tubes remove up to 0.2 micromoles of CO2 per tube per minute. The rate of removal increases with photon flux up to at least 225 micromoles m-2 s-1 (PPF); peak rates of 0.35 micromoles of CO2 per tube per minute have been achieved with Chlorella vulgaris. These rates correspond to between 120 and 210 micromoles of CO2 removed per square meter of projected area per minute.

Non-NASA Center↗
NASA NTRS2017

Validation of SMAP Surface Soil Moisture Products with Core Validation Sites

The NASA Soil Moisture Active Passive (SMAP) mission has utilized a set of core validation sites as the primary methodology in assessing the soil moisture retrieval algorithm performance. Those sites provide well calibrated in situ soil moisture measurements within SMAP product grid pixels for diverse conditions and locations.The estimation of the average soil moisture within the SMAP product grid pixels based on in situ measurements is more reliable when location specific calibration of the sensors has been performed and there is adequate replication over the spatial domain, with an up-scaling function based on analysis using independent estimates of the soil moisture distribution. SMAP fulfilled these requirements through a collaborative CalVal Partner program.This paper presents the results from 34 candidate core validation sites for the first eleven months of the SMAP mission. As a result of the screening of the sites prior to the availability of SMAP data, out of the 34 candidate sites 18 sites fulfilled all the requirements at one of the resolution scales (at least). The rest of the sites are used as secondary information in algorithm evaluation. The results indicate that the SMAP radiometer-based soil moisture data product meets its expected performance of 0.04 cu m/cu m volumetric soil moisture (unbiased root mean square error); the combined radar-radiometer product is close to its expected performance of 0.04 cu m/cu m, and the radar-based product meets its target accuracy of 0.06 cu m/cu m (the lengths of the combined and radar-based products are truncated to about 10 weeks because of the SMAP radar failure). Upon completing the intensive CalVal phase of the mission the SMAP project will continue to enhance the products in the primary and extended geographic domains, in co-operation with the CalVal Partners, by continuing the comparisons over the existing core validation sites and inclusion of candidate sites that can address shortcomings.

SMAP↗
NASA NTRS1986

Design of adaptive control systems by means of self-adjusting transversal filters

The design of closed-loop adaptive control systems based on nonparametric identification was addressed. Implementation is by self-adjusting Least Mean Square (LMS) transversal filters. The design concept is Model Reference Adaptive Control (MRAC). Major issues are to preserve the linearity of the error equations of each LMS filter, and to prevent estimation bias that is due to process or measurement noise, thus providing necessary conditions for the convergence and stability of the control system. The controlled element is assumed to be asymptotically stable and minimum phase. Because of the nonparametric Finite Impulse Response (FIR) estimates provided by the LMS filters, a-priori information on the plant model is needed only in broad terms. Following a survey of control system configurations and filter design considerations, system implementation is shown here in Single Input Single Output (SISO) format which is readily extendable to multivariable forms. In extensive computer simulation studies the controlled element is represented by a second-order system with widely varying damping, natural frequency, and relative degree.

Merhav, S. J.↗
NASA NTRS1993

Comparison of TOPEX/Poseidon orbit determination solutions obtained by the Goddard Space Flight Center Flight Dynamics Division and Precision Orbit Determination Teams

Orbit determination results are obtained by the Goddard Space Flight Center (GSFC) Flight Dynamics Division (FDD) using the Goddard Trajectory Determination System (GTDS) and a real-time extended Kalman filter estimation system to process Tracking Data and Relay Satellite (TDRS) System (TDRSS) measurements in support of the Ocean Topography Experiment (TOPEX)/Poseidon spacecraft navigation and health and safety operations. GTDS is the operational orbit determination system used by the FDD, and the extended Kalman fliter was implemented in an analysis prototype system, the Real-Time Orbit Determination System/Enhanced (RTOD/E). The Precision Orbit Determination (POD) team within the GSFC Space Geodesy Branch generates an independent set of high-accuracy trajectories to support the TOPEX/Poseidon scientific data. These latter solutions use the Geodynamics (GEODYN) orbit determination system with laser ranging tracking data. The TOPEX/Poseidon trajectories were estimated for the October 22 - November 1, 1992, timeframe, for which the latest preliminary POD results were available. Independent assessments were made of the consistencies of solutions produced by the batch and sequential methods. The batch cases were assessed using overlap comparisons, while the sequential cases were assessed with covariances and the first measurement residuals. The batch least-squares and forward-filtered RTOD/E orbit solutions were compared with the definitive POD orbit solutions. The solution differences were generally less than 10 meters (m) for the batch least squares and less than 18 m for the sequential estimation solutions. The differences among the POD, GTDS, and RTOD/E solutions can be traced to differences in modeling and tracking data types, which are being analyzed in detail.

Doll, C.↗
NASA NTRS1996

Dependence of the Normalized Radar Cross Section of Water Waves on Bragg Wavelength-Wind Speed Sensitivity

Measurements of the normalized radar cross section (sigma(sup o)) made by the YSCAT ultrawideband scatterometer during an extended deployment on the Canada Centre for Inland Waters(CCIW) Research Tower located at Lake Ontario are analyzed and compared with anemometer wind measurements to study the sensitivity of (sigma(sup o)) to the wind speed as a function of the Bragg wavelength. This paper concentrates on upwind and downwind azimuth angles in the wind speed range of 4.5-12 m/s. While YSCAT collected measurements of sigma(sup o) at a variety of frequencies and incidence angles, this paper focuses on frequencies of 2.0, 3.05, 5.30, 10.02, and 14.0 GHz and incidence angles within the Bragg regime, 30-50 deg. Adopting a power law model to describe the relationship between sigma(sup o) and wind speed, both wind speed exponents and upwind/downwind (u/d) ratios of sigma(sup o) are found using least squares linear regression. The analysis of the wind speed exponents and u/d ratios show that shorter Bragg wavelengths (Lambda less than 4 cm) are the most sensitive to wind speed and direction. Additionally, vertical polarization (V-pol) sigma(sup o) is shown to be more sensitive to wind speed than horizontal polarization (H-pol) sigma(sup o), while the H-pol u/d ratio is larger than the V-pol u/d ratio.

Long, David G.↗
NASA NTRS1975

A quantitative magnetospheric model derived from spacecraft magnetometer data

The model is derived by making least squares fits to magnetic field measurements from four Imp satellites. It includes four sets of coefficients, representing different degrees of magnetic disturbance as determined by the range of Kp values. The data are fit to a power series expansion in the solar magnetic coordinates and the solar wind-dipole tilt angle, and thus the effects of seasonal north-south asymmetries are contained. The expansion is divergence-free, but unlike the usual scalar potential expansion, the model contains a nonzero curl representing currents distributed within the magnetosphere. The latitude at the earth separating open polar cap field lines from field lines closing on the day side is about 5 deg lower than that determined by previous theoretically derived models. At times of high Kp, additional high-latitude field lines extend back into the tail. Near solstice, the separation latitude can be as low as 75 deg in the winter hemisphere. The average northward component of the external field is much smaller than that predicted by theoretical models; this finding indicates the important effects of distributed currents in the magnetosphere.

Mead, G. D.↗
NASA NTRS1993

A solution scheme for the Euler equations based on a multi-dimensional wave model

A scheme for the solution of scalar advection on an unstructured mesh has been developed, tested, and extended to the Euler equations. The scheme preserves a linear function exactly, and yields nearly monotone results. The flux function associated with the Euler scheme is based on a discrete 'wave model' for the system of equations. The wave model decomposes the solution gradient at a location into shear waves, entropy waves and acoustic waves and calculates the speeds, strengths and directions associated with the waves. The approach differs from typical flux-difference splitting schemes in that the waves are not assumed to propagate normal to the faces of the control volumes; directions of propagation of the waves are instead computed from solution-gradient information. Results are shown for three test cases, and two different wave models. The results are compared to those from other approaches, including MUSCL and Galerkin least squares schemes.

Powell, Kenneth G.↗
NASA NTRS2022

The Spectral-Timing Analysis of Cygnus X-1 With Insight-HXMT

Cygnus X-1, as the first discovered black hole binary, is a key source for understanding the mechanisms of state transitions, and the scenarios of accretion in extreme gravity fields. We present a spectral-timing analysis of observations taken with the Insight-HXMT mission, focusing on the spectral-state dependent timing properties in the broad energy range of 1–150 keV, thus extending previous RXTE-based studies to both lower and higher energies. Our main results are the following: a) We successfully use a simple empirical model to fit all spectra, confirming that the reflection component is stronger in the soft state than in the hard state; b) The evolution of the total fractional root mean square (rms) depends on the selected energy band and the spectral shape, which is a direct result of the evolution of the power spectral densities (PSDs); c) In the hard/intermediate state, we see clear short-term variability features and a positive correlation between central frequencies of the variability components and the soft photon index Γ(sub 1), also at energies above 15 keV. The power spectrum is dominated by red noise in the soft state instead. These behaviors can be traced to at least 90 keV; d) The coherence and the phase-lag spectra show different behaviors dependent on different spectral shapes.

Binary X-rays↗
NASA NTRS2005

Experimental and Analytical Determination of the Geometric Far Field for Round Jets

An investigation was conducted at the NASA Glenn Research Center using a set of three round jets operating under unheated subsonic conditions to address the question: "How close is too close?" Although sound sources are distributed at various distances throughout a jet plume downstream of the nozzle exit, at great distances from the nozzle the sound will appear to emanate from a point and the inverse-square law can be properly applied. Examination of normalized sound spectra at different distances from a jet, from experiments and from computational tools, established the required minimum distance for valid far-field measurements of the sound from subsonic round jets. Experimental data were acquired in the Aeroacoustic Propulsion Laboratory at the NASA Glenn Research Center. The WIND computer program solved the Reynolds-Averaged Navier-Stokes equations for aerodynamic computations; the MGBK jet-noise prediction computer code was used to predict the sound pressure levels. Results from both the experiments and the analytical exercises indicated that while the shortest measurement arc (with radius approximately 8 nozzle diameters) was already in the geometric far field for high-frequency sound (Strouhal number >5), low-frequency sound (Strouhal number <0.2) reached the geometric far field at a measurement radius of at least 50 nozzle diameters because of its extended source distribution.

Koch, L. Danielle↗
NASA NTRS2005

Adaptive Filtering Using Recurrent Neural Networks

A method for adaptive (or, optionally, nonadaptive) filtering has been developed for estimating the states of complex process systems (e.g., chemical plants, factories, or manufacturing processes at some level of abstraction) from time series of measurements of system inputs and outputs. The method is based partly on the fundamental principles of the Kalman filter and partly on the use of recurrent neural networks. The standard Kalman filter involves an assumption of linearity of the mathematical model used to describe a process system. The extended Kalman filter accommodates a nonlinear process model but still requires linearization about the state estimate. Both the standard and extended Kalman filters involve the often unrealistic assumption that process and measurement noise are zero-mean, Gaussian, and white. In contrast, the present method does not involve any assumptions of linearity of process models or of the nature of process noise; on the contrary, few (if any) assumptions are made about process models, noise models, or the parameters of such models. In this regard, the method can be characterized as one of nonlinear, nonparametric filtering. The method exploits the unique ability of neural networks to approximate nonlinear functions. In a given case, the process model is limited mainly by limitations of the approximation ability of the neural networks chosen for that case. Moreover, despite the lack of assumptions regarding process noise, the method yields minimum- variance filters. In that they do not require statistical models of noise, the neural- network-based state filters of this method are comparable to conventional nonlinear least-squares estimators.

Parlos, Alexander G.↗
NASA NTRS2002

Pulsed Electromagnetic Acceleration of Plasma: A Review

Much have been learned in the acceleration mechanisms involved in accelerating a plasma electromagnetically in the laboratory over the last 40 years since the early review by Winston Bostik of 1963, but the accumulated understanding is very much scattered throughout the literature. This literature extends back at least to the early sixties and includes Rosenbluth's snowplow model, discussions by Ralph Lovberg, Colgate's boundary-layer model of a current sheet, many papers from the activity at Columbia by Robert Gross and his colleagues, and the relevant, 1-D unsteady descriptions developed from the U. of Maryland theta-pinch studies. Recent progress on the understanding of the pulsed penetration of magnetic fields into collisionless or nearly collisionless plasmas are also be reviewed. Somewhat more recently, we have the two-dimensional, unsteady results in the collisional regime associated with so-called wall-instability in large radius pinch discharges and also in coaxial plasma guns (e.g., Plasma Flow Switch). Among other things, for example, we have the phenomenon of a high- density plasma discharge propagating in a cooaxial gun as an apparently straight sheet (vs paraboloid) because mass re-distribution (on a microsecond timescale) compensates for the 1/r- squared variation of magnetic pressure. We will attempt to collate some of this vast material and bring some coherence tc the development of the subject.

Thio, Y. C. Francis↗
NASA NTRS2002

Computational Thermochemistry of Jet Fuels and Rocket Propellants

The design of new high-energy density molecules as candidates for jet and rocket fuels is an important goal of modern chemical thermodynamics. The NASA Glenn Research Center is home to a database of thermodynamic data for over 2000 compounds related to this goal, in the form of least-squares fits of heat capacities, enthalpies, and entropies as functions of temperature over the range of 300 - 6000 K. The chemical equilibrium with applications (CEA) program written and maintained by researchers at NASA Glenn over the last fifty years, makes use of this database for modeling the performance of potential rocket propellants. During its long history, the NASA Glenn database has been developed based on experimental results and data published in the scientific literature such as the standard JANAF tables. The recent development of efficient computational techniques based on quantum chemical methods provides an alternative source of information for expansion of such databases. For example, it is now possible to model dissociation or combustion reactions of small molecules to high accuracy using techniques such as coupled cluster theory or density functional theory. Unfortunately, the current applicability of reliable computational models is limited to relatively small molecules containing only around a dozen (non-hydrogen) atoms. We propose to extend the applicability of coupled cluster theory- often referred to as the 'gold standard' of quantum chemical methods- to molecules containing 30-50 non-hydrogen atoms. The centerpiece of this work is the concept of local correlation, in which the description of the electron interactions- known as electron correlation effects- are reduced to only their most important localized components. Such an advance has the potential to greatly expand the current reach of computational thermochemistry and thus to have a significant impact on the theoretical study of jet and rocket propellants.

Crawford, T. Daniel↗
NASA NTRS2007

Sensitivity Testing of the NSTAR Ion Thruster

During the Extended Life Test of the DS1 flight spare ion thruster, the engine was subjected to sensitvity testing in order to characterize the macroscopic dependence of discharge chamber sensitivity to a +\-3% vatiation in main flow, cathode flow and beam current, and to +\5% variation in beam and accelerator voltage, was determined for the minimum- (THO), half- (TH8) and full power (TH15) throttle levels. For each power level investigared, 16 high/low operating conditions were chosen to vary the flows, beam current, and grid voltages in in a matrix that mapped out the entire parameter space. The matrix of data generated was used to determine the partial derivative or senitivity of the dependent parameters--discharge voltage, discharge current, discharge loss, double-to-single-ion current ratio, and neutralizer-keeper voltage--to the variation in the independent parameters--main flow, cathode flow, beam current, and beam voltage. The sensititivities of each dependent parameter with respect to each independent parameter were determined using a least-square fit routine. Variation in these sensitivities with thruster runtime was recorded over the duration of the ELT, to detemine if discharge performance changed with thruster wear. Several key findings have been ascertained from the sensitivity testing. Discharge operation is most sensitve to changes in cathode flow and to a lesser degree main flow. The data also confirms that for the NSTAR configuration plasma production is limited by primary electron input due to the fixed neutral population. Key sensitivities along with their change with thruster wear (operating time) will be presented. In addition double ion content measurements with an ExB probe will also be presented to illustrate beam ion production and content sensitivity to the discharge chamber operating parameteres.

electric propulsion↗
NASA NTRS2005

Constraining Oceanic Sources and Sinks of CO2 by Inverse Modeling: A First Step Toward Assimilation of In Situ and Remote Data

The primary goal of this collaborative project, lead by Nicholas Gruber at UCLA, is to better quantify the air-sea fluxes of pre-industrial and anthropogenic carbon dioxide CO2 for the global ocean with a newly developed inverse method. The method involves using ocean general circulation models (OGCMs) to characterize the subsurface distribution of tracers injected for different regions of the ocean surface. By linearly combining a series of these basis (or Greens) functions, one should be able to recreate any arbitrary ocean tracer field. In our case, we are using the global database of ocean dissolved inorganic carbon (DIC) from the WOCE/JGOFS global CO2 survey from the 1990s. The observed DIC field has been partitioned into pre-industrial and anthropogenic components using the empirical C* technique. We then solve for an optimal combination of the regional basis functions that best match the two sets observations in a least-squares sense with singular value decomposition. This in turn provides solutions for the regional air-sea CO2 fluxes (pre-industrial and anthropogenic) that are independent of more traditional estimates derived from surface water disequilibria in the partial pressure of CO2 and wind speed gas exchange parameterizations. We are apply the method to a suite of nine different OGCMs to test the robustness of the derived surface fluxes fields to differences in ocean circulation. Together with a collegue Keith Lindsay at NCAR, S. Doney's role in this project was to create the appropriate steady-state and temporally evolving basis functions from the NCAR CSM-1 ocean circulation model. Because of a slow start across the project in defining a common set of experimental protocols, we requested and were given a one year, no cost extension. Both the steady-state and temporally evolving basis functions were delivered to the UCLA group for further analysis. Additional information is included in the original extended abstract.

Doney, Scott↗