Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Ground observations of a space laser for the assessment of its in-orbit performance

Open Access Open Access
  • check crossmark updates
    Check for updates

Abstract

The wind mission Aeolus of the European Space Agency was a groundbreaking achievement for Earth observation. Between 2018 and 2023, the space-borne lidar instrument ALADIN onboard the Aeolus satellite measured atmospheric wind profiles with global coverage, which contributed to improving the accuracy of numerical weather prediction. The precision of the wind observations, however, declined over the course of the mission due to a progressive loss of the atmospheric backscatter signal. The analysis of the root cause was supported by the Pierre Auger Observatory in Argentina whose fluorescence detector registered the ultraviolet laser pulses emitted from the instrument in space, thereby offering an estimation of the laser energy at the exit of the instrument for several days in 2019, 2020, and 2021. The reconstruction of the laser beam not only allowed for an independent assessment of the Aeolus performance, but also helped to improve the accuracy in the determination of the laser beam’s ground track on single pulse level. The results presented in this paper set a precedent for the monitoring of space lasers by ground-based telescopes and open new possibilities for the calibration of cosmic-ray observatories.

Published by Optica Publishing Group under the terms of the Creative Commons Attribution 4.0 License. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI.

1. INTRODUCTION

The Aeolus mission by the European Space Agency (ESA) was launched on 22 August 2018 and successfully completed on 30 April 2023, exceeding its planned three-year life time by 18 months, before its satellite reentry into the Earth atmosphere on 28 July 2023. Its single payload, the Atmospheric Laser Doppler Instrument (ALADIN), was the first European lidar and the first Doppler wind lidar in space [1,2]. As an Earth Explorer mission, the satellite was designed as a technology demonstrator for future operational wind lidar space missions. However, Aeolus wind measurements were already operationally used for several years by numerous weather services, including the German Weather Service [3], Météo France [4], and the European Centre for Medium-Range Weather Forecasts [5].

The success of the Aeolus mission was achieved despite the slow but steady degradation of the atmospheric return signal of ALADIN, which reduced the precision of the wind observations. The root cause of the signal decline could not be identified unambiguously in the first year after the loss started in summer 2019. In particular, the question whether the loss occurred on the emit path or the receive path of the instrument or both could not be clarified beyond doubt. During this phase, a collaboration between the Aeolus team, including ESA and the Aeolus Data Innovation and Science Cluster (DISC), and scientists from the Pierre Auger Observatory was established when the laser beam from ALADIN was serendipitously registered by the telescopes of the Observatory in Argentina. The cooperation resulted in an unprecedented ground-based energy measurement of a space-borne laser after 25 years of space lidar missions [610]. It was the first time that a laser beam from a space-borne lidar was directly observed by the ground-based Auger Observatory, as its telescopes are optimized to detect fluorescence light from air showers and, therefore, only sensitive to ultra-violet (UV) wavelengths. The previous space missions operated with green and near-infrared laser wavelengths except for the two-week period of the Lidar In-space Technology Experiment (LITE) in 1994 with an additional UV laser beam [6]. The independent assessment of emitted laser energy by the Auger Observatory verified the hypothesis that the ALADIN signal is already degraded before the laser beam is transmitted to the atmosphere, and therefore strongly supported the root cause analysis. In addition, the reconstruction of the laser beam revealed a geometrical offset between its assumed and actual ground track, thereby pointing out an error in the retrieval of the geolocation within the Aeolus ground processor.

The benefit of the Auger observations for the monitoring of the Aeolus performance and the verification of the measurement ground track will be discussed in this paper. The ALADIN instrument and the Pierre Auger Observatory are described in Section 2.A and Section 2.B, respectively. The methods for the detection and reconstruction of the laser beam as well as for the determination of the laser energy are elaborated in Supplement 1. Section 3 covers the results of the measured ground track (Section 3.A) and energy of the laser beam (Section 3.B), while in Section 3.C the ground-based energy estimates are compared with the signal evolution that was measured with the ALADIN detectors in-orbit. The paper concludes with an outlook to future space laser missions in Section 4.

2. METHODS

A. Atmospheric Laser Doppler Instrument

The Doppler wind lidar instrument ALADIN was the single payload of the Aeolus satellite that circled the Earth on a sun-synchronous orbit at about 320 km altitude with a repeat cycle of seven days. It was composed of a pulsed frequency-stabilized UV laser transmitter, transmit-receive optics (TRO), a Cassegrain-type telescope in monostatic configuration (where the signal emission and reception are realized via the same telescope), and a dual-channel receiver that was sensitive for both molecular and particle backscatter from clouds and aerosols [1]. The wind measurement principle of ALADIN relied on detecting frequency differences between the emitted and the backscattered laser pulses. Due to the Doppler effect, the frequency of the outgoing pulse was shifted upon backscattering from particles and molecules that moved with the ambient wind. In order to observe a significant fraction of the Doppler frequency shift from the horizontal wind speed along its line-of-sight (LOS), the telescope was pointed with an angle of 35° off nadir (37.5° to 37.7° at the geolocation of the intersection with the surface due to the Earth’s curvature).

The design of ALADIN is illustrated in the simplified schematic in Fig. 1. The instrument was equipped with two fully redundant laser transmitters, referred to as flight models (FMs) A and B, which could be switched by means of a flip-flop mechanism (FFM). The transmitters were realized as diode-pumped injection-seeded Nd:YAG lasers in a master-oscillator power-amplifier configuration with subsequent frequency-tripling to generate narrowband and frequency-stable emission at a UV wavelength of 354.8 nm. Detailed information on the design and performance of the ALADIN laser transmitters and the challenges during the development stage can be found in Refs. [1117].

 figure: Fig. 1.

Fig. 1. Simplified optical layout of the ALADIN instrument onboard Aeolus. HR: highly-reflective mirror; FFM: flip-flop mechanism; BS: beam splitter; QWP: quarter-wave plate; FS: field stop.

Download Full Size | PDF

Before exiting the laser transmitter and entering the TRO, a small portion of the UV light was detected on a photodiode for monitoring purposes. When being guided through the TRO, the laser beam passed a quarter-wave plate that ensures circular polarization and a beam expander (magnification 3.4) that increased its diameter from 6.2 mm to about 21 mm. The nanosecond laser pulses (repetition rate 50.5 Hz) were transmitted to the atmosphere by the telescope whose primary mirror has a diameter of 1.5 m. The telescope with its magnification of 41.7 expanded the laser beam to a diameter of about 0.9 m and a divergence of about 18–20 µrad (86% encircled energy) in the atmosphere. Thus, the footprint of the laser beam was around 8 m on the Earth’s surface from a satellite altitude of 320 km and an off-nadir angle of 37.7°, resulting in a range of 414 km. The backscattered signal that was collected by the same telescope was directed to the optical receiver after passing through a field stop (FS) with a diameter of only 88 µm to limit the receiver field of view (FoV) to about 18 µrad, thereby diminishing the influence of the solar background radiation and accounting for the high angular sensitivity of the spectrometers.

It should be noted that the instrument and laser divergence were designed such that the laser divergence was smaller than the receiver FoV, but the as-built performance of the laser and the beam-expanding optics resulted in a slightly larger laser divergence of 20 µrad. The resulting loss at the field stop was accounted for in the receive path transmission budget.

The ALADIN receiver consisted of two complementary channels to derive the Doppler frequency shift from both the narrowband (FWHM $\approx 50\,\,\rm MHz$) Mie backscatter from particles like clouds and aerosols and from the broadband (FWHM $\approx 3.8\,\,\rm GHz$ at 355 nm and 293 K) Rayleigh-Brillouin molecular backscatter [1,2,1719]. Accumulation charge-coupled-device detectors were used in the two channels to finally register the Mie and Rayleigh backscatter signals.

B. Pierre Auger Observatory

The Pierre Auger Observatory is located in the Argentinian Pampa, at around ${35.2^ \circ}{\rm S}$, ${69.3^ \circ}\rm W$. It was designed to detect ultrahigh-energy cosmic rays and to accurately measure their properties, such as energy, arrival direction, and mass composition. While the origin of cosmic rays of lower energies can be attributed to the sun (mainly protons, $E\lesssim{10^9}\,\,\rm eV$) or sources within our Galaxy like supernova remnants (nuclei up to iron with atomic number $Z$, $E\lesssim Z \times {10^{15}}\,\,\rm eV$), cosmic rays with energies above ${\sim}{10^{18}}\,\,\rm eV$ originate from outside of our Galaxy; see, e.g., Ref. [20] for a recent review on ultrahigh-energy cosmic rays. The mechanisms that are able to accelerate particles to such energies are, however, still under debate [21]. These cosmic rays with highest energies have an extremely low flux, which makes a direct detection in space unfeasible. Instead, the atmosphere is used as a detection volume. As the cosmic rays enter the atmosphere, they interact with the nuclei of the air and create extensive air showers of secondary particles, carrying the information of the primary particle.

The Pierre Auger Observatory started taking data in 2004. It uses a hybrid detection method for the measurement of air showers, which is realized by the combination of the fluorescence detector (FD) and the surface detector (SD). The secondary shower particles that reach the ground are measured using the SD. It consists of 1660 tanks, each filled with 12,000 liters of ultra-pure water, covering an area of $3000\,\,\rm km^{2}$ [22]. Each station contains three photomultiplier tubes (PMTs), which are able to detect Cherenkov light emitted by the traversing shower particles in the tank. As the shower particles pass through the atmosphere, they excite nitrogen molecules in air, which, in turn, emit fluorescence light with wavelengths between about 280 and 430 nm; see, e.g., Ref. [23]. This light is measured by the FD [24], consisting of four sites. Each site contains six telescopes with a FoV of approximately ${30}^\circ \times {30}^\circ$; one site thus covers 180° horizontally. The sites are at the periphery of the SD array, overlooking the full area. The telescopes are equipped with a wavelength filter for background reduction and the full optical setup enables the detection of light for wavelengths between 310 and 410 nm. The $13\;{{\rm m}^2}$ spherical mirror of a telescope projects the shower image onto a camera, consisting of 440 PMTs as pixels, arranged in a honeycomb pattern. The arrangement of those components is illustrated in Fig. 2.

 figure: Fig. 2.

Fig. 2. Structure of one FD telescope of the Pierre Auger Observatory with some important features highlighted. Taken from Ref. [22].

Download Full Size | PDF

The FD operates exclusively during dark nights, particularly excluding periods around the full moon. Additionally, it requires clear atmospheric conditions to avoid scattering and attenuation of the light from the air shower by clouds. This results in a duty cycle of about 15% throughout a year [24].

Since the FD is capable of detecting pulses of faint UV light, a measurement of the Aeolus laser using the telescopes of the Pierre Auger Observatory was possible, as illustrated in Fig. 3. Due to the nature of the sun-synchronous dusk-dawn orbit of Aeolus, it always passed close to the sunrise or sunset over any given point on the surface of the Earth. This limited the opportunities for measurements with the FD. Additional to the aforementioned limitations by the moon cycle, a measurement could only take place if the satellite passage time fell within the astronomical night, which occurred only during the southern-hemisphere winter months, i.e., between May and August. Furthermore, the orbit had a seven-day repeat cycle. Therefore, a measurement could only be taken once per week. Overall, this resulted in up to six observations of the Aeolus laser per year, when clouds were not preventing the visibility.

 figure: Fig. 3.

Fig. 3. Geometry of the Aeolus laser beam being detected by one of the four fluorescence detectors of the Pierre Auger Observatory.

Download Full Size | PDF

The Aeolus laser was visible in the FD due to the scattering of the laser beams with air molecules as well as aerosol and cloud particles in the atmosphere. Similar to the fluorescence light of showers, this light could be measured by the telescopes, creating an image of the laser beam in the camera constituting a line of pixels. This served as the basis for the geometrical reconstruction of the laser beam. If the beam was seen by only one telescope, the geometry could be obtained only by a so-called monocular reconstruction. In a next step, the energy of the laser beam at the output of the ALADIN telescope was estimated by considering the molecular and particle scattering of the light when passing through the atmosphere under consideration of optical losses along the path. A detailed description of the reconstruction of the ALADIN laser beam and the estimation of the laser energy is provided in Supplement 1.

The first Auger measurement of the Aeolus laser occurred in June 2019. Over the next four years measurements could be conducted on 16 occasions. For the purpose of this analysis, the best nights without clouds or large aerosol contamination were selected. This resulted in three overpasses—one in each year between 2019 and 2021, specifically on 3 August 2019, 27 June 2020, and 17 July 2021. No overpass in 2022 met all our requirements. While the Aeolus laser was seen by several of the telescopes as it passed over the array (this is especially the case for the 2021 overpass; see also Fig. 4), for the energy reconstruction that is presented in the following sections we only used data from one set of telescopes (the northernmost FD site) that observed all three laser overpasses discussed in this paper. In this way, we minimized the calibration uncertainty between the individual telescopes and ensured an optimal comparability between the different data sets.

 figure: Fig. 4.

Fig. 4. Measured ground tracks for an altitude of 1400 m for three sample Aeolus overpasses in the years 2019, 2020, and 2021. Each FD site has a 180° FoV and is oriented to overlook the enclosed area. Consequently, the events in 2019 and 2020 were observed primarily by the northern and southern sites, as the eastern site is situated too far away. Additionally, a fraction of the 2019 events was visible to the western FD site, where the southernmost events appeared at the edge of its FoV.

Download Full Size | PDF

3. RESULTS

After applying the aforementioned methods for reconstructing the laser beam geometry, the ground tracks formed by the impacting laser beams can be obtained. It should be noted that while the geometry of the laser beams is fully determined by the fits described in Supplement 1, a calculation of the true impact point on the ground would require a detailed surface elevation model to be included, which does not apply to these calculations. Instead, the ground tracks can be defined for a certain altitude; thus the ground position of an event can be propagated along the reconstructed axis to match this altitude, to ensure comparability of the reconstructed positions.

A. Ground Track

Figure 4 shows the reconstructed ground tracks for an altitude of 1400 m (WGS84) for three sample overpasses in the years 2019, 2020, and 2021. Notable is the change of the laser ground position for the year 2021 further to the east after an adjustment of the satellite orbit, which was obtained via a dedicated maneuver to allow for better measurements of the laser beam by the Observatory. The passages during the former two years happened at larger distances and thus fewer events and shorter tracks were observed.

To study the accuracy of these reconstructions, simulated laser events were employed. Using the reconstructed event positions as in Fig. 4, a straight line can be fitted through these positions. On this line then many equidistant positions can be chosen that act as locations of the simulated laser beams. The inclination of the simulated laser shots is likewise obtained from the geometrical reconstruction of measured laser shots. In this manner, a realistic scenario can be simulated. In these simulations the light flux and camera response are calculated, so that the same reconstruction mechanisms can be used as for measured data, allowing for a direct comparison. The majority of measured events (extending beyond the $1\sigma$ range) are reconstructed within 200 m of the simulated position.

When comparing the Aeolus track position determined by the Auger Observatory with that in the Aeolus data products, a horizontal offset of about 0.075° longitude (6.8 km) became evident. Figure 5 shows an example of an overpass on 17 July 2021, but the offset was observed systematically for all overpasses. In order to facilitate a comparison between the two methods, a reference altitude was selected for which the positions are evaluated. The tracks presented in Fig. 5 were compared at an altitude of 10 km.

 figure: Fig. 5.

Fig. 5. Comparison between the geolocation derived with different Aeolus data processor versions before (L1A v7.11.2, orange line) and after (L1A v7.12, purple line) the fix of the error and measurements by the Auger Observatory (blue dots) at 10 km altitude for the overpass on 17 July 2021. The offset amounts to 6.8 km between Auger and v7.11.2, and 0.8 km between Auger and v7.12.

Download Full Size | PDF

The offset could be traced back to an error in the Aeolus Level 1A (L1A) processor for the calculation of the geolocation of the Aeolus observations, which was fixed in version 7.12 (Baseline 14). The L1A processor [25] uses the Earth Observation Mission Customer Furnished Item (CFI) Software [26] provided by ESA to calculate the longitude/latitude position and altitude values. The relevant CFI routines use for their calculation the time and an identifier stating if the time is UTC (Coordinated Universal Time), GPS (Global Position System), or TAI (International Atomic Time). In two places a wrong combination of time and identifier had been provided to the CFI routines, which led to a slightly wrong calculation of the longitude and latitude values. However, as the attitude angles of the platform, which are relevant for the wind retrieval, were not affected by the timing issue, there was no impact on the accuracy of the wind measurements.

It should be noted here that this geolocation error of about 6.8 km was not identified by other ground-based or airborne validation measurements of vertical profiles for wind speed and aerosol/cloud optical properties. Typically, a horizontal co-location requirement of 50–150 km was used for the validation of Aeolus products, and the horizontal change of atmospheric properties would not be significant enough on scales of only a few km to be detectable with the Aeolus horizontal resolution. Also, comparisons of the Aeolus wind profiles with outputs from numerical weather-prediction models did not reveal this horizontal geolocation error, as the typical grid size of those models is of the order of 3–10 km, but typical horizontal resolution is about four to eight times this grid spacing [35].

Another independent method for the validation of the geolocation including altitude would be the use of a digital elevation model (DEM) as demonstrated for airborne measurements [27,28]. Here the derived altitude of the surface return from the lidar is compared with the altitude from the DEM for this specific location. By using DEMs with horizontal resolution of the order of 100 m at locations with strongly varying altitudes (e.g., mountain regions) one could potentially identify errors in the retrieval of the horizontal geolocation even with the challenge of seven free parameters (three coordinates, three angles, one timing parameter) in the underlying equations.

 figure: Fig. 6.

Fig. 6. Reconstructed energies for three sample Aeolus overpasses in 2019, 2020, and 2021. The average energy per overpass is marked by the dashed line. The energy is given in mJ by the top axis and in PeV (${10^{15}}\,\,\rm eV$) by the bottom axis.

Download Full Size | PDF

The Auger measurements of the Aeolus ground-track positions also allow the calculation of the accuracy and precision of the Aeolus pointing. Auger measurements are available for each Aeolus pulse. However, the Aeolus ground track is reported in the Aeolus L1A product only on measurement level (mean over 30 pulses). To be able to compare the two, the Aeolus pulse times are calculated by using the laser pulse frequency of 50.5 Hz and the measurement centroid time from the L1A product. The ground track geolocations are then calculated by performing a time-based interpolation of the longitude and latitude values reported in the L1A product with respect to the calculated pulse times. These times, longitude, and latitude values are then compared with the results from Auger to assess the pointing accuracy and precision of Aeolus.

For 17 July 2021 and the Baseline 14 data shown in Fig. 5, a pointing accuracy of 0.06 km along track and 0.82 km across track was determined. The precision ($2\sigma$) is 1.28 km along track and 0.93 km across track. These values are a composite of Aeolus pointing errors, interpolation errors, and Auger measurement errors and, thus, are upper limits for the Aeolus pointing accuracy and precision. Nevertheless, they are well within the Aeolus mission requirements of 2.0 km ($2\sigma$) for the horizontal geolocation defined on observation level (i.e., mean over 600 pulses).

B. Laser Energy

As described above, the measured laser data can also be used to reconstruct the laser beam energy. An event-wise energy reconstruction of each laser transmission can be obtained. Analyzing the same three sample nights as for the geometry reconstruction from the years 2019, 2020, and 2021, the three energy distributions shown in Fig. 6 are obtained. The histograms were scaled to show the relative instead of the absolute number of events per bin for the sake of better visibility. As mentioned above, the years differ in the total number of measured events due to the changing laser-telescope distance. Most notable is the decreasing energy across the years, which is discussed further below.

Again, the simulated laser shots were used to evaluate the accuracy of the energy reconstruction. With the same simulation setup as for the ground track geometry study, the energies of simulated events were reconstructed and compared to the simulation input. Two simulation studies were made, one with the track position of the 2019/2020 case and a simulated energy of 34 mJ and one with the 2021 position and an energy of 20 mJ. A small systematic bias of the reconstructed energies towards lower values was found. With slight differences depending on the track positions and between different FD sites, the bias between simulated and reconstructed energy is up to 3.7%. The reconstructed energy values for the three sample overpasses are corrected for the bias and are

$$\begin{split}{E(2019)}&={ 33.1_{- 0.8}^{+ 1.9}\;{\rm mJ} ,}\\{E(2020)}&={ 23.7_{- 0.6}^{+ 1.7}\;{\rm mJ} ,}\\{E(2021)}&={ 17.3_{- 0.4}^{+ 0.9}\;{\rm mJ}.}\end{split}$$

These values are obtained using the absolute calibration of the fluorescence telescopes [29,30] and are fully corrected for atmospheric attenuation occurring after the laser beam exited the telescope, as implied by Eq. (S3) in Supplement 1. The quoted uncertainties are statistical uncertainties and contain contributions from the standard error of the mean energy of the events of one overpass and the statistical uncertainty of the aerosol optical depth at the time of observation.

In addition to these statistical uncertainties, the systematic energy-scale uncertainty of the fluorescence telescopes also needs to be taken into account, which amounts to 13%. It is dominated by the uncertainties of the absolute calibration and aerosol optical depth. Furthermore, it contains contributions from the light-collection efficiency, reconstruction bias, molecular atmosphere, multiple scattering, and the long-term calibration stability [31,32]. This energy-scale uncertainty is correlated across the individual years, i.e., the energy determined in each year could be under- or overestimated by the same amount within this uncertainty. Therefore, in the relative time evolution discussed in the following this does not add to the uncertainty, but instead is cancelled out when formulating the energy ratios.

C. Aeolus Signal Evolution

The reconstructed energies from the Auger Observatory allowed for an independent assessment of the in-orbit instrument performance of ALADIN between 2019 and 2021. A few weeks before the first Auger measurement in August 2019, a switchover to the redundant FM-B laser was performed after the emitted energy of the nominal FM-A laser had decreased significantly from 65 to 41 mJ between December 2018 and May 2019 [16], as recorded by the monitoring photodiode located behind the laser transmitter. The energy degradation was later traced back to a gradual misalignment of the master oscillator caused by a thermal drift of the laser bench, which resulted in a non-optimal temperature set point of the master oscillator depending on the emitted laser frequency. After its switch-on in late June 2019, the FM-B laser delivered an initial energy of 67 mJ and significantly slower power degradation, ensuring a higher signal-to-noise ratio (SNR) of the atmospheric backscatter return and, hence, lower random error for the wind observations. However, despite the stable laser energy, which was even increased to more than 100 mJ by several laser adjustments, the atmospheric signal levels decreased by more than 70% over the three years following the switchover, leading to a degrading precision of the Aeolus data with the random error of the Rayleigh wind results increasing from about 5 m/s to more than 8 m/s. In October/November 2022, the instrument was switched back to the FM-A laser whose performance could be optimized based on the knowledge gained with the second laser during the nearly 40 months of its operation. As a result, stable output energy above 50 mJ was achieved until the end of the nominal Aeolus operations on 30 April 2023, which was followed by a period of seven weeks for dedicated instrument tests.

 figure: Fig. 7.

Fig. 7. Signal evolution of the Aeolus instrument: the purple curve represents the laser energy as measured at the output of the respective transmitter [flight models (FMs) A and B]. The red curve denotes the atmospheric return signal that is detected on the Rayleigh receiver under clear-air conditions. The black dots indicate the Auger measurements with error bars that describe only their statistical errors (see text).

Download Full Size | PDF

Figure 7 depicts the signal evolution of ALADIN over the course of the Aeolus mission. The laser energy that is measured on a photodiode at the output of the transmitter is shown in purple, while the atmospheric return signal detected on the receiver is represented in red. The latter is obtained from the Rayleigh channel after selecting observations under clear-air conditions to ensure almost pure molecular backscatter signals. The energy levels estimated from the Auger laser-beam reconstruction are plotted as black dots. To compare the development of the signal energies measured at the different locations along the laser-beam path (output of the laser transmitter, Auger reconstruction of the satellite-telescope output, receiver), the signal levels in Fig. 7 are normalized to the respective values from 3 August 2019, i.e., the date of the first Auger measurement. As a consequence, the systematic calibration uncertainty of the normalized Auger energy estimates (13%) cancels out, leaving only the statistical error given in Eq. (1). Also, the uncertainty of the first Auger measurement (normalized to one) propagates into the uncertainty of the subsequent relative energy points.

The progression of the reconstructed energies from the Auger Observatory is in good agreement with the decline in the atmospheric signal levels registered at the ALADIN detectors between 2019 and 2021. This result implies that the signal loss observed during the operation of the FM-B laser occurred along the emit path of the instrument, i.e., between the laser output and the telescope. The first Auger estimates in July/August 2021 strongly supported the root-cause analysis of the signal loss and ultimately led to the decision to switch back to the FM-A laser in November 2022. Despite the much lower output energy compared to $\text{FM-B}$ at the end of its operation (101 mJ), the switchover to FM-A increased the atmospheric signal by a factor of 2.2, to a level corresponding to the laser energy of the first FM-A period in early 2019. The signal loss experienced during the FM-B period was thus fully recovered, confirming that it had occurred on the optics behind the laser output, more specifically on the optics that are unique to the FM-B, most likely within the relay optics, including the FFM, which guide the redundant FM-B laser beam onto the nominal optical path (see Fig. 1). The actual loss mechanism is the subject of ongoing studies as of the writing of the paper. Laser-induced contamination and laser-induced damage of the affected optics are currently assessed to be the most probable causes.

The fact that the energy decrease observed by the Auger Observatory in 2020 (${-}28\%$) and 2021 (${-}48\%$) with respect to the reference measurement in 2019 is smaller than the signal loss registered on the ALADIN detectors (${-}34\%$ and ${-}53\%$) points to an additional loss mechanism in the receive path, most likely clipping of the atmospheric return signal at the field stop.

It is also notable that the absolute laser energy in 2019 reconstructed by the Auger Observatory ($\approx 33\,\,\rm mJ$) is lower than what would be expected at the Aeolus telescope output (44 mJ) when considering the specified emit-path transmission of around 0.704 and correction factors for the laser energy. This points to additional transmission loss in the emit path of ALADIN. An initial loss of about a factor of two compared to pre-launch simulations was already observed after launch for the atmospheric-backscatter signals of the Rayleigh channel [33].

4. OUTLOOK

The methods presented in this study have proven to be highly valuable for the analysis and validation of satellite lidar parameters. Furthermore, the measurements of the satellite laser beams by the Pierre Auger Observatory offer some promising applications that extend beyond the scope of this paper. One possible future application is the use of the measured laser light to gain knowledge about the state of the atmosphere above the Observatory. Since the satellite laser passes over the entire Observatory within a few seconds, this provides a wide range of measurements from many different distances and angles. Assuming a constant laser energy for the duration of the passage, this can be used to obtain information about the aerosol load in the atmosphere. Similar to the algorithm described above, a likelihood fit can be devised that includes both the laser energy and the variables of a parametric aerosol model as adjustable quantities. The potential of using the complete measured data of one overpass for such a combined fit of energy and aerosol parameters is being investigated. While this method of aerosol determination is not suitable to serve as a regular monitoring device due to the scarce availability of Aeolus overpasses, it could provide an opportunity for individual cross-checks with the atmospheric monitoring devices that are employed at the Pierre Auger Observatory. A first analysis of this type was presented in Ref. [34].

Another exciting prospect are future space lidar missions like EarthCARE, a satellite with the goal of measuring global cloud and aerosol profiles, which will include, among other instruments, a 355 nm lidar [35]. In addition, Aeolus will have a follow-up mission, Aeolus-2, which is under development [36] and planned to be launched in the early 2030s. These future space lidar missions are not only intriguing because measurements like the ones shown in this paper can be repeated and the data set of observatory-measured satellite-laser events can be extended, but also due to a different geometry of the satellite orbits. Due to its dawn-dusk orbit with 6 am/pm equator crossing times Aeolus was only visible in the observatory during a few months in the southern-hemisphere winter. A similar potential measurement in the northern hemisphere, more specifically at the Telescope Array Project in Utah, United States [37], was investigated, but found not to be possible due to the passage time of Aeolus [38]. For the mentioned upcoming EarthCARE mission, however, thanks to its 2 am/pm orbit, the possibility exists that the satellite laser will be seen in both the Telescope Array and the Pierre Auger Observatory within a short time frame and throughout the year. This would allow for the unique opportunity of directly comparing the energy calibration of both experiments and to significantly reduce existing uncertainties in the cosmic-ray flux and anisotropy [39,40].

Funding

European Space Agency (4000126336/18/I-BG); United Nations Educational, Scientific and Cultural Organization (UNESCO); European Particle Physics Latin American Network; Marie Curie-IRSES/EPLANET; Grainger Foundation; National Science Foundation (0450696); U.S. Department of Energy (DE-AC02-07CH11359, DE-FG02-99ER41107, DE-FR02-04ER41300, DE-SC0011689); María de Maeztu Unit of Excellence (MDM-2016-0692); Red Nacional Temática de Astropartículas (FPA2015-68783-REDT); Federación Española de Enfermedades Raras; Junta de Andalucía (P18-FR-4314, SOMM17/6104/UGR); Xunta de Galicia (ED431C 2017/07); Ministerio de Economía y Competitividad (FPA2017-85114-P, PID2019-104676GB-C32); Javna Agencija za Raziskovalno Dejavnost RS (I0-0033, N1-0111, P1-0031, P1-0385); Romanian National Core Program LAPLAS VII (PN 23 21 01 02, PN-III-P1-1.1-TE-2021-0924/TE57/2022); CNCS-UEFISCDI (30N/2023); Ministerul Cercetării, Inovării şi Digitalizării; Fundação para a Ciência e a Tecnologia; Narodowe Centrum Nauki (2016/22/M/ST9/00198, 2016/23/B/ST9/01635, 2020/39/B/ST9/01398, 2022/45/B/ST9/02163); Ministerstwo Edukacji i Nauki (2022/WK/12, DIR/WK/2018/11); Dutch national e-infrastructure with the support of SURF Cooperative; Nederlandse Organisatie voor Wetenschappelijk Onderzoek; Ministerie van Onderwijs, Cultuur en Wetenschap; Dirección General de Asuntos del Personal Académico, Universidad Nacional Autónoma de México; Universidad Nacional Autónoma de México; Consejo Nacional de Ciencia y Tecnología (167733); Big Data and Quantum Computing, NextGenerationEU (CN_00000013); ICSC Centro Nazionale di Ricerca in High Performance Computing; Ministero degli Affari Esteri e della Cooperazione Internazionale; CETEMPS Center of Excellence; Ministero dell’Università e della Ricerca; Istituto Nazionale di Astrofisica; Istituto Nazionale di Fisica Nucleare Sezione di Padova; Ministerium für Wissenschaft, Forschung und Kunst Baden-Württemberg; Ministerium für Kultur und Wissenschaft des Landes Nordrhein-Westfalen; Helmholtz-Gemeinschaft Deutscher Forschungszentren (HGF); Helmholtz Alliance for Astroparticle Physics; Ministerium für Finanzen Baden-Württemberg; Deutsche Forschungsgemeinschaft; Bundesministerium für Bildung und Forschung; Institut Lagrange de Paris (ANR-11-IDEX-0004-02, LABEX ANR-10-LABX-63); Département Sciences de l’Univers (SDU-INSU/CNRS); Institut National de Physique Nucléaire et de Physique des Particules; Conseil Régional Ile-de-France; Centre National de la Recherche Scientifique; Centre de Calcul (IN2P3/CNRS); Czech Republic (CZ.02.1.01/0.0/0.0/16_013/0001402, CZ.02.1.01/0.0/0.0/17_049/0008422, CZ.02.1.01/0.0/0.0/18_046/0016010, LM2015038, LM2018102, MSMT CR LTT18004); Ministério da Ciência, Tecnologia, Inovações e Comunicações; Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro; Financiadora de Estudos e Projetos; Conselho Nacional de Desenvolvimento Científico e Tecnológico; Fonds Wetenschappelijk Onderzoek; Fonds De La Recherche Scientifique - FNRS; Australian Research Council; NDM Holdings and Valle Las Leñas; Gobierno de la Provincia de Mendoza; Municipalidad de Malargüe; Consejo Nacional de Investigaciones Científicas y Técnicas; Agencia Nacional de Promoción Cientfica y Tecnológica; Comisión Nacional de Energía Atómica, Gobierno de Argentina; Fundação de Amparo à Pesquisa do Estado de São Paulo (1999/05404-3, 2010/07359-6, 2019/10151-2); FP7 People: Marie-Curie Actions (101107047).

Acknowledgment

The successful installation, commissioning, and operation of the Pierre Auger Observatory would not have been possible without the strong commitment and effort from the technical and administrative staff in Malargüe.

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Supplemental document

See Supplement 1 for supporting content.

The Pierre Auger Collaboration

The Pierre Auger Collaboration

A. Abdul Halim,13 P. Abreu,73 M. Aglietta,55,53 I. Allekotte,1 K. Almeida Cheminant,71 A. Almela,7,12 R. Aloisio,46,47 J. Alvarez-Muñiz,80 J. Ammerman Yebra,80 G.A. Anastasi,55,53 L. Anchordoqui,87 B. Andrada,7 S. Andringa,73 Anukriti,77 L. Apollonio,60,50 C. Aramo,51 P.R. Araújo Ferreira,43 E. Arnone,64,53 J. C. Arteaga Velázquez,68 P. Assis,73 G. Avila,11 E. Avocone,58,47 A.M. Badescu,76 A. Bakalova,33 F. Barbato,46,47 A. Bartz Mocellin,86 J.A. Bellido,13,70 C. Berat,37 M.E. Bertaina,64,53 G. Bhatta,71 M. Bianciotto,64,53 P.L. Biermann,i V. Binet,5 K. Bismark,40,7 T. Bister,81,82 J. Biteau,38,b J. Blazek,33 C. Bleve,37 J. Blümer,42 M. Boháčová,33 D. Boncioli,58,47 C. Bonifazi,8,27 L. Bonneau Arbeletche,22 N. Borodai,71 J. Brack,k P.G. Brichetto Orchera,7 F.L. Briechle,43 A. Bueno,79 S. Buitink,15 M. Buscemi,48,62 A. Bwembya,81,82 M. Büsken,40,7 K.S. Caballero-Mora,67 S. Cabana-Freire,80 L. Caccianiga,60,50 R. Caruso,59,48 A. Castellina,55,53 F. Catalani,19 G. Cataldi,49 L. Cazon,80 M. Cerda,10 A. Cermenati,46,47 J.A. Chinellato,22 J. Chudoba,33 L. Chytka,34 R.W. Clay,13 A.C. Cobos Cerutti,6 R. Colalillo,61,51 A. Coleman,91 M.R. Coluccia,49 R. Conceição,73 A. Condorelli,38 G. Consolati,50,56 M. Conte,57,49 F. Convenga,58,47 D. Correia dos Santos,29 P.J. Costa,73 C.E. Covault,85 M. Cristinziani,45 C.S. Cruz Sanchez,3 S. Dasso,4,2 K. Daumiller,42 B.R. Dawson,13 R.M. de Almeida,29 J. de Jesús,7,42 S.J. de Jong,81,82 J.R.T. de Mello Neto,27,28 I. De Mitri,46,47 J. de Oliveira,18 D. de Oliveira Franco,22 F. de Palma,57,49 V. de Souza,20 B.P. de Souza de Errico,27 E. De Vito,57,49 A. Del Popolo,59,48 O. Deligny,35 N. Denner,33 L. Deval,42,7 A. di Matteo,53 M. Dobre,74 C. Dobrigkeit,22 J.C. D’Olivo,69 L.M. Domingues Mendes,73 Q. Dorosti,45 J.C. dos Anjos,16 R.C. dos Anjos,26 J. Ebr,33 F. Ellwanger,42 M. Emam,81,82 R. Engel,40,42 I. Epicoco,57,49 M. Erdmann,43 A. Etchegoyen,7,12 C. Evoli,46,47 H. Falcke,81,83,82 J. Farmer,90 G. Farrar,89 A.C. Fauth,22 N. Fazzini,f F. Feldbusch,41 F. Fenu,42,e A. Fernandes,73 B. Fick,88 J.M. Figueira,7 A. Filipčič,78,77 T. Fitoussi,42 B. Flaggs,91 T. Fodran,81 T. Fujii,90,g A. Fuster,7,12 C. Galea,81 C. Galelli,60,50 B. Garca,6 C. Gaudu,39 H. Gemmeke,41 F. Gesualdi,7,42 A. Gherghel-Lascu,74 P.L. Ghia,35 U. Giaccari,49 J. Glombitza,43,h F. Gobbi,10 F. Gollan,7 G. Golup,1 J.P. Gongora,11 J.M. González,1 N. González,7 I. Goos,1 A. Gorgi,55,53 M. Gottowik,80 T.D. Grubb,13 F. Guarino,61,51 G.P. Guedes,23 E. Guido,45 M. Gómez Berisso,1 P.F. Gómez Vitale,11 D. Góra,71 L. Gülzow,42 S. Hahn,40 P. Hamal,33 M.R. Hampel,7 P. Hansen,3 D. Harari,1 V.M. Harvey,13 A. Haungs,42 T. Hebbeker,43 C. Hojvat,f P. Horvath,34 M. Hrabovský,34 T. Huege,42,15 J.R. Hörandel,81,82 A. Insolia,59,48 P.G. Isar,75 P. Janecek,33 J.A. Johnsen,86 J. Jurysek,33 K.H. Kampert,39 B. Keilhauer,42 A. Khakurdikar,81 V.V. Kizakke Covilakam,7,42 H.O. Klages,42 M. Kleifges,41 F. Knapp,40 N. Kunka,41 J. Köhler,42 B.L. Lago,17 N. Langner,43 M.A. Leigui de Oliveira,25 Y. Lema-Capeans,80 A. Letessier-Selvon,36 I. Lhenry-Yvon,35 L. Lopes,73 L. Lu,92 Q. Luce,40 J.P. Lundquist,77 A. Machado Payeras,22 M. Majercakova,33 D. Mandat,33 B.C. Manning,13 P. Mantsch,f S. Marafico,35 F.M. Mariani,60,50 A.G. Mariazzi,3 I.C. Mariş,14 G. Marsella,62,48 D. Martello,57,49 S. Martinelli,42,7 M.A. Martins,80 O. Martínez Bravo,65 H.J. Mathes,42 J. Matthews,a G. Matthiae,63,52 E. Mayotte,86,39 S. Mayotte,86 P.O. Mazur,f G. Medina-Tanco,69 J. Meinert,39 D. Melo,7 A. Menshikov,41 C. Merx,42 S. Michal,34 M.I. Micheletti,5 L. Miramonti,60,50 S. Mollerach,1 F. Montanet,37 L. Morejon,39 C. Morello,55,53 K. Mulrey,81,82 R. Mussa,53 W.M. Namasaka,39 S. Negi,33 L. Nellen,69 K. Nguyen,88 G. Nicora,9 M. Niechciol,45 D. Nitz,88 D. Nosek,32 V. Novotny,32 L. Nožka,34 A. Nucita,57,49 L.A. Núñez,31 C. Oliveira,20 M. Palatka,33 J. Pallotta,9 S. Panja,33 G. Parente,80 T. Paulsen,39 J. Pawlowsky,39 M. Pech,33 R. Pelayo,66 L.A.S. Pereira,24 E.E. Pereira Martins,40,7 J. Perez Armand,21 L. Perrone,57,49 S. Petrera,46,47 C. Petrucci,58,47 T. Pierog,42 M. Pimenta,73 M. Platino,7 B. Pont,81 M. Pothast,82,81 M. Pourmohammad Shahvar,62,48 P. Privitera,90 M. Prouza,33 A. Puyleart,88 C. Pérez Bertolli,7,42 J. Pękala,71 S. Querchfeld,39 J. Rautenberg,39 D. Ravignani,7 J.V. Reginatto Akim,22 M. Reininghaus,40 J. Ridky,33 F. Riehn,80 M. Risse,45 V. Rizi,58,47 W. Rodrigues de Carvalho,81 E. Rodriguez,7,42 J. Rodriguez Rojo,11 M.J. Roncoroni,7 S. Rossoni,44 M. Roth,42 E. Roulet,1 A.C. Rovero,4 P. Ruehl,45 A. Saftoiu,74 M. Saharan,81 F. Salamida,58,47 H. Salazar,65 G. Salina,52 J.D. Sanabria Gomez,31 E.M. Santos,21 E. Santos,33 F. Sarazin,86 R. Sarmento,73 R. Sato,11 P. Savina,92 V. Scherini,57,49 H. Schieler,42 M. Schimassek,35 M. Schimp,39 D. Schmidt,42 O. Scholten,15,j H. Schoorlemmer,81,82 P. Schovánek,33 F.G. Schröder,91,42 J. Schulte,43 T. Schulz,42 C.M. Schäfer,40 S.J. Sciutto,3 M. Scornavacche,7,42 A. Segreto,54,48 S. Sehgal,39 S.U. Shivashankara,77 G. Sigl,44 G. Silli,7 O. Sima,74,c K. Simkova,15 F. Simon,41 R. Smau,74 R. Šmída,90 P. Sommers,l J.F. Soriano,87 R. Squartini,10 M. Stadelmaier,33 S. Stanič,77 J. Stasielak,71 P. Stassi,37 M. Straub,43 S. Strähnz,40 T. Suomijärvi,38 A.D. Supanitsky,7 Z. Svozilikova,33 Z. Szadkowski,72 F. Sánchez,7 F. Tairli,13 A. Tapia,30 C. Taricco,64,53 C. Timmermans,82,81 O. Tkachenko,42 P. Tobiska,33 C.J. Todero Peixoto,19 B. Tomé,73 Z. Torrès,37 A. Travaini,10 P. Travnicek,33 C. Trimarelli,58,47 M. Tueros,3 M. Unger,42 L. Vaclavek,34 M. Vacula,34 J.F. Valdés Galicia,69 L. Valore,61,51 E. Varela,65 D. Veberič,42 C. Ventura,28 I.D. Vergara Quispe,3 V. Verzi,52 J. Vicha,33 J. Vink,84 J. Vlastimil,33 S. Vorobiov,77 A. Vásquez-Ramírez,31 C. Watanabe,27 A.A. Watson,d A. Weindl,42 L. Wiencke,86 H. Wilczyński,71 D. Wittkowski,39 B. Wundheiler,7 B. Yue,39 A. Yushkov,33 O. Zapparrata,14 E. Zas,80 D. Zavrtanik,77,78 M. Zavrtanik78,77

1 Centro Atómico Bariloche and Instituto Balseiro (CNEA-UNCuyo-CONICET), San Carlos de Bariloche, Argentina

2 Departamento de Física and Departamento de Ciencias de la Atmósfera y los Océanos, FCEyN, Universidad de Buenos Aires and CONICET, Buenos Aires, Argentina

3 IFLP, Universidad Nacional de La Plata and CONICET, La Plata, Argentina

4 Instituto de Astronomía y Física del Espacio (IAFE, CONICET-UBA), Buenos Aires, Argentina

5 Instituto de Física de Rosario (IFIR) – CONICET/U.N.R. and Facultad de Ciencias Bioquímicas y Farmacéuticas U.N.R., Rosario, Argentina

6 Instituto de Tecnologas en Detección y Astropartículas (CNEA, CONICET, UNSAM), and Universidad Tecnológica Nacional – Facultad Regional Mendoza (CONICET/CNEA), Mendoza, Argentina

7 Instituto de Tecnologías en Detección y Astropartículas (CNEA, CONICET, UNSAM), Buenos Aires, Argentina

8 International Center of Advanced Studies and Instituto de Ciencias Físicas, ECyT-UNSAM and CONICET, Campus Miguelete – San Martín, Buenos Aires, Argentina

9 Laboratorio Atmósfera – Departamento de Investigaciones en Láseres y sus Aplicaciones – UNIDEF (CITEDEF-CONICET), Argentina

10 Observatorio Pierre Auger, Malargüe, Argentina

11 Observatorio Pierre Auger and Comisión Nacional de Energía Atómica, Malargüe, Argentina

12 Universidad Tecnológica Nacional – Facultad Regional Buenos Aires, Buenos Aires, Argentina

13 University of Adelaide, Adelaide, S.A., Australia

14 Université Libre de Bruxelles (ULB), Brussels, Belgium

15 Vrije Universiteit Brussels, Brussels, Belgium

16 Centro Brasileiro de Pesquisas Fisicas, Rio de Janeiro, RJ, Brazil

17 Centro Federal de Educação Tecnológica Celso Suckow da Fonseca, Petropolis, Brazil

18 Instituto Federal de Educação, Ciência e Tecnologia do Rio de Janeiro (IFRJ), Brazil

19 Universidade de São Paulo, Escola de Engenharia de Lorena, Lorena, SP, Brazil

20 Universidade de São Paulo, Instituto de Física de São Carlos, São Carlos, SP, Brazil

21 Universidade de São Paulo, Instituto de Física, São Paulo, SP, Brazil

22 Universidade Estadual de Campinas, IFGW, Campinas, SP, Brazil

23 Universidade Estadual de Feira de Santana, Feira de Santana, Brazil

24 Universidade Federal de Campina Grande, Centro de Ciencias e Tecnologia, Campina Grande, Brazil

25 Universidade Federal do ABC, Santo André, SP, Brazil

26 Universidade Federal do Paraná, Setor Palotina, Palotina, Brazil

27 Universidade Federal do Rio de Janeiro, Instituto de Física, Rio de Janeiro, RJ, Brazil

28 Universidade Federal do Rio de Janeiro (UFRJ), Observatório do Valongo, Rio de Janeiro, RJ, Brazil

29 Universidade Federal Fluminense, EEIMVR, Volta Redonda, RJ, Brazil

30 Universidad de Medellín, Medellín, Colombia

31 Universidad Industrial de Santander, Bucaramanga, Colombia

32 Charles University, Faculty of Mathematics and Physics, Institute of Particle and Nuclear Physics, Prague, Czech Republic

33 Institute of Physics of the Czech Academy of Sciences, Prague, Czech Republic

34 Palacky University, Olomouc, Czech Republic

35 CNRS/IN2P3, IJCLab, Université Paris-Saclay, Orsay, France

36 Laboratoire de Physique Nucléaire et de Hautes Energies (LPNHE), Sorbonne Université, Université de Paris, CNRS-IN2P3, Paris, France

37 Univ. Grenoble Alpes, CNRS, Grenoble Institute of Engineering Univ. Grenoble Alpes, LPSC-IN2P3, 38000 Grenoble, France

38 Université Paris-Saclay, CNRS/IN2P3, IJCLab, Orsay, France

39 Bergische Universität Wuppertal, Department of Physics, Wuppertal, Germany

40 Karlsruhe Institute of Technology (KIT), Institute for Experimental Particle Physics, Karlsruhe, Germany

41 Karlsruhe Institute of Technology (KIT), Institut für Prozessdatenverarbeitung und Elektronik, Karlsruhe, Germany

42 Karlsruhe Institute of Technology (KIT), Institute for Astroparticle Physics, Karlsruhe, Germany

43 RWTH Aachen University, III. Physikalisches Institut A, Aachen, Germany

44 Universität Hamburg, II. Institut für Theoretische Physik, Hamburg, Germany

45 Universität Siegen, Department Physik – Experimentelle Teilchenphysik, Siegen, Germany

46 Gran Sasso Science Institute, L’Aquila, Italy

47 INFN Laboratori Nazionali del Gran Sasso, Assergi (L’Aquila), Italy

48 INFN, Sezione di Catania, Catania, Italy

49 INFN, Sezione di Lecce, Lecce, Italy

50 INFN, Sezione di Milano, Milano, Italy

51 INFN, Sezione di Napoli, Napoli, Italy

52 INFN, Sezione di Roma “Tor Vergata”, Roma, Italy

53 INFN, Sezione di Torino, Torino, Italy

54 Istituto di Astrofisica Spaziale e Fisica Cosmica di Palermo (INAF), Palermo, Italy

55 Osservatorio Astrofisico di Torino (INAF), Torino, Italy

56 Politecnico di Milano, Dipartimento di Scienze e Tecnologie Aerospaziali , Milano, Italy

57 Università del Salento, Dipartimento di Matematica e Fisica “E. De Giorgi”, Lecce, Italy

58 Università dell’Aquila, Dipartimento di Scienze Fisiche e Chimiche, L’Aquila, Italy

59 Università di Catania, Dipartimento di Fisica e Astronomia “Ettore Majorana”, Catania, Italy

60 Università di Milano, Dipartimento di Fisica, Milano, Italy

61 Università di Napoli “Federico II”, Dipartimento di Fisica “Ettore Pancini”, Napoli, Italy

62 Università di Palermo, Dipartimento di Fisica e Chimica ”E. Segrè”, Palermo, Italy

63 Università di Roma “Tor Vergata”, Dipartimento di Fisica, Roma, Italy

64 Università Torino, Dipartimento di Fisica, Torino, Italy

65 Benemérita Universidad Autónoma de Puebla, Puebla, México

66 Unidad Profesional Interdisciplinaria en Ingeniería y Tecnologías Avanzadas del Instituto Politécnico Nacional (UPIITA-IPN), México, D.F., México

67 Universidad Autónoma de Chiapas, Tuxtla Gutiérrez, Chiapas, México

68 Universidad Michoacana de San Nicolás de Hidalgo, Morelia, Michoacán, México

69 Universidad Nacional Autónoma de México, México, D.F., México

70 Universidad Nacional de San Agustin de Arequipa, Facultad de Ciencias Naturales y Formales, Arequipa, Peru

71 Institute of Nuclear Physics PAN, Krakow, Poland

72 University of Łódź, Faculty of High-Energy Astrophysics, Łódź, Poland

73 Laboratório de Instrumentação e Física Experimental de Partículas – LIP and Instituto Superior Técnico – IST, Universidade de Lisboa – UL, Lisboa, Portugal

74 “Horia Hulubei” National Institute for Physics and Nuclear Engineering, Bucharest-Magurele, Romania

75 Institute of Space Science, Bucharest-Magurele, Romania

76 University Politehnica of Bucharest, Bucharest, Romania

77 Center for Astrophysics and Cosmology (CAC), University of Nova Gorica, Nova Gorica, Slovenia

78 Experimental Particle Physics Department, J. Stefan Institute, Ljubljana, Slovenia

79 Universidad de Granada and C.A.F.P.E., Granada, Spain

80 Instituto Galego de Física de Altas Enerxías (IGFAE), Universidade de Santiago de Compostela, Santiago de Compostela, Spain

81 IMAPP, Radboud University Nijmegen, Nijmegen, The Netherlands

82 Nationaal Instituut voor Kernfysica en Hoge Energie Fysica (NIKHEF), Science Park, Amsterdam, The Netherlands

83 Stichting Astronomisch Onderzoek in Nederland (ASTRON), Dwingeloo, The Netherlands

84 Universiteit van Amsterdam, Faculty of Science, Amsterdam, The Netherlands

85 Case Western Reserve University, Cleveland, OH, USA

86 Colorado School of Mines, Golden, CO, USA

87 Department of Physics and Astronomy, Lehman College, City University of New York, Bronx, NY, USA

88 Michigan Technological University, Houghton, MI, USA

89 New York University, New York, NY, USA

90 University of Chicago, Enrico Fermi Institute, Chicago, IL, USA

91 University of Delaware, Department of Physics and Astronomy, Bartol Research Institute, Newark, DE, USA

92 University of Wisconsin-Madison, Department of Physics and WIPAC, Madison, WI, USA

a Louisiana State University, Baton Rouge, LA, USA

b Institut universitaire de France (IUF), France

c also at University of Bucharest, Physics Department, Bucharest, Romania

d School of Physics and Astronomy, University of Leeds, Leeds, United Kingdom

e now at Agenzia Spaziale Italiana (ASI). Via del Politecnico 00133, Roma, Italy

f Fermi National Accelerator Laboratory, Fermilab, Batavia, IL, USA

g now at Graduate School of Science, Osaka Metropolitan University, Osaka, Japan

h now at ECAP, Erlangen, Germany

i Max-Planck-Institut für Radioastronomie, Bonn, Germany

j also at Kapteyn Institute, University of Groningen, Groningen, The Netherlands

k Colorado State University, Fort Collins, CO, USA

l Pennsylvania State University, University Park, PA, USA

REFERENCES

1. European Space Agency, “ADM-Aeolus science report,” ESA SP-1311 (2008).

2. O. Reitebuch, “The spaceborne wind lidar mission ADM-Aeolus,” in Atmospheric Physics: Background–Methods–Trends, U. Schumann, ed. (Springer, 2012), pp. 815–827.

3. A. Martin, M. Weissmann, and A. Cress, “Investigation of links between dynamical scenarios and particularly high impact of Aeolus on numerical weather prediction (NWP) forecasts,” Weather Clim. Dyn. 4, 249–264 (2023). [CrossRef]  

4. V. Pourret, M. Šavli, J.-F. Mahfouf, et al., “Operational assimilation of Aeolus winds in the Météo-France global NWP model ARPEGE,” Q. J. R. Meteorolog. Soc. 148, 2652–2671 (2022). [CrossRef]  

5. M. P. Rennie, L. Isaksen, F. Weiler, et al., “The impact of Aeolus wind retrievals on ECMWF global weather forecasts,” Q. J. R. Meteorolog. Soc. 147, 3555–3586 (2021). [CrossRef]  

6. D. M. Winker, R. H. Couch, and M. P. McCormick, “An overview of LITE: NASA’s lidar in-space technology experiment,” Proc. IEEE 84, 164–180 (1996). [CrossRef]  

7. J. B. Abshire, X. Sun, H. Riris, et al., “Geoscience laser altimeter system (GLAS) on the ICESat Mission: on-orbit measurement performance,” Geophys. Res. Lett. 32, 405 (2005). [CrossRef]  

8. D. Winker, M. Vaughan, and B. Hunt, “The CALIPSO mission and initial results from CALIOP,” Proc. SPIE 6409, 640902 (2006). [CrossRef]  

9. M. J. McGill, J. E. Yorks, V. S. Scott, et al., “The cloud-aerosol transport system (CATS): a technology demonstration on the International Space Station,” Proc. SPIE 9612, 96120A (2015). [CrossRef]  

10. A. J. Martino, J. Cavanaugh, A. Gibbons, et al., “ICESat-2/ATLAS at 4 years: instrument performance and projected life,” Proc. SPIE 12512, 125120A (2023). [CrossRef]  

11. A. Cosentino, A. Mondello, A. Sapia, et al., “High energy, single frequency, tunable laser source operating in burst mode for space based lidar applications,” Proc. SPIE 10568, 1056817 (2017). [CrossRef]  

12. L. Mondin and P. Bravetti, “Aeolus high energy UV laser wavelength measurement and frequency stability analysis,” Proc. SPIE 10563, 105633B (2017). [CrossRef]  

13. L. Mondin, F. F. Era, H. Kheyrandish, et al., “Caveats and pitfalls of high energy UV laser operation on ground,” Proc. SPIE 10562, 105621U (2017). [CrossRef]  

14. W. Riede, P. Allenspacher, A. Ciapponi, et al., “Qualification of laser optics for upcoming ESA space laser missions,” Proc. SPIE 10564, 105642K (2018). [CrossRef]  

15. N. Bartels, P. Allenspacher, W. Riede, et al., “Removal of laser-induced contamination on ALADIN laser optics by UV/ozone cleaning,” Proc. SPIE 11173, 111731C (2019). [CrossRef]  

16. O. Lux, D. Wernham, P. Bravetti, et al., “High-power and frequency-stable ultraviolet laser performance in space for the wind lidar on Aeolus,” Opt. Lett. 45, 1443–1446 (2020). [CrossRef]  

17. O. Lux, C. Lemmerz, F. Weiler, et al., “ALADIN laser frequency stability and its impact on the Aeolus wind error,” Atmos. Meas. Tech. 14, 6305–6333 (2021). [CrossRef]  

18. B. Witschas, M. O. Vieitez, E.-J. van Duijn, et al., “Spontaneous Rayleigh–Brillouin scattering of ultraviolet light in nitrogen, dry air, and moist air,” Appl. Opt. 49, 4217–4227 (2010). [CrossRef]  

19. B. Witschas, “Analytical model for Rayleigh–Brillouin line shapes in air,” Appl. Opt. 50, 267–270 (2011). [CrossRef]  

20. L. A. Anchordoqui, “Ultra-high-energy cosmic rays,” Phys. Rep. 801,1–93 (2019). [CrossRef]  

21. R. A. Batista, J. Biteau, M. Bustamante, et al., “Open questions in cosmic-ray research at ultrahigh energies,” Front. Astron. Space Sci. 6, 23 (2019). [CrossRef]  

22. A. Aab, P. Abreu, M. Aglietta, et al., “The Pierre Auger cosmic ray observatory,” Nucl. Instrum. Methods Phys. Res., Sect. A 798, 172–213 (2015). [CrossRef]  

23. M. Ave, M. Bohacova, B. Buonomo, et al., “Spectrally resolved pressure dependence measurements of air fluorescence emission with AIRFLY,” Nucl. Instrum. Meth. A 597, 41–45 (2008). [CrossRef]  

24. J. Abraham, P. Abreu, M. Aglietta, et al., “The fluorescence detector of the Pierre Auger Observatory,” Nucl. Instrum. Methods Phys. Res., Sect. A 620, 227–251 (2010). [CrossRef]  

25. O. Reitebuch, D. Huber, and I. Nikolaus, “ADM-aeolus algorithm theoretical basis document (ATBD) level-1B products, Version: 4.4,” 2018, https://earth.esa.int/eogateway/documents/20142/37627/Aeolus-L1B-Algorithm-ATBD.pdf.

26. European Space Agency, “Earth observation mission CFI software, Version: 4.12,” 2013, https://eop-cfi.esa.int/index.php/mission-cfi-software/eocfi-software.

27. F. Weiler, “Bias correction using ground echoes for the airborne demonstrator of the wind lidar on the ADM-Aeolus mission,” Master’s thesis (University of Innsbruck, 2017).

28. A. Amediek and M. Wirth, “Pointing verification method for spaceborne lidars,” Remote Sens. 9, 56 (2017). [CrossRef]  

29. J. T. Brack, R. Meyhandan, G. J. Hofman, et al., “Absolute photometric calibration of large aperture optical systems,” Astropart. Phys. 20, 653–659 (2004). [CrossRef]  

30. J. T. Brack, R. Cope, A. Dorofeev, et al., “Absolute calibration of a large-diameter light source,” JINST 8, P05014 (2013). [CrossRef]  

31. B. Dawson, “The energy scale of the Pierre Auger observatory,” in PoS International Cosmic Ray Conference (ICRC2019) (2020), p. 231.

32. V. Verzi, “The energy scale of the Pierre Auger observatory,” in Proceedings 33rd International Cosmic Ray Conference (2013), p. 376.

33. O. Reitebuch, C. Lemmerz, O. Lux, et al., “Initial assessment of the performance of the first wind lidar in space on Aeolus,” EPJ Web Conf. 237, 01010 (2020). [CrossRef]  

34. P. Abreu, M. Aglietta, J. M. Albury, et al., “Satellite data for atmospheric monitoring at the Pierre Auger observatory,” in PoS International Cosmic Ray Conference (ICRC2021) (2021), p. 235.

35. A. J. Illingworth, H. W. Barker, A. Beljaars, et al., “The EarthCARE Satellite: The next step forward in global measurements of clouds, aerosols, precipitation, and radiation,” Bull. Am. Meteorol. Soc. 96, 1311–1332 (2015). [CrossRef]  

36. A. Heliere, D. Wernham, G. Mason, et al., “Status of Aeolus-2 mission pre-development activities,” Proc. SPIE 12777, 1277709 (2023). [CrossRef]  

37. H. Kawai, S. Yoshida, H. Yoshii, et al., “Telescope array experiment,” Nucl. Phys. B Proc. Suppl. 175–176, 221–226 (2008). [CrossRef]  

38. F. Knapp, “Analysis of laser shots of the Aeolus satellite observed with the fluorescence telescopes of the Pierre Auger observatory,” Master’s thesis (Karlsruher Institut für Technologie, Institut für Astroteilchenphysik, 2021).

39. A. A. Halim, P. Abreu, M. Aglietta, et al., “Update on the searches for anisotropies in UHECR arrival directions with the Pierre Auger Observatory and the Telescope Array,” in PoS International Cosmic Ray Conference (ICRC2023) (2023), p. 521.

40. T. Tsunesada, “Measurement of UHECR energy spectrum with the Pierre Auger observatory and the telescope array,” in PoS International Cosmic Ray Conference (ICRC2023) (2023), p. 406.

Supplementary Material (1)

NameDescription
Supplement 1       Supplemental document

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (7)

Fig. 1.
Fig. 1. Simplified optical layout of the ALADIN instrument onboard Aeolus. HR: highly-reflective mirror; FFM: flip-flop mechanism; BS: beam splitter; QWP: quarter-wave plate; FS: field stop.
Fig. 2.
Fig. 2. Structure of one FD telescope of the Pierre Auger Observatory with some important features highlighted. Taken from Ref. [22].
Fig. 3.
Fig. 3. Geometry of the Aeolus laser beam being detected by one of the four fluorescence detectors of the Pierre Auger Observatory.
Fig. 4.
Fig. 4. Measured ground tracks for an altitude of 1400 m for three sample Aeolus overpasses in the years 2019, 2020, and 2021. Each FD site has a 180° FoV and is oriented to overlook the enclosed area. Consequently, the events in 2019 and 2020 were observed primarily by the northern and southern sites, as the eastern site is situated too far away. Additionally, a fraction of the 2019 events was visible to the western FD site, where the southernmost events appeared at the edge of its FoV.
Fig. 5.
Fig. 5. Comparison between the geolocation derived with different Aeolus data processor versions before (L1A v7.11.2, orange line) and after (L1A v7.12, purple line) the fix of the error and measurements by the Auger Observatory (blue dots) at 10 km altitude for the overpass on 17 July 2021. The offset amounts to 6.8 km between Auger and v7.11.2, and 0.8 km between Auger and v7.12.
Fig. 6.
Fig. 6. Reconstructed energies for three sample Aeolus overpasses in 2019, 2020, and 2021. The average energy per overpass is marked by the dashed line. The energy is given in mJ by the top axis and in PeV ( ${10^{15}}\,\,\rm eV$ ) by the bottom axis.
Fig. 7.
Fig. 7. Signal evolution of the Aeolus instrument: the purple curve represents the laser energy as measured at the output of the respective transmitter [flight models (FMs) A and B]. The red curve denotes the atmospheric return signal that is detected on the Rayleigh receiver under clear-air conditions. The black dots indicate the Auger measurements with error bars that describe only their statistical errors (see text).

Equations (1)

Equations on this page are rendered with MathJax. Learn more.

E ( 2019 ) = 33.1 0.8 + 1.9 m J , E ( 2020 ) = 23.7 0.6 + 1.7 m J , E ( 2021 ) = 17.3 0.4 + 0.9 m J .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.