1. Introduction
The Aircraft Meteorological Data Relay (AMDAR) program of the World Meteorological Organization (WMO) was initiated to employ existing meteorological equipment and data links [Aircraft Communications Addressing and Reporting System (ACARS)] on commercial jet aircraft to acquire upper-air measurements in real time (WMO AMDAR Panel 2007). AMDAR data contain at least time, position, altitude, temperature, and wind speed and direction. A few aircraft additionally carry prototype humidity sensors (Helms et al. 2009), but mixing ratio is rarely reported.
AMDAR data are extensively used as input to numerical weather prediction models (Moninger et al. 2003). Operational quality control is essentially ensured by monitoring the differences between AMDAR reports and the model background (see Drüe et al. 2008, for more details). The differences are typically so small (Cardinali et al. 2003) that forecast centers put AMDAR data into the highest of all confidence classes (Schraff and Hess 2002).
The robustness of the operational measurement systems used is without doubt, due to the strict regulations applied (e.g., see FAA 1968, 1995). Nevertheless, the allowed tolerances are defined by flight engineering needs, not by the requirements for meteorological applications.
In addition to instrument uncertainties, it is known that airborne measurements can have significant errors caused, for example, by variations in speed or load or uneven orientation during maneuvers (e.g., see Drüe and Heinemann 2001; Cooper and Baumgartner 1989; Painting 2003). Such systematic errors are mostly well correlated with flight dynamic parameters (see, e.g., Vörsmann 1990; Nacass 1992). By comparing AMDAR reports to other measurements or numerical forecast models, it was found that the resulting errors of AMDAR reports usually consist of a smaller random and a larger systematic contribution (Drüe et al. 2008, 2010; Ballish and Kumar 2008). The systematic contribution was shown to be associated with aircraft type and height (Drüe et al. 2008; Ballish and Kumar 2006) as well as with season (Ballish and Kumar 2008; Drüe et al. 2010).
During these previous investigations, indications of further systematic dependencies of AMDAR data errors on the flight dynamic state also have been found. From flight testing experience (Stoliker et al. 2005), the following flight dynamic parameters have been identified as candidates for being most influential (A. Hoff 2007, personal communication): heading, sideslip, pitch and roll angles, airspeed (or Mach number), and load (“weight over delta”; e.g., see Stoliker et al. 2005). AMDAR reports, however, do not contain such information, except for a flag indicating that the roll angle exceeds a threshold of 5° (Painting 2003). However, not all participating airlines set this flag accordingly.
The objective of the current investigation, hence, is to develop a method to retrieve the actual trajectory and as many flight dynamic parameters as possible from AMDAR data. This paper is organized in the following way: After this introduction in section 1, previous efforts and available input data are described in section 2. Section 3 summarizes the algorithm and section 4 presents the verification results. The summary and conclusions are given in section 5.
2. Data and models
a. AMDAR data
1) Data flow
On board a participating aircraft (a “unit”), the AMDAR software gathers meteorological values and compiles and transmits messages via ACARS to the airline operating the unit (WMO AMDAR Panel 2007; Painting 2003). The airline converts each single report into one of the two data formats defined by the WMO: either the FM-42 (human readable) or FM-94 Binary Universal Form for the Representation (BUFR) format (Painting 2003). The converted data are transferred to the national weather service that distributes them worldwide. In this paper, we use AMDAR data from the German weather service Deutscher Wetterdienst (DWD) observation database.
To ensure the privacy of the airlines, all information on flight origin, destination, and aircraft registration is removed and most flight numbers are replaced by anonymous unit identifiers (Painting 2003). Hence, neither aircraft registration nor type or airline can be determined from AMDAR data.
The data centers of some participating weather services, however, maintain lists indicating registration, airline, and type of each unit. Although aircraft type may be available for quality control and research, aircraft registration and airline are never disseminated. Nevertheless, information on both the type and operating airline of most units was available for the present study, resulting from ongoing research cooperation with DWD. Because the publication of this information is not permitted, no unit numbers will be specified in the following and specific airlines will referenced by randomly chosen single letters “A”–“K.”
2) AMDAR software
Although some airlines have their own software, there are two standards established (Stickland 2002): The first is ACARS Teledyne Aircraft Conditioning and Management System (ACMS) AMDAR software (AAA), which was developed from the original WMO specifications and is freely available. The second is ARINC-620 Meteorological Report, which is the industry standard issued by Aeronautical Radio Incorporated (ARINC).
Both standards minimize the number of digits transmitted via ACARS to save telecommunication costs and ensure timely transmission. One result is that both standards strip seconds from both time and the latitude–longitude position of the report (Painting 2003). After this, AMDAR data never contain seconds, although the format of distributed data has room for them (Painting 2003).
Another result is that if a high initial data rate is chosen after takeoff, then no actual position information is transmitted by ARINC-620 software, in order to squeeze more reports in one transmission block. The position distributed with such reports is generated on the ground by linear interpolation. Or, in case of at least one airline, the position of the first data point is simply repeated. We will refer to such ground-generated position values as “faked positions” hereafter. Faked positions can be marked by setting a flag in the FM-94 format, but apparently not every airline sets the flag accordingly. In the case of increased roll angle a flag in FM-94 format can be set, and in the case of increased roll or pitch angle “unsteady flight” can be reported as flight phase, instead of “ascent” or “descent” (Painting 2003) but again, apparently not every airline uses such flagging.
As a consequence of such limited time and position resolution and unreliable quality flags, it is almost impossible to know the actual trajectory, or any flight dynamic parameters, from AMDAR data alone.
3) Collection strategy
AMDAR data are collected following various schemes, at fixed time intervals, at fixed geographic positions, or at fixed pressure levels during ascent and descent (AMDAR Panel 2010, personal communication). The European Regional AMDAR Program (E-AMDAR) mostly collects data at a high rate during selected ascents or descents, yielding vertical atmospheric profiles (Pastre 2000), while the U.S. national program and others collect data continuously throughout all flight phases.
The need to develop an error model for level flight, however, seems less urgent because en route flight is usually straight, even, at cruise speed, and in clean configuration [with flaps fully retracted; see Stengel (2004)]. Conditions are nominal with respect to sensor calibrations for this phase of flight and the expected errors are minimal (Stickney et al. 1994). Our study hence focuses on ascent and descent.
In ascent and descent, the WMO specification allows for observations at fixed time intervals or at fixed pressure levels (Painting 2003). Observations should be denser near the ground to better resolve the atmospheric boundary layer (see Fig. 1). To determine the data collection schemes in use, the number of data in the first 12 min was counted for every profile taken over Europe in January 2008. The results were compared to the number of data per minute calculated for all possible combinations of t1…5 and p1,3. It was found that most aircraft of the same airline share an identical scheme and timer values. Most of the schemes detected to be typical for a fleet were confirmed by AMDAR Panel members (2007, 2009 personal communication). However, a fraction of most fleets do not use the typical scheme, probably because of incoherent software or leased aircraft.

Definition of flight phases, timer values t1–t5, and pressure intervals p1 and p3 used in the text.
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1

Definition of flight phases, timer values t1–t5, and pressure intervals p1 and p3 used in the text.
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1
Definition of flight phases, timer values t1–t5, and pressure intervals p1 and p3 used in the text.
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1
b. BADA model
Base of Aircraft Data (BADA) is an aircraft performance model (APM) developed by the EUROCONTROL Experimental Center (EEC) mainly for trajectory simulation and prediction by air traffic management (ATM) facilities. It contains performance characteristics and operating procedure data for a total of 294 supported (or equivalent) types EEC (2009b).


The BADA model for each aircraft is based on an operations performance model and airline procedure model (EEC 2009a). The operations performance model defines aircraft drag, thrust, and fuel consumption depending on altitude, mass, speed, and aerodynamic configuration (i.e., flap setting). The airline procedure model defines the typical aircraft operation characteristics, which include takeoff and landing speeds, standard cruising speed, and the acceleration or deceleration profile during ascent or descent, respectively, which includes the usual choice of the aerodynamic configuration. The empirical coefficients of each BADA model are determined by the coupled multivariate nonlinear fitting of model trajectories to measured aircraft positions (see EEC 2009a, for details).
The model allows a number of input parameters to be set, but when they were varied in reasonable ranges it turned out that for the purposes of this study it is sufficient to use the precalculated performance table data (PTD) provided within BADA.
In the BADA documentation, there is no clear distinction between altitude, that is, vertical distance to the sea level, and height, that is, vertical distance to the surface. Because the present study mainly focuses on (the lower part of) ascent and descent, the BADA vertical coordinate is associated to height in the following.
c. Trajectory models
Recorded data from AMDAR units cannot be obtained because the flight data recorder (FDR) stores only a limited time history and can only be read out after dismantling it from the aircraft (Lufthansa Technik 2009, personal communication).
In Drüe et al. (2008), for each of three consecutive data points a circle was fitted and the direction of the tangent at the center point was used as heading. However, results are not acceptable if the spacing of the actual positions where AMDAR data are taken is on the order of or below the resolution of position.
Air traffic control (ATC) radar tracks aircraft with a time resolution of several seconds, including the height transmitted by the aircraft transponder (Duke 2009). However, AMDAR reports cannot be associated with radar tracks because of the missing aircraft registration and flight number.
ATC uses trajectory prediction or a so-called trajectory solver for planning the flow of aircraft through control sectors and to schedule arrival times (Erzberger and Tobias 1986; Slattery and Zhao 1997). Although such models are optimized for prediction of the vertical trajectory (i.e., height versus time), they are potentially usable for purposes of the present study. Using such models is, however, prevented by the fact that they are factually not available to researchers outside the organizations that develop them.
3. Algorithm description
a. Data preparation
For each profile, the AMDAR data are retrieved from a database and sorted by increasing time and decreasing (in ascent) or increasing (in descent) static pressure, the aircraft type is read from the list of units and the corresponding performance coefficients are extracted from BADA.
To facilitate geometric calculations, AMDAR positions are projected onto an x–y plane by an oblique stereographic projection using the Conformal Map Function (CMAPF) version 2.0 routines (Taylor 1997). For each profile, the position of the lowest data point is chosen as the projection center.
During approach, aircraft often have to interrupt their descent at certain altitudes. If they undershoot such an altitude, then a slight climb might be required. If in such cases, multiple AMDAR reports are generated within the same minute, then sorting them by ascending static pressure does not yield the correct order of the data points. Hence, for each possible permutation, the algorithm calculates the distance needed to travel from the last position in the preceding minute j − 1 along all of the data points in minute j to the first position in the following minute j + 1, and the permutation with the shortest distance traveled is chosen.
b. Scheme recognition
Detection of a known collection scheme not only easily yields report seconds (a more general method will be described in section 3), but also provides information about positions being faked. Therefore, a method was developed that tests if one of the known schemes applies. It evaluates the first 6 min of an ascent or the last 6 min of a descent, respectively, and calculates the ARINC timer values, if needed.
The success of this procedure was tested by calculating the probability of detection (POD), which gives the percentage of cases in which the detected scheme matches the scheme that is typical for the airline operating that particular unit. Table 1 summarizes results for the AMDAR profiles taken in January 2008 over Europe. Except for airlines F, H, and K (which have reported only a small number of profiles), the POD is always well above 50%. The average POD is even 80% in ascent and almost 95% in descent.
Probability of detection of the data collection schemes for (top) ascent and (bottom) descent. Collection schemes: “?” unknown, “1pt” one-phased timed, “2pt” two-phased timed, “pl” pressure levels, or “A” two-phased timed starting at full minute. Data represent the E-AMDAR profiles taken over Europe in January 2008.


For example, in case of airline B about 7% of the fleet are aircraft with different avionics (Giraud 2010) that typically implement pressure-based reporting (AMDAR Panel 2009, personal communication). Hence, POD cannot reach values of more than roughly 93%. With respect to such general conditions, the calculated POD values seem rather satisfying.
c. Trajectory fitting
1) First guess
In case a particular reporting scheme is detected, all faked data point positions are flagged. In the case of an unknown two-phased timed or pressure-triggered ascent scheme, data taken at a higher rate than every 15 s are assumed to be faked, except for the first data point.
From the remaining data, a first guess of the trajectory is determined. For every data point i, heading Φpp(i) is estimated from the forward difference of the data position, and the values are smoothed using a running window of a 3-min flight time. From the estimated heading, the amount |ΔΦpp(i)| of the turn angle at each data point is calculated. In case the sense of turn is different from the last data point or the amount of the turn angle exhibits a local minimum, the value is set to zero to distinguish closely following turns. Each group of one or more consecutive data points where |ΔΦpp(i)| exceeds a threshold value of 1.5° is regarded as a candidate for a turn. If more than five candidates are detected, then the candidates with the highest peak values are selected. If less than five candidates are detected, additional candidates are added at the center of the longest distance between two candidates until five candidates are reached.
Based on the first guess, the aircraft trajectory is fitted. The fit is done first horizontally and then vertically, since Slattery and Zhao (1997) have demonstrated that this procedure represents a calculation cost–effective way to achieve sufficiently precise results. As in Slattery and Zhao (1997), the wind is also neglected here so that airspeed equals groundspeed for trajectory fitting.
2) Horizontal
For every profile, up to six raw trajectories are composed from two templates: straight flight and a turn with a 30° bank angle (EEC 2009a). If a profile contains less than 13 (nonfaked) data points, the number of segments is limited to N/2 − 1, in which N is the number of data points. For each of these raw trajectories, a minimal set of parameters is optimized to minimize a cost function. The trajectory that yields the lowest cost function, that is, the best fit, is chosen as the horizontal solution.





Distance to a 1-min part (bold gray line) of an exemplary trajectory (dashed gray line) as function of the position. Contours give distance in kilometers.
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1

Distance to a 1-min part (bold gray line) of an exemplary trajectory (dashed gray line) as function of the position. Contours give distance in kilometers.
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1
Distance to a 1-min part (bold gray line) of an exemplary trajectory (dashed gray line) as function of the position. Contours give distance in kilometers.
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1
Respectively, di and de are the distance between the trajectory start position and the first data point (subscript i) and the trajectory end position and the last data point (subscript e).






Starting from the first guess, the trajectory parameters are optimized by Powell’s quadratic convergent method (Press et al. 1986) until J reaches a minimum, allowing a remaining fractional tolerance of 10−5.
The weight factors in Eqs. (4)–(6) were determined empirically. For each type of data collection scheme, between 2 and 12 sample profiles were selected from AMDAR data take over Europe in January and July 2008. The samples ranged in complexity from a straight departure to an approach that was interrupted in a holding pattern. The weights were varied until a trajectory was found for all of the samples that seemed reasonable by visual inspection (see Figs. 3 and 4).

Example of a fitted 2-turn trajectory. Data are taken form the descent of an Airbus A319 to Frankfurt airport (EDDF) on 5 Jul 2006. Data point positions reported by AMDAR (black squares), the fitted trajectory (dashed line), and the fitted data point positions (open squares). The three bold black lines mark the runways of Frankfurt airport.
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1

Example of a fitted 2-turn trajectory. Data are taken form the descent of an Airbus A319 to Frankfurt airport (EDDF) on 5 Jul 2006. Data point positions reported by AMDAR (black squares), the fitted trajectory (dashed line), and the fitted data point positions (open squares). The three bold black lines mark the runways of Frankfurt airport.
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1
Example of a fitted 2-turn trajectory. Data are taken form the descent of an Airbus A319 to Frankfurt airport (EDDF) on 5 Jul 2006. Data point positions reported by AMDAR (black squares), the fitted trajectory (dashed line), and the fitted data point positions (open squares). The three bold black lines mark the runways of Frankfurt airport.
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1

Example of an ascent fitted to data from Lufthansa A340–600 “Leverkusen” departing from Chennai airport (VOMM) on 13 Aug 2008. Data are taken from the verification dataset described in section 4. Presentation as in Fig. 3, but some data points represent faked positions (gray squares).
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1

Example of an ascent fitted to data from Lufthansa A340–600 “Leverkusen” departing from Chennai airport (VOMM) on 13 Aug 2008. Data are taken from the verification dataset described in section 4. Presentation as in Fig. 3, but some data points represent faked positions (gray squares).
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1
Example of an ascent fitted to data from Lufthansa A340–600 “Leverkusen” departing from Chennai airport (VOMM) on 13 Aug 2008. Data are taken from the verification dataset described in section 4. Presentation as in Fig. 3, but some data points represent faked positions (gray squares).
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1
3) Vertical
Because altitude is measured and reported with good accuracy, a vertical trajectory fit is only needed if the reporting scheme is pressure based or unknown. If no known scheme reproduces the observed number of data per minute, then a vertical trajectory fit is done.
(i) Known collection scheme
Applying a known scheme is strikingly difficult, because single data points may be missing, as explained in the previous section. Hence, the algorithm applies a three-pass method: First, the sequence of reporting times is modeled by assuming that the first data point is taken at the full minute. Second, the second of the first data point is varied until the modeled number of data in every minute is larger than or equal to that measured and the sum of the differences is minimal. In the third pass, the seconds of data for minutes with less data points than that modeled are estimated by interpolating time according to altitude and choosing the nearest modeled time.
(ii) Unknown collection scheme


After the seconds are estimated, the data point positions are recalculated by integrating traveled distance along the horizontal trajectory (see Figs. 5 and 6).

Vertical profile of the descent shown in Fig. 3. The altitudes of the AMDAR reports are marked (black bars); each bar fills the one minute of the report. The fitted trajectory (dashed line) and the fitted data point positions (open squares).
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1

Vertical profile of the descent shown in Fig. 3. The altitudes of the AMDAR reports are marked (black bars); each bar fills the one minute of the report. The fitted trajectory (dashed line) and the fitted data point positions (open squares).
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1
Vertical profile of the descent shown in Fig. 3. The altitudes of the AMDAR reports are marked (black bars); each bar fills the one minute of the report. The fitted trajectory (dashed line) and the fitted data point positions (open squares).
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1

Vertical profile of the ascent shown in Fig. 4. Same presentation as in Fig. 5.
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1

Vertical profile of the ascent shown in Fig. 4. Same presentation as in Fig. 5.
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1
Vertical profile of the ascent shown in Fig. 4. Same presentation as in Fig. 5.
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1
d. Flight dynamic state
1) Speed






In the previous calculations, the wind was neglected, but to include the wind, it is assumed at this point that the speed values calculated from the fitted trajectory are ground speed. The track angle Ψ(i), that is, the direction of the motion relative to the ground, is calculated from the geometry of the fitted trajectory. Subtracting the wind vector contained in the AMDAR data then yields the airspeed vector (see Drüe et al. 2008), that is, airspeed υas(i) and heading Φ(i) of the aircraft at data point i are calculated.
The vertical velocity of the aircraft is calculated from the differences of altitude and time to the nearest data point that is more than 5 s away. This minimum distance was introduced to reduce numerical noise in the case of very closely spaced data.
2) Aircraft mass








Therefore, Eq. (13) is evaluated for all data points of a profile and the lower 85% percentile value of the calculated masses is taken as estimation for mac.
3) Orientation
The aircraft roll angle is set to zero on straight segments and equal to the predefined roll angle in turns.




















Typical civil jet aircraft use a moderate sweep (Jenkinson et al. 1999) angle between 25°, for example, on Airbus A320, and 37°, on Boeing 747 (Morichon 2006). Because Λ is not included in the BADA model, the algorithm uses a fixed value of 30°. This simplification causes an uncertainty of about 7% in CLα, which corresponds to roughly only ±0.1° in angle of attack.


4. Verification
To verify the results of the algorithm, datasets are needed that were taken by typical commercial jet aircraft and contain attitude angles as well as atmospheric variables.
The Global Aircraft Dataset (GADS; Tenenbaum 1991) collected by the State University of New York does contain such information, but unfortunately the project does not collect data during ascent and descent (J. Tenenbaum 2010, personal communication).
a. CARIBIC
The most useful dataset for this study is taken by the Civil Aircraft for the Regular Investigation of the Atmosphere Based on an Instrument Container (CARIBIC) project, which is carried out to study and monitor important chemical and physical processes in the earth’s atmosphere using long-distance passenger flights (Brenninkmeijer et al. 2007).
CARIBIC regularly deploys an airfreight container with automated scientific instruments to the Lufthansa A340–600, called “Leverkusen.” This aircraft is equipped with an air and particle (aerosol) inlet underneath the fuselage that is permanently installed. Along with the data of its own instrumentation, CARIBIC records navigation data from the aircraft data bus with 1-s resolution. We used the navigation data from all 22 CARIBIC flights in 2008 for validation of the algorithm. Unfortunately, AMDAR data from this aircraft were not available for all flights.
Therefore a simulator program was developed that mimics onboard software conforming to the ARINC-620 meteorological reporting standard, including data smoothing (Painting 2003). The algorithm described in section 3 was then applied to the simulated AMDAR data and the results were compared to the measured values for each data point. It should be noted that this procedure also ensures that the verification is performed with data that were not used to develop the algorithm.
Although some assumptions made in section 3 seem quite crude, the agreement is mostly very good. The average horizontal distance of the measured position to the reconstructed position is just 460 m and the mean difference between the measured and reconstructed altitude is 7.6 m. True airspeed is reconstructed with a mean difference of only 0.4 m s−1. Linear regression yields a slope of 0.944 and an intercept of 1.3 m s−1. The scatter, however, is greater than that of position and the regression coefficient is only 0.68 (see Table 2). The vertical speed, however, is rather well reconstructed with a slope of 0.964, an intercept of 0.3 m s−1, and a regression coefficient of 0.94.
Regression results of estimated vs measured values for selected variables: icept denotes the regression intercept, r2 denotes the correlation coefficient, and THdg and AOA stand for true heading and AOA, respectively.


The estimated aircraft mass (Fig. 7) is 276 t on average, with 330 t after takeoff and 222 t before landing. Because the true mass values are not known, a verification of this estimate is not possible. Nevertheless the values obtained seem quite reasonable: for 32 out of 44 profiles, the estimated mass is inside the permitted weight limits, departing aircraft are heavier than arriving ones, and the average fuel burn is 90 t, which seems realistic (Jenkinson et al. 1999).

Aircraft mass estimated from the validation dataset. Bars indicate the total number of profiles per 25-t bin. The gray portion of each bar gives the number of ascents. The bandwidth of permitted takeoff mass for A340–600 is marked by a white plot background.
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1

Aircraft mass estimated from the validation dataset. Bars indicate the total number of profiles per 25-t bin. The gray portion of each bar gives the number of ascents. The bandwidth of permitted takeoff mass for A340–600 is marked by a white plot background.
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1
Aircraft mass estimated from the validation dataset. Bars indicate the total number of profiles per 25-t bin. The gray portion of each bar gives the number of ascents. The bandwidth of permitted takeoff mass for A340–600 is marked by a white plot background.
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1
The good agreement of reconstructed and measured values promises reasonable estimation of the aircraft orientation. And indeed, true heading and pitch angle are estimated with a slope near unity, an almost zero intercept, and regression coefficients above 0.65 (see Fig. 8). The scatter of the angle of attack is rather large but still the slope is at 0.65 and the intercept is near zero. Roll angle, however, is not met very well, probably because it is prescribed as a fixed value, which is apparently less realistic than that estimated in literature (e.g., EEC 2009a; Jenkinson et al. 1999).

Scatterplot of estimated (AMDAR) vs measured (CARIBIC) angles. (left) True heading and (right) pitch angle. The line of identity (dashed line) and the linear fit to the data (solid line).
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1

Scatterplot of estimated (AMDAR) vs measured (CARIBIC) angles. (left) True heading and (right) pitch angle. The line of identity (dashed line) and the linear fit to the data (solid line).
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1
Scatterplot of estimated (AMDAR) vs measured (CARIBIC) angles. (left) True heading and (right) pitch angle. The line of identity (dashed line) and the linear fit to the data (solid line).
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1
b. ATRA
A second, much smaller dataset was supplied by Deutsches Zentrum für Luft- und Raumfahrt (DLR) from its Advanced Technology Research Aircraft (ATRA), which is an A320–200 modified for aviation research (Renouard-Vallet et al. 2010). Because this aircraft is rather new and mostly operates using nonstandard procedures (H. Becker, DLR, 2010, personal communication), data from only four flights could be used.
As from the CARIBIC dataset, AMDAR data were generated by the ARINC-620 simulator and used as input for the algorithm described in section 3. Because the data supplied by DLR do not contain meteorological values, wind and temperature data of the nearest operational radiosonde were used.
Compared to the CARIBIC dataset, the results are even better with respect to correlation (Table 2). The slope is even closer to unity, except for pitch. However, no measured values of the angle of attack are available and none of the generated AMDAR data points is in a turn. Hence, the ATRA dataset cannot further support the verification of roll and attack angle.
5. Summary and conclusions
To enable development of an error model for AMDAR weather reports from commercial aircraft, information about the flight dynamic state of the aircraft is needed. Because such information is neither included in the data reported nor can be substituted from other sources, a method was developed to estimate flight dynamic state variables only from AMDAR data and the aircraft type.
AMDAR data do not contain seconds, neither in position nor in time. Hence, estimating time and position of each report with greater precision is a primary goal of the algorithm. The only additional information needed is the aircraft type, which is not disclosed to the public but is available to weather data centers.
The algorithm first tries to determine the data collection scheme in order to determine (presumably) faked positions and (time) seconds of the data points. The detection was shown to exhibit a POD of 80% in ascent and 94% in descent.
The most probable horizontal trajectory is determined by consecutively fitting trajectory templates to the nonfaked AMDAR data positions. In case no known data collection scheme was detected, the seconds (time) are determined by fitting a scaled BADA vertical speed profile to the data point altitudes. From the fitted positions and times, speed, mass, and orientation are estimated applying the BADA total energy model, assuming thrust minus drag (T − D) is always at its BADA reference value.
To validate the algorithm, two rare datasets of flight dynamic recordings from commercially used jet aircraft types were used. On these data, it could be shown that the new algorithm estimates the actual position with good accuracy, that is, a position with less than 0.5 km and altitude with less than 8-m average deviation. The estimated mass, although no reference data are available, appears rather reasonable and mostly keeps in the permitted mass limits. Pitch, heading, and angle of attack are estimated with rather good agreement, but angle of attack exhibits a huge scatter. The output roll angle is not a good estimation because it is prescribed as a fixed value, which is a general assumption made by BADA and others, which might be too simple.
The computation time for one profile on a standard desktop computer is about 1.5 s. This value, however, overestimates the computational cost of the presented algorithm. First, the code is experimental and offers wide potential for optimization. Second, once an error model is developed and the most relevant parameters are known, either the model could be simplified to yield only these or measurements of these parameters could be included in a future version of AMDAR reports.
Although the envisioned error model is not yet developed, the implications of the above results on the expected outcome may be estimated: From the Investigation of the Greenland Boundary Layer over Summit (IGLOS) field experiment (Drüe and Heinemann 2003), the pitch dependency of the true air temperature determined by a Dornier 228 standard temperature probe was estimated to amount roughly to 0.03 K (°)−1. Hypothesizing that the same value applies for an Airbus A319 allows estimation of the error by neglecting a nonzero pitch (Fig. 9). The value is significantly nonzero at all levels (except one) and amounts about half of the typical systematic deviation of AMDAR temperature (Drüe et al. 2008). Such a correction would be capable of reducing systematic AMDAR errors while increasing the random error only rather slightly.

Hypothetical temperature error by neglecting pitch. Data are taken from the descent shown in Fig. 3.
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1

Hypothetical temperature error by neglecting pitch. Data are taken from the descent shown in Fig. 3.
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1
Hypothetical temperature error by neglecting pitch. Data are taken from the descent shown in Fig. 3.
Citation: Journal of Atmospheric and Oceanic Technology 28, 7; 10.1175/2011JTECHA1499.1
Although the scatter of some variables may by large, the flight dynamic state is estimated with good accuracy, on average. I hence conclude that the algorithm presented meets the objectives of the present study and may serve as basis for the envisioned error model or similar applications.
Acknowledgments
AMDAR data were supplied by Deutscher Wetterdienst within a continued cooperation on AMDAR data quality and application. The E-AMDAR technical coordinator, the WMO AMDAR technical coordinator, the Global Systems Division (GSD) of the Earth System Research Laboratory (ESRL) at the U.S. National Oceanic and Atmospheric Administration (NOAA) and the National Centers for Environmental Prediction (NCEP) of the U.S. National Weather Service (NWS) are thankfully acknowledged for supplying AMDAR unit metadata. Airbus A340 data were kindly supplied by the CARIBIC project through the Max-Planck-Institute for Chemistry, Atmospheric Chemistry Department, Mainz, Germany. Airbus A320 (ATRA) data were kindly supplied by DLR, Braunschweig, Germany.
REFERENCES
Ballish, B., and Kumar K. , 2006: Comparison of aircraft and radiosonde temperature biases at NCEP. Preprints, 10th Symp. on Integrated Observing and Assimilation Systems for the Atmosphere, Oceans, and Land Surface (IOAS-AOLS), Boston, MA, Amer. Meteor. Soc., 3.5. [Available online at http://ams.confex.com/ams/pdfpapers/103076.pdf.]
Ballish, B., and Kumar V. K. , 2008: Systematic differences in aircraft and radiosonde temperatures. Bull. Amer. Meteor. Soc., 89, 1689–1707.
Brenninkmeijer, C. A. M., and Coauthors, 2007: Civil aircraft for the regular investigation of the atmosphere based on an instrumented container: The new CARIBIC system. Atmos. Chem. Phys., 7, 4953–4976.
Cardinali, C., Isaksen L. , and Andersson E. , 2003: Use and impact of automated aircraft data in a global 4DVAR data assimilation system. Mon. Wea. Rev., 131, 1865–1877.
Cooper, W. A., and Baumgartner D. , 1989: Meeting review: Workshop on airborne instrumentation. NCAR Tech. Note NCAR/TN-330+PROC, 1–71.
Drüe, C., and Heinemann G. , 2001: Airborne investigation of arctic boundary-layer fronts over the marginal ice zone of the Davis Strait. Bound.-Layer Meteor., 101, 261–292.
Drüe, C., and Heinemann G. , 2003: Investigation of the Greenland atmospheric boundary layer over summit 2002 (IGLOS): Field phase report. Alfred Wegener Institute for Polar and Marine Research, Reports on Polar and Marine Research 447, 88 pp. [Available online at http://epic.awi.de/Publications/BerPolarforsch2003447.pdf.]
Drüe, C., Frey W. , Hoff A. , and Hauf T. , 2008: Aircraft type-specific errors in AMDAR weather reports from commercial aircraft. Quart. J. Roy. Meteor. Soc., 134, 229–239.
Drüe, C., Hauf T. , and Hoff A. , 2010: Comparison of boundary-layer profiles and layer detection by AMDAR and WRT/RASS at Frankfurt airport. Bound.-Layer Meteor., 135, 407–432.
Duke, G., 2009: Air Traffic Control. Midland, 112 pp.
EEC, 2009a: Base of Aircraft Data (BADA) aircraft performance modelling report. Eurocontrol Experimental Centre Technical/Scientific Rep. 2009-009, 68 pp. [Available online at http://www.eurocontrol.int/eec/gallery/content/public/document/eec/report/2009/009_BADA_Aircraft_performance_modelling.pdf.]
EEC, 2009b: User manual for the Base of Aircraft Data (BADA) revision 3.7. Eurocontrol Experimental Centre Technical/Scientific Rep. 2009-003, 102 pp. [Available online at http://www.eurocontrol.int/eec/gallery/content/public/document/eec/report/2009/003_BADA_3_7_User_manual.pdf.]
Erzberger, H., and Tobias L. , 1986: A time-based concept for terminal-area traffic management. National Aeronuatics and Space Administration Ames Research Center Tech. Memo. 88243, 19 pp.
FAA, 1968: Maximum allowable airspeed indicator systems. Aircraft Certification Service, Federal Aviation Administration, Department of Transportation, Technical Standard Order Vol. TSO-C46a, 4 pp.
FAA, 1995: Temperature instruments. Aircraft Certification Service, Federal Aviation Administration, Department of Transportation, Technical Standard Order Vol. TSO-C43c, 4 pp.
Giraud, N., 2010: Airfleets.net - The Book. Airfleets.net, 168 pp.
Helms, D., Johnston K. L. , Sanger G. , Taubvurtzel B. , Petersen R. A. , Homans A. , and Hoff A. , 2009: Testing and deployment of the water vapor sensing system II. Preprints, 25th Conf. on International Interactive Information and Processing Systems (IIPS) for Meteorology, Oceanography, and Hydrology, Phoenix, AZ, Amer. Meteor. Soc., 4A.5. [Available online at http://ams.confex.com/ams/89annual/techprogram/paper_149759.htm.]
Jenkinson, L. R., Simpkin P. , and Rhodes D. , 1999: Civil Jet Aircraft Design. Butterworth-Heinemann, 436 pp.
Moninger, W. R., Mamrosh R. D. , and Pauley P. M. , 2003: Automated meteorological reports from commercial aircraft. Bull. Amer. Meteor. Soc., 84, 203–216.
Morichon, L., 2006: Selected statistics in aircraft design. Diploma thesis, Department of Vehicle Technology and Aircraft Construction (Fahrzeugtechnik und Flugzeugbau), Hamburg University of Applied Sciences, 62 pp.
Nacass, P. L., 1992: Theoretical errors on airborne measurements of static pressure, impact temperature, air flow angle, and air flow speed. NCAR Tech. Note 385+STR, 62 pp.
Painting, J. D., 2003: AMDAR reference manual. World Meteorological Organization Tech. Rep. WMO-No. 958, 84 pp.
Pastre, C., 2000: What is new in our programmes? EUMETNET News, Vol. 6, EUMETNET, Paris, France, 13–15.
Press, W. H., Teukolsky S. A. , Vetterling W. T. , and Flannery B. P. , 1986: Numerical Recipes in Fortran. Cambridge University Press, 848 pp.
Renouard-Vallet, G., Saballus M. , Schmithals G. , Schirmer J. , Kallo J. , and Friedrich K. A. , 2010: Improving the environmental impact of civil aircraft by fuel cell technology: Concepts and technological progress. Energy Environ. Sci., 3, 1458–1468.
Schraff, C. and Hess R. , 2002: Datenassimilation für das LM (Data assimilation for the LM). Promet, 27 (3/4), 156–164.
Slattery, R., and Zhao Y. , 1997: Trajectory synthesis for air traffic automation. J. Guidance Control Dyn., 20, 232–238.
Stengel, R. F., 2004: Flight Dynamics. Princeton University Press, 845 pp.
Stickland, J., 2002: AMDAR status at December 2001. Expert Team on Observational Data Requirements and Redesign of the Global Observing System Fourth Session, World Meteorological Organization Report to Commission for Basic Systems Open Programme Area Group on Integrated Observing Systems, 9 pp.
Stickney, T. M., Shedlov M. W. , and Thompson D. I. , 1994: Goodrich total temperature sensors. Goodrich Corporation Tech. Rep. 5755, 32 pp.
Stoliker, F., and Coauthors, 2005: Introduction to flight test engineering. RTO AGARDograph 300 Flight Test Techniques Series RTO-AG-300 AC/323(SCI-FT3)TP/74, Vol. 14, Research and Technology Organisation of the North Atlantic Treaty Organization, 456 pp.
Taylor, A. D., 1997: Conformal map transformations for meteorological modelers. Comput. Geosci., 23, 63–75.
Tenenbaum, J., 1991: Jet stream winds: Comparisons of analyses with independent aircraft data over southwest Asia. Wea. Forecasting, 6, 320–336.
Vörsmann, P., 1990: Meteopod, an airborne system for measurements of mean wind, turbulence, and other meteorological parameters. Onde Electr., 70, 31–38.
WMO AMDAR Panel, 2007: The international AMDAR program. World Meteorological Organization Information Flyer, 12 pp. [Available online at http://www.wmo.int/amdar/Publications/Final%20Production%20AMDAR%20Flyer.pdf.]