# Impact of Lock-In Time Constant on Remote Monitoring of Trace Gas in the Atmospheric Column Using Laser Heterodyne Radiometer (LHR)

Fengjiao Shen, Gaoxuan Wang, Zhengyue Xue, Tu Tan, Zhensong Cao, Xiaoming Gao, Weidong Chen

► To cite this version:

Fengjiao Shen, Gaoxuan Wang, Zhengyue Xue, Tu Tan, Zhensong Cao, et al., Impact of Lock-In Time Constant on Remote Monitoring of Trace Gas in the Atmospheric Column Using Laser Heterodyne Radiometer (LHR). Remote Sensing, 2022, 14 (12), pp.2923. 10.3390/rs14122923. hal-04291001

HAL Id: hal-04291001

https://ulco.hal.science/hal-04291001v1

Submitted on 22 Nov 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

Technical Note

# Impact of Lock-In Time Constant on Remote Monitoring of Trace Gas in the Atmospheric Column Using Laser Heterodyne Radiometer (LHR)

Fengjiao Shen<sup>1,2</sup>, Gaoxuan Wang<sup>2</sup>, Zhengyue Xue<sup>3</sup>, Tu Tan<sup>3</sup>, Zhensong Cao<sup>3</sup>, Xiaoming Gao<sup>3</sup>
and Weidong Chen<sup>2,*</sup><sup>†</sup>

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Shen, F.; Wang, G.; Xue, Z.; Tan, T.; Cao, Z.; Gao, X.; Chen, W. Impact of Lock-In Time Constant on Remote Monitoring of Trace Gas in the Atmospheric Column Using Laser Heterodyne Radiometer (LHR). *Remote Sens.* **2022**, *14*, 2923. https://doi.org/10.3390/rs14122923

Academic Editors: Filomena Romano and Elisabetta Ricciardelli

Received: 17 May 2022

Accepted: 16 June 2022

Published: 18 June 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** The time constant selected for lock-in amplification (LIA) has a crucial impact on observed line shapes in laser heterodyne spectroscopy, in particular in the case of ground-based remote monitoring of trace gas in the atmospheric column using laser heterodyne radiometer (LHR). Conventional simulation could not allow validation of LHR spectra measured in a real and complex atmospheric environment exhibiting large temporal and spatial variability (humidity, temperature, pressure, etc) that impact significantly the measured LHR spectra profiles. High-precision spectral measurement is thus crucial to avoid any spectral distortion resulting from the measurement. In this paper, the impact of LIA time constant on spectral line shape is investigated for LHR operating in continuous laser tuning mode, based on analysis of laboratory heterodyne spectra, in terms of signal-to-noise ratio (SNR), line width broadening, absorption depth and line shift. With respect to the given frequency scanning speed in continuous mode and to the halfwidth of the absorption feature to scan, a reasonable scanning time $\Delta T_{scan}$, the time needed for scanning laser frequency through the halfwidth $\Delta v_{HWHM}$ of the absorption line, equal to or longer than 14 times of the LIA time constant $\tau$ is concluded in order to efficiently reduce the noise while without significant shift and distortion of the line shape. Experimental validation was carried out using a laser heterodyne absorption spectroscopy approach in the laboratory. Four different combinations of time constants $\tau$ and scanning time $\Delta T_{scan}$ were used to record heterodyne spectra of a CH$_4$ absorption line near 1242.00 $cm^{-1}$ in continuous laser tuning mode. An optimal combination of a scanning time of 137 ms with a time constant of 1 ms was obtained. This optimal combination was used for ground-based measurements of CH$_4$ and N$_2$O in the atmospheric column by LHR. The extracted LHR spectrum is in good agreement with a referenced TCCON (Total Carbon Column Observing Network) FT-IR (Fourier-transform infrared) spectrum.

**Keywords:** lock-in amplification time constant; laser heterodyne spectrum; continuous tuning of laser frequency; high-precision; line shape; laser heterodyne radiometer

## 1. Introduction

Since the last decade, there has been a growing development of laser heterodyne radiometers (LHR) [1–5] for ground-based remote sensing of greenhouse gases (GHGs, such as nitrous oxide (N$_2$O), methane (CH$_4$), water vapor (H$_2$O), carbon dioxide (CO$_2$) and ozone (O$_3$) in the atmospheric column. Though the vertical air transport in the troposphere is smaller than in the longitudinal direction [6] but is nonetheless of critical importance for the study of pollutant transport, climate change, global environmental parameters

<footnote n="1">¹ School of Advanced Manufacturing Engineering, Hefei University, Hefei 230601, China; shenfj@hfuu.edu.cn</footnote>

<footnote n="2">² Laboratoire de Physico-Chimie de l'Atmosphère, Université du Littoral Côte d'Opale, 59140 Dunkerque, France; gaoxuanwang@zju.edu.cn</footnote>

<footnote n="3">³ Anhui Institute of Optics and Fine Mechanics, Chinese Academy of Sciences, Hefei 230031, China; xzy1106@mail.ustc.edu.cn (Z.X.); tantu@aiofm.ac.cn (T.T.); zscao@aiofm.ac.cn (Z.C.); xmgao@aiofm.ac.cn (X.G.)</footnote>

<footnote>* Correspondence: chen@univ-littoral.fr; Tel.: +33-328658264</footnote>
and anthropogenic effects on atmospheric chemistry [7,8]. Vertical concentration profiles of key atmospheric trace gases are essential in validating chemical transport models, satellite observations, and provide evidence for environmental assessment and policy formulation [9,10]. LHR exhibits its unique advantages in comparison with conventionally used Fourier-transform spectrometer (FTS) for ground-based measurements of vertical concentration distribution, which include high spectral resolution (<10⁻³ cm⁻¹, mainly determined by the selected electronic filter bandwidths), high sensitivity (within a factor of ~2 of the quantum noise limit) [2], high spatial resolution owing to a very small field of view (~0.4 mrad, compared to ~2.4 mrad of the conventionally used FTS) and cost-effective compact instrumental dimension.

In an LHR system, a lock-in amplifier (LIA) is usually used to recover weak heterodyne signals buried in noise [11,12] through phase-sensitive detection (PSD) of signals and rejection of noises by a low-pass filter (whose cut-off frequency is determined by the LIA time constant τ) [13]. Meanwhile, the used LIA time constant might significantly impact the line shape of the recorded LHR spectrum [3]. The impact of the LIA time constant on laser heterodyne spectrum measurements was discussed for LHR operating in point-by-point tuning mode of laser frequency, where a time interval (between two sampling points) longer than ~6 times of the LIA time constant τ should be used to avoid deformation of LHR spectrum line shape [4,5].

In this paper, we report on an investigation on the impact of LIA time constant τ on LHR spectrum line shape in continuous tuning mode of the laser frequency. Continuous tuning mode is particularly interesting in fast extraction of high-precision line shape of LHR spectrum compared with point-by-point tuning mode [14]. Relationship between scanning time ΔT<sub>scan</sub> (the time needed for scanning laser frequency through the halfwidth Δν<sub>HWHM</sub> of the absorption line, which depends upon the frequency scanning speed ν<sub>sc</sub>) and the LIA time constant τ in continuous tuning mode of laser frequency has been experimentally investigated, which was then validated by laser heterodyne absorption experiment performed on the LHR measurement of the CH<sub>4</sub> absorption spectrum.

**2. Materials and Methods**

2.1. *Principle of Laser Heterodyne Radiometer*

The principle of LHR is schematically shown in Figure 1. In an LHR, the electromag-
netic field of the sunlight E<sub>S</sub>(t) = A<sub>S0</sub>(ω)cos(ωt) at frequency ω is mixed in a photomixer
with an electromagnetic field E<sub>L</sub>(t) = A<sub>L0</sub>(ω<sub>L</sub>)cos(ω<sub>L</sub>t) at frequency ω<sub>L</sub> from a local os-
cillator (LO). A beat note at an intermediate frequency (IF, |ω-ω<sub>L</sub>|), produced from the
photomixing, is amplified by a low-noise amplifier and filtered by a band-pass filter. This
heterodyne signal is then square-law detected and fed into an LIA, the output of the LIA is
sampled by a laptop for further data processing. The heterodyne absorption spectrum of
atmospheric species can be recovered by scanning the LO frequency across the absorption
line (leading to a change in intensity of the resulting radio frequency (RF) beat signal). The
shape of the ground-recovered atmospheric absorption spectrum contains information
on the vertical distributions of the target absorbers. It is possible to retrieve trace gas
abundances at different altitudes through a retrieval algorithm based on an appropriate
inversion model [3–5].
<figure><img src="image_5.png" /><figcaption>Figure 1. Basic principle of a typical laser heterodyne radiometer (LHR). RF: Radio frequency; IF: Intermediate frequency; LIA: Lock-in amplifier; LO: Local oscillator.</figcaption></figure>

## 2.2. Analytical Model

In continuous tuning mode, deformation of the LHR spectrum line shape is mainly introduced in the process of signal demodulation of the LIA, in which the whole LHR spectrum passes through a low-pass filter of the LIA, the cut-off frequency and order of the low-pass filter may have different degrees of influence on the LHR spectrum.

The cut-off frequency (fc) of the low-pass (LP) filter is driven by the LIA time constant $τ$ (in [s]) as $f_c = 1/(2πτ)$ [15]. The time constant $τ$ is a very critical parameter for the use of LIA, which impacts noise suppression and the line shape of the measured LHR spectrum. The longer time constant $τ$ can significantly filter out noise, while it also results in a blue-shift of the spectral line center (shifting to the frequency scan direction/decrease in wavelength) and then distorts the line shape (becoming broadened and asymmetric) [15,16].

The frequency shift $δν$ can be approximately expressed by [17]:

$$δν = \frac{dν}{dt}τ = v_{sc}τ \quad (1)$$

where $dν/dt = v_{sc}$ is the frequency scanning speed in [cm$^{-1}$/s]. In continuous periodical scanning, the frequency scanning speed is the ratio of the frequency scanning range $Δv_{scan}$ to the scanning period T (in [s]) that is the reciprocal of the scanning frequency.

In order to avoid deformation of the LHR spectrum line shape, the LIA time constant $τ$ has to be appropriately chosen with respect to the time $ΔT_{scan}$. Simulations of the impact of $τ$ have been performed in order to investigate the optimal $τ$ with respect to the experimentally applied frequency scanning speed. The effective transfer function of the LP filter (Figure 2) can be expressed by Equation (2).

<figure><img src="image_6.png" /><figcaption>Figure 2. Scheme of a LP filter. R is the resistor and C is the capacitor.</figcaption></figure>

According to Kirchhoff’s Laws and the definition of capacitance:

$$V_{in}(t) - V_{out}(t) = RC \frac{dV_{out}(t)}{dt} \quad (2)$$
where R is the resistance and C is the capacitance of the LP filter. For simplicity, assuming that samples of the input and output are taken at evenly spaced points in time separated by Δt, Equation (2) can be expressed in terms of Vout at time k, Vin at time k and dVout of Vout at time k and k−1, respectively:

$$V_{\text{in}}(k) - V_{\text{out}}(k) = RC \cdot \frac{V_{\text{out}}(k) - V_{\text{out}}(k-1)}{\Delta t} \quad (3)$$

Rearranging terms in Equation (3) gives a recurrence relation:

$$V_{\text{out}}(k) = V_{\text{in}}(k) \cdot \left( \frac{\Delta t}{RC+\Delta t} \right) + V_{\text{out}}(k-1) \cdot \left( \frac{RC}{RC+\Delta t} \right) \quad (4)$$

The first term on the right side of Equation (4) is the contribution from the input, and the second term is the inertia from the previous output. This differential model of the LP filter described in Equation (4) provides a way to determine the output data based on the input and the previous output values.

Based on Equation (4), the effect of the LP filter for frequency scan in continuous tuning mode is simulated by using a pseudocode algorithm with the help of Matlab. In the simulation, Vin is a Hitran simulation spectrum of CH₄ around 1233.45 cm⁻¹ with an HWHM line width of 0.014 cm⁻¹, and Vout is the simulation output resulting from different experimental parameters: scanning period T, scanning speed vsc, scanning time ΔTscan (ΔvHWHM/vsc) over a spectral range of Δvscan = 1 cm⁻¹, and time constant τ. The simulation results are overviewed in Table 1 and presented in Figure 3.

<table><caption>**Table 1.** Simulation of LP filter effect on absorption line shape with different τ. HWHM: Half Width at Half Maximum.</caption><tbody><tr><td>Halfwidth Δv<sub>HWHM</sub> (cm<sup>-1</sup>)</td><td>0.014</td><td>0.014</td><td>0.014</td><td>0.014</td><td>0.014</td><td>0.014</td><td>0.014</td><td>0.014</td><td>0.014</td><td>0.014</td></tr><tr><td>Scanning range Δv<sub>scan</sub> (cm<sup>-1</sup>)</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><td>Scanning period T (s)</td><td>0.3</td><td>0.5</td><td>1</td><td>1.5</td><td>2</td><td>2.5</td><td>3</td><td>5</td><td>10</td><td>30</td></tr><tr><td>Scan frequency v (Hz)</td><td>3.33</td><td>2</td><td>1</td><td>0.67</td><td>0.5</td><td>0.4</td><td>0.33</td><td>0.2</td><td>0.1</td><td>0.033</td></tr><tr><td>Scanning speed v<sub>sc</sub> (cm<sup>-1</sup>/s) = Δv<sub>scan</sub>/T</td><td>3.33</td><td>2</td><td>1</td><td>0.67</td><td>0.5</td><td>0.4</td><td>0.33</td><td>0.2</td><td>0.1</td><td>0.033</td></tr><tr><td>Scanning time ΔT<sub>scan</sub> (ms)</td><td>4.17</td><td>6.95</td><td>13.9</td><td>20.85</td><td>27.8</td><td>34.75</td><td>41.7</td><td>69.5</td><td>139</td><td>417</td></tr><tr><td rowspan="6">τ (ms)</td><td>0.1</td><td>0.1</td><td>0.1</td><td>0.1</td><td>0.1</td><td>0.3</td><td>0.3</td><td>0.3</td><td>1</td><td>1</td></tr><tr><td>0.3</td><td>0.3</td><td>0.3</td><td>0.3</td><td>0.3</td><td>1</td><td>1</td><td>1</td><td>3</td><td>3</td></tr><tr><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>3</td><td>3</td><td>3</td><td>10</td><td>10</td></tr><tr><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>10</td><td>10</td><td>10</td><td>30</td><td>30</td></tr><tr><td>10</td><td>10</td><td>10</td><td>10</td><td>10</td><td>30</td><td>30</td><td>30</td><td>100</td><td>100</td></tr><tr><td>30</td><td>30</td><td>30</td><td>30</td><td>30</td><td>-</td><td>-</td><td>-</td><td>-</td><td>300</td></tr><tr><td></td><td>-</td><td>-</td><td>100</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr></tbody></table>

As can be seen in Figure 3, when the scanning time ΔTscan is 9–23 times of time constant τ, the output spectrum (red) is almost consistent with the input spectrum (black) and the noise can be efficiently rejected without causing obvious shift and distortion of the line shape. The scanning time ΔTscan (Table 1 in **blue bold**) as a function of the optimum time constant τ (Table 1 in **red bold**) is plotted in Figure 4. The mathematical relationship between the time constant τ and the scanning time ΔTscan is obtained from a linear fit ($r^2 = 99.25$):

$$\Delta T_{\text{scan}} = \Delta v_{\text{HWHM}}/v_{\text{sc}} \geq 14 \cdot \tau \quad (5)$$
<figure><img src="image_7.png" /><figcaption>**Figure 3.** Simulation of output spectra from the LP filter using different parameters (see Table 1), in comparison to the Hitran simulation of input spectrum of CH₄ near 1233.45 cm⁻¹.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Plot of $\Delta T_{\text{scan}}$, vs. time constant $\tau$, associated with a linear fit.</figcaption></figure>

Based on Equation (5), with respect to the given frequency scanning speed in the continuous mode and to the halfwidth of the spectrum to scan, the reasonable scanning time $\Delta T_{\text{scan}}$ should be equal to or longer than 14 times of the LIA time constant $\tau$ in order to efficiently reduce the noise without significant shift and distortion of the line shape.

## 2.3. *Experimental System Description*

An experimental investigation has been performed on measurements of heterodyne spectra under different $\tau$ and $\nu_{sc}$ combinations, based on a blackbody simulating sunlight emission.

The developed proof of concept (PoC) LHR prototype operating in the mid-infrared (mid-IR) near 8 µm is schematically shown in Figure 5. As a typical LHR, it includes a sunlight collection system associated with a sun-tracker, a local oscillator (tunable laser source), a photomixer to generate the beating signal at RF, an RF receiver, a signal modulation and demodulation module, as well as a frequency metrology module (in our case, direct absorption of the target molecule in a 12.5 cm long single-pass reference cell was used). To characterize and optimize the performance of the developed PoC LHR, a stable blackbody (BB) source was used as sunlight. Heterodyne measurements of CH₄ absorption spectra were performed in the laboratory using a 12.5 cm single-pass absorption cell.

Radiation from the blackbody (BB, MODELS 67030, Newport Inc. Florida, USA) at 1000 °C was collected and collimated with off-axis parabolic mirror (OAPM) 1 (f1 = 120 mm), then focused and directed by OAPM 2 (f2 = 75 mm) towards a 12.5 cm long single-pass cell (Figure 5a) that was filled with CH₄ mixture in air at atmospheric pressure. The BB beam was filtered by an optical filter in germanium in order to let only the transmission of BB radiation at wavelengths longer than 7 µm and shorter than 12 µm in our case to minimize source shot noise [18]. The beam emerging from the filter was amplitude modulated at 1 kHz using a chopper (Model MC1F10HP, Thorlabs, Inc. Newton, NJ, USA) located at the real focal point of OAPM2 and then directed to a 5%R/95%T beam splitter (BS, in CaF₂). The BB beam (at $\omega$) was superimposed with a LO laser beam (at $\omega_L$) from an external cavity quantum cascade laser (EC-QCL, TLS-41000-MHF, Daylight solutions Inc. San Diego, CA, USA) [19] which was beam size-expanded by two OAPMs (f3 = 12.7 mm and f4 = 38.1 mm), intensity-controlled and polarization-adjusted using a wideband IR polarizer (2–12 µm, PGC-5, Innpho. Inc. Verona, NJ, USA). The combined beams were focused onto a VIGO photomixer (PVI-4TE-10.6, VIGO System S.A. Ozarow Mazowiecki, Poland) by OAPM 5 (f5 = 12.7 mm). The beat note signal at $|\omega-\omega_L|$ in the RF domain, generated from the photomixer, entered the RF receiver through a 27–33 MHz band-pass (BP) filter (Model 15542, Mini-Circuits®) followed by 2-stage low-noise amplifiers (Model ZFL-500LN+, Mini-Circuits®) with a gain of 24 dB each. After that, the amplified beat signal was connected to a square-law Schottky diode (Model 8472B, Keysight Inc. Colorado
Springs, CO, USA) for measurement of the RF signal power that is proportional to the square of the input beat signal amplitude. The output signal from the square-law detector was demodulated with a lock-in amplifier (LIA, Model 5110, AMETEK Inc. Berwyn, PA, USA with time constant range: 1 ms–300 s.). A National Instruments data acquisition card (PCI 6251, NI Inc. Austin, TX, USA) was used to digitalize the output signal from the LIA via a custom-written LabVIEW program. The data were then transferred to a laptop for further data processing.

<figure><img src="image_9.png" /><figcaption>**Figure 5.** (a) Schematic of the developed LHR prototype. BB: Blackbody; EC-QCL: external cavity quantum cascade laser; OAPM: Off-axis parabolic mirror; BS: Beam splitter; BP Filter: Band-pass filter; DAQ: Data acquisition card. (b) Picture of the mid-IR LHR receiver (without sun-tracker part). Red dot line: sunlight beam path; orange dot line: blackbody beam path; yellow dot line: He-Ne laser beam path; blue dot line: EC-QCL beam path; green dot line: combined LO-BB beams to the photomixer.</figcaption></figure>
One part of the EC-QCL beam (transmitted by the BS) is directed to either a 12.5 cm long custom-made Fabry-Pérot etalon (with a free spectral range, FSR = 0.0333 cm⁻¹) for relative frequency calibration or to a single-pass cell filled with target gas sample for determination of the EC-QCL absolute frequency.

The LO laser frequency was scanned across an appropriate CH₄ absorption line (around 1233.46 cm⁻¹) by applying a sine-wave (with an amplitude of +3 VPP and an offset of +1.6 V_D (limited by the EC-QCL structure) from a function generator to the PZT (Piezoelectric ceramic transducer) element of the EC-QCL. Heterodyne spectra of CH₄ absorption of the incident BB radiation (Figure 6, black) were recorded using four combinations of different τ and vsc (Table 2) and fitted to a Lorentz profile (Figure 6, blue and red).

<table><caption>**Table 2.** Used time constants at a given scanning speed of the LO frequency vsc. (Note: the spectral half width ΔvHWHM of the used CH₄ line around 1233.46 cm⁻¹ is 0.0538 cm⁻¹ [20] at ambient pressure).</caption><tr><td>Scanning Frequency (mHz)</td><td>20</td><td>50</td><td>100</td><td>250</td></tr><tr><td>Scanning period T (s)</td><td>25</td><td>10</td><td>5</td><td>2</td></tr><tr><td>Scanning range Δv<sub>scan</sub> (cm<sup>-1</sup>)</td><td>0.787</td><td>0.787</td><td>0.787</td><td>0.787</td></tr><tr><td>halfwidth Δv<sub>HWHM</sub> (cm<sup>-1</sup>)</td><td>0.0538</td><td>0.0538</td><td>0.0538</td><td>0.0538</td></tr><tr><td>Scanning speed v<sub>sc</sub> (cm<sup>-1</sup>/s) = Δv<sub>scan</sub>/T</td><td>0.03148</td><td>0.0787</td><td>0.1574</td><td>0.3935</td></tr><tr><td>Scanning time ΔT<sub>scan</sub> (ms) = Δv<sub>HWHM</sub>/v<sub>sc</sub></td><td>1709</td><td>684</td><td>342</td><td>137</td></tr><tr><td rowspan="4">τ (ms)</td><td>100</td><td>30</td><td>10</td><td>1</td></tr><tr><td>300</td><td>100</td><td>30</td><td>10</td></tr><tr><td>1000</td><td>300</td><td>100</td><td>30</td></tr><tr><td>-</td><td>1000</td><td>300</td><td>-</td></tr></table>
<figure><img src="image_10.png" /><figcaption>**Figure 6.** BB-based heterodyne spectra of CH₄ absorption at different scanning speed ($v_{sc}$) (i.e., different given scanning time $\Delta T_{scan}$) combined with different time constants ($\tau$), respectively. Line shape fits were performed using the Lorentz profile (red or blue).</figcaption></figure>
# 3. Results

## 3.1. Experimental Results

Based on the analysis of the heterodyne absorption spectra shown in Figure 6 (black), signal-to-noise ratio (SNR), line width broadening ($\Delta v_b$), absorption depth, line shift ($\delta v$), and noise level ($1\sigma$) of these spectra were analyzed and summarized in Table 3, the optimal time constants ($\tau$) at given scanning speeds of LO frequency ($v_{sc}$) and scanning times ($\Delta T_{scan}$) were obtained from the current work, as shown in Table 3 in bold and Figure 6 in red.

<table><caption>**Table 3.** Analysis of BB-based heterodyne spectra of CH$_4$ absorption at different scanning speeds of the LO frequency ($v_{sc}$) (i.e., different given scanning time $\Delta T_{scan}$) with different time constants ($\tau$). $1\sigma$: Standard deviation.</caption><thead><tr><th>$v_{sc}$ (cm<sup>-1</sup>/s)</th><th>$\Delta T_{scan}$ (ms)</th><th>$\tau$ (ms)</th><th>SNR</th><th>$\Delta v_b$ (cm<sup>-1</sup>)</th><th>Absorption Depth</th><th>$\delta v$ (cm<sup>-1</sup>)</th><th>Fit Residual ($1\sigma$)</th></tr></thead><tbody><tr><td rowspan="3">0.03148</td><td rowspan="3">1709</td><td>100</td><td>5</td><td>0.0545</td><td>0.0126</td><td>0.0032</td><td>0.0026</td></tr><tr><td><strong>300</strong></td><td><strong>8</strong></td><td><strong>0.0656</strong></td><td><strong>0.0114</strong></td><td><strong>0.0095</strong></td><td><strong>0.0015</strong></td></tr><tr><td>1000</td><td>8</td><td>0.0711</td><td>0.0062</td><td>0.0325</td><td>0.0008</td></tr><tr><td rowspan="3">0.0787</td><td rowspan="3">684</td><td>30</td><td>5</td><td>0.0548</td><td>0.01258</td><td>0.0031</td><td>0.00230</td></tr><tr><td><strong>100</strong></td><td><strong>6</strong></td><td><strong>0.0640</strong></td><td><strong>0.01244</strong></td><td><strong>0.0081</strong></td><td><strong>0.00210</strong></td></tr><tr><td>300</td><td>6</td><td>0.0722</td><td>0.01154</td><td>0.0243</td><td>0.00180</td></tr><tr><td rowspan="5">0.1574</td><td rowspan="5">342</td><td>1000</td><td>6</td><td>0.1222</td><td>0.00570</td><td>0.0821</td><td>0.00089</td></tr><tr><td>10</td><td>5</td><td>0.0550</td><td>0.01251</td><td>0.0016</td><td>0.00230</td></tr><tr><td>30</td><td>6</td><td>0.0627</td><td>0.01204</td><td>0.0048</td><td>0.00209</td></tr><tr><td><strong>100</strong></td><td><strong>7</strong></td><td><strong>0.0672</strong></td><td><strong>0.01101</strong></td><td><strong>0.0181</strong></td><td><strong>0.00156</strong></td></tr><tr><td>300</td><td>7</td><td>0.0812</td><td>0.00824</td><td>0.0490</td><td>0.00115</td></tr><tr><td rowspan="3">0.3935</td><td rowspan="3">137</td><td><strong>1</strong></td><td><strong>7</strong></td><td><strong>0.0541</strong></td><td><strong>0.01261</strong></td><td><strong>0.0004</strong></td><td><strong>0.00174</strong></td></tr><tr><td>10</td><td>9</td><td>0.0599</td><td>0.01181</td><td>0.0041</td><td>0.00129</td></tr><tr><td>30</td><td>8</td><td>0.0733</td><td>0.00736</td><td>0.0120</td><td>0.00096</td></tr></tbody></table>

As can be seen, the selection of the LIA time constant $\tau$ should be well matched to the scanning speed of the LO frequency as well as the scanning time. When the time constant is too small, the noise could not be efficiently removed. While a longer (than the optimal one) time constant would result in reduced absorption depth, a shift in absorption center frequency and distortion of absorption line shape.

Comparing the above four optimal combinations for heterodyne measurements, it is seen that the combination of the scanning speed of the LO frequency of 0.3935 cm⁻¹/s (the scanning time of 137 ms) and the LIA time constant of 1 ms could provide a relatively low noise level and the highest absorption depth without an obvious shift in absorption center frequency and distortion of the line shape. This relatively fast scanning time $\Delta T_{scan}$ of 137 ms makes it easier to acquire a fast extraction of the high-precision line shape of LHR spectrum. This combination was therefore selected for our further LHR measurements.

## 3.2. Field Measurement

Using the optimal combination of the scanning time of 137 ms and the time constant of 1 ms, the developed LHR instrument was deployed to a field campaign in Dunkerque (51.05°N/2.34°E) for ground-based remote sensing of greenhouse gases in the atmospheric column.

The LO laser frequency was scanned around 1262.00 cm⁻¹ to probe absorption lines of N₂O (at 1261.9871 cm⁻¹) and CH₄ (at 1261.6450 and 1262.2249 cm⁻¹). An RF filter of 225–270 MHz pass-band was used to remove 1/f noise as well as other extraneous noises
outside of the bandwidth, and this pass-band resulted in a ~0.003 cm⁻¹ double-sideband spectral resolution for the remote measurements.

The measured CH₄ and N₂O heterodyne spectrum in the atmospheric column is presented in Figure 7 (blue curve), accompanied by a referenced TCCON (Total Carbon Column Observing Network) FT-IR (Fourier-transform infrared) spectrum (Figure 7 (red curve)) measured using a Bruker IFS 125HR FTS (resolution: 0.02 cm⁻¹) of the TCCON observation network [21] on the QualAir platform of University of Paris 6. As can be seen, the LHR spectrum, with higher SNR, is in good agreement with the referenced FT-IR spectrum. Due to interference of nearby water vapor absorption that changes with time and space [22], wings on either side of the LHR spectrum are sometimes slightly different from those of the referenced FT-IR spectrum. The overall spectral line shapes of both two spectra are well consistent. In comparison to FTS, LHR is more suitable for field deployment owing to its room-temperature operation capacity, compact size, and lightweight (0.6 × 0.9 m²/25 kg for the present LHR vs. 2 × 4 m²/700 kg for a traditional FTS).

<figure><img src="image_11.png" /><figcaption>**Figure 7.** Atmospheric absorption spectrum of CH₄ and N₂O recorded with a Bruker IFS 125HR Fourier-transform spectrometer (FTS) in Paris (in red), with our LHR in Dunkerque (in blue).</figcaption></figure>

**4. Discussion**

Study and experimental investigation of the impact of LIA time constant on observed LHR line shape in continuous tuning laser frequency mode are carried out. Based on the present work, the scanning time $\Delta T_{\text{scan}}$ of laser frequency across the halfwidth of the absorption line $\Delta\nu_{\text{HWHM}}$, the scanning speed of the LO laser frequency $\nu_{\text{sc}}$ and the LIA time constant $\tau$ should be appropriately matched such that $\Delta T_{\text{scan}} = \Delta\nu_{\text{HWHM}}/\nu_{\text{sc}} \ge 14\tau$ in order to obtain high-precision heterodyne spectrum (with minimum line shift and line shape distortion). This relationship was validated by laser heterodyne measurements of CH₄ absorption spectra in the laboratory using four combinations of different time constants $\tau$ and scanning time $\Delta T_{\text{scan}}$ in continuous tuning mode. The optimal combination of the scanning time of 137 ms and the time constant of 1 ms, which can efficiently reduce the noise without introducing obvious line shift and distortion, was therefore selected for ground-based measurements of CH₄ and N₂O in the atmospheric column. The extracted LHR spectrum is in good agreement with a referenced TCCON FT-IR spectrum.

**Author Contributions:** Conceptualization, F.S., W.C. and Z.C.; methodology, F.S., W.C. and T.T.; software, F.S., G.W., T.T. and Z.X.; validation, F.S. and W.C.; formal analysis, F.S.; investigation, F.S.; resources, F.S.; data curation, F.S.; writing—original draft preparation, F.S.; writing—review and editing, F.S.; visualization, F.S.; supervision, W.C. and X.G.; project administration, F.S. and W.C.; funding acquisition, F.S. and W.C. All authors have read and agreed to the published version of the manuscript.
**Funding:** This research was funded by the Talent fund of Hefei University under Grant, grant number 22040421002, University Natural Science Research Program of Anhui Province, grant number KJ2021A0979, Key disciplines of Hefei University: Information and Communication Engineering, grant number 2018xk03, French national research agency (ANR) under the CaPPA, grant number ANR-10-LABX-005 and the CPER ECRIN contracts.

**Data Availability Statement:** Not applicable.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Lu, X.; Huang, Y.; Wu, P.; Liu, D.; Ma, H.; Wang, G.; Cao, Z. Distributed Feedback Interband Cascade Laser Based Laser Heterodyne Radiometer for Column Density of HDO and CH4 Measurements at Dunhuang, Northwest of China. *Remote Sens.* **2022**, *14*, 1489. [CrossRef]

2. Shen, F.J.; Wang, G.X.; Wang, J.J.; Tan, T.; Wang, G.S.; Jeseck, P.; Té, Y.V.; Gao, X.M.; Chen, W.D. Transportable mid-infrared laser heterodyne radiometer operating in the shot-noise dominated regime. *Opt. Lett.* **2021**, *46*, 3171–3174. [CrossRef]

3. Lu, X.J.; Cao, Z.S.; Tan, T.; Huang, Y.B.; Gao, X.M.; Rao, R.Z. Instrument line shape function of laser heterodyne spectrometer. *Acta Phys. Sin.* **2019**, *68*, 064208. [CrossRef]

4. Tsai, T.; Rose, R.A.; Weidmann, D.; Wysocki, G. Atmospheric vertical profiles of O₃, N₂O, CH₄, CCl₂F₂, and H₂O retrieved from external-cavity quantum-cascade laser heterodyne radiometer measurements. *Appl. Opt.* **2012**, *51*, 8779–8792. [CrossRef] [PubMed]

5. Wang, J.J.; Wang, G.S.; Tan, T.; Zhu, G.; Sun, C.Y.; Cao, Z.S.; Chen, W.D.; Gao, X.M. Mid-infrared laser heterodyne radiometer (LHR) based on a 3.53 µm room-temperature interband cascade laser. *Opt. Exp.* **2019**, 27, 9610–9619. [CrossRef] [PubMed]

6. Brasseur, G.P.; Jacob, D.J. *Modeling of Atmospheric Chemistry*; Cambridge University Press: Cambridge, UK, 2017; pp. 84–204.

7. Ehhalt, D.H. Concentrations and distributions of atmospheric trace gases. *Ber. Bunsenges. Phys. Chem.* **1992**, *96*, 229–240. [CrossRef]

8. Tsai, T. External cavity quantum cascade lasers for spectroscopic applications. Doctoral Dissertation, Princeton University, Princeton, NJ, USA, 2012.

9. Mahieu, E.; Duchatelet, P.; Demoulin, P.; Walker, K.A.; Dupuy, E.; Froidevaux, L.; Randall, C.; Catoire, V.; Strong, K.; Boone, C.D.; et al. Validation of ACE-FTS v2.2 measurements of HCl, HF, CCl₃F and CCl₂F₂ using space-, balloon- and ground-based instrument observations. *Atmos. Chem. Phys.* **2008**, *8*, 6199–6221. [CrossRef]

10. Marquart, S.; Ponater, M.; Mager, F.; Sausen, R. Future development of contrail cover, optical depth, and radiative forcing: Impacts of increasing air traffic and climate change. *J. Clim.* **2003**, *16*, 2890–2904. [CrossRef]

11. Scofield, J.H. Frequency domain description of a lock-in amplifier. *Am. J. Phys.* **1994**, *62*, 129–133. [CrossRef]

12. Wang, G.X.; Rebou, S.; Choque, J.; Fertein, E.; Chen, W.D. Circular regression-based dual-phase lock-in amplifier for coherent recovery of weak signal. *Sensors* **2017**, *17*, 2615. [CrossRef] [PubMed]

13. SRS. About Lock-in Amplifiers. Available online: http://www.thinksrs.com/downloads/PDFs/ApplicationNotes/AboutLIAs.pdf (accessed on 16 May 2015).

14. Rohart, F.; Mejri, S.; Sow, P.; Tokunaga, S.K.; Chardonnet, C.; Darquié, B.; Dinesan, H.; Fasci, E.; Castrillo, A.; Gianfrani, L. Absorption line shape recovery beyond the detection bandwidth limit: Application to the precision spectroscopic measurement of the Boltzmann constant. *Phys. Rev. A Am. Phys. Soc.* **2014**, *90*, 042506–042630. [CrossRef]

15. Pelant, I.; Valenta, J. Luminescence Spectroscopy of Semiconductors. Oxford Scholarship Online. 2012. Available online: https://oxford.universitypressscholarship.com/view/10.1093/acprof:oso/9780199588336.001.0001/acprof-9780199588336 (accessed on 1 May 2012).

16. Brodersen, S. Noise filter effects in recording spectroscopes. *J. Opt. Soc. Am.* **1953**, *43*, 1216–1220. [CrossRef]

17. Demtröder, W. *Laser Spectroscopy: Basic Concepts and Instrumentation*; Springer: Berlin/Heidelberg, Germany, 1996.

18. Kostiuk, T.; Mumma, M.J. Remote sensing by IR heterodyne spectroscopy. *Appl. Opt.* **1983**, *22*, 2644–2654. [CrossRef] [PubMed]

19. Shen, F.J.; Akil, J.; Wang, G.X.; Poupin, C.; Cousin, R.; Siffert, S.; Fertein, E.; Ba, T.; Chen, W.D. Real-time monitoring of N₂O production in a catalytic reaction process using mid-infrared quantum cascade laser. *J. Quant. Spectrosc. Radiat. Transf.* **2018**, *221*, 1–7. [CrossRef]

20. Gordon, I.E.; Rothman, L.S.; Hill, C.; Kochanov, R.V.; Tan, Y.; Bernath, P.F.; Birk, M.; Boudon, V.; Campargue, A.; Chance, K.V.; et al. The HITRAN2016 molecular spectroscopic database. *J. Quant. Spectrosc. Rad. Transf.* **2017**, *203*, 3–69. [CrossRef]

21. Té, Y.; Jeseck, P.; Franco, B.; Mahieu, E.; Jones, N.P.; Paton-Walsh, C.; Griffith, D.W.T.; Buchholz, R.R.; Hadji-Lazaro, J.; Hurtmans, D.; et al. Seasonal variability of surface and column carbon monoxide over the megacity Paris, high-altitude Jungfraujoch and Southern Hemispheric Wollongong stations. *Atmos. Chem. Phys.* **2016**, *16*, 10911–10925. [CrossRef]

22. Kämpfer, N. *Monitoring Atmospheric Water Vapour: Ground-Based Remote Sensing and In Situ Methods*; Springer: New York, NY, USA, 2013; Volume 13, pp. 1–7.