the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
The EarthCARE mission – science and system overview
Tobias Wehr
Georgios Tzeremes
Kotska Wallace
Hirotaka Nakatsuka
Yuichi Ohno
Rob Koopman
Stephanie Rusli
Maki Kikuchi
Michael Eisinger
Toshiyuki Tanaka
Masatoshi Taga
Patrick Deghaye
Eichi Tomita
Dirk Bernaerts
The Earth Cloud, Aerosol and Radiation Explorer (EarthCARE) is a satellite mission implemented by the European Space Agency (ESA), in cooperation with the Japan Aerospace Exploration Agency (JAXA), to measure global profiles of aerosols, clouds and precipitation properties together with radiative fluxes and derived heating rates. The simultaneous measurements of the vertical structure and horizontal distribution of cloud and aerosol fields, together with outgoing radiation, will be used in particular to evaluate their representation in weather forecasting and climate models and to improve our understanding of cloud and aerosol radiative impact and feedback mechanisms. To achieve the objective, the goal is that a retrieved scene with footprint size of 10 km × 10 km is measured with sufficiently high resolution that the atmospheric vertical profile of short-wave (solar) and long-wave (thermal) flux can be reconstructed with an accuracy of 10 W m−2 at the top of the atmosphere.
To optimise the performance of the two active instruments, the platform will fly at a relatively low altitude of 393 km, with an equatorial revisit time of 25 d. The scientific payload consists of four instruments: an atmospheric lidar, a cloud-profiling radar with Doppler capability, a multi-spectral imager and a broadband radiometer. Co-located measurements from these instruments are processed in the ground segment, which produces and distributes a wide range of science data products. As well as the Level 1 (L1) product of each instrument, a large number of multiple-instrument L2 products have been developed, in both Europe and Japan, benefiting from the data synergy. An end-to-end simulator and several test scenes have been developed that simulate EarthCARE observations and provide a development and test environment for L1 and L2 processors.
Within this paper the EarthCARE observational requirements are addressed. An overview is given of the space segment with a detailed description of the four science instruments, demonstrating how the observational requirements will be met. Furthermore, the elements of the space segment and ground segment that are relevant for science data users are described and the data products are introduced.
- Article
(10255 KB) - Full-text XML
- BibTeX
- EndNote
The Earth Cloud, Aerosol and Radiation Explorer (EarthCARE) will provide global profiles of clouds, aerosols and precipitation along with co-located radiative flux measurements. Atmospheric microphysical properties and associated radiative fluxes will be used to evaluate the representation of aerosols, clouds and precipitation in weather forecast and climate models and help improve their parameterisation schemes.
The EarthCARE satellite includes four scientific instruments: an atmospheric lidar, a Doppler cloud radar, a multi-spectral imager and a broadband radiometer. The instruments have been characterised, calibrated, tested and integrated on the spacecraft. The lidar operates at UV wavelength and is equipped with a high-spectral-resolution receiver and depolarisation measurement channel. It will measure vertical profiles of aerosols and thin clouds and allow for classification of aerosol types. The highly sensitive 94 GHz (W-band) cloud radar is equipped with Doppler measurement capability and will provide measurements of clouds and precipitation with a sensitivity partly overlapping with the lidar but penetrating through or deep into the cloud, far beyond where the lidar signal attenuates. Furthermore, it will measure the vertical velocities of cloud particles. The imager has four solar and three thermal channels and will provide across-track swath information on clouds and aerosols. It will also be used to construct retrieved 3D cloud–aerosol–precipitation scenes for radiative transfer calculations. The Broad-Band Radiometer measures solar and thermal radiances in three fixed viewing directions along the flight track so that accurate top-of-atmosphere flux estimates can be derived. EarthCARE will fly in a sun-synchronic low-Earth orbit at a relatively low altitude in order to maximise the performance of the active instruments.
The EarthCARE-retrieved geophysical data products will include, for example, target classification, the vertical profile of microphysical properties of ice, mixed and liquid clouds, particle fall speed, precipitation parameters, and aerosol type. By applying 1D and 3D radiative transfer models to those retrieved data, heating rates and radiative fluxes will be calculated. These calculations will be evaluated against the flux estimates derived from EarthCARE Broad-Band Radiometer measurements. The scientific data products have been defined, and their retrieval algorithms have been developed. The performance of the retrieval algorithms has been evaluated. Other papers in this special issue will describe the details of these data products and their retrieval algorithms.
The EarthCARE Payload Data Ground Segment (PDGS) produces data products up to the geophysical retrievals co-located at the satellite instruments' fields of view (Level 2 data products). (Higher processing levels and applications, such as data product exploitation for climate/weather model improvements, are out of the scope of this paper.)
This paper presents the EarthCARE mission background, the satellite with its scientific payload and the ground segment. The paper is intended as background information for the primary objective of this special issue, namely, the description of the theoretical bases and validation of the various geophysical retrieval algorithms and data products developed in Europe, Japan and Canada for EarthCARE. A follow-up publication is under preparation that describes the overall data-processing chain and provides an overview of the data products (Eisinger et al., 2023).
Section 2 addresses the overarching scientific needs that justify the mission. The required observations are discussed in Sect. 3. A brief system overview is given in Sect. 4. The satellite instruments are described in Sect. 5. The orbit is described in Sect. 6. The mission phases are described in Sect. 7. The mission ground segment, including a summary of the geophysical data products, is presented in Sect. 8.
EarthCARE will acquire key observables needed to better understand the role of cloud, aerosol and radiation in climate. The mission was selected in 2000 for pre-feasibility study, in 2001 for Phase A study and in 2004 for implementation. Understanding of the role of cloud radiative feedback and the direct and indirect radiative feedback of aerosols, as well as their representation in weather in climate models, has significantly improved since that time. Yet, high-resolution and vertically resolved observations of clouds and aerosols, co-located with measurements of reflected solar and emitted thermal radiation from Earth, remain critically important parameters to further the understanding of cloud and aerosol radiative impact and their feedback mechanisms in a warming climate.
The observations of the EarthCARE mission will address these science questions directly. The mission will measure (1) vertically resolved profiles of clouds, aerosols and precipitation at the satellite nadir track; (2) 3D cloud–aerosol scenes extending across track; and (3) co-located solar and thermal broadband radiances and fluxes.
Climate feedbacks from clouds are significant and by far the most inconsistent across climate models. The IPCC Sixth Assessment Report (IPCC, 2021) assumes a positive net cloud feedback in response to global warming with high confidence. Major advances in the understanding of cloud processes, compared to the Fifth Assessment Report, have decreased the uncertainty range in the cloud feedback by about 50 % and assume a net cloud feedback of +0.42 (−0.10 to 0.94) W m−2 K−1. Still, clouds remain the largest contribution to overall uncertainty in climate feedbacks. The IPCC report also addresses aerosol–cloud interaction and points out that “CMIP6 models generally represent more processes that drive aerosol–cloud interactions than the previous generation of climate models, but there is only medium confidence that those enhancements improve their fitness-for-purpose of simulating radiative forcing of aerosol–cloud interactions.” Finally – related to the objectives of EarthCARE – the report acknowledges that CMIP6 models still have deficiencies in simulating precipitation patterns, particularly in the tropical ocean.
The largest uncertainty of general circulation models (GCMs) in the prediction of climate sensitivity to CO2 increase is related to cloud feedback due to the diversity of cloud formation processes and the feedback mechanisms depending on cloud macrophysical and microphysical properties such as type, thickness, phase, location and height (Sherwood et al., 2020). In particular, height and vertical extent of cloud over the globe, which are essential for cloud feedback evaluation, cannot be satisfactorily determined by passive space sensors, let alone ground-based and airborne observations, but require the use of space-borne active, range-resolving, lidar and radar systems. Regarding detection of cloud feedback to a warming climate, Chepfer et al. (2014) call for a long-term lidar record in order to detect changes in clouds; similarly Takahashi et al. (2019) pointed out that, using W-band radar, long data records extending into the 2030s are needed to reliably detect trends in cloud heights. Vaillant de Guélis et al. (2018) stressed that traditional satellite cloud records poorly constrain cloud feedback due to the lack of vertical cloud profile information and, therefore, cloud lidars are required to reduce climate sensitivity uncertainties. Substantial advances in the understanding of interactions of cloud and climate were achieved with CloudSat and CALIPSO but will have to be continued in order to allow trend detection (Winker et al., 2017).
Illingworth et al. (2015) provided an overview of the EarthCARE mission and its scientific goals. As addressed in the paper, the objective of the mission is to provide observations of global cloud and aerosol profiles together with radiation observations in order to provide data for the improvement of clouds and aerosol parameterisation in climate and weather models. Required observations include global vertical profiles of natural and anthropogenic aerosol, vertical profiles of ice and liquid cloud, precipitation, and derived vertical heating rates. EarthCARE's observed top-of-atmosphere radiances and fluxes will both contribute to Earth radiation budget studies and, more importantly, be used in verification of the calculated heating rate and flux profiles calculated from retrieved cloud–aerosol scenes. Beyond these original EarthCARE mission objectives related to climate and numerical weather prediction (NWP) model improvements, it is now envisaged by ECMWF to assimilate EarthCARE lidar and radar observations operationally into the Integrated Forecasting System in anticipation of a positive impact on the forecast accuracy (Janisková and Fielding, 2020).
EarthCARE will both extend the record of the former space-borne lidar and radar observations of CloudSat and CALIPSO and provide increased performance of such observations, by using a lidar with high spectral resolution and two polarisation channels and an accurately co-located W-band radar with Doppler capability – which will be the first in space – and with an increased sensitivity of about 5 dB compared to CloudSat. With the help of a push-broom imager, 3D cloud–aerosol scenes will be retrieved, and heating rate and radiative flux profiles will be calculated and compared to co-located top-of-atmosphere fluxes derived from the Broad-Band Radiometer measurements.
EarthCARE also builds on the success of the CloudSat and CALIPSO missions of the A-Train (Stephens et al., 2018). Launched in 2006, both satellites have exceeded their design lifetime by far and delivered unprecedented measurements of cloud and aerosol profiles and precipitation. In comparison to the pre-CloudSat/pre-CALIPSO era, when only passive space-borne cloud observations were available, the understanding of cloud net forcing has been significantly improved through provision of true measured cloud heights and optical thicknesses derived from measured ice and water paths and the identification of multi-layered cloud systems. Stephens et al. (2018) provide an overview of how CloudSat–CALIPSO measurements provide a new understanding of clouds and their processes and radiative effects. Winker et al. (2013) describe how CALIPSO observations were used to characterise the global 3D distributions of aerosol and their seasonal and inter-annual variations, providing unprecedented data sets for comparisons to aerosol models.
EarthCARE will extend and advance these observations with a radar with Doppler capability to give information on convective motions as well as ice precipitation and rainfall speeds, leading to improved drizzle, rainfall and snowfall rates. The additional sensitivity compared with CloudSat will enable it to better detect thin ice clouds and much more low-level stratus and stratocumulus. The lidar will operate in the UV and is equipped with a high-spectral-resolution receiver to discriminate between particulate and molecular backscatter, thus providing direct measurements of the extinction profile of clouds and aerosol plus a polarisation channel for advanced classification of aerosol type and ice particle characteristics.
Earth Radiation Budget missions are also relevant heritage for EarthCARE so that measured cloud and aerosol profiles can be utilised in the context of co-located radiative flux measurements. Namely, CERES (Wielicki et al., 1996) observations have been used to evaluate radiative transfer properties derived from CloudSat and CALIPSO cloud and aerosol retrievals (Kato et al., 2011). EarthCARE will fly its own broadband radiometer to derive top-of-atmosphere radiative fluxes (Velázquez-Blázquez et al., 2023b), allowing for similar studies.
The EarthCARE Mission Requirements Document (MRD; Wehr, 2006) was written after the completion of the feasibility study (Mission Phase A). The MRD defines the scientific mission objectives and observational requirements as a reference for the system (including satellite, instruments, ground segment, data products) design and implementation.
The MRD prescribes that EarthCARE shall meet these mission objectives by simultaneously measuring the vertical structure and horizontal distribution of cloud and aerosol fields together with outgoing radiation over all climate zones. It shall measure properties of aerosol layers (occurrences, profile of extinction, boundary layer height, aerosol type); clouds (3D cloud boundaries, occurrence of ice/liquid/super-cooled layers); vertical profiles of ice and liquid water content along with effective ice particles and shape or effective droplet size, respectively; vertical velocities (cloud convective motion and ice sedimentation); and drizzle and rainfall rates, co-located with top-of-atmosphere solar and thermal fluxes.
The key accuracy requirement provided in the MRD stipulates that a retrieved (cloud–aerosol–precipitation) scene with a footprint size of 10 km × 10 km shall be sufficiently accurate for its atmospheric vertical profile of short-wave (solar) and long-wave (thermal) flux to be reconstructed with an accuracy of 10 W m−2 at the top of the atmosphere.
Therefore, the geophysical data products retrieved by EarthCARE must provide, firstly, a comprehensive set of cloud and aerosol vertical properties for each observed scene. Secondly, these scenes must be extended in the across-track dimension in order to establish 3D scenes sufficiently wide for the application of – depending on the scene complexity – 1D or 3D Monte Carlo radiative transfer calculations providing vertical profiles of heating rates and (broadband) solar and thermal fluxes, co-located with the retrieved cloud–precipitation and aerosol profiles. The so-derived fluxes (and radiances) at the top of the atmosphere are, thirdly, compared to the measurements of the EarthCARE Broad-Band Radiometer for verification of the cloud–aerosol retrievals and consistency of the radiative transfer calculations. The required spatial resolution for cloud–precipitation–aerosol profiles along track is typically 1 to several kilometres and the vertical resolution 100 to a few hundred metres, depending on the target. The radiative properties are evaluated for scenes of 100 km2 size.
Table 1 lists the MRD observational requirements for clouds and aerosols. The requirements are derived from the 10 W m−2 requirement, with the exception of the vertical velocity requirement. For this, the MRD requires the measurement of vertical motion within clouds, in particular for characterising convection, estimating sedimentation velocity of ice particles in cirrus, quantifying drizzle fluxes in stratocumulus and estimating heavier rain rates. The MRD requires retrieved cloud products on a horizontal grid of 10 km, while individual measurements and retrievals should be made on a 1 km grid.
Global sampling is required at a fixed orbit local time in order to avoid aliasing of diurnal and seasonal variations. The choice of the orbit equatorial crossing time is a compromise depending on both instrument performances and diurnal variations in atmospheric properties and radiation fields. Afternoon orbits are generally preferred, as convection (over land) is mainly initiated in the early afternoon. Furthermore, near-noon retrievals of top-of-atmosphere broadband short-wave fluxes are most representative of the diurnal average, and relative errors in short-wave radiance measurements and flux retrievals are smallest when the short-wave signal is largest. In order to minimise the effect of sunglint, the orbit must not be too close to noon. The MRD requires an Equator crossing time between 13:15 and 14:00 MLST.
In order to fully exploit the synergistic value of the observations of the EarthCARE instruments, their footprints must be co-located. The MRD requires that three instruments used for cloud, precipitation and aerosol retrievals are co-located within 350 m (goal 200 m) and with the radiation measurements within 1 km.
The data latency for Level 1 data products (Sect. 8.2) has to be compatible with monitoring and data assimilation into NWP models. When the MRD was written in 2006, advances in NWP modelling were anticipated that would allow for the assimilation of EarthCARE data, which has been confirmed in the mean time (Janisková and Fielding, 2020). For climatological studies and evaluation of climate and NWP models, delayed data latencies of several days to a month might be adequate; however, a fast data availability for Level 1 and Level 2 data would still be very important for fast detection of any possible data degradation so that corrective actions can be taken.
The EarthCARE mission has been developed in full compliance with its observational requirements. The satellite and instrument performances have been verified and confirmed. A sophisticated data-processing system has been developed, and the comprehensive data products and their retrieval algorithms – presented in this special issue of this journal – are fully compliant with the EarthCARE observational requirements.
The EarthCARE satellite is a three-axis, stabilised, custom-built, carbon-fibre-reinforced polymer (CFRP) platform, optimised for minimum mass and high stability. The streamlined shape of the spacecraft, with its trailing solar panel, minimises its cross section in order to reduce residual atmospheric drag at its relatively low flight altitude of 393 km (nominal mean spherical altitude). The low altitude has been selected in order to optimise the performance of the two active instruments, the lidar and the radar. The orbit is sun-synchronous, with a descending-node crossing time of 14:00 MLST, an inclination of 97∘ and revisit time of 25 d (389 orbits). The deployable solar array delivers a power of 1710 W (average at end of life) for a nominal average power requirement of 1670 W and use of lithium-ion batteries with a capacity of 324 A h. The total mass of the satellite is 2350 kg, including 313 kg of hydrazine on-board propellant, sufficient for 3 years of operation plus 1 year reserve.
The satellite includes four scientific instruments. The ATmospheric LIDar (ATLID) is a cloud–aerosol high-spectral-resolution lidar (HSRL) operating at 355 nm wavelength. It is equipped with a receiver that separates molecular from particulate backscatter and depolarised total backscatter. It will provide measurements of the vertical profiles of aerosol and (thin) clouds. The Cloud Profiling Radar (CPR) is a 94 GHz (W-band) cloud radar with Doppler capability that will provide cloud profiles, rain estimates and particle vertical velocity. The Multi-Spectral Imager (MSI), with four solar and three thermal channels, provides cloud and aerosol observations extended in the across-track direction. The MSI is essential for creating 3D cloud–aerosol scenes out of the lidar- and radar-derived vertical profiles in order to calculate radiative properties of these scenes, including heating rates and top-of-atmosphere radiances and fluxes, which will be compared to those measured by the fourth EarthCARE instrument, the Broad-Band Radiometer (BBR). The satellite, including the instruments ATLID, MSI and BBR, is developed by ESA, with Airbus Defence and Space GmbH, Germany, as the prime contractor, while the CPR is developed by JAXA and the National Institute of Information and Communications Technology (NICT) in Japan.
The satellite command and control uses S-band up- or downlinks at the ground station in Kiruna (Sweden). An X-band system is used for science data communication using both Kiruna and Inuvik (Canada) ground stations. The science data packages are sent to the Payload Data Ground Segment (PDGS) at the ESA Centre for Earth Observation (ESRIN), Frascati, for the production of Level 0 data. Subsequent Level 1 and Level 2 operational processing is performed at JAXA for CPR Level 1 and JAXA Level 2 data products and at ESA for ATLID, MSI and BBR Level 1 and ESA Level 2 data products. Both ESA and JAXA Level 1 and Level 2 data products will be made available to the global user community.
The sophisticated data-processing system is described by Eisinger et al. (2023). The geophysical data products and their retrieval algorithms are described in the papers of this special issue.
The satellite science payload consists of four instruments: the ATLID, the CPR, the MSI and the BBR. Figure 1 shows an artist's impression of the satellite and indicates the location of the four instruments. Figure 2 shows a schematic sketch of the viewing geometry of the instruments. The satellite instruments have been previously described, in particular by Wallace et al. (2014, 2016) and Heliere et al. (2017). The CPR has been described previously by Nakatsuka et al. (2012). However, meanwhile, some details of the instruments have changed and their calibration and performance verifications have been completed.
5.1 The ATmospheric LIDar – ATLID
The ATLID is a high-spectral-resolution (HSR) atmospheric backscatter light detection and ranging (lidar) instrument (HSRL) for the detection of cloud boundaries and profiling of optically thin clouds and aerosols. The instrument has been developed by Airbus Defence and Space (France).
The ATLID is used for the retrieval of a lidar feature mask (Nishizawa et al., 2023; van Zadelhoff et al., 2023a) and, subsequently, cloud and aerosol profile retrievals (Donovan et al., 2022; Nishizawa et al., 2023; Okamoto et al., 2023; Wandinger et al., 2023) and, in synergy with MSI data, to derive cloud top height/column properties (Haarig et al., 2023) and aerosol properties (Kudo et al., 2016). The ATLID is furthermore used in synergy with the CPR and MSI to derive synergistic target classification and atmospheric profile products (Irbah et al., 2023; Mason et al., 2022; Okamoto et al., 2023).
The ATLID is designed as a self-standing instrument, reducing the mechanical coupling of instrument–platform interfaces and allowing for better flexibility in the satellite integration sequence. The transceiver architecture selected for this instrument is bistatic, meaning that the two redundant emission chains are separated from the receiver chain. Consequently, there is a single optical element exposed to vacuum that is illuminated by the laser. This approach was selected in order to minimise the risk of laser-induced contamination (LIC) of optical surfaces. Figure 3 shows a schematic sketch of the bistatic design, including the two laser sources. The ATLID weighs 555 kg, consumes 538 W during operation and generates 753 kB s−1 of data during operation.
In order to ensure the longevity of the mission (minimum 2.5 years of operations after commissioning), there are two completely independent, cold redundant emission chains. Each emission chain is manufactured by Leonardo (Pomezia, Italy) and henceforth referred to as a transmitter assembly (TxA). Each TxA is comprised of a reference laser head (RLH), developed by TESAT (Backnang, Germany), seeding via optical fibre a power laser head (PLH), plus the associated transmitter laser electronics (TLE). A major evolution with respect to the ALADIN (Cosentino et al., 2012) laser of the Aeolus mission (Straume et al., 2020), the ATLID PLH is a sealed and pressurised enclosure at 1.2 atm dry air (pressurisation verified for over 10 years of space operations). This approach, in addition to careful selection and screening of included materials as well as extensive bake-out of all optical assemblies, ensures a mitigation of LIC (Wernham et al., 2010). A 2-year dedicated test campaign (accelerated test extrapolated to lifetime duration via contaminant heating and higher pulse repetition frequency) demonstrated the advantages of this approach. The PLH design has a master oscillator power amplifier (MOPA) laser architecture. The seeded master oscillator (MO) is an end-pumped Nd:YAG rod, pumped at 808 nm by collimated high-power laser diode stacks. The MO is equipped with 1+1 cold redundant diodes stacks to compensate for diode ageing losses and is isostatically mounted on a dedicated thick optical plate that is totally decoupled from all thermal sources. The MO is significantly derated (operates at one-third of its peak power) in order to provide a margin for possible energy recovery actions, as well as to extend the lifetime of its optical coatings. Special care was taken to ensure the entire laser is operating at low optical fluence, as coating longevity is significantly affected by high fluence. The pulse generated by the MO is amplified by a double-pass pump amplifier before passing to a frequency-tripling conversion crystal that is set to a wavelength of 355 nm. Another change to the design, with respect to the ALADIN laser, is the UV section. The beam coming out of the third harmonic generator is expanded immediately, to ensure very low fluence in the UV section, and the number of optical elements in UV are kept to a minimum (and partly enclosed with an open cover to minimise the field of view of possible contaminants). Finally, the laser unit is equipped with five photodetectors that allow precise knowledge of the emitted power (measurement accuracy 2 %) to permit control of the emitted frequency of the laser. The output from this laser consists of linearly polarised laser pulses, less than 35 ns pulse length, at a pulse repetition rate of 51 Hz and with a 31-to-35 mJ pulse energy. The output beam size is approximately 8×9 mm2 at the laser exit, with a divergence of 36 µrad after the × 6.7 beam expander. It points 3∘ backwards of nadir in order to minimise specular reflection on ice clouds. Figure 4 plots the UV energy during the PLH flight spare lifetime test performed by Leonardo.
Figure 5 shows schematically the optical path of the ATLID receiver. Light generated by the TxA and backscattered in the atmosphere is collected by an afocal Cassegrain telescope with a diameter of 62 cm. The ATLID receiver separately measures the co-polar and cross-polar backscattered signals. The co-polar signal is further separated, using a high-spectral-resolution etalon (HSRE) developed by Thales Alenia Space (Switzerland), into a spectrally broadened molecular backscatter (“Rayleigh”) signal channel and a particulate backscatter channel (usually physically inaccurately referred to as “Mie”), according to the principle of HSRL. Therefore, the ATLID has three receiver channels, referred to as the cross-polar channel, and the Rayleigh and Mie channels (which are both co-polar).
Figures 6 and 7 show the Fabry–Pérot etalon design and its Mie channel HSR module spectral proportion. The separated signals from the three channels are coupled into fibres and routed to three identical and separate memory charge-coupled devices (MCCDs). The MCCD, developed by Teledyne e2v (UK), is able to measure single-photon events to meet the worst-case radiometric performance requirements. The selected design provides a high response, together with extremely low noise thanks to on-chip storage of the echo samples, which allows delayed read-out at very low pixel frequency (typically below 50 kHz). Combined with an innovative read-out stage and sampling technique, the detection chain provides an extremely low read-out noise (<2 e− rms per sample). With respect to the hot pixel issues observed in the Aeolus mission (Weiler et al., 2021), ATLID pixel level CCD behaviour has been characterised under flight-like operational conditions and the dark current values do not present the same level of defect, with no hot pixels observed; this is due to a different CCD control sequence, which reduces the time that the charge resides in the device by a factor of 10. Each detector achieves a vertical resolution of 103 m in the vertical range from the surface up to 20 km altitude and of about 500 m in the vertical range of 20 km up to 40 km. The along-track sampling distance is 140 m; however, in order to increase the signal-to-noise performance, it is envisaged to integrate two consecutive lidar pixels, leading to an effective along-track spatial resolution of about 280 m.
In order to align the emission beam to the telescope field of view (FoV), the instrument implements a co-alignment control loop. Each PLH contains a beam-steering mechanism (BSM) that is controlled by the co-alignment sensor (CAS) in the receiver chain (a beam splitter divides off approximately 7 % of the incoming signal towards the co-alignment detector prior to the polarisation separation). This active alignment approach allows us to correct for depressurisation, structural as well as thermo-elastic deformations that could be expected after launch or during orbit. The MCCD sensor used in the CAS is identical to the three other ATLID channels but operates uncooled; subsequently its performance is degraded.
The ATLID instrument at launch is fully compliant with the performance requirements for the current orbit selection. The instrument was designed with significant performance margins in order to remain compliant even in a worst-case end-of-mission scenario, including the design case high orbit of up to 425 km. Table 2 summarises the most important performance requirements, as well as the measured beginning of life and extrapolated worst case under end-of-life conditions at the nominal operational orbit. The variations in the two lasers are caused by the accuracy of measuring the emitted power. In this aspect the ATLID project benefited significantly from the lessons learned from ALADIN, with improved monitoring setup and the optics after the UV monitoring photodiode kept to a minimum. Figures 8 and 9 provide a more detailed plot of the channel retrieval accuracies versus latitude during summer as well as winter conditions.
In order to achieve the above-mentioned performance conditions, a number of instrument calibrations take place systematically on orbit. Some of them are adapted with every pulse, while other calibrations require the instrument to switch into a specific mode (dark current calibration requires the laser to be switched off). The accuracy of each derived input signal is not a product of one calibration but a result of several in-flight calibrations as well as an a priori good knowledge of the instrument characteristics and capabilities. The calibration requirements for absolute calibration accuracies are derived from the periodic measurement of total noise in darkness, the pulse-to-pulse background subtraction accuracy and the systematic calculation of lidar constants on both channels, as well as the spectral cross-talk knowledge accuracy. Furthermore, each type of calibration performed requires a specific observation scene as well as certain conditions to be met (e.g. preference of night measurement, avoidance of South Atlantic Anomaly, specific cloud formations, avoidance of ocean overflights). Figure 3 summarises the most critical calibration requirements as well as the spectral and polarisation cross-talk channel knowledge.
The pre-launch testing and calibration of the ATLID has been described by do Carmo et al. (2021).
The ATLID calibration algorithms (Level 1 processing) and Level 1b data products are described by Eisinger et al. (2023). The Level 1 data product consists of geolocated attenuated particulate (Mie) backscatter, attenuated molecular (Rayleigh) backscatter and attenuated cross-polar backscatter.
5.2 The Cloud Profiling Radar – CPR
The EarthCARE CPR is a 94 GHz pulse radar which measures altitude and Doppler velocity from the received echo. It is the world’s first space-borne cloud radar that estimates vertical profiles of the cloud and its vertical motion. Its minimum sensitivity is better than −35 dBZ with a Doppler velocity measurement accuracy of better than 1.3 m s−1 for specification and also 1.0 m s−1 for target performance. To achieve those performance criteria, the CPR is equipped with a large reflector with a diameter of 2.5 m. The instrument vertical range of observation will be selected via three modes, which select between altitude range window settings of 20, 18 or 16 km. The selection is made depending on the latitude. The lowest measurement altitude also extends to −1 km in order to permit the use of surface backscatter. The instrument has been jointly developed by JAXA and NICT. It is a 270 kg instrument with 316 W power consumption and generates a science data rate of up to 270 kB s−1. With the 2.5 m diameter antenna deployed, its dimensions are m.
The CPR Level 1 product contains radar reflectivity and Doppler velocity detected from cloud and precipitation particles. Those products are used to derive cloud mask and cloud and precipitation properties (Okamoto et al., 2023; Kollias et al., 2022; Mroz et al., 2023). In synergistic use with the ATLID and MSI, those data are also used to classify atmospheric targets (Irbah et al., 2023) and to derive those profiles that are important for radiative transfer calculation (Mason et al., 2022).
The external view of the EarthCARE CPR and major components' layout inside the platform are shown in Fig. 10. Also a simple block diagram of the CPR is shown in Fig. 11. The CPR consists of a TRS (transmitter receiver subsystem) and SPU (signal processor unit), QOF (quasi-optical feeder) with sub-reflector and MREF (main reflector). The TRS mainly consists of an LPT (low-power transmitter), HPT (high-power transmitter) and RCV (receiver): the 94 GHz pulse signal is generated by the LPT, the signal is amplified by the HPT to more than 1.5 kW using an EIK (extended-interaction klystron) and the signal is transmitted towards the QOF. The RCV receives echoes from clouds and ground through the QOF and down-converts these to an intermediate frequency (IF), which is 60 MHz, and the IF signal is transferred to the SPU. The TRS is a fully redundant system; in particular, the HPT is linked with both nominal and redundant sides via a magic T, which is a 3 dB hybrid waveguide coupler used for power division to both sides of the HPT. The leak signal from the QOF during the transmitting of the signal pulses is used as reference for the Doppler capability. The noise source (hot and ambient) within the RCV, which can be switched from the antenna (observation) receiving channel, is used to calibrate the receiver system. The QOF has two functions: primary feeds to the MREF and separation of the Tx (transmitting) and Rx (receiving) signals. Input signal to the QOF, which is linearly polarised, is transformed to circular polarisation by means of an FSP (frequency-selective polariser) system in the QOF. A LHCP (left-hand circular polarisation) signal is transmitted from the CPR to the Earth, and a RHCP (right-hand circular polarisation) signal is returned to the CPR. The diameter of the MREF aperture is 2.5 m. This size was derived from the restriction of the fairing size of the assumed rocket launcher and the need to maximise the antenna gain and therefore decrease the estimated error in Doppler velocity retrieval. Manufacturing a large antenna reflector of 2.5 m was challenging for the frequency band, and thus an accurate and light-weighted antenna was demanded for space-borne millimetre-wave radars.
The frequency of radar signal is 94 GHz, which is the same as the CloudSat CPR (Stephens et al., 2008), but the sensitivity of the radar is much better due to its lower orbit and larger antenna aperture size. The minimum radar reflectivity of the EarthCARE CPR is −35 dBZ defined at the top of the atmosphere (20 km) under the condition of 10 km horizontal integration. The observation range of the Doppler velocity is ±10 m s−1, with an accuracy of at least 1.3 m s−1 for specification and also 1.0 m s−1 for target performance in the ideal-case beam-pointing error (bias and harmonic elements) removal achieved in orbit, for cloud echoes of more than −19 dBZ when they are integrated over a 10 km horizontal distance. The transmit pulse width is 3.3 µs corresponding to the vertical resolution of 500 m, which is also the same as the CloudSat CPR. However, the received echo is over-sampled at 100 m, which is more precise than that of CloudSat (250 m). The vertical observation range can be selected from three different range window settings of 20, 18 or 16 km considering the satellite latitude, and the lowest measurement altitude extends to −1 km (below the ground surface), allowing the utilisation of surface backscatter. Minimising the range window would be needed to achieve a higher PRF (pulse repetition frequency) for more accurate Doppler measurement. A variable PRF function has been implemented for altering the PRF with satellite altitude and latitude. The instantaneous footprint size of the CPR is about 750 m at ground level under the assumption of a maximum geodetic altitude of 450 km. Horizontally averaged data along track are produced with the minimum pixel size of 500 m for a fine signal-to-noise ratio (SNR) and reducing data volume. The CPR obtains several types of calibration data to ensure and maintain performance. Specifically, the CPR is designed to acquire transmit power by using the QOF's power monitor capability. Receiver performance (e.g. noise figure – NF) is determined with the noise source’s signal in the RCV. During internal calibration mode, linearity and bias of the logarithmic amplifier in the SPU are acquired utilising the internal IF signal source, a step attenuator and the data acquired by a terminated logarithmic amplifier.
The overall performance of the CPR will be evaluated using an ARC (active radar calibrator) on the ground for which the CPR is configured to external calibration mode. The ARC is used to confirm the CPR transmitter performance (ARC in receiver mode) and receiver performance (ARC in transmitter mode) as well as overall performance (ARC in transponder mode with returning signal delay to avoid the contamination with ground echo). In addition, radar signal acquired by the ARC during the ARC fly-over will provide the antenna pattern and information on pulse length.
For the sea surface calibration (sea surface calibration mode), the EarthCARE satellite will perform a roll manoeuvre at regular intervals (e.g. once a month). Sigma 0 for various incident angles are useful to evaluate the CPR performance, since the normalised radar cross section of the sea surface echo (namely, sigma 0) has clear incident angle dependency and its shape is dependent on sea surface wind. In addition, other natural targets will be considered for the CPR Doppler performance evaluation.
In CPR Level 1 processing, various CPR calibrations (explained above) are conducted to retrieve calibrated radar reflectivity and Doppler velocity, which are used for higher-level product processing. In order to derive radar reflectivity, received echo power is calibrated with calibration load in the CPR receiver. Then, the calibrated received echo power is transformed to the radar reflectivity with use of noise level estimation and CPR transmit power estimation. For the Doppler velocity product, pulse-pair covariance is calibrated with transmit radio frequency (RF) phase change between inter-pulses (almost negligible) and the satellite velocity contamination in the CPR beam direction. Then, Doppler velocity is derived from calibrated pulse-pair covariance. CPR Level 1 processing also includes surface echo bin detection, estimation of the normalised radar cross section of the surface and calculation of the Doppler velocity of the surface. Those products are useful to validate geometric accuracy, radiometric accuracy and Doppler product accuracy.
The CPR Level 1b product contains the radar reflectivity and vertical Doppler velocity (Eisinger et al., 2023).
5.3 The Multi-Spectral Imager – MSI
The MSI will provide scene context information left and right of the satellite nadir track where the active instruments, the ATLID and CPR, are measuring vertical cloud, aerosol and precipitation profiles. MSI observations are used to derive “MSI-only” cloud and aerosol data products across the swath (Docter et al., 2023; Hünerbein et al., 2023, 2022; Nakajima et al., 2019). The MSI nadir pixel is used in the retrieval of synergistic products with the ATLID (Haarig et al., 2023) or the ATLID and CPR (Mason et al., 2022; Okamoto et al., 2023). Across-track data of selected MSI channels are furthermore used for the creation of 3D scenes (Qu et al., 2023), which are needed for 1D and 3D radiative transfer calculations to calculate radiative properties (Cole et al., 2022; Oikawa et al., 2018; Okata et al., 2017; Yamauchi et al., 2023) and perform closure assessment (Barker et al., 2023). The MSI is also used to improve the unfiltering of BBR-measured radiances (Velázquez-Blázquez et al., 2023a).
The MSI is a nadir-pointing push-broom imager with a swath width of 150 km. The swath is shifted sideways and extends 115 km to the left and 35 km to the right side of nadir in order to minimise the number of pixels affected by sunglint. (The satellite flies southwards on the day side of the orbit and crosses the Equator in the early afternoon.) The pixel sampling size is 500 m. The instrument has two cameras, a VNS (visible (VIS), near-infrared (NIR), short-wave infrared (SWIR)) camera with four “solar channels” and a TIR (thermal infrared) camera with three “thermal infrared channels”.
Figure 12 shows an illustration of the MSI optical bench. The two apertures of the VNS camera and the TIR aperture are pointing towards the Earth. The TIR cold-space view is used for calibration of the TIR, and the sun view is used for the VNS calibration. The front-end electronics (FEE) are shared between the two cameras. (They are not visible in the illustration as they are located behind the TIR camera. Also, the instrument control unit (ICU) is not shown as it is enclosed within the satellite.) The optical layout of the VNS camera is shown in Fig. 13 and the TIR in Fig. 14. Table 4 describes the MSI channel definition and the radiometric performance.
The instrument has a power consumption of 60 W and a total mass of 50 kg, which is 40 kg for the optical bench module and 10 kg for the ICU. The maximum data rate is 611 kB s−1. SSTL (UK) is the instrument prime contractor and manufacturer of the TIR camera and FEE. The VNS camera is built by TNO (the Netherlands). The ICU is made by TAS (UK), the VNS detectors by Xenics (Belgium), the TIR microbolometer array by ULIS (France) and the blackbody by ABSL (UK).
The VNS camera. This camera images light from the scene onto four different focal planes (one per channel). Each is equipped with a linear photodiode array (silicon detectors) with 512 pixels of 25 µm. The target scene illuminates 360 pixels. Bands VIS, NIR and SWIR1 share a common 4.7 mm aperture and are separated by a system of dichroics and mirrors. The SWIR2 band has a larger, separate 10.4 mm aperture. Its corresponding focal plane is stabilised at 235 K via a passive cooling system. The field of view is 11.5∘ at for the VIS, NIR and SWIR1 channels and at for SWIR2. The optical components include glasses such as NLaSF31, ZnS and Ge in the SWIR channels. For the calibration, the camera view can be switched to a solar viewing port, where a pair of Spectralon diffusers are illuminated by the sun. For the protection of the camera and for the acquisition of the dark signal, the aperture of the instrument can be closed.
A standard InGaAs detector is used for the SWIR1 channel. The SWIR2 channel uses an extended-wavelength InGaAs detector, which offers better performance at high-temperature operations than MCT (mercury cadmium telluride). The SWIR2 detector is cooled to 235 K via a dedicated radiator. Each diode array includes two read-out circuits, which allows the use of odd–even redundancy. The use of the same read-out electronics for each detector array allows for a single design serving all four channels, which simplifies the FEE. Multiple detector exposures are accumulated within a single ground line for an improved SNR. The pixels at the edges of each array, which are not illuminated, are used for monitoring the detector thermal leakage current.
The VNS calibration unit consists of a carousel with two openings, one nadir pointing and one sun pointing. Inside the unit, a rotating disc can be positioned so that the nadir view is open for Earth observations (nominal operation) while the sun acquisition port is closed. Dark signal acquisition will be performed when both nadir and sun ports are closed. For solar calibration, the nadir is closed and the sunlight illuminates a pair of quasi-volume diffusers (QVDs), one for the VIS–NIR–SWIR1 aperture and one for the SWIR2 aperture. Two pairs of QVDs are installed that use the same optical chain to the detectors as is used in nominal operation such that any change in the optics chain transmission can be detected and corrected in the ground processor. The first QVD pair is used for a daily radiometric calibration of VNS channels against the known scene provided by the sun. The second QVD pair is used monthly in order to monitor ageing of the QVD pair that is used daily, which can be caused by contamination and solar UV exposure due to the pair's more frequent use. The solar calibration is performed over the south polar region while the satellite is flying over the shadowed Earth with the sun still within the line of sight. The solar aperture is baffled in order to avoid light scattered by the satellite or the solar panels.
The VNS calibration is based on
where Ωg(t) is the BSDF (bidirectional scattering distribution function) of the Earth scene at time t, SE(t) is the Earth signal at time t, SD(t) is the signal from the short-term calibration diffuser at time t, SD(t′) is the signal from the short-term calibration diffuser at time t′, SM(t′) is the signal from the long-term diffuser at time t′ and ΩM(t0) is the BSDF of the long-term diffuser characterised on the ground.
The measured VNS channel filter functions are shown in Fig. 15. The VNS channel central wavelengths have a small dependency on the viewing angle, caused by bandpass filters on curved optical surfaces (lenses), resulting in an angular swath dependency of the central wavelength as shown in Fig. 16. We refer to this effect as the spectral misalignment effect (SMILE) (although this is not identical to some previous uses of the term “smile” in the literature, e.g. as discussed by Fisher et al., 1998). This effect, as well as a slight dependence of the point spread function on the viewing angle, is considered in the MSI Level 2 retrievals.
The TIR camera. This camera uses a two-part imaging system, consisting of a telescope and a relay, to focus the three channels onto a single area-array detector. A rotating mirror (calibration mirror) allows the line of sight to be redirected for calibration purpose onto an internal, warm blackbody and to deep space. This calibration is done once per orbit.
Figure 14 illustrates the TIR optical system, where the calibration mirror reflects the incoming signal (from Earth or blackbody or cold space) onto lens L1 that forms the aperture, and the elements L1, M1 and L2 create an image of the scene at long focal length. The beam then enters the rectangular filter and dichroic optical components (FDOC) where it is split by two dichroic beam splitters (D1, D2) into three beams, which are reflected onto a common image plane and pass through their respective filters that define the response for each of the channels. An intermediate focus is formed within the FDOC so that the filters define a field stop. The system is close to tele-centric in the region of the dichroics and filters. The filters are tilted up to 10∘ in order to control stray reflections such that they all systematically reflect radiation from an internal reference inside the inner shroud, thereby providing control of the out-of-band radiation directed by the filters to the detectors. The detectors form an image of both the filters as well as the mounting frame, which acts as an internal reference body for the temperature of the rear optics structure so that it can be corrected for by subtraction from the scene signal.
The beams are re-imaged onto the microbolometer array with substantial de-magnification, with the final image at approximately with a ± 11.3∘ field of view. The rear optics are contained in a thermally isolated enclosure, temperature-controlled to better than 100 mK. The optics include zinc selenide, germanium and gallium arsenide elements. Athermalisation is achieved by structural means; an Invar metering rod holds the detector in a position fixed with respect to the optics lens barrel over a 50 K temperature range. The optic surface shapes include a number of aspherical components in order to meet the imaging requirements. The optical system operates near the diffraction limit. It uses time delay integration (TDI) to improve the signal-to-noise ratio: signals from a predetermined number of rows (N) are co-added from N successive frames, which effectively multiplies the total ground signal from each ground pixel by the TDI factor, N. MSI TIR uses a TDI factor of 19, applied on board, which permits an improvement in the SNR by a factor of almost 4.5.
The detector is a ULIS microbolometer area-array detector of 385 × 288 pixels at 35 µm pitch. The device is mounted in a flat package with an integral thermoelectric cooler (TEC) for thermal control.
For reference area subtraction, non-scene pixels are also imaged onto defined areas of the detector and are called reference blocks, which view the internal reference body (whose temperature follows that of the rear optics structure). Signal from these areas is used in the ground processor in a column-correction step to compensate for noise associated with the different read-out circuits in the detector and the correction of small thermal drifts over time.
The TIR is calibrated according to
where RE is radiance from the Earth scene, RB is the radiance from the internal calibration blackbody, SE is the signal recorded from Earth scene, SC is the signal recorded from cold-space observation and SB is the signal recorded from internal blackbody observation.
The TIR-measured filter functions are shown in Fig. 17.
The MSI Level 1b product contains the calibrated radiances for the solar channels (VNS channels) and the brightness temperatures for the thermal channels (TIR channels). The MSI Level 1c product contains the calibrated radiances and brightness temperatures remapped on a reference channel. The Level 1 products are described by Eisinger et al. (2023).
5.4 The Broad-Band Radiometer – BBR
The core objective of the EarthCARE mission is to study how aerosols, clouds and precipitation impact radiation. While the ATLID, CPR and MSI are used to derive profiles and even 3D aerosol–cloud–precipitation scenes, the BBR measures the co-located reflected solar and emitted thermal radiation. From the unfiltered BBR radiance measurements (Velázquez-Blázquez et al., 2023a), broadband fluxes at the top of the atmosphere are derived (Velázquez-Blázquez et al., 2023b). A key activity making use of these radiances and fluxes is the assessment of radiative closure, as described by Barker et al. (2023), where radiative transfer calculation (1D for every scene plus 3D Monte Carlo calculations for as many scenes as computationally affordable) is used on the aerosol–cloud–precipitation scenes retrieved from the ATLID, CPR and MSI to calculate heating rates and top-of-atmosphere fluxes and/or radiances, which are then compared to those measured by the BBR. The guiding performance requirement for the BBR is that top-of-atmosphere solar and thermal flux can be derived with an accuracy of 10 W m−2 for a scene size of 10 km × 10 km. Note that the performance requirements of the BBR are defined for a ground pixel size of 10 km × 10 km; however, the technical implementation of the instrument allows for a variable integrated pixel size.
For a completely homogeneous scene, the conversion of nadir radiance to flux is straightforward, but in the general case of inhomogeneous scenes, the scene should ideally be observed simultaneously from all directions. As this is not possible with a single satellite instrument, the BBR uses instead three fixed fields of view, one forward pointing along the ground track, one nadir and one backward pointing, so that each scene is observed from three directions. Velázquez-Blázquez et al. (2023b) describe how angular distribution models are used to estimate the flux from these observations.
The BBR development was led by Thales Alenia Space (UK), and the optics unit was built by STFC RAL Space (UK). The BBR has three fixed telescopes: pointed at nadir, in the forward direction at 50∘ relative to nadir and in backward direction at 50∘ relative to nadir (leading to a zenith angle at the respective fore and aft pixel of 54–55∘; see Fig. 2). Each view is mapped onto a microbolometer linear array detector. For each of the three telescopes, the instrument measures the total spectrally integrated radiation from 0.25 to >50 µm, alternating with the observation of the short-wave-only part (0.25 to 4 µm) by means of periodically switching a short-wave filter into the field of view. The individual pixel size results from the detector array resolution (across-track direction) and the total-wave and short-wave sampling speed; the instantaneous field of view is nominally 0.64 km × 0.64 km at nadir and 1.876 km (along track) × 1.069 km (across track) for the two slanted views. In the data processing on the ground, these pixels are integrated to (nominally) 10 km × 10 km pixels, as indicated in Fig. 22. The long-wave (LW) component is derived by subtracting the short-wave (SW) from the total-wave (TW) measurement.
The instrument requirements have all been defined with a nominal pixel size of 10 km × 10 km. However, the instrument's concept of high-resolution sampling of the total wave and short wave allows for flexibility in the integrated scene size. It is useful for some applications to integrate scenes of approximately the same area of 100 km2 – in order to maintain the signal-to-noise performance – but in a narrower strip around the nadir with accordingly longer along-track dimension, as is done by Qu et al. (2023), who construct 5 km × 21 km scenes for radiative transfer applications. (Note a 21 km scene length is used for this application, rather than 20 km, because of the 7 km along-track periods of the joint retrieval grid, defined in the Level 1d data product X-JSG; Eisinger et al., 2023.)
The BBR radiance product BM-RAD therefore includes the nominal pixel size of 10 km × 10 km as well as high-resolution data that the user can use to integrate to other desired pixel sizes.
It is important to advise users interested in the BBR radiances not to use the BBR Level 1 products (B-NOM, B-SNG) that have not been unfiltered because in these products the instrument effects have not yet been fully corrected for. Instead, the use of the BM-RAD product (Velázquez-Blázquez et al., 2023a), which is unfiltered, is advised. The so-called unfiltering (removal of the instrument effects) is part of the BM-RAD (Level 2) processing, where MSI scene data are (optionally) used for improved results.
BBR instrument. The TW channel is defined at 0.25 to >50 µm and the SW channel at <0.25 to 4.0 µm. The measured filter response is shown in Fig. 20. The dynamic range is 0 to 550 W m−2 sr−1 for TW and 0 to 450 W m−2 sr−1 for SW. Table 5 lists the radiometric performance against a set of reference scenes pre-defined for maximum and minimum SW and LW signals. The radiometric stability is 0.5 % for both channels.
The three views (fore, nadir and aft) are implemented using three identical fixed-mirror, dedicated telescopes imaging onto linear array detectors. In order to produce a uniform system point spread function (PSF), achieved by summing of pixels, and to minimise diffraction, each focus is deliberately de-focused to create a triangular pixel response. For the oblique (fore and aft) views less de-focus is applied, as otherwise spectral resolution would be decreased at the increased range of the targets. The detectors used are 30-pixel, linear microbolometer arrays from INO (Canada) (Proulx et al., 2017), onto which a special gold–black coating has been applied and laser-etched for the 0.1 mm × 0.1 mm pixels in order to increase their sensitivity to thermal infrared (TIR). The BBR's radiator is located approximately at the middle of the optics unit. The telescopes are mounted onto a baseplate that is parallel to that radiator.
The linear detector arrays correspond to an area on the ground of approximately 0.648 km (along track) and ± 10.2 km across track, for the nadir view, and approximately 1.88 km by ± 16 km for the slanted views. The ground sampling distance between one TW and the next TW sample (or between SW and SW, respectively) is 0.8 km. The reasons for the large width of the slanted-view ground track are, firstly, the use of identical telescopes and detector arrays for all three views and, secondly, allowing for compensation of a shifting field of view due to the Earth's curvature without having to apply any telescope steering by simply selecting the pixels relevant to the desired scene accordingly. Figure 21 shows the acquisition sequence of individual SW and TW measurements. Figure 22 indicates the resulting PSF of an integrated 10 km × 10 km pixel.
A chopper drum with four apertures – two of which house 2 mm thick, curved quartz filters – rotates continuously around the telescopes, cycling the signal onto the detectors through SW–drum skin–TW–drum skin, etc., with a chopping scheme based on GERB heritage (Harries et al., 2005). The drift dimension during detector integration is much smaller than the detector footprint. The nominal duration of each measurement view, 25 ms, is not long enough for the detector to reach a stable voltage, and the digitised waveform of the measured view is processed to predict the asymptote of the exponential rise in the detected signal. The chopper drum rotates continuously, and its speed is programmable. Around the chopper drum sits the calibration target drum, which houses hot and cold blackbodies, a diffuser, three fold mirrors to view the diffuser from each telescope, and the telescope baffles. The calibration target drum is rotated in steps as required. The chopper drum and calibration target drum mechanisms are nested and located on the opposite side to the telescope assembly baseplate. Both mechanisms operate on dry-lubricated bearings, the chopper bearing with a self-lubricating (solid composite PGM-HT) cage and the calibration bearing with a lead bronze cage. Separate mounted boxes are located on the opposite side of the radiator and house the signal conditioning electronics and the visible calibration (VisCal) system. The VisCal system houses a complex arrangement of mirrors; internal baffling; and an external illumination baffle that directs solar illumination onto the diffuser, which is mounted on the calibration target drum.
Figure 18 shows an exploded schematic of the instrument. Figure 19 shows photographs of some key components of the instrument.
The BBR will perform sampling such that the views from each of the three observation angles are spatially coincident, slightly separated in time by just over a minute between acquisitions. This requires the instrument to accommodate the residual shifts in scene position of the three fields of view that arise from Earth rotation between measurements after the implementation of yaw steering. For each telescope view a 10 km × 10 km scene point spread function (PSF) is synthesised by summation of individual pixel PSF responses along and across track. Pixels are selected to form scenes at 1 km intervals along track. To reduce the effect of overlapping radiances (from pixels outside of the nominal footprint) the PSFs of edge pixels are weighted accordingly. The match of the scenes, for the two channels and the three view directions, is controlled by the integrated energy (IE) requirements, Eq. (4). The IE is the ratio of the radiance measured by the instrument in a directional view i and a spectral channel j over a squared area of size d centred around a position (x0,y0) on Earth to the radiance measured by the instrument from the entire large and spatially uniform scene, according to
where i stands for the instrument along-track directional view (nadir, forward, backward), j stands for the spectral channel (SW, TW), Rij(λ) is the instrument spectral response function in directional view i and spectral channel j, and PSFij is the instrument point spread function in directional view i and spectral channel j. The IE of a ground scene in the SW and TW, respectively, over a square target of size ΔX, centred around the system PSF barycentre of the nadir view, considering the total radiance of the reference scene, is required to be identical within 5 % (goal 2 %) for the three directional views, according to the integrated energy requirement of the nominal 10 km × 10 km pixel:
where k stands for the forward and the backward view.
BBR calibration. BBR calibration relies on the use of blackbodies. For LW calibration, gain and offsets in the measured voltage signals are removed using hot and cold blackbodies that are autonomously rotated into the field of view every 88 s. This calibration frequency is decided by the temperature stability of the chopper drum skin and can be modified, if necessary, on orbit. The four instrument blackbodies are operated in hot redundancy and have been characterised on the ground against a blackbody that is traceable to radiometric standards. Every 6 months on orbit, the blackbody temperatures are swapped, allowing for a check on instrument linearity by recording data over a range of radiance.
For SW calibration, LW gain is transferable to the SW, with knowledge of the filter spectral response that has been characterised on the ground. SW gain is monitored in flight via measurements made with the VisCal system to detect changes in instrument response due to, for example, ageing. The wall of each telescope baffle incorporates a set of three monitor photodiodes in order to monitor in turn ageing of the VisCal optical path. The solar calibration is performed over approximately 30 orbits every 2 months.
The PSFs of the pixels in the microbolometer arrays are also characterised on the ground in order to formulate the algorithm for summing pixels to form the 10 km synthesised scene-PSFs.
The BBR Level 1b product (Eisinger et al., 2023) contains in particular calibrated irradiances in various spatial resolutions. Instrument effects have not been removed from the Level 1b product. The so-called unfiltering process takes place as part of the Level 2 processing; therefore, the science data user should use the BM-RAD Level 2 data product.
The orbit of the EarthCARE satellite is sun-synchronous, with a descending-node crossing time of 14:00 MLST ± 5 min at the Equator, an inclination of ≈ 97∘ and revisit time of 25 d (389 orbits). Table 7 shows the key parameters of the EarthCARE orbit.
Immediately after launch, a record of the orbits completed by the satellite is maintained through the use of an orbit number, starting with the first orbit and increasing incrementally. Orbit numbers are used for various data headers and file naming (e.g. for science data) to provide an unambiguous reference to the orbit within which the data were collected. The start and end of a particular orbit, and therefore the increment of the orbit count, are taken to be the ascending-node crossing (ANX) (i.e. one complete EarthCARE orbit starts from one ANX and ends at the next ANX). Orbit number 1 for the EarthCARE mission begins at the first ANX after launcher separation.
The orbit parameters according to Table 7 are actively maintained, in particular, to compensate for atmospheric drag. On-board thrusters are used for orbit maintenance manoeuvres, so-called “Δv” manoeuvres. (This name hints at the use of the satellite on-board thrusters which will be used to achieve a difference, Δ, in the satellite velocity, v.)
The baseline Δv manoeuvre (so-called 0∘ in-plane manoeuvre, meaning the spacecraft is accelerated along its flight vector in the orbital plane) is used for ground-track and attitude maintenance through semi-major axis and eccentricity control. Depending on the level of solar activity, it must be carried out (around the Equator) every 10 d (high solar activity) or, in the best case of low solar activity, every 30 d. The full operation takes 2000 s outside nominal operation while the instruments are switched down for safety, with a maximum period of 600 s operations of the thrusters. As an alternative, small Δv manoeuvres might be used, which are of much shorter maximum duration, allowing the spacecraft to stay in nominal pointing mode, and do not require the instruments to be switched off. In this case, the attitude disturbance lasts only a few seconds while the instrument science data acquisition continues. This operation lasts in total about 500 s, with a thruster burn duration of approximately 3 s, leading to a total out-of-specification pointing period of about 5 s, which is flagged in the data products accordingly. The small Δv manoeuvres would be grouped together on consecutive orbits during approximately 1 week with many small Δv manoeuvres (10 to 90 pairs per week). Both baseline and small Δv manoeuvres will be tested during the commissioning phase (see Sect. 7) in order to identify the preferred baseline for routine orbit maintenance operations.
There are also two other manoeuvres: (1) the first is a 90∘, out-of-plane manoeuvre (meaning the spacecraft is accelerated perpendicularly to its flight vector in the orbital plane) which is used to control inclination and mean local solar time. It can only be performed in ascending node (eclipse), within ± 1 month of equinox. There are two strategies proposed, either a single, large manoeuvre at mission midpoint or smaller manoeuvres at every equinox which use more fuel but slightly increase the accuracy to which the mean local solar time is maintained as well as reduce the size of each individual manoeuvre. (2) The second is a 180∘ out-of-plane retrograde manoeuvre, in order to lower the orbit, meaning the satellite velocity is reduced while keeping the flight vector in the orbital plane. As EarthCARE has no forward-pointing thrusters, the satellite will be rotated by 180∘ so that it is flying backwards for the duration of the manoeuvre. The only foreseen possible use of this manoeuvre is for a collision avoidance emergency.
The 0∘ in-plane Δv manoeuvres will maintain the satellite orbit so that its ground track will stay within the so-called dead band, which is the required accuracy of the predictable on-ground orbit track and has been set for EarthCARE to ± 25 km with respect to the nominal orbit nadir track. This means that a predicted overpass over a reference ground location has an uncertainty of ± 25 km.
The satellite altitude will be maintained within ± 2 km (worst-case altitude variation is ± 2.5 km). The exact orbit parameters, however, are always available in the Level 1 and Level 2 data products.
The ESA Flight Operations Segment (FOS; see Sect. 8.1) generates daily a new predicted-orbit file that covers the next 7 d, listing orbit state vectors. FOS also produces daily a reconstituted orbit file that contains the last 24 h of actual orbit state vectors. Operators of ground-based observatories or research aircraft can use these data for their planning of measurement activities coinciding with an EarthCARE overflight, for example, for the purpose of collecting ground-based data complementary to the EarthCARE satellite observations or for EarthCARE product validation.
The mission will be broken down into the following operational phases:
-
Launch and early orbit phase (LEOP). This phase includes launch, S-band communication acquisition, orbit acquisition, and deployment of the solar panel and CPR antenna.
-
Commissioning phase (also referred to as Phase E1). This phase includes a platform functional check (satellite basic functions and health verification), ground segment acquisition (including X-band communication), instrument switch-on and functional/health verification, in-orbit verification of Level 1 and Level 2 performance, specific operations for the calibration/validation campaign, and Level 1 and Level 2 processor maintenance and evolution activities. The validation of the satellite and its instruments' performance will start during the commissioning phase. The commissioning phase is expected to be completed 6 months after launch.
-
Operational phase (also referred to as Phase E2). This phase includes science data acquisition, data (re)processing and distribution, the continuation of data product monitoring and validation, Level 2 processor maintenance and evolution, and satellite and instrument calibration maintenance. This phase is expected to last at least 2.5 years, with a sufficient consumables margin for 1 additional year. With the start of this phase, Level 1 and Level 2 data products are expected to become publicly available.
-
End of mission (also referred to Phase F). During this phase, the satellite will no longer be operational; however, ground segment activities, such as reprocessing, may still continue.
8.1 Ground segment overview
The EarthCARE System consists of the EarthCARE satellite (space segment), the ground segment at ESA and JAXA, the interfaces to the EarthCARE users, and ECMWF as the provider of weather forecast data (auxiliary data product X-MET). ECMWF also receives EarthCARE data for quality monitoring and assimilation.
The EarthCARE ground segment is composed of the following:
-
The ESA Flight Operations Segment (FOS). This is located in the ESA European Space Operations Centre (ESOC), Darmstadt, Germany. It is in particular responsible for control and manoeuvring of the satellite, commanding and health monitoring of the satellite and the instruments, mission planning, orbit (control, determination, prediction), on-board software management, and dissemination and archiving of housekeeping telemetry (HK/TM) data.
-
The ESA PDGS. This is located in ESA ESRIN, Frascati, Italy. It is in particular responsible for the acquisition of science data; the generation of Level 0, Level 1 (except for CPR) and ESA Level 2 data products; the archiving of data products and auxiliary data; user services; the distribution of Level 1 and Level 2 data products to users; and the monitoring of the performance of the instruments, as well as the overall mission performance, the provision of diagnostic data to the algorithm developers and validation teams, the maintenance of the instruments' configuration settings, the calibration of the instruments, and the mission planning.
-
The JAXA Satellite Applications and Operations Center (SAOC). This is located in Tsukuba, Japan, and is in particular responsible for the reception of CPR Level 0 data and related .aux files received from PDGS and the production of the CPR Level 1 data product, support to ESA for CPR operations, CPR performance monitoring, CPR calibration and instrument parameter updates as required, the reception or processing of other Level 1 (and Level 2) products from ESA, user services, and the distribution of Level 1 and Level 2 data products to JAXA users.
The S-band station, at Kiruna, Sweden, will downlink the telemetry (TM) data from the satellite and uplink the telecommand (TC) data. The two X-band stations for transmission of the science data will be located in Kiruna, Sweden, and Inuvik, Canada.
Figure 23 illustrates the elements and data flow of the overall EarthCARE ground segment structure. Auxiliary meteorological data are provided by ECMWF.
Figure 24 shows the components of JAXA’s EarthCARE mission operation system. It consists of a data-processing system (DPS), a data transfer and management system (E-XING), a consolidated data dissemination system (G-Portal), and the third-generation JAXA Supercomputer System (JSS3) used for product reprocessing. DPS receives CPR L0, L1 of ESA sensors and other data necessary for data processing from ESA, and it generates CPR L1b and higher-level products. E-XING transfers relevant data between DPS and facilities inside and outside of JAXA, and it archives the data in the dedicated storage system. G-Portal disseminates JAXA and ESA products on the Internet, providing users with functions of catalogue search and data retrieval. The reprocessing of products after an algorithm update is performed on the Linux server cluster deployed in JSS3.
Data latency. The EarthCARE ground segment has been designed for fast data processing to prevent data accumulation and backlog. All instruments' Level 1b/1c and Level 2a data products will be made available to users within 24 h of sensing, Level 2b products will be made available to users within 48 h of sensing and any contingency case will be recovered within 5 d.
8.2 Data products
All data products generated in the frame of the EarthCARE mission are categorised per product level. The definition of each product level has been specifically tuned for the EarthCARE mission and can be found in Table 8.
L1 products of an instrument of EarthCARE will be developed by the agency responsible for the instrument; namely L1 products of the ATLID, MSI and BBR are developed by ESA and the CPR L1 product is produced by JAXA. As for L2 products, both agencies develop algorithms independently, although continuous exchange of information is being conducted between Japan and Europe through the Joint Algorithm Development Endeavor (JADE) under the framework of the Joint Mission Advisory Group (JMAG). The production model and L2 data products of ESA and JAXA are described by Eisinger et al. (2023).
Users can acquire products of both agencies from the websites of both agencies; the ESA website provides JAXA products in addition to ESA products, and the JAXA G-Portal website distributes ESA and JAXA products as well. Therefore, users can access all the products in a single place.
Eisinger et al. (2023)Qu et al. (2022)Donovan et al. (2023)Wandinger et al. (2022)Mason et al. (2023)Hashino et al. (2013, 2016)Roh et al. (2023)L0 data products consist of re-sorted science packages of the individual instruments transferred from the satellite to the ground segment. They are not available to users and therefore not described. The ESA PDGS processes the L0 data products of the ATLID, MSI and BBR into the calibrated instrument Level 1b (L1b) data products. The JAXA PDGS does the same for the CPR data products. Level 1a data products are not produced.
The L1b products are calibrated and geolocated instrument data in scientific units and consist of the data products C-NOM, A-NOM, M-NOM, B-NOM and B-SNG, plus additional calibration products, which are generally not of interest to scientific users and therefore not described here. M-RGR is a re-gridded MSI L1c product. (Note that users interested in BBR radiances are advised to use BM-RAD instead of B-NOM or B-SNG.) As the only L1d data product, the X-JSG data product contains the “joint standard grid” used by many of the L2a and L2b processors as a retrieval grid. This is required in particular for synergistic retrievals as every instrument has its own characteristic spatial sampling, while for the geophysical product retrieval, mapping of the various input data onto a joint grid is a prerequisite. The L1b, L1c and L1d data products are described by Eisinger et al. (2023).
The L2 data products include a comprehensive range of geophysical parameters related to aerosols, clouds, precipitation and radiation. Figures 25 and 26 provide an overview of the ESA and JAXA L2a and L2b data products containing retrieved aerosol, cloud, precipitation and radiation parameters.
A complete list of all L2a and L2b data products with references to the corresponding publications describing the products, retrieval algorithms and their validation is given by Eisinger et al. (2023).
8.3 Level 2 retrievals' supporting science
For preparation and verification of the scientific processors that produce the L2a and L2b data products, a number of supporting scientific activities have been carried out. In particular, the development of an EarthCARE end-to-end simulator and suitable test scenes was required. Furthermore, an aerosol model as a basis for aerosol typing was developed and a study was carried out to inter-compare the microphysical assumptions used in the cloud retrievals in the different frequency/wavelength domains applicable to EarthCARE observations.
To this end, the EarthCARE end-to-end simulator (E3SIM) has been developed to simulate EarthCARE observations and provide a development and test environment for L1 and L2 processors (Eisinger et al., 2023).
For the purpose of algorithm development and testing, test scenes have been generated using the Canadian GEM model (Qu et al., 2022). These test scenes have been processed by E3SIM to generate simulated CPR, ATLID, MSI and BBR L1b data products (Donovan et al., 2023).
To stimulate algorithm development in Japan, JAXA has been developing the Joint Simulator for Satellite Sensors (Joint-Simulator), which is a suite of software that simulates EarthCARE observations using numerical model data as input (Hashino et al., 2013, 2016). The basic structure of the Joint-Simulator is inherited from the Satellite Data Simulator Unit (SDSU; Masunaga et al., 2010) and NASA Goddard SDSU (Matsui et al., 2013). The EarthCARE sensors are simulated by the EarthCARE Active Sensor Simulator, EASE (Okamoto et al., 2003, 2007, 2008; Nishizawa et al., 2008), for the Doppler CPR and ATLID; by Rstar-6b Nakajima and Tanaka (1986, 1988) for the MSI; and by MSTRN-X (Sekiguchi and Nakajima, 2008) for the BBR. The main purpose of the Joint-Simulator is to simulate the signals from EarthCARE’s four sensors within the model for validating the aerosol and cloud fields simulated by a cloud-resolving model (CRM) and general circulation model (GCM) after the launch. The Joint-Simulator has also been used with a secondary purpose for algorithm evaluation before the launch. Here, the Joint-Simulator is applied to a Japanese global cloud-resolving model, the Nonhydrostatic ICosahedral Atmospheric Model (NICAM), and produces EarthCARE Level 1 test products that are then used for preparing EarthCARE algorithms after the launch. For example, Hagihara et al. (2021) studied Doppler velocity simulations from the Joint-Simulator and showed the standard deviation of random error in the cirrus clouds case for the Doppler simulation was decreased to 0.5 m s−1 for −10 dBZe after 10 km horizontal integration. Hagihara et al. (2022) showed another study on a global assessment of Doppler velocity errors of the EarthCARE CPR using NICAM and the Joint-simulator. The MSI observation data are expected to be unique where their centre wavelength has a dependency on their position in the across-track observation swath (SMILE effect; see Sect. 5.3). Wang et al. (2022) investigated the impact of this MSI SMILE effect on the cloud properties using the simulated L1b data by the Joint-Simulator.
For aerosol retrievals, a common baseline for the development, evaluation and implementation of EarthCARE algorithms is required to ensure consistency of different aerosol products from the multi-instrument platform and to facilitate the conformity to specifications of broadband optical properties needed for the EarthCARE radiative-closure efforts. This hybrid end-to-end aerosol classification model is described by Wandinger et al. (2022).
A range of L2a and L2b data products are available with retrieved aerosol, cloud and precipitation quantities from one or several of the EarthCARE instruments, exploiting the fact that the CPR, ATLID and MSI operate in wavelength or frequency regions that are very different from each other.
A total of 14 ESA L2a and 10 L2b products are based on individually optimised retrieval algorithms for each one of the EarthCARE instruments or combination of instruments for synergistic exploitation. Different instruments have different and specific sensitivities to the various atmospheric parameters of interest, and retrieval algorithms must use specific assumptions, in particular, related to microphysics and optical properties. Mason et al. (2023) provide an inter-comparison overview of the different cloud, aerosol and precipitation retrievals of the European/Canadian L2a and L2b data products.
After the commissioning phase, expected to last 6 months after launch, there will be full and open access to the EarthCARE L1 and L2 data sets, free of charge, provided online via the ESA and JAXA websites. Access to data is granted upon electronic online registration, with the user accepting the terms and conditions. The terms and conditions cover all types of use. No specific proposal describing the use of the data needs to be submitted. During the commissioning phase, data access is limited to the accepted participants in the validation programme.
The EarthCARE mission is implemented by ESA in cooperation with JAXA to acquire key observables needed to better understand the role of clouds, aerosol and radiation in climate. The satellite platform and three instruments, the cloud–aerosol lidar instrument ATLID, the imager MSI and the radiometer BBR, are provided by ESA; the cloud Doppler radar CPR is provided by JAXA. All instruments have been characterised, and their required performances have been confirmed. The satellite will fly in a sun-synchronous polar orbit with a descending node at 14:00 MLST mean local solar time at the Equator. The EarthCARE ground segment will operationally produce a large number of data products, including comprehensive sets of geophysical data related to aerosol, cloud, precipitation and vertical motion of cloud particles, together with corresponding heating profiles, broadband radiative fluxes and radiances. The mission will be launched in 2024.
This overview paper refers generally to software code for Level 2 processors and for atmospheric models that have been used to prepare data and simulations. The references for software repositories can be found in each of the dedicated papers that are contained within this special issue on EarthCARE Level 2 algorithms and data products.
The EarthCARE Level 2 demonstration products from simulated scenes are available at https://doi.org/10.5281/zenodo.7117115 (van Zadelhoff et al., 2023b).
DB, the ESA EarthCARE Project Manager, and ET, the JAXA EarthCARE CPR Project Manager, provided overall guidance on the technical and programmatic context. KW contributed in particular to the MSI and BBR instrument and performance descriptions. GT contributed to the ATLID instrument and performance description. HN and YO contributed to the CPR instrument and performance description. PD contributed to the description of the ESA ground segment. MT contributed to the description of the JAXA ground segment. ME ensured consistency with (Eisinger et al., 2023). RK, SR, MK and TT contributed to the science background and data product overview. TW and TK prepared the manuscript with contributions from all co-authors.
The contact author has declared that none of the authors has any competing interests.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This article is part of the special issue “EarthCARE Level 2 algorithms and data products”. It is not associated with a conference.
We would like to thank our industrial partners working on the satellite and payload, including the satellite prime, Airbus Defence and Space (Germany); the base platform responsible, Airbus Defence and Space (UK); and the instrument teams at NEC (Japan), SSTL (UK), TNO (the Netherlands), Airbus Defence and Space (France), Selex Galileo (Italy), STFC RAL Space (UK), Thales Alenia Space (UK), and Leonardo (Italy). Sincere thanks also go to the members of the Joint Mission Advisory Group that advises ESA and JAXA on science aspects of the mission and to the European, Canadian and Japanese members of the Level 2 development teams and ECMWF, as well as to the principal investigators and project scientists of NASA missions CloudSat and CALIPSO and the CERES instrument for their long-standing support and advice. We would also like to thank our colleagues and members of the ESA and JAXA teams and NICT.
A special acknowledgement is appended at the final edit, following the sudden and unexpected passing of ESA’s EarthCARE mission scientist and this paper’s principal author, Tobias Wehr. Tobias was eagerly looking forward to EarthCARE’s launch, a mission to which he had dedicated a not inconsiderable span of his career. His support for the science community, collaborative approach and enthusiasm for the mission science will not be forgotten and nor will the tragedy that he will not get to see the fruition of his work.
This paper was edited by Ulla Wandinger and reviewed by two anonymous referees.
Barker, H. W., Cole, J. N. S., Qu, Z., Villefranque, N., and Shephard, M.: Radiative closure assessment of retrieved cloud and aerosol properties for the EarthCARE mission: the ACMB-DF product, EGUSphere, in preparation, 2023. a, b
Chepfer, H., Noel, V., Winker, D., and Chiriaco, M.: Where and when will we observe cloud changes due to climate warming?, Geophys. Res. Lett., 41, 8387–8395, https://doi.org/10.1002/2014GL061792, 2014. a
Cole, J. N. S., Barker, H. W., Qu, Z., Villefranque, N., and Shephard, M. W.: Broadband Radiative Quantities for the EarthCARE Mission: The ACM-COM and ACM-RT Products, Atmos. Meas. Tech. Discuss. [preprint], https://doi.org/10.5194/amt-2022-304, in review, 2022. a, b
Cosentino, A., D'Ottavi, A., Sapia, A., and Suetta, E.: Spaceborne lasers development for ALADIN and ATLID instruments, in: 2012 IEEE International Geoscience and Remote Sensing Symposium, 5673–5676, https://doi.org/10.1109/IGARSS.2012.6352324, 2012. a
do Carmo, J. P., de Villele, G., Wallace, K., Lefebvre, A., Ghose, K., Kanitz, T., Chassat, F., Corselle, B., Belhadj, T., and Bravetti, P.: ATmospheric LIDar (ATLID): Pre-Launch Testing and Calibration of the European Space Agency Instrument That Will Measure Aerosols and Thin Clouds in the Atmosphere, Atmosphere, 12, 76, https://doi.org/10.3390/atmos12010076, 2021. a
Docter, N., Preusker, R., Filipitsch, F., Kritten, L., Schmidt, F., and Fischer, J.: Aerosol optical depth retrieval from the EarthCARE multi-spectral imager: the M-AOT product, EGUsphere, 2023, 1–31, https://doi.org/10.5194/egusphere-2023-150, 2023. a
Donovan, D., van Zadelhoff, G.-J., and Wang, P.: The ATLID L2a profile processor (A-AER, A-EBD, A-TC and A-ICE products), Atmos. Meas. Tech., in preparation, 2022. a
Donovan, D. P., Kollias, P., Velázquez Blázquez, A., and van Zadelhoff, G.-J.: The Generation of EarthCARE L1 Test Data sets Using Atmospheric Model Data Sets, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2023-384, 2023. a, b
Eisinger, M., Wehr, T., Kubota, T., Bernaerts, D., Wallace, K., and Marnas, F.: The EarthCARE Mission – Science Data Processing Chain Overview, EGUSphere, in preparation, 2023. a, b, c, d, e, f, g, h, i, j, k, l, m
Fisher, J., Baumback, M. M., Bowles, J. H., Grossmann, J. M., and Antoniades, J. A.: Comparison of low-cost hyperspectral sensors, in: Imaging Spectrometry IV, edited by: Descour, M. R. and Shen, S. S., vol. 3438, 23–30, International Society for Optics and Photonics, SPIE, https://doi.org/10.1117/12.328112, 1998. a
Haarig, M., Hünerbein, A., Wandinger, U., Docter, N., Bley, S., Donovan, D., and van Zadelhoff, G.-J.: Cloud top heights and aerosol columnar properties from combined EarthCARE lidar and imager observations: the AM-CTH and AM-ACD products, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2023-327, 2023. a, b
Hagihara, Y., Ohno, Y., Horie, H., Roh, W., Satoh, M., Kubota, T., and Oki, R.: Assessments of Doppler Velocity Errors of EarthCARE Cloud Profiling Radar Using Global Cloud System Resolving Simulations: Effects of Doppler Broadening and Folding, IEEE Trans. Geosci. Remote Sens., 60, 1–9, 2021. a
Hagihara, Y., Ohno, Y., Horie, H., Roh, W., Satoh, M., and Kubota, T.: Global evaluation of Doppler velocity errors of EarthCARE Cloud Profiling Radar using global storm-resolving simulation, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2022-1255, 2022. a
Harries, J. E., Russell, J. E., Hanafin, J. A., Brindley, H., Futyan, J., Rufus, J., Kellock, S., Matthews, G., Wrigley, R., Last, A., Mueller, J., Mossavati, R., Ashmall, J., Sawyer, E., Parker, D., Caldwell, M., Allan, P. M., Smith, A., Bates, M. J., Coan, B., Stewart, B. C., Lepine, D. R., Cornwall, L. A., Ricketts, D. R. C. M. J., Drummond, D., Smart, D., Cutler, R., Dewitte, S., Clerbaux, N., Gonzalez, L., Ipe, A., Bertrand, C., Joukoff, A., Crommelynck, D., Nelms, N., LLewellyn-Jones, D. T., Butcher, G., Smith, G. L., Szewczyk, Z. P., Mlynczak, P. E., Slingo, A., Allan, R. P., and Ringer, M. A.: The Geostationary Earth Radiation Budget Project, B. Am. Meteorol. Soc., 86, 945–960, https://doi.org/10.1175/BAMS-86-7-945, 2005. a
Hashino, T., Satoh, M., Hagihara, Y., Kubota, T., Matsui, T., Nasuno, T., and Okamoto, H.: Evaluating cloud microphysics from NICAM against CloudSat and CALIPSO, J. Geophys. Res.-Atmos., 118, 7273–7292, 2013. a, b
Hashino, T., Satoh, M., Hagihara, Y., Kato, S., Kubota, T., Matsui, T., Nasuno, T., Okamoto, H., and Sekiguchi, M.: Evaluating Arctic cloud radiative effects simulated by NICAM with A-train, J. Geophys. Res.-Atmos., 121, 7041–7063, 2016. a, b
Heliere, A., Wallace, K., do Cormo, J. P., Eisinger, M., and Lefebvre, A.: EarthCARE instruments description, European Space Agency, https://earth.esa.int/eogateway/documents/20142/37627/EarthCARE-instrument-descriptions.pdf (last access: 3 July 2023), 2017. a
Hünerbein, A., Bley, S., Horn, S., Deneke, H., and Walther, A.: Cloud mask algorithm from the EarthCARE multi-spectral imager: the M-CM products, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2022-1240, 2022. a
Hünerbein, A., Bley, S., Deneke, H., Meirink, J. F., van Zadelhoff, G.-J., and Walther, A.: Cloud optical and physical properties retrieval from EarthCARE multi-spectral imager: the M-COP products, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2023-305, 2023. a
Illingworth, A. J., Barker, H. W., Beljaars, A., Ceccaldi, M., Chepfer, H., Clerbaux, N., Cole, J., Delanoë, J., Domenech, C., Donovan, D. P., Fukuda, S., Hirakata, M., Hogan, R. J., Huenerbein, A., Kollias, P., Kubota, T., Nakajima, T., Nakajima, T. Y., Nishizawa, T., Ohno, Y., Okamoto, H., Oki, R., Sato, K., Satoh, M., Shephard, M. W., Velázquez-Blázquez, A., Wandinger, U., Wehr, T., and van Zadelhoff, G.-J.: The EarthCARE Satellite: The Next Step Forward in Global Measurements of Clouds, Aerosols, Precipitation, and Radiation, B. Am. Meteorol. Soc., 96, 1311–1332, https://doi.org/10.1175/BAMS-D-12-00227.1, 2015. a
IPCC 2021: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Masson-Delmotte, V., Zhai, P., Pirani, A., Connors, S. L., Péan, C., Berger, S., Caud, N., Chen, Y., Goldfarb, L., Gomis, M. I., Huang, M., Leitzell, K., Lonnoy, E., Matthews, J. B. R., Maycock, T. K., Waterfield, T., Yelekçi, O., Yu, R., and Zhou, B.: Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, https://doi.org/10.1017/9781009157896, in press, 2021. a
Irbah, A., Delanoë, J., van Zadelhoff, G.-J., Donovan, D. P., Kollias, P., Puigdomènech Treserras, B., Mason, S., Hogan, R. J., and Tatarevic, A.: The classification of atmospheric hydrometeors and aerosols from the EarthCARE radar and lidar: the A-TC, C-TC and AC-TC products, Atmos. Meas. Tech., 16, 2795–2820, https://doi.org/10.5194/amt-16-2795-2023, 2023. a, b
Janisková, M. and Fielding, M. D.: Direct 4D-Var assimilation of space-borne cloud radar and lidar observations. Part II: Impact on analysis and subsequent forecast, Q. J. Roy. Meteorol. Soc., 146, 3900–3916, https://doi.org/10.1002/qj.3879, 2020. a, b
Kato, S., Rose, F. G., Sun-Mack, S., Miller, W. F., Chen, Y., Rutan, D. A., Stephens, G. L., Loeb, N. G., Minnis, P., Wielicki, B. A., Winker, D. M., Charlock, T. P., Stackhouse, P. W., Xu, K.-M., and Collins, W. D.: Improvements of top-of-atmosphere and surface irradiance computations with CALIPSO-, CloudSat-, and MODIS-derived cloud and aerosol properties, J. Geophys. Res., 116, D19209, https://doi.org/10.1029/2011JD016050, 2011. a
Kollias, P., Puidgomènech Treserras, B., Battaglia, A., Borque, P., and Tatarevic, A.: Processing reflectivity and Doppler velocity from EarthCARE’s cloud profiling radar: the C-FMR, C-CD and C-APC products, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2022-1284, 2022. a
Kudo, R., Nishizawa, T., and Aoyagi, T.: Vertical profiles of aerosol optical properties and the solar heating rate estimated by combining sky radiometer and lidar measurements, Atmos. Meas. Tech., 9, 3223–3243, https://doi.org/10.5194/amt-9-3223-2016, 2016. a
Mason, S. L., Hogan, R. J., Bozzo, A., and Pounder, N. L.: A unified synergistic retrieval of clouds, aerosols and precipitation from EarthCARE: the ACM-CAP product, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2022-1195, 2022. a, b, c
Mason, S., Hogan, R. J., Donovan, D., van Zadelhoff, G.-J., Kollias, P., Treserras, B. P., Qu, Z., Cole, J. N. S., Hünerbein, A., and Doctor, N.: An intercomparison of EarthCARE cloud, aerosol and precipitation retrieval products, EGUSphere, in preparation, 2023. a, b
Masunaga, H., Matsui, T., Tao, W.-K., Hou, A. Y., Kummerow, C. D., Nakajima, T., Bauer, P., Olson, W. S., Sekiguchi, M., and Nakajima, T. Y.: Satellite data simulator unit: A multisensor, multispectral satellite simulator package, B. Am. Meteorol. Soc., 91, 1625–1632, 2010. a
Matsui, T., Iguchi, T., Li, X., Han, M., Tao, W.-K., Petersen, W., L'Ecuyer, T., Meneghini, R., Olson, W., Kummerow, C. D., Hou, A. Y., Schwaller, M. R., Stocker, E. F., and Kwiatkowski, J.: GPM satellite simulator over ground validation sites, B. Am. Meteorol. Soc., 94, 1653–1660, https://doi.org/10.1175/BAMS-D-12-00160.1, 2013. a
Mroz, K., Treserras, B. P., Battaglia, A., Kollias, P., Tatarevic, A., and Tridon, F.: Cloud and precipitation microphysical retrievals from the EarthCARE Cloud Profiling Radar: the C-CLD product, Atmos. Meas. Tech., 16, 2865–2888, https://doi.org/10.5194/amt-16-2865-2023, 2023. a
Nakajima, T. and Tanaka, M.: Matrix formulations for the transfer of solar radiation in a plane-parallel scattering atmosphere, J. Quant. Spectrosc. Ra., 35, 13–21, 1986. a
Nakajima, T. and Tanaka, M.: Algorithms for radiative intensity calculations in moderately thick atmospheres using a truncation approximation, J. Quant. Spectrosc. Ra., 40, 51–69, 1988. a
Nakajima, T. Y., Ishida, H., Nagao, T. M., Hori, M., Letu, H., Higuchi, R., Tamaru, N., Imoto, N., and Yamazaki, A.: Theoretical basis of the algorithms and early phase results of the GCOM-C (Shikisai) SGLI cloud products, Prog. Earth Planet. Sci., 6, 1–25, 2019. a
Nakatsuka, H., Kimura, T., Seki, Y., Kadosaki, G., Iide, Y., Okada, K., Yamaguchi, J., Takahashi, N., Ohno, Y., Horie, H., and Sato, K.: Design and development status of the EarthCARE Cloud Profiling Radar, in: 2012 IEEE International Geoscience and Remote Sensing Symposium, 2415–2418, https://doi.org/10.1109/IGARSS.2012.6351004, 2012. a
Nishizawa, T., Sugimoto, N., Matsui, I., Shimizu, A., Tatarov, B., and Okamoto, H.: Algorithm to retrieve aerosol optical properties from high-spectral-resolution lidar and polarization Mie-scattering lidar measurements, IEEE T. Geosci. Remote Sens., 46, 4094–4103, 2008. a
Nishizawa, T., Kudo, R., Oikawa, E., Higurashi, A., Yoshitaka, J., and Okamoto, H.: Retrieval of aerosol optical properties from the EarthCARE high-spectral resolution lidar measurements, Atmos. Meas. Tech., in preparation, 2023. a, b
Oikawa, E., Nakajima, T., and Winker, D.: An evaluation of the shortwave direct aerosol radiative forcing using CALIOP and MODIS observations, J. Geophys. Res.-Atmos., 123, 1211–1233, 2018. a
Okamoto, H., Iwasaki, S., Yasui, M., Horie, H., Kuroiwa, H., and Kumagai, H.: An algorithm for retrieval of cloud microphysics using 95-GHz cloud radar and lidar, J. Geophys. Res.-Atmos., 108, D7, https://doi.org/10.1029/2001JD001225, 2003. a
Okamoto, H., Nishizawa, T., Takemura, T., Kumagai, H., Kuroiwa, H., Sugimoto, N., Matsui, I., Shimizu, A., Emori, S., Kamei, A., and Nakajima, T.: Vertical cloud structure observed from shipborne radar and lidar: Midlatitude case study during the MR01/K02 cruise of the research vessel Mirai, J. Geophys. Res.-Atmos., 112, D8, https://doi.org/10.1029/2006JD007628, 2007. a
Okamoto, H., Nishizawa, T., Takemura, T., Sato, K., Kumagai, H., Ohno, Y., Sugimoto, N., Shimizu, A., Matsui, I., and Nakajima, T.: Vertical cloud properties in the tropical western Pacific Ocean: Validation of the CCSR/NIES/FRCGC GCM by shipborne radar and lidar, J. Geophys. Res.-Atmos., 113, D24, https://doi.org/10.1029/2008JD009812, 2008. a
Okamoto, H., Sato, K., Oikawa, E., Ishimoto, H., Ohno, Y., Horie, H., Hagihara, Y., Nishizawa, T., Kudo, R., Higurashi, A., Jin, Y., Nakajima, T. Y., Wang, M., Roh, W., Satoh, M., Suzuki, K., Kubota, T., Yamauchi, A., Sekiguchi, M., and Nagao, T. M.: New perspectives of clouds, radiation and dynamics from EarthCARE observation, Atmos. Meas. Tech., in preparation, 2023. a, b, c, d
Okata, M., Nakajima, T., Suzuki, K., Inoue, T., Nakajima, T., and Okamoto, H.: A study on radiative transfer effects in 3-D cloudy atmosphere using satellite data, J. Geophys. Res.-Atmos., 122, 443–468, 2017. a
Proulx, C., Allard, M., T. Pope, B. T., Williamson, F., Julien, C., Larouche, C., Delderfield, J., and Parker, D.: Performance characterization of the EarthCARE BBR Detectors, Proc. SPIE 10565, International Conference on Space Optics – ICSO 2010, 1056508 (20 November 2017), https://doi.org/10.1117/12.2309219, 2017. a
Qu, Z., Donovan, D. P., Barker, H. W., Cole, J. N. S., Shephard, M. W., and Huijnen, V.: Numerical Model Generation of Test Frames for Pre-launch Studies of EarthCARE’s Retrieval Algorithms and Data Management System, Atmos. Meas. Tech. Discuss. [preprint], https://doi.org/10.5194/amt-2022-300, in review, 2022. a, b
Qu, Z., Barker, H. W., Cole, J. N. S., and Shephard, M. W.: Across-track extension of retrieved cloud and aerosol properties for the EarthCARE mission: the ACMB-3D product, Atmos. Meas. Tech., 16, 2319–2331, https://doi.org/10.5194/amt-16-2319-2023, 2023. a, b
Roh, W., Satoh, M., Hagihara, Y., Ohno, Y., and Horie, H.: Investigation of the performances of the EarthCARE CPR using NICAM and the ground data in Japan, EGUSphere [preprint], https://doi.org/10.5194/amt-16-3331-2023, 2023. a
Sekiguchi, M. and Nakajima, T.: A k-distribution-based radiation code and its computational optimization for an atmospheric general circulation model, J. Quant. Spectrosc. Ra., 109, 2779–2793, 2008. a
Sherwood, S. C., Webb, M. J., Annan, J. D., Armour, K. C., Forster, P. M., Hargreaves, J. C., Hegerl, G., Klein, S. A., Marvel, K. D., Rohling, E. J., Watanabe, M., Andrews, T., Braconnot, P., Bretherton, C. S., Foster, G. L., Hausfather, Z., von der Heydt, A. S., Knutti, R., Mauritsen, T., Norris, J. R., Proistosescu, C., Rugenstein, M., Schmidt, G. A., Tokarska, K. B., and Zelinka, M. D.: An Assessment of Earth's Climate Sensitivity Using Multiple Lines of Evidence, Rev. Geophys., 58, e2019RG000678, https://doi.org/10.1029/2019RG000678, 2020. a
Stephens, G., Winker, D., Pelon, J., Trepte, C., Vane, D., Yuhas, C., L’Ecuyer, T., and Lebsock, M.: CloudSat and CALIPSO within the A-Train: Ten Years of Actively Observing the Earth System, B. Am. Meteorol. Soc., 99, 569–581, https://doi.org/10.1175/BAMS-D-16-0324.1, 2018. a, b
Stephens, G. L., Vane, D. G., Tanelli, S., Im, E., Durden, S., Rokey, M., Reinke, D., Partain, P., Mace, G. G., Austin, R., L'Ecuyer, T., Haynes, J., Lebsock, M., Suzuki, K., Waliser, D., Wu, D., Kay, J., Gettelman, A., Wang, Z., and Marchand, R.: CloudSat mission: Performance and early science after the first year of operation, J. Geophys. Res.-Atmos., 113, D8, https://doi.org/10.1029/2008JD009982, 2008. a
Straume, A., Rennie, M., Isaksen, L., de Kloe, J., Marseille, G.-J., Stoffelen, A., Flament, T., Stieglitz, H., Dabas, A., Huber, D., Reitebuch, O., Lemmerz, C., Lux, O., Marksteiner, U., Weiler, F., Witschas, B., Meringer, M., Schmidt, K., Nikolaus, I., Geiss, A., Flamant, P., Kanitz, T., Wernham, D., von Bismarck, J., Bley, S., Fehr, T., Floberghagen, R., and Parinello, T.: ESA's Space-Based Doppler Wind Lidar Mission Aeolus – First Wind and Aerosol Product Assessment Results, EPJ Web Conf., 237, 01007, https://doi.org/10.1051/epjconf/202023701007, 2020. a
Takahashi, H., Lebsock, M. D., Richardson, M., Marchand, R., and Kay, J. E.: When Will Spaceborne Cloud Radar Detect Upward Shifts in Cloud Heights?, J. Geophys. Res.-Atmos., 124, 7270–7285, https://doi.org/10.1029/2018JD030242, 2019. a
Vaillant de Guélis, T., Chepfer, H., Guzman, R., Bonazzola, M., Winker, D. M., and Noel, V.: Space lidar observations constrain longwave cloud feedback, Sci. Rep., 8, 1–8, https://doi.org/10.1038/s41598-018-34943-1, 2018. a
van Zadelhoff, G.-J., Donovan, D. P., and Wang, P.: Detection of aerosol and cloud features for the EarthCARE lidar ATLID: the A-FM product, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2023-145, 2023a. a
van Zadelhoff, G.-J., Barker, H. W., Baudrez, E., Bley, S., Clerbaux, N., Cole, J. N. S., de Kloe, J., Docter, N., Domenech, C., Donovan, D. P., Dufresne, J.-L., Eisinger, M., Fischer, J., García-Marañón, R., Haarig, M., Hogan, R. J., Hünerbein, A., Kollias, P., Koopman, R., Madenach, N., Mason, S. L., Preusker, R., Puigdomènech Treserras, B., Qu, Z., Ruiz-Saldaña, M., Shephard, M., Velázquez-Blazquez, A., Villefranque, N., Wandinger, U., Wang, P., and Wehr, T.: EarthCARE level-2 demonstration products from simulated scenes, Zenodo [data set], https://doi.org/10.5281/zenodo.7117115, 2023b. a
Velázquez-Blázquez, A., Baudrez, E., Clerbaux, N., and Domenech, C.: Unfiltering of the EarthCARE Broadband Radiometer (BBR) observations: the BM-RAD product, Atmos. Meas. Tech., in preparation, 2023a. a, b, c
Velázquez-Blázquez, A., Baudrez, E., Clerbaux, N., Domenech, C., Marañón, R. G., and Madenach, N.: Towards instantaneous top-of-atmosphere fluxes from EarthCARE: The BMA-FLX product, EGUSphere, in preparation, 2023b. a, b, c
Wallace, K., Perez-Albinana, A., Lemanczyk, J., Heliere, A., Wehr, T., Eisinger, M., Lefebvre, A., Nakatsuka, H., and Tomita, E.: The EarthCARE satellite payload, in: Sensors, Systems, and Next-Generation Satellites XVIII, edited by: Meynart, R., Neeck, S. P., and Shimoda, H., vol. 9241, p. 92410F, International Society for Optics and Photonics, SPIE, https://doi.org/10.1117/12.2067208, 2014. a
Wallace, K., Hélière, A., Lefebvre, A., Eisinger, M., and Wehr, T.: Status of ESA's EarthCARE mission, passive instruments payload, in: Earth Observing Systems XXI, edited by: Butler, J. J., Xiong, X. J., and Gu, X., vol. 9972, p. 997214, International Society for Optics and Photonics, SPIE, https://doi.org/10.1117/12.2236498, 2016. a
Wandinger, U., Floutsi, A. A., Baars, H., Haarig, M., Ansmann, A., Hünerbein, A., Docter, N., Donovan, D., van Zadelhoff, G.-J., Mason, S., and Cole, J.: HETEAC – The Hybrid End-To-End Aerosol Classification model for EarthCARE, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2022-1241, 2022. a, b
Wandinger, U., Haarig, M., Baars, H., Donovan, D., and van Zadelhoff, G.-J.: Cloud top heights and aerosol layer properties from EarthCARE lidar observations: the A-CTH and A-ALD products, EGUSphere [preprint], https://doi.org/10.5194/egusphere-2023-748, 2023. a
Wang, M., Nakajima, T. Y., Roh, W., Satoh, M., Suzuki, K., Kubota, T., and Yoshida, M.: Evaluation of the smile effect on the Earth Clouds, Aerosols and Radiation Explorer (EarthCARE)/Multi-Spectral Imager (MSI) cloud product, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2022-736, 2022. a
Wehr, T. (Ed.): EarthCARE Mission Requirements Document, Earth and Mission Science Division, European Space Agency, https://doi.org/10.5270/esa.earthcare-mrd.2006, 2006. a
Weiler, F., Kanitz, T., Wernham, D., Rennie, M., Huber, D., Schillinger, M., Saint-Pe, O., Bell, R., Parrinello, T., and Reitebuch, O.: Characterization of dark current signal measurements of the ACCDs used on board the Aeolus satellite, Atmos. Meas. Tech., 14, 5153–5177, https://doi.org/10.5194/amt-14-5153-2021, 2021. a
Wernham, D., Alves, J., Pettazzi, F., and Tighe, A. P.: Laser-induced contamination mitigation on the ALADIN laser for ADM-Aeolus, in: Laser-Induced Damage in Optical Materials: 2010, edited by: Exarhos, G. J., Gruzdev, V. E., Menapace, J. A., Ristau, D., and Soileau, M. J., vol. 7842, 394 – 405, International Society for Optics and Photonics, SPIE, https://doi.org/10.1117/12.867268, 2010. a
Wielicki, B. A., Barkstrom, B. R., Harrison, E. F., LeeIII, R. B., Smith, G. L., and Cooper, J. E.: Clouds and the Earth's Radiant Energy System (CERES):An Earth Observing System Experiment, B. Am. Meteorol. Soc., 77, 853–868, 1996. a
Winker, D., Chepfer, H., Noel, V., and Cai, X.: Observational Constraints on Cloud Feedbacks: The Role of Active Satellite Sensors, Surv. Geophys., 38, 1483–1508, https://doi.org/10.1007/s10712-017-9452-0, 2017. a
Winker, D. M., Tackett, J. L., Getzewich, B. J., Liu, Z., Vaughan, M. A., and Rogers, R. R.: The global 3-D distribution of tropospheric aerosols as characterized by CALIOP, Atmos. Chem. Phys., 13, 3345–3361, https://doi.org/10.5194/acp-13-3345-2013, 2013. a
Yamauchi, A., Suzuki, K., Oikawa, E., Sekiguchi, M., and Nagao, T. M.: Description and validation of the Japanese algorithms for radiative flux and heating rate products with all four EarthCARE instruments, EGUSphere, in preparation, 2023. a
- Abstract
- Introduction
- Mission science objectives
- Observational requirements
- EarthCARE system overview
- Satellite instruments
- EarthCARE orbit
- Mission phases
- Ground segment and data products
- Conclusions
- Code availability
- Data availability
- Author contributions
- Competing interests
- Disclaimer
- Special issue statement
- Acknowledgements
- Review statement
- References
- Abstract
- Introduction
- Mission science objectives
- Observational requirements
- EarthCARE system overview
- Satellite instruments
- EarthCARE orbit
- Mission phases
- Ground segment and data products
- Conclusions
- Code availability
- Data availability
- Author contributions
- Competing interests
- Disclaimer
- Special issue statement
- Acknowledgements
- Review statement
- References