ExoMars 痕量气体轨道飞行器:One Martian Year of Science ExoMars 痕量气体轨道飞行器:一个火星年的科学
要点:
介绍了 ExoMars 痕量气体轨道飞行器上大气化学套件仪器对重力波的观测结果
观察到的波活动、势能、动量通量和波阻力的全球分布与模型预测非常吻合 The observed global distribution of wave activity, potential energy, momentum flux, and wave resistance is in good agreement with model predictions
Figure 1. Spectroscopy of CO_(2)\mathrm{CO}_{2} and H_(2)O\mathrm{H}_{2} \mathrm{O} absorption in the diffraction order #223 of ACS-MIR (panel A) and an example of retrieved atmospheric temperature and density vertical profiles (panel B). (a) Transmission spectra measured at tangent altitudes of 150, 90 and 29 km (red dots) on a background of the best-fitted models (black solid lines); (b) Vertical profiles of temperature (left) and atmospheric number density (right) derived from the MCD (black dashed line), from ACS-MIR (red dots), and from ACS-NIR (Fedorova et al., 2020) (green dots). Error bars for the temperature values express 1- sigma\sigma uncertainties of the retrievals. 图 1. ACS-MIR 的衍射级数#223 中 CO_(2)\mathrm{CO}_{2} 和 H_(2)O\mathrm{H}_{2} \mathrm{O} 的光谱吸收(面板 A)以及检索的大气温度和密度垂直剖面的示例(面板 B)。(a) 在 150、90 和 29 公里的切线高度测得的传输光谱(红点),背景为最佳拟合模型(黑色实线);(b) 从 MCD(黑色虚线)、ACS-MIR(红点)和 ACS-NIR(Fedorova 等,2020)(绿色点)得出的温度(左)和大气数密度(右)的垂直剖面。温度值的误差条表示检索的 1- sigma\sigma 不确定性。
covers 0.15-0.3 mum0.15-0.3 \mu \mathrm{~m} range. The instrument’s resolving power is lambda//Delta lambda∼25,000\lambda / \Delta \lambda \sim 25,000 and the signal-to-noise ratio varies between 1,000 and 10,000. The vertical resolution of MIR depends on the integration time ( ∼2s\sim 2 \mathrm{~s} per image) and ranges from 0.5 to 2.5 km . The transmission is obtained by division of the solar spectrum passed through the atmosphere to the reference one, which is measured above the altitude of 200 km , where the absorption by the atmosphere is negligible. 覆盖 0.15-0.3 mum0.15-0.3 \mu \mathrm{~m} 范围。仪器的分辨率为 lambda//Delta lambda∼25,000\lambda / \Delta \lambda \sim 25,000 ,信噪比在 1,000 到 10,000 之间变化。MIR 的垂直分辨率取决于积分时间( ∼2s\sim 2 \mathrm{~s} 每幅图像)并在 0.5 到 2.5 公里之间变化。传输是通过将穿过大气的太阳光谱与在 200 公里以上测量的参考光谱进行除法获得的,在该高度,大气的吸收可以忽略不计。
In this study, we use the 2.66-2.68 mum2.66-2.68 \mu \mathrm{~m} portion of the spectrum from the grating position #4, the echelle diffraction order #223, which includes a wing of the 2.7 mumCO_(2)2.7 \mu \mathrm{~m} \mathrm{CO}_{2} absorption band (Figure 1a). 在本研究中,我们使用来自光栅位置#4 的 2.66-2.68 mum2.66-2.68 \mu \mathrm{~m} 光谱部分,埃谢尔衍射级数#223,其中包括 2.7 mumCO_(2)2.7 \mu \mathrm{~m} \mathrm{CO}_{2} 吸收带的一部分(图 1a)。
Strong absorption lines of CO_(2)\mathrm{CO}_{2} allow for retrieving temperature and density in the Martian atmosphere with good sensitivity. 强吸收线的 CO_(2)\mathrm{CO}_{2} 允许以良好的灵敏度获取火星大气中的温度和密度。
3. Methods 3. 方法
3.1. Retrieval of Temperature Profiles 3.1. 温度剖面的检索
The retrieval scheme consists of several iterations. On the first step, we retrieve temperature and pressure from the rotational structure of CO_(2)\mathrm{CO}_{2} absorption bands in spectral intervals without H_(2)O\mathrm{H}_{2} \mathrm{O} lines (see Figure 1a). A priori altitude profiles of T(z)T(z) and p(z)p(z) as well as one of the CO_(2)VMR\mathrm{CO}_{2} \mathrm{VMR}, are taken from the Mars Climate Database (MCD) for a specified occultation in MY34 (Millour et al., 2018). On the second step, we simultaneously retrieve temperature and CO_(2)\mathrm{CO}_{2} concentration, while the pressure profile is kept constant assuming the hydrostatic equilibrium p_(hyd)(z)=p_(0)(z_(0))exp[-int_(z_(0))^(z)(g(z^('))M(z^(')))/(RT(z^(')))dz^(')]p_{h y d}(z)=p_{0}\left(z_{0}\right) \exp \left[-\int_{z_{0}}^{z} \frac{g\left(z^{\prime}\right) M\left(z^{\prime}\right)}{R T\left(z^{\prime}\right)} d z^{\prime}\right], where gg is the acceleration of gravity, MM is the atmospheric molar mass and RR is the gas constant. The reference pressure p_(0)p_{0} is chosen at an altitude z_(0)z_{0}, usually around 30-50km30-50 \mathrm{~km}, where uncertainties of the fitting are smallest. We repeat the second step 5-7 times until the profiles reach convergence. In each iteration, we apply the Tikhonov regularization (Tikhonov & Arsenin, 1977) for the temperature and concentration altitude profiles with a smoothing coefficient less than 5 km . It defines the shortest wavelength to 5-6km5-6 \mathrm{~km} when analyzing vertical wavy structures. The third step focuses only on CO_(2)\mathrm{CO}_{2} and H_(2)O\mathrm{H}_{2} \mathrm{O} concentration retrievals over the entire wavenumber range in 检索方案由几个迭代组成。在第一步中,我们从 CO_(2)\mathrm{CO}_{2} 吸收带的旋转结构中检索温度和压力,光谱区间内没有 H_(2)O\mathrm{H}_{2} \mathrm{O} 线(见图 1a)。 T(z)T(z) 和 p(z)p(z) 的先验高度剖面以及 CO_(2)VMR\mathrm{CO}_{2} \mathrm{VMR} 之一,均取自火星气候数据库(MCD),用于指定的 MY34(Millour 等,2018)掩蔽。在第二步中,我们同时检索温度和 CO_(2)\mathrm{CO}_{2} 浓度,同时保持压力剖面不变,假设静水平衡 p_(hyd)(z)=p_(0)(z_(0))exp[-int_(z_(0))^(z)(g(z^('))M(z^(')))/(RT(z^(')))dz^(')]p_{h y d}(z)=p_{0}\left(z_{0}\right) \exp \left[-\int_{z_{0}}^{z} \frac{g\left(z^{\prime}\right) M\left(z^{\prime}\right)}{R T\left(z^{\prime}\right)} d z^{\prime}\right] ,其中 gg 是重力加速度, MM 是大气摩尔质量, RR 是气体常数。参考压力 p_(0)p_{0} 选择在高度 z_(0)z_{0} ,通常在 30-50km30-50 \mathrm{~km} 附近,此时拟合的不确定性最小。我们重复第二步 5-7 次,直到剖面达到收敛。在每次迭代中,我们对温度和浓度高度剖面应用 Tikhonov 正则化(Tikhonov & Arsenin,1977),平滑系数小于 5 公里。它定义了在分析垂直波状结构时的最短波长 5-6km5-6 \mathrm{~km} 。第三步仅关注整个波数范围内的 CO_(2)\mathrm{CO}_{2} 和 H_(2)O\mathrm{H}_{2} \mathrm{O} 浓度检索
order #223 (Figure 1a) using the p(z)p(z) and T(z)T(z) profiles already found. This step is not a subject of the present paper. 订单 #223 (图 1a) 使用已找到的 p(z)p(z) 和 T(z)T(z) 配置文件。此步骤不是本论文的主题。
A similar fitting procedure, including the hydrostatic approximation, has been used in the work by Fedorova et al. (2020) (proprietary code) and Alday et al. (2019) (the NEMESIS code (Irwin et al., 2008),) in their retrievals of temperature and pressure from the ACS data. We validated our atmospheric temperature and number density profiles with simultaneous and collocated occultation measurements by ACS-NIR (Fedorova et al., 2020). An example comparison is presented in Figure 1b. A weaker CO_(2)\mathrm{CO}_{2} absorption band at 1.58 mum\mu \mathrm{m} measured by NIR allows for detection up to 110-120km110-120 \mathrm{~km}, or the density of ∼10^(12)cm^(-3)\sim 10^{12} \mathrm{~cm}^{-3}, while the band at 2.7 mum2.7 \mu \mathrm{~m} observed by MIR is measurable up to 160-170km160-170 \mathrm{~km}, or ∼10^(9)cm^(-3)\sim 10^{9} \mathrm{~cm}^{-3}. The lowermost altitude of the temperature profile retrieval is conditioned by the aerosol opacity and by the saturation of the CO_(2)\mathrm{CO}_{2} absorption lines. 一种类似的拟合程序,包括静水近似,已在 Fedorova 等人(2020 年)(专有代码)和 Alday 等人(2019 年)(NEMESIS 代码(Irwin 等,2008 年))的工作中使用,用于从 ACS 数据中提取温度和压力。我们通过 ACS-NIR 的同时和重叠的掩星测量(Fedorova 等,2020 年)验证了我们的气温和数密度剖面。图 1b 中展示了一个比较示例。NIR 测得的 1.58 的较弱吸收带允许检测到高达 110-120km110-120 \mathrm{~km} ,或 ∼10^(12)cm^(-3)\sim 10^{12} \mathrm{~cm}^{-3} 的密度,而 MIR 观察到的 2.7 mum2.7 \mu \mathrm{~m} 的带可测量高达 160-170km160-170 \mathrm{~km} ,或 ∼10^(9)cm^(-3)\sim 10^{9} \mathrm{~cm}^{-3} 。温度剖面提取的最低高度受气溶胶不透明度和 CO_(2)\mathrm{CO}_{2} 吸收线的饱和度的限制。
Each temperature value in a vertical profile was retrieved by fitting a modeled transmission spectrum J_("mod ")J_{\text {mod }} to the measured one J_("mes ")J_{\text {mes }} at a specified altitude. We model the spectra by the Beer-Lambert law 每个垂直剖面中的温度值是通过将模型传输光谱拟合到指定高度的测量光谱中获取的。我们通过比尔-朗伯定律对光谱进行建模。
where n(z)n(z) are gaseous concentrations, sigma(T,p)\sigma(T, p) are absorption cross-sections of CO_(2)\mathrm{CO}_{2} and H_(2)O\mathrm{H}_{2} \mathrm{O} correspondingly for specific temperature T(z)T(z) and pressure p(z)p(z) at an altitude zz, and tau_(a)\tau_{a} is aerosol slant opacity. A transfer between the linear [cm^(-2)]\left[\mathrm{cm}^{-2}\right] and the volume [cm^(-3)]\left[\mathrm{cm}^{-3}\right] concentrations is performed using the well-known “onion-peeling” method with the numeric integration over all altitude layers z_(i)z_{i} above the ii th one. Molecular cross-sections are calculated line-by-line on a basis of the HITRAN2016 database (Gordon et al., 2017) considering pressure-broadening coefficients of the H_(2)O\mathrm{H}_{2} \mathrm{O} lines suitable for a CO_(2)\mathrm{CO}_{2}-rich atmosphere ( Ga mache et al., 2016) and self-broadening in the case of CO_(2)\mathrm{CO}_{2}. Then we convolve the modeled spectrum by the previously determined instrument line shape (ILS) using wavenumber calibrations (see details in Alday et al., 2019). The fitting procedure is conducted by minimizing the “chi-square” function chi^(2)=sum_(i)A^(2)(v_(i))\chi^{2}=\sum_{i} A^{2}\left(v_{i}\right), A(v_(i))=[J_("mod ")(v_(i))-J_("mes ")(v_(i))]//delta JA\left(v_{i}\right)=\left[J_{\text {mod }}\left(v_{i}\right)-J_{\text {mes }}\left(v_{i}\right)\right] / \delta J, where delta J\delta J are transmittance uncertainties, and the sum is taken over all considered spectral points (pixels). Our optimization algorithm to search for the chi^(2)\chi^{2} minimum is based on partial derivatives of the Jacobian matrix del A//del X\partial A / \partial X (Marquardt, 1963), where XX is a vector of free parameters, that is, temperature, CO_(2)\mathrm{CO}_{2} concentration, H_(2)O\mathrm{H}_{2} \mathrm{O} mixing ratio, and aerosol slant opacity. Here, a significant contribution to the Jacobian comes from the rotational absorption lines, which are strongly sensitive to the temperature variability in the spectral range of interest. 其中 n(z)n(z) 是气体浓度, sigma(T,p)\sigma(T, p) 是 CO_(2)\mathrm{CO}_{2} 的吸收截面, H_(2)O\mathrm{H}_{2} \mathrm{O} 对应于特定温度 T(z)T(z) 和压力 p(z)p(z) 在一个高度 zz ,并且 tau_(a)\tau_{a} 是气溶胶斜率不透明度。通过使用著名的“洋葱剥皮”方法,利用数值积分在所有高度层 [cm^(-2)]\left[\mathrm{cm}^{-2}\right] 之上进行线性 [cm^(-3)]\left[\mathrm{cm}^{-3}\right] 和体积 z_(i)z_{i} 浓度之间的转换 ii 。分子截面是基于 HITRAN2016 数据库(Gordon 等,2017)逐行计算的,考虑了适合 H_(2)O\mathrm{H}_{2} \mathrm{O} 的 CO_(2)\mathrm{CO}_{2} -富气氛的压力展宽系数(Ga mache 等,2016)以及在 CO_(2)\mathrm{CO}_{2} 的情况下的自展宽。然后,我们通过先前确定的仪器线形(ILS)对建模的光谱进行卷积,使用波数校准(详见 Alday 等,2019)。 拟合过程是通过最小化“卡方”函数进行的,其中是透射不确定性,求和是在所有考虑的光谱点(像素)上进行的。我们搜索最小值的优化算法基于雅可比矩阵的偏导数,其中是自由参数的向量,即温度、浓度、混合比和气溶胶斜率不透明度。在这里,雅可比矩阵的一个重要贡献来自于旋转吸收线,这些吸收线对感兴趣光谱范围内的温度变化非常敏感。
3.2. Derivation of Wave Disturbances 3.2. 波动扰动的推导
Gravity wave-induced perturbations of temperature T^(')T^{\prime} are sought by separating the mean, or background profile bar(T)(z)\bar{T}(z) from the measured one T(z)T(z) : 重力波引起的温度扰动 T^(')T^{\prime} 通过将平均值或背景剖面 bar(T)(z)\bar{T}(z) 与测量值分离来寻找 T(z)T(z) :
T^(')=T- bar(T),T^{\prime}=T-\bar{T},
where the bar denotes an appropriate averaging. Generally, it implies averaging over wave phases, or spatial and temporal scales that are larger than the periods and wavelengths of contributing GW harmonics. In the case of almost instantaneous (with respect to the periods of GWs) occultation profiles, only separation in vertical scales is possible. 其中的条形表示适当的平均。一般来说,这意味着在波相位或空间和时间尺度上进行平均,这些尺度大于贡献的重力波谐波的周期和波长。在几乎瞬时(相对于重力波的周期)遮蔽轮廓的情况下,仅可能在垂直尺度上进行分离。
John and Kumar (2013) and Ehard et al. (2015) reviewed several common methods of the partition of measured temperature and/or density profiles into the “mean” and wave components. They work well if a clear separation in vertical wavelengths does exist between GWs and large-scale motions belonging to the background. This is not always the case in the Martian atmosphere, because vertical scales of disturbances associated with tides, planetary waves, and other motions may overlap with those due to GWs. It is desirable to retain the former in the background, but one still has to set a vertical scale Lambda_(z)\Lambda_{z} that separates GWs from the larger-scale features. In the following, we assumed Lambda_(z)=30km\Lambda_{z}=30 \mathrm{~km}. This value may lead to an overestimation of the retrieved wave activity by including non-GW perturbations, but at least no large-scale GW components are missed. Concerning the short-wavelength part of the spectrum, the limited vertical resolution favors detection of larger-scale waves, leaving out harmonics with smaller scales unobserved. Alexander (1998) 约翰和库马尔(2013)以及埃哈德等人(2015)回顾了几种将测量的温度和/或密度剖面划分为“均值”和波动成分的常见方法。如果存在 GWs 与属于背景的大尺度运动之间的清晰垂直波长分离,它们的效果很好。但在火星大气中,这种情况并不总是成立,因为与潮汐、行星波和其他运动相关的扰动的垂直尺度可能与 GWs 的尺度重叠。希望在背景中保留前者,但仍然需要设定一个垂直尺度 Lambda_(z)\Lambda_{z} ,以将 GWs 与更大尺度特征分开。在接下来的部分中,我们假设 Lambda_(z)=30km\Lambda_{z}=30 \mathrm{~km} 。这个值可能会导致通过包括非 GW 扰动而高估检索到的波动活动,但至少不会遗漏大型 GW 成分。关于光谱的短波长部分,有限的垂直分辨率有利于检测更大尺度的波动,未观察到较小尺度的谐波。亚历山大(1998)
Figure 2. Separation of the observed temperature into the mean and wave components for two characteristic profiles: dominated by large vertical-scale (orbit 2892n1, upper row) and small-scale disturbances (orbit 3251n1, lower row). Left column is for the mean temperature bar(T)(z)\bar{T}(z), the right one is for the relative perturbations T^(')(z)// bar(T)(z)T^{\prime}(z) / \bar{T}(z) (in percent). The legend describes the applied methods. Red dashed lines correspond to the Fourier decomposition, green and yellow lines are for the sliding polynomial fit with 2-km2-\mathrm{km} and 11-km shift steps, correspondingly, and the blue lines are for the 7-th order polynomial fit. The observed temperature profiles are given with the solid black lines. Shaded area denotes the uncertainty of the measurements. 图 2. 将观测温度分离为均值和波动分量的两个特征剖面:由大垂直尺度(轨道 2892n1,上排)和小尺度扰动(轨道 3251n1,下排)主导。左列为均温 bar(T)(z)\bar{T}(z) ,右列为相对扰动 T^(')(z)// bar(T)(z)T^{\prime}(z) / \bar{T}(z) (以百分比表示)。图例描述了所应用的方法。红色虚线对应于傅里叶分解,绿色和黄色线分别为滑动多项式拟合,步长为 2-km2-\mathrm{km} 和 11 公里,蓝色线为 7 阶多项式拟合。观测温度剖面用实线黑色表示。阴影区域表示测量的不确定性。
has quantified this “observational filter” and pointed out that some large-scale harmonics refracted by the mean wind beyond the lowest resolution may be missing in observations as well. 已经量化了这种“观察滤波器”,并指出一些被平均风折射的超出最低分辨率的大规模谐波在观测中可能也会缺失。
We explore three methods: spectral filtering, sliding least-square polynomial fit and high-order polynomial fit. The former two have been discussed in relation to lidar and space-based measurements in the atmosphere of Earth (Ehard et al., 2015; John & Kumar, 2013, and the references therein), while the latter was applied to profiles obtained in the terrestrial (e.g., Spiga et al., 2008) and Martian atmosphere (Jesch et al., 2019; Terada et al., 2017; Yiğit et al., 2015). Since the ACS data are distributed irregularly over the altitude, they were first interpolated (oversampled) to an evenly spaced 500-m500-\mathrm{m} grid. We used only the temperature data with errors <= 10K\leq 10 \mathrm{~K}. Spectral filtering was performed using Fourier decomposition within sliding 60 - km intervals ( +-30km\pm 30 \mathrm{~km} around each point), and zero-order Fourier coefficients were used to calculate the background temperature. The examples are shown in Figure 2 for two characteristic profiles T(z)T(z). They visibly differ: the profile in Figure 2a (orbit 2892n1) contains large-scale disturbances, while the one in Figure 2c (orbit 3251n1) comprises mostly smaller-scale fluctuations. This method yields smooth mean 我们探讨了三种方法:谱滤波、滑动最小二乘多项式拟合和高阶多项式拟合。前两者已在与激光雷达和地球大气中的空间测量相关的研究中讨论过(Ehard et al., 2015; John & Kumar, 2013,以及其中的参考文献),而后者则应用于在地球(例如,Spiga et al., 2008)和火星大气中获得的剖面(Jesch et al., 2019; Terada et al., 2017; Yiğit et al., 2015)。由于 ACS 数据在高度上分布不规则,因此首先对其进行了插值(过采样),以形成均匀间隔的 500-m500-\mathrm{m} 网格。我们仅使用了误差为 <= 10K\leq 10 \mathrm{~K} 的温度数据。谱滤波是通过在滑动 60 公里区间内进行傅里叶分解来执行的( +-30km\pm 30 \mathrm{~km} 围绕每个点),并使用零阶傅里叶系数计算背景温度。图 2 展示了两个特征剖面的示例 T(z)T(z) 。它们明显不同:图 2a(轨道 2892n1)中的剖面包含大规模扰动,而图 2c(轨道 3251n1)中的剖面主要由小规模波动组成。该方法产生平滑的平均值。
temperature profiles and, as a result, large deviations from the mean (Figures 2 b and 2 d ). This is in particular obvious below 60 km and in the upper part of the domain (panels b and d). 温度剖面,因此,与均值存在较大偏差(图 2b 和 2d)。这在 60 公里以下和区域的上部(面板 b 和 d)尤为明显。
For the sliding polynomial fit, we used a procedure described in the work of Whiteway and Carswell (1995). The background profiles are obtained by fitting cubic polynomials within the 60-km60-\mathrm{km} sliding intervals. Ob-\mathrm{Ob}- servational errors were used as weights that assign a significance to the measurements at each altitude. At first, the intervals were shifted up from the bottom to top by a certain distance (shown in Figures 2a and 2c for 2 and 11 km ), and then the procedure was repeated for the downward shifts starting from the top. The overlapping values of fits from each range were then averaged. Thus obtained profiles were then smoothed using a moving average. At the bottom of the profiles, we had to decrease the width of the sliding windows due to large spurious variations in fitted polynomials and in order to make most of the observational data. The upper and lower 4 km of thus obtained profiles have to be excluded anyway, because of the poor behavior of fitting polynomials, which cannot be averaged with counterparts from other sliding windows. This method occasionally produces disturbances oscillating not around zero. To correct for these numerical biases, we perform detrending by applying the Theil-Sen estimator (Sen, 1968; Theil, 1950) and fitting a linear function to the perturbation profile. The Theil-Sen estimator is a robust method, which is used for determining the linear regression taking the median of the slopes of all lines that can be drawn through the given data set. The linear function is then subtracted from the profile to obtain the corrected temperature. 对于滑动多项式拟合,我们使用了 Whiteway 和 Carswell(1995)工作中描述的程序。背景轮廓是通过在滑动区间内拟合三次多项式获得的。观测误差被用作权重,为每个高度的测量分配重要性。起初,区间从底部向顶部移动一定距离(在图 2a 和 2c 中显示为 2 和 11 公里),然后从顶部开始重复向下移动的过程。然后,将每个范围的拟合重叠值进行平均。这样获得的轮廓随后使用移动平均进行平滑。在轮廓的底部,由于拟合多项式中存在较大的虚假变化,我们必须减小滑动窗口的宽度,以便利用大部分观测数据。由于拟合多项式的表现不佳,上下 4 公里的轮廓必须被排除,因为它们无法与其他滑动窗口的对应值进行平均。该方法偶尔会产生不围绕零波动的干扰。 为了纠正这些数值偏差,我们通过应用 Theil-Sen 估计量(Sen, 1968; Theil, 1950)进行去趋势处理,并拟合一个线性函数到扰动轮廓上。Theil-Sen 估计量是一种稳健的方法,用于确定线性回归,取所有可以通过给定数据集绘制的线的斜率的中位数。然后,从轮廓中减去线性函数,以获得校正后的温度。
The results for the sliding polynomial fit are plotted in Figure 2 for the 2 and 11 km shift steps with green and yellow lines, correspondingly. It is seen that they are very close and, thus, the background and disturbances depend on the sliding step to a minor degree. The method shows some useful features in comparison with spectral filtering. The fitted mean curves in the regions of large-scale disturbances (Case 1) follow the observed temperature profiles closer (Figure 2a) and are smoother where small-scale structure dominates (Case 2) (Figure 2c, between 70 and 130 km ). This produces smaller wave amplitudes in Case 1, and reveals more wavy structures in Case 2. Especially plausible results are in the bottom of the profiles, where GWs are expected to have smaller amplitudes (due to larger density). 滑动多项式拟合的结果在图 2 中绘制,2 公里和 11 公里的移动步长分别用绿色和黄色线表示。可以看出,它们非常接近,因此背景和干扰在很小程度上依赖于滑动步长。与光谱滤波相比,该方法显示出一些有用的特征。在大尺度干扰区域(案例 1)中,拟合的平均曲线更接近观察到的温度剖面(图 2a),而在小尺度结构主导的区域(案例 2)(图 2c,70 到 130 公里之间)则更平滑。这在案例 1 中产生了较小的波幅,并在案例 2 中揭示了更多的波动结构。特别可信的结果出现在剖面的底部,预计重力波的幅度较小(由于密度较大)。
We next explored the technique of fitting higher-order polynomials in the entire interval of heights. In particular, the seventh-order polynomial fit, which was previously used for extracting GWs on Mars (Jesch et al., 2019; Yiğit et al., 2015), produces most plausible results. They are presented in Figure 2 with dashed and dotted blue lines. It is immediately seen that thus obtained wave disturbances are in a very good agreement with those derived by the sliding polynomial fit method, especially for profiles containing small-scale features (Figure 2d). For profiles dominated by large-scale perturbations, the agreement is also good in terms of the determined vertical structure of the wave, although the magnitudes are often exaggerated (Figure 2 b ). The weak point of the method is that it occasionally produces spurious disturbances near the edges of the vertical domain with vertical gradients of the mean temperature directed opposite to the measured profiles. After careful consideration of the three methods applied to the available measurements, we selected the sliding third-order polynomial fit as the most appropriate and robust. 我们接下来探讨了在整个高度区间内拟合高阶多项式的技术。特别是,第七阶多项式拟合,之前用于提取火星上的重力波(Jesch et al., 2019; Yiğit et al., 2015),产生了最可信的结果。它们以虚线和点线蓝色线条的形式呈现在图 2 中。可以立即看出,这样获得的波动与通过滑动多项式拟合方法得出的结果非常一致,特别是对于包含小尺度特征的剖面(图 2d)。对于以大尺度扰动为主的剖面,尽管幅度往往被夸大,但在确定波的垂直结构方面也有良好的一致性(图 2b)。该方法的弱点在于,它偶尔会在垂直域的边缘产生虚假扰动,且平均温度的垂直梯度与测量剖面相反。在对三种应用于可用测量的方法进行仔细考虑后,我们选择了滑动三阶多项式拟合作为最合适和稳健的方法。
3.3. Wave Activity and Potential Energy 3.3. 波动活动和潜在能量
The GW field is often characterized by the magnitude of fluctuations |T^(')|=( bar(T^('2)))^(1//2)\left|T^{\prime}\right|=\left(\overline{T^{\prime 2}}\right)^{1 / 2} and wave potential energy (per unit mass) GW 场通常以波动的幅度 |T^(')|=( bar(T^('2)))^(1//2)\left|T^{\prime}\right|=\left(\overline{T^{\prime 2}}\right)^{1 / 2} 和波动势能(每单位质量)
gg is the acceleration of gravity and c_(p)c_{p} is the specific heat capacity at constant pressure. The amplitude of the wave packet at a given height | T^(')(z)∣T^{\prime}(z) \mid (hereafter called “wave activity”) represents an envelope of the measured profile T^(')(z)T^{\prime}(z). We calculated it by performing Fourier decomposition in each sliding 60-km60-\mathrm{km} vertical gg 是重力加速度, c_(p)c_{p} 是定压下的比热容。给定高度处波包的振幅 | T^(')(z)∣T^{\prime}(z) \mid (以下简称“波动活动”)表示测量轮廓的包络 T^(')(z)T^{\prime}(z) 。我们通过在每个滑动 60-km60-\mathrm{km} 垂直 进行傅里叶分解来计算它。
Figure 3. Wave activity | T^(')∣T^{\prime} \mid (left column) and potential energy (per unit mass, right column) for the same as in Figure 2 representative profiles. Dashed blue lines indicate quantities calculated for the entire spectrum, dashed red lines are for accounting two longest harmonics only. Shaded areas denote observational errors. 图 3. 波动活动 | T^(')∣T^{\prime} \mid (左列)和潜在能量(每单位质量,右列)与图 2 中代表性剖面相同。虚线蓝线表示为整个谱计算的量,虚线红线仅用于考虑两个最长的谐波。阴影区域表示观测误差。
interval and, based on Parceval’s identity, summing up contributions of all harmonics. Examples of thus obtained envelopes and potential energy for the same selected profiles as in Section 3.2 are presented in Figure 3. 基于帕尔塞瓦尔定理的间隔,并汇总所有谐波的贡献。图 3 中展示了与第 3.2 节中选择的相同轮廓所获得的包络和潜在能量的示例。
Blue and red dashed lines denote the quantities calculated from the entire spectrum and by accounting for contributions of only two largest harmonics. It is seen that the neglect of shorter-scale harmonics, as was occasionally done in analyses of satellite observations (e.g., Ern et al., 2004), introduces little error to the estimated GW activity. However, the neglect of short-scale harmonics may lead to a noticeable underestimation of wave potential energy (cf. Figures 3b and 3d). 蓝色和红色虚线表示从整个光谱计算的量,以及仅考虑两个最大谐波的贡献。可以看出,忽略短尺度谐波(如在卫星观测分析中偶尔进行的那样,参见 Ern 等,2004)对估计的重力波活动引入的误差很小。然而,忽略短尺度谐波可能会导致波动势能的明显低估(参见图 3b 和 3d)。
3.4. Momentum Flux and Momentum Deposition 3.4. 动量通量和动量沉积
Another useful characteristic of the GW field is the vertical flux of horizontal momentum, or “momentum flux” for brevity, F=(F_(x),F_(y),0)=rho_(0)( bar(u^(')w^(')), bar(v^(')w^(')),0)\mathbf{F}=\left(F_{x}, F_{y}, 0\right)=\rho_{0}\left(\overline{u^{\prime} w^{\prime}}, \overline{v^{\prime} w^{\prime}}, 0\right), where rho_(0)\rho_{0} is the mean density and ( {:u^('),v^('),w^('))\left.u^{\prime}, v^{\prime}, w^{\prime}\right) are the components of wave-induced perturbations of wind velocity u^(')\mathbf{u}^{\prime} along with the two horizontal and the vertical axis, correspondingly. Momentum flux is constant for conservatively propagating waves. Breaking/dissipating GWs deposit their momentum to the mean flow, thus inducing an acceleration or deceleration (depending on the sign) of the horizontal flow. GW 场的另一个有用特征是水平动量的垂直通量,或简而言之的“动量通量”, F=(F_(x),F_(y),0)=rho_(0)( bar(u^(')w^(')), bar(v^(')w^(')),0)\mathbf{F}=\left(F_{x}, F_{y}, 0\right)=\rho_{0}\left(\overline{u^{\prime} w^{\prime}}, \overline{v^{\prime} w^{\prime}}, 0\right) ,其中 rho_(0)\rho_{0} 是平均密度,( {:u^('),v^('),w^('))\left.u^{\prime}, v^{\prime}, w^{\prime}\right) 是风速波动的分量 u^(')\mathbf{u}^{\prime} ,分别沿着两个水平轴和垂直轴。动量通量对于保守传播的波是恒定的。破裂/耗散的重力波将其动量传递给平均流,从而引起水平流的加速或减速(取决于符号)。
The direction of the flux cannot be determined from the occultation measurements, however total (or absolute) momentum fluxes for a harmonic F_(k,m)=sqrt(F_(x,k,m)^(2)+F_(y,k,m)^(2))F_{k, m}=\sqrt{F_{x, k, m}^{2}+F_{y, k, m}^{2}} can be estimated (e.g., Ern et al., 2004, sect. 4): 通量的方向无法从掩星测量中确定,但可以估计谐波的总(或绝对)动量通量 F_(k,m)=sqrt(F_(x,k,m)^(2)+F_(y,k,m)^(2))F_{k, m}=\sqrt{F_{x, k, m}^{2}+F_{y, k, m}^{2}} (例如,Ern 等,2004 年,第 4 节):
where k_(h)k_{h} and mm are the horizontal and vertical wavenumbers, correspondingly, and |T_(k,m)^(')|\left|T_{k, m}^{\prime}\right| is the amplitude. The latter two are found from the Fourier decomposition, whereas k_(h)k_{h} cannot be derived from our measurements. 其中 k_(h)k_{h} 和 mm 分别是水平和垂直波数, |T_(k,m)^(')|\left|T_{k, m}^{\prime}\right| 是振幅。后两个是通过傅里叶分解得到的,而 k_(h)k_{h} 不能从我们的测量中推导出来。
The total flux FF is the sum of contributions of individual harmonics F=sum_(m)F_(k,m)F=\sum_{m} F_{k, m}. Since the horizontal wavenumber k_(h)k_{h} cannot be obtained from the measurements, it, therefore, serves as a scaling factor for the derived profiles of FF and momentum forcing (5). The densest atmospheric footprint at a target point in occultation geometry is 400-500km400-500 \mathrm{~km} horizontally, depending on the height. This constrains the upper limit for unresolved wavelengths. In our calculations, we assumed a representative horizontal wavelength lambda_(h)=2pi//k_(h)=300km\lambda_{h}=2 \pi / k_{h}=300 \mathrm{~km}, the value typically used in numerical general circulation models (Yiğit et al., 2018), and allowing for more direct comparison with simulations. The results for two representative profiles, same as in Figures 2 and 3, are given in Figure 4. 总通量 FF 是各个谐波贡献的总和 F=sum_(m)F_(k,m)F=\sum_{m} F_{k, m} 。由于无法从测量中获得水平波数,因此它作为推导的 FF 和动量强迫(5)的缩放因子。目标点在掩蔽几何中的最密集大气足迹是 400-500km400-500 \mathrm{~km} 水平的,取决于高度。这限制了未解析波长的上限。在我们的计算中,我们假设一个代表性的水平波长 lambda_(h)=2pi//k_(h)=300km\lambda_{h}=2 \pi / k_{h}=300 \mathrm{~km} ,这是数值一般环流模型中通常使用的值(Yiğit 等,2018),并允许与模拟进行更直接的比较。两个代表性剖面的结果,与图 2 和图 3 相同,见图 4。
To demonstrate the sensitivity of the calculations to the used parameters of the technique, we plotted with different colors the profiles of momentum fluxes (per unit mass) F//rho_(0)F / \rho_{0} and GW momentum deposition, that is, wave drag aa obtained from the full spectrum and taking account of only two major harmonics. In addition, the results are shown for the interval shifts 2 and 7 km . It is immediately seen that these details play little role, and the calculations of fluxes and wave drag are very robust when the measured temperature profile is dominated by large-scale features (Figure 4, the upper row). It is different for profiles containing smaller vertical-scale disturbances (Figure 4, the lower row): their neglect leads to an underestimation of the fluxes and wave drag, and the smaller vertical shifts reveal finer structure associated with dissipation of individual spectral harmonics. 为了演示计算对所用技术参数的敏感性,我们用不同的颜色绘制了动量通量(每单位质量)和重力波动量沉积,即波阻力的剖面,这些数据是从全谱中获得的,并仅考虑了两个主要谐波。此外,结果显示在 2 公里和 7 公里的间隔变化中。可以立即看出,这些细节作用不大,当测量的温度剖面主要由大尺度特征主导时,通量和波阻力的计算非常稳健(图 4,上排)。对于包含较小垂直尺度扰动的剖面(图 4,下排),情况则不同:忽略它们会导致通量和波阻力的低估,而较小的垂直变化揭示了与单个谱谐波耗散相关的更精细结构。
4. Results and Discussions 4. 结果与讨论
4.1. Case Study 4.1. 案例研究
Spectral analysis of the obtained set of profiles (described in the next subsection) has demonstrated greater contribution of larger-scale disturbances in all cases. However, each individual profile was unique. Two examples with and without small vertical-scale components have been presented above. We next consider a case with a relatively broad spectrum of wave-like perturbations with large amplitudes (about twice as large as those in orbit 3251 n 1 ). The retrieved temperature for the orbit 4926 n 1 along with the fitted background profile are plotted in Figure 5a. The envelope in Figure 5b clearly shows that the amplitude gradually ceases its exponential growth with height and becomes nearly constant above ∼110km\sim 110 \mathrm{~km}. 所获得的剖面的光谱分析(在下一小节中描述)表明,在所有情况下,大尺度扰动的贡献更大。然而,每个单独的剖面都是独特的。上面展示了两个例子,一个有小的垂直尺度分量,一个没有。接下来,我们考虑一个具有相对宽广波动谱的案例,波动幅度较大(约为轨道 3251 n 1 的两倍)。轨道 4926 n 1 的检索温度以及拟合的背景剖面绘制在图 5a 中。图 5b 中的包络线清楚地显示,幅度随着高度的增加逐渐停止其指数增长,并在 ∼110km\sim 110 \mathrm{~km} .以上几乎保持不变。
The reason for this so-called wave “saturation” can be seen from the behavior of the squared Brunt-Väisälä frequency N^(2)(z)N^{2}(z) (Figure 5c, black). N^(2)N^{2} calculated from the background profiles (Figure 5c) remains relatively constant with height (up to about 120 km ) suggesting convective stability of the mean state. N^(2)N^{2} from the original profiles (see Figure 5c, red-dashed) shows large swings associated with temperature disturbances. Near 110km,N^(2)110 \mathrm{~km}, N^{2} drops almost to zero as the result of the temperature gradient (associated with a large amplitude of the disturbances) approaching the adiabatic lapse rate. Enhanced wave dissipation due to a combination of physical processes (Yiğit et al., 2018) in the vicinity of the convective instability severely limits the GW amplitude, leading to the decrease of the momentum flux above this altitude and peaking of the mean flow acceleration (Figure 5d) at almost 2,000ms^(-1)sol^(-1)2,000 \mathrm{~m} \mathrm{~s}^{-1} \mathrm{sol}^{-1}. In the analyzed data set, such large numbers are not common and occur only occasionally. Application of a smaller vertical shift of sliding intervals shows finer structure of the GW momentum flux and drag, but do not significantly modify the magnitudes. 这种所谓的波“饱和”现象的原因可以从平方布伦特-维萨拉频率的行为中看出(图 5c,黑色)。从背景剖面(图 5c)计算得出的频率在高度上保持相对恒定(最高约 120 公里),这表明平均状态的对流稳定性。来自原始剖面(见图 5c,红色虚线)显示出与温度扰动相关的大幅波动。由于温度梯度(与扰动的大幅度相关)接近绝热递减率,几乎降到零。由于物理过程的组合(Yiğit 等,2018)在对流不稳定性附近增强了波的耗散,严重限制了重力波的幅度,导致该高度以上的动量通量减少,并且平均流加速度在几乎达到峰值(图 5d)。在分析的数据集中,这样的大数字并不常见,仅偶尔出现。 较小的垂直位移滑动区间的应用显示了重力波动量通量和阻力的更精细结构,但并未显著改变其大小。
Figure 4. Absolute momentum flux (per unit mass) and the momentum forcing for two representative profiles (orbits 2892n1 and 3251n1, upper and lower rows, correspondingly). The legend describes the profiles calculated using the full spectrum and only two major harmonics along with sliding interval steps 2 and 7 km . 图 4. 绝对动量通量(每单位质量)和两个代表性剖面的动量强迫(轨道 2892n1 和 3251n1,上下行,分别)。图例描述了使用全谱和仅使用两个主要谐波计算的剖面,以及滑动间隔步长为 2 和 7 公里。
4.2. Spatial Distribution of Gravity Wave Activity 4.2. 重力波活动的空间分布
In this section, we use the data obtained by the ACS instrument in MY34, at solar longitudes from L_(s)=164^(@)L_{s}=164^{\circ} to 354^(@)354^{\circ}. The data set contains altogether 144 occultation profiles: 84 in the northern hemisphere and 60 in the southern one. The latitude-solar longitude coverage is shown in Figure 6 with red and blue dots representing morning and evening occultation measurements, correspondingly. The longitudinal orbit coverage was fairly uniform, and is not discussed here. 在本节中,我们使用在 MY34 由 ACS 仪器获得的数据,太阳经度范围从 L_(s)=164^(@)L_{s}=164^{\circ} 到 354^(@)354^{\circ} 。数据集总共包含 144 个掩星剖面:北半球 84 个,南半球 60 个。纬度-太阳经度覆盖范围如图 6 所示,红点和蓝点分别代表早晨和晚上的掩星测量。经度轨道覆盖相当均匀,这里不作讨论。
A significant portion of observations were made during the global dust storm of MY34, which started between L_(s)=185^(@)L_{s}=185^{\circ} and 190^(@)190^{\circ}, attained its maximum around L_(s)=220^(@)L_{s}=220^{\circ}, and gradually decreased until L_(s)~~290^(@)L_{s} \approx 290^{\circ}. A regional storm occurred at the end of MY34 between approximately L_(s)=325^(@)L_{s}=325^{\circ} and 345^(@)345^{\circ}. Figure 7 presents latitude-altitude distribution of the derived GW parameters averaged over the entire period of observations depicted in Figure 6. 在 MY34 全球沙尘暴期间,观察的一个重要部分是在 L_(s)=185^(@)L_{s}=185^{\circ} 和 190^(@)190^{\circ} 之间进行的,达到最大值是在 L_(s)=220^(@)L_{s}=220^{\circ} ,并逐渐减少直到 L_(s)~~290^(@)L_{s} \approx 290^{\circ} 。在 MY34 结束时,发生了一场区域性风暴,大约在 L_(s)=325^(@)L_{s}=325^{\circ} 和 345^(@)345^{\circ} 之间。图 7 展示了在图 6 所示的整个观察期间平均的派生 GW 参数的纬度-高度分布。
It shows that the mean amplitude of GW-induced temperature fluctuations ( |T^(')|\left|T^{\prime}\right|, Figure 7a) grows with height reaching up to ∼10K\sim 10 \mathrm{~K} near the top of the domain. At higher altitudes ( 170-220km170-220 \mathrm{~km} ), the in situ measurements with Neutral Gas and Ion Mass Spectrometer (NGIMS) on board MAVEN revealed even larger GW magnitudes over the same time (Leelavathi et al., 2020; Yiğit et al., 2021). The latitudinal structure of the GW activity in the mesosphere and lower thermosphere is not uniform. For comparison, we overplotted the zonal wind simulated with the Max Planck Institute (MAOAM) MGCM https://mars.mipt.ru/ for 这表明 GW 引起的温度波动的平均幅度( |T^(')|\left|T^{\prime}\right| ,图 7a)随着高度的增加而增长,达到 ∼10K\sim 10 \mathrm{~K} ,接近区域的顶部。在更高的高度( 170-220km170-220 \mathrm{~km} ),搭载在 MAVEN 上的中性气体和离子质谱仪(NGIMS)的原位测量显示,在相同时间内 GW 的幅度甚至更大(Leelavathi 等,2020;Yiğit 等,2021)。中层大气和下热层中 GW 活动的纬向结构并不均匀。为了比较,我们叠加了由马克斯·普朗克研究所(MAOAM)MGCM 模拟的纬向风https://mars.mipt.ru/。
MY34 and averaged over the same interval of L_(s)L_{s} as in the observations. The wind distribution varied during this time from the equinoctial to solstitial and back to the equinoctial types. The result reflects the largest contribution of the prograde and retrograde jets during the perihelion solstice. It is seen that the regions with large wave amplitudes encircle the upper edges of two midlatitude jets. This is the result of intensive filtering of individual harmonics by strong background winds. For the wave potential energy, which is a quadratic function of wave amplitudes, this pattern is even more obvious (Figure 7b). MY34 和在观察中相同时间间隔的平均值。此期间的风分布从春分型变化到夏至型,再回到春分型。结果反映了在近日点夏至期间,顺行和逆行喷流的最大贡献。可以看到,大波幅区域环绕着两个中纬度喷流的上缘。这是强背景风对个别谐波进行强烈过滤的结果。对于波动势能,它是波幅的二次函数,这种模式更加明显(图 7b)。
Figure 7 c shows that GW momentum fluxes reach local maxima near the mesopause ( 100-125km100-125 \mathrm{~km} ) giving evidence of very intensive wave breakdown/dissipation in this region. The peaks of the associated momentum deposition approximately coincide (Figure 7d). They too wrap around the edges of the jets in the middle atmosphere. It is noteworthy that such distribution of the GW drag is 图 7c 显示,重力波动量通量在中间大气的中层边缘附近达到局部最大值( 100-125km100-125 \mathrm{~km} ),这表明该区域存在非常强烈的波动破裂/耗散。相关动量沉积的峰值大致重合(图 7d)。它们也环绕在中层大气的喷流边缘。值得注意的是,重力波阻力的这种分布是
Figure 6. Latitude-solar longitude ( L_(s)L_{s} ) distribution of the ACS MIR occultation profiles used in this study. Morning and evening measurements are shown in red and blue, correspondingly. 图 6. 本研究中使用的 ACS MIR 掩蔽剖面的纬度-太阳经度( L_(s)L_{s} )分布。早晨和晚上的测量分别用红色和蓝色表示。
very similar to that predicted by a Martian GCM (Medvedev et al., 2011, Figures 3 and 7) for the solstice and equinox, respectively, and represents the first (to the best of our knowledge) observational validation of the model predictions. The magnitudes of the GW drag, although defined up to the constant k_(h)k_{h}, agree with the simulations (using a similar k_(h)k_{h} ) as well. 与火星 GCM(Medvedev 等,2011 年,图 3 和图 7)对至日和春分的预测非常相似,并且代表了我们所知的模型预测的首次观测验证。GW 拖曳的大小,尽管定义为常数 k_(h)k_{h} ,也与模拟结果(使用类似的 k_(h)k_{h} )一致。
4.3. Amplitude Dependence on Mean Temperature and BruntVäisälä Frequency 4.3. 振幅对平均温度和布伦特-维萨拉频率的依赖
In situ measurements with NGIMS on board MAVEN showed a clear anti-correlation between relative density fluctuations in the upper thermosphere and the ambient temperature (England et al., 2017; Terada et al., 2017; Vals et al., 2019; Yiğit et al., 2015). It was linked to convective instability as a dominant mechanism that limits growth of GW 在 MAVEN 上搭载的 NGIMS 进行的原位测量显示,上层热层中相对密度波动与环境温度之间存在明显的反相关关系(England 等,2017;Terada 等,2017;Vals 等,2019;Yiğit 等,2015)。这与对流不稳定性相关,作为限制重力波增长的主要机制。
where |u^(')|\left|u^{\prime}\right| is the amplitude of fluctuations of horizontal velocity in the wave, cc is its horizontal phase velocity and bar(u)\bar{u} is the background wind. When |u^(')|\left|u^{\prime}\right| approaches |c- bar(u)||c-\bar{u}|, increasing dissipation limits |u^(')|\left|u^{\prime}\right| thus that the ratio |u^(')|//|c- bar(u)|\left|u^{\prime}\right| /|c-\bar{u}| becomes constant. The linear convective instability threshold demands a unit ratio, however observations suggested a ratio of 0.7 (Fritts et al., 1988; Equation 2), and the theoretical consideration of the nonlinear diffusion mechanism yielded 1//sqrt2~~0.7071 / \sqrt{2} \approx 0.707 (Medvedev & Klaassen, 2000, Section 7). Regardless of the precise number, Equation 7 establishes proportionality between the amplitude of relative temperature/density perturbations and squared mean Brunt-Väisälä frequency under the saturation condition. Near the exobase, where the majority of NGIMS/MAVEN observations were taken, the vertical gradient d bar(T)//dzd \bar{T} / d z is small and can be neglected in Equation 4, thus giving the inverse proportionality of relative perturbation amplitudes and bar(T)\bar{T}. 其中 |u^(')|\left|u^{\prime}\right| 是波中水平速度波动的幅度, cc 是其水平相速度, bar(u)\bar{u} 是背景风。当 |u^(')|\left|u^{\prime}\right| 接近 |c- bar(u)||c-\bar{u}| 时,增加的耗散限制了 |u^(')|\left|u^{\prime}\right| ,使得比率 |u^(')|//|c- bar(u)|\left|u^{\prime}\right| /|c-\bar{u}| 变得恒定。线性对流不稳定性阈值要求单位比率,然而观察结果表明比率为 0.7(Fritts 等,1988;方程 2),而非线性扩散机制的理论考虑得出了 1//sqrt2~~0.7071 / \sqrt{2} \approx 0.707 (Medvedev & Klaassen,2000,第 7 节)。无论精确数字如何,方程 7 确立了相对温度/密度扰动幅度与饱和条件下的平方平均 Brunt-Väisälä 频率之间的比例关系。在外层大气基底附近,大多数 NGIMS/MAVEN 观测数据是在此处获得的,垂直梯度 d bar(T)//dzd \bar{T} / d z 较小,可以在方程 4 中忽略,从而给出了相对扰动幅度与 bar(T)\bar{T} 的反比例关系。
ACS/TGO occultation data cover altitudes below the exobase, where d bar(T)//dzd \bar{T} / d z can no longer be neglected. Therefore, we plotted in Figure 8a the amplitudes of relative temperature perturbations for all orbits as functions of N^(2)N^{2}. ACS/TGO 消失数据覆盖了低于外逸层的高度,在这里 d bar(T)//dzd \bar{T} / d z 不再可以忽略。因此,我们在图 8a 中绘制了所有轨道的相对温度扰动幅度作为 N^(2)N^{2} . 的函数。
It is seen that red and blue dots corresponding to morning and evening measurements show no clear dependence on N^(2)N^{2} at all altitudes. To explore this further, we over-plotted the linear regression of the form |T^(')|// bar(T)=alphaN^(2)+beta\left|T^{\prime}\right| / \bar{T}=\alpha N^{2}+\beta and put the values of alpha\alpha and beta\beta in the legend. The coefficients alpha\alpha are far less than those expected from Equation 7, that is, several tens or hundreds, depending on the characteristic vertical wavenumber mm. The distinction between morning and evening amplitudes is also insignificant, except above 100 km , where morning values are slightly larger. 可以看出,早晨和晚上的测量对应的红点和蓝点在所有高度上都没有明显的依赖关系。 N^(2)N^{2} 为了进一步探讨这一点,我们叠加了形式为 |T^(')|// bar(T)=alphaN^(2)+beta\left|T^{\prime}\right| / \bar{T}=\alpha N^{2}+\beta 的线性回归,并在图例中放置了 alpha\alpha 和 beta\beta 的值。 alpha\alpha 系数 mm 远低于方程 7 所预期的值,即几十或几百,具体取决于特征垂直波数。{{6}} 早晨和晚上的幅度差异也不显著,除了在 100 公里以上,早晨的值略大一些。
Figure 8. Amplitudes of relative temperature disturbances as functions of the squared Brunt-Väisälä frequency (a) and mean temperature (b) at different heights. Red and blue dots are for the morning and evening measurements, correspondingly. Linear regressions of the form (a)। T^(')∣// bar(T)=alphaN^(2)+betaT^{\prime} \mid / \bar{T}=\alpha N^{2}+\beta and (b) |T^(')|// bar(T)=alpha bar(T)+beta\left|T^{\prime}\right| / \bar{T}=\alpha \bar{T}+\beta are shown with thin solid lines, and the values of the respective coefficients are given in the legends. 图 8. 不同高度下相对温度扰动的幅度与平方布伦特-维萨拉频率(a)和平均温度(b)的关系。红点和蓝点分别代表早晨和晚上的测量。线性回归形式为(a)。 T^(')∣// bar(T)=alphaN^(2)+betaT^{\prime} \mid / \bar{T}=\alpha N^{2}+\beta 和(b) |T^(')|// bar(T)=alpha bar(T)+beta\left|T^{\prime}\right| / \bar{T}=\alpha \bar{T}+\beta 以细实线表示,相应系数的值在图例中给出。
Figure 8 b presents the dependencies of amplitudes of relative temperature disturbances as functions of the mean temperature. They are nearly uniform. Although regression coefficients show a weak negative trends at all altitudes, their magnitudes are much smaller than to those observed previously (of the order of 0.5-1) near the exobase (England et al., 2017; Terada et al., 2017; Vals et al., 2019; Yiğit et al., 2015). A similar lack of correlation between GW amplitudes and atmospheric temperature was found from TGO aerobraking measurements at altitudes between 100 and 130 km (Jesch et al., 2019, Figure 12). The atmospheric drag data were collected between L_(s)=332^(@)L_{s}=332^{\circ} of MY33 and L_(s)=132^(@)L_{s}=132^{\circ} of MY34. The ACS observations after the aerobraking cover the dusty second half of MY34. Thus, the absence of correlation between GW amplitudes and the background temperature in the lower thermosphere appear to be independent of the season and dust conditions. In the upper thermosphere (Leelavathi et al., 2020, Figure 10d), found a positive correlation during the same second half of MY34, instead of a clear negative correlation over the first (“non-dusty”) half of the year. Our results in the adjacent region (around 140 km ) show no visible change, neither strong negative trend previously found in the MAVEN/NGIMS observations, nor indication of a positive trend. This means that convective instability may not be the main mechanism responsible for damping GWs in the thermosphere, at least during dust storms. 图 8b 展示了相对温度扰动幅度与平均温度的依赖关系。它们几乎是均匀的。尽管回归系数在所有高度上显示出微弱的负趋势,但其幅度远小于之前观察到的(约为 0.5-1)在外气层附近(England 等,2017;Terada 等,2017;Vals 等,2019;Yiğit 等,2015)。在 100 到 130 公里的高度上,从 TGO 气动制动测量中发现 GW 幅度与大气温度之间的相关性也很弱(Jesch 等,2019,图 12)。大气阻力数据是在 MY33 的 L_(s)=332^(@)L_{s}=332^{\circ} 和 L_(s)=132^(@)L_{s}=132^{\circ} 的 MY34 之间收集的。气动制动后的 ACS 观测覆盖了 MY34 的多尘后半段。因此,GW 幅度与下热层背景温度之间缺乏相关性似乎与季节和尘埃条件无关。在上热层(Leelavathi 等,2020,图 10d),在 MY34 的同一后半段发现了正相关,而不是在第一半段(“无尘”)的明显负相关。 我们在相邻区域(约 140 公里)的结果显示没有明显变化,既没有在 MAVEN/NGIMS 观测中发现的强负趋势,也没有正趋势的迹象。这意味着对流不稳定性可能不是在热层中抑制重力波的主要机制,至少在沙尘暴期间不是。
5. Summary and Conclusions 5. 摘要和结论
We have presented the results of gravity wave (GW) retrievals obtained from the Atmospheric Chemistry Suite instrument on board the ExoMars Trace Gas Orbiter (ACS/TGO), which observed solar occultation spectra. GW disturbances are derived from the vertical temperature profiles retrieved from one of the three instrument channels - the mid-infrared ACS/MIR. The uniqueness of the data is that they continuously cover a broad range of altitudes from the Martian troposphere to the thermosphere ( 20-160km20-160 \mathrm{~km} ) and have a relatively high ( 0.5-2.5km0.5-2.5 \mathrm{~km} ) vertical resolution. 我们展示了从火星探测器(ExoMars Trace Gas Orbiter,ACS/TGO)上的大气化学套件仪器获得的重力波(GW)检索结果,该仪器观察了太阳掩蔽光谱。GW 扰动是从三个仪器通道之一的垂直温度剖面中提取的 - 中红外 ACS/MIR。数据的独特之处在于它们持续覆盖从火星对流层到热层的广泛高度范围( 20-160km20-160 \mathrm{~km} ),并具有相对较高的( 0.5-2.5km0.5-2.5 \mathrm{~km} )垂直分辨率。
Several techniques of separating GW components from the background temperature have been studied. The sliding-window least squares polynomial fitting method have demonstrated to be the most robust and 已经研究了几种将 GW 成分与背景温度分离的技术。滑动窗口最小二乘多项式拟合方法已被证明是最稳健的。
effective. The procedure was applied to 144 measurements collected over the second half of MY34 to derive vertical profiles of GW disturbances as well as further wave characteristics: amplitude, wave potential energy, absolute vertical flux of horizontal momentum and absolute momentum forcing produced by breaking/ dissipating GWs (“GW drag”). The main results are listed below. 有效的。该程序应用于在 MY34 下半年收集的 144 个测量数据,以推导出 GW 扰动的垂直剖面以及进一步的波动特征:振幅、波动势能、水平动量的绝对垂直通量和由破裂/消散的 GW 产生的绝对动量强迫(“GW 拖曳”)。主要结果如下。
Amplitudes of GW-induced temperature fluctuations, generally, grow with height, while breaking/saturation processes often limit the wave amplitude growth at higher altitudes. Based on a half-year average, wave amplitudes are around 8-14K8-14 \mathrm{~K} near the mesopause, and often exceed these values in individual profiles. 引起的温度波动的幅度通常随着高度的增加而增加,而破裂/饱和过程往往限制了高海拔地区波幅的增长。基于半年的平均值,波幅在中间层附近约为 8-14K8-14 \mathrm{~K} ,并且在个别剖面中通常超过这些值。
The mesopause ( 100-120km100-120 \mathrm{~km} ) is the region of the strongest GW breaking/dissipation, which is evidenced by a local maximum of momentum fluxes and their vertical divergence, that is, GW drag. Similarly, a large GW drag of hundreds of ms^(-1)\mathrm{m} \mathrm{s}^{-1} sol ^(-1)^{-1} in the mesopause region has been demonstrated by MGCMs (e.g., Yiğit et al., 2018). 中间层顶 ( 100-120km100-120 \mathrm{~km} ) 是重力波破裂/耗散最强的区域,这通过动量通量及其垂直散度的局部最大值得以证明,即重力波拖曳。同样,MGCMs(例如,Yiğit 等,2018)已证明中间层顶区域的重力波拖曳可达数百 ms^(-1)\mathrm{m} \mathrm{s}^{-1} 太阳 ^(-1)^{-1} 。
The spatial (altitude-latitude) distribution of the wave drag also agrees well with modeling results (e.g., Medvedev et al., 2011). This is the first direct observational validation of model predictions. 波阻力的空间(高度-纬度)分布与模型结果(例如,Medvedev 等,2011)也很好地一致。这是对模型预测的第一次直接观测验证。
We did not find positive correlation between amplitudes of relative temperature perturbations and the Brunt-Väisälä frequency at all heights. This correlation is a more general formulation of the anti-correlation found near the exobase (England et al., 2017; Terada et al., 2017; Vals et al., 2019; Yiğit et al., 2015) that accounts for vertically varying mean temperature. 我们没有发现相对温度扰动幅度与所有高度的布伦特-维萨拉频率之间的正相关。这种相关性是对在外层大气基线附近发现的反相关性(England et al., 2017; Terada et al., 2017; Vals et al., 2019; Yiğit et al., 2015)的更一般的表述,考虑了垂直变化的平均温度。
The presented GW activity retrievals extending from the middle troposphere to the thermosphere, as derived from the ExoMars data, highlight the role of atmospheric gravity waves as a whole atmosphere phenomenon on Mars. Mars’ thin and windy atmosphere favors strong gravity wave generation, thus an accurate characterization of gravity waves is absolutely essential for a better understanding of the Martian climate (Yiğit & Medvedev, 2019). 从 ExoMars 数据得出的中层对流层到热层的 GW 活动检索,突显了大气重力波作为火星整体大气现象的作用。火星稀薄而多风的大气有利于强重力波的产生,因此准确表征重力波对于更好地理解火星气候至关重要(Yiğit & Medvedev, 2019)。
Data Availability Statement 数据可用性声明
The ACS data are available from ESA’s Planetary Science Archive at https://archives.esac.esa.int/psa/#!Table\\%20View//ACS=\backslash \% 20 \mathrm{View} / \mathrm{ACS}= instrument. The data for latitude-altitude cross-sections of GW characteristics such as amplitudes, potential energy, fluxes, accelerations and zonal wind are available at the Mendeley database (Starichenko, 2021). ACS 数据可从 ESA 的行星科学档案馆获取,网址为 https://archives.esac.esa.int/psa/#!Table\\%20View//ACS=\backslash \% 20 \mathrm{View} / \mathrm{ACS}= 仪器。关于 GW 特征的纬度-高度截面数据,如振幅、势能、通量、加速度和纬向风,可在 Mendeley 数据库中获取(Starichenko, 2021)。
Acknowledgments 致谢
ExoMars is the space mission of ESA and Roscosmos. The ACS experiment is led by IKI Space Research Institute in Moscow. Science operations of ACS are funded by Roscosmos and ESA. The analyses of temperature profiles and wave structures at IKI are funded by the grant #20-42-09035 of the Russian Science Foundation. ExoMars 是欧洲航天局和俄罗斯航天局的太空任务。ACS 实验由莫斯科的 IKI 空间研究所领导。ACS 的科学操作由俄罗斯航天局和欧洲航天局资助。IKI 对温度剖面和波动结构的分析由俄罗斯科学基金会的资助#20-42-09035 提供资金。
References 参考文献
Alday, J., Wilson, C. F., Irwin, P. G. J., Olsen, K. S., Baggio, L., Montmessin, F., et al. (2019). Oxygen isotopic ratios in Martian water vapour observed by ACS MIR on board the ExoMars Trace Gas Orbiter. Astronomy & Astrophysics, 630, A91. https://doi. org/10.1051/0004-6361/201936234 Alday, J., Wilson, C. F., Irwin, P. G. J., Olsen, K. S., Baggio, L., Montmessin, F., et al. (2019). 氧同位素比在火星水蒸气中通过 ExoMars 追踪气体轨道器上的 ACS MIR 观察到。天文学与天体物理学, 630, A91. https://doi. org/10.1051/0004-6361/201936234
Alexander, M. J. (1998). Interpretations of observed climatological patterns in stratospheric gravity wave variance. Journal of Geophysical Research, 103(D8), 8627-8640. https://doi.org/10.1029/97jd03325 亚历山大,M. J. (1998)。对平流层重力波方差中观察到的气候模式的解释。《地球物理研究杂志》,103(D8),8627-8640。https://doi.org/10.1029/97jd03325
Ando, H., Imamura, T., & Tsuda, T. (2012). Vertical wavenumber spectra of gravity waves in the Martian atmosphere obtained from Mars Global Surveyor radio occultation data. Journal of the Atmospheric Sciences, 69, 2906-2912. https://doi.org/10.1175/JAS-D-11-0339.1 安藤,H.,今村,T.,& 津田,T.(2012)。从火星全球探测器无线电掩星数据获得的火星大气中重力波的垂直波数谱。《大气科学杂志》,69,2906-2912。https://doi.org/10.1175/JAS-D-11-0339.1
Creasey, J. E., Forbes, J. M., & Hinson, D. P. (2006a). Global and seasonal distribution of gravity wave activity in Mars’ lower atmosphere derived from MGS radio occultation data. Geophysical Research Letters, 33. https://doi.org/10.1029/2005GL024037 Creasey, J. E., Forbes, J. M., & Hinson, D. P. (2006a). 从 MGS 无线电掩蔽数据得出的火星下层大气中重力波活动的全球和季节分布。地球物理研究快报, 33. https://doi.org/10.1029/2005GL024037
Creasey, J. E., Forbes, J. M., & Keating, G. M. (2006b). Density variability at scales typical of gravity waves observed in Mars’ thermosphere by the MGS accelerometer. Geophysical Research Letters, 33(L22814). https://doi.org/10.1029/2006gl027583 Creasey, J. E., Forbes, J. M., & Keating, G. M. (2006b). 火星热层中通过 MGS 加速度计观察到的重力波典型尺度下的密度变异性. 地球物理研究快报, 33(L22814). https://doi.org/10.1029/2006gl027583
Ehard, B., Kaifler, B., Kaifler, N., & Rapp, M. (2015). Evaluation of methods for gravity wave extraction from middle-atmospheric lidar temperature measurements. Atmospheric Measurement Techniques, 8(11), 4645-4655. https://doi.org/10.5194/amt-8-4645-2015 Ehard, B., Kaifler, B., Kaifler, N., & Rapp, M. (2015). 从中层大气激光雷达温度测量中提取重力波的方法评估. 大气测量技术, 8(11), 4645-4655. https://doi.org/10.5194/amt-8-4645-2015
England, S. L., Liu, G., Yiğit, E., Mahaffy, P. R., Elrod, M., Benna, M., et al. (2017). MAVEN NGIMS observations of atmospheric gravity waves in the Martian thermosphere. Journal of Geophysical Research: Space Physics, 122, 2310-2335. https://doi.org/10.1002/2016JA023475 英国,S. L.,刘,G.,Yiğit,E.,Mahaffy,P. R.,Elrod,M.,Benna,M.,等(2017)。MAVEN NGIMS 对火星热层大气重力波的观测。《地球物理研究杂志:空间物理》,122,2310-2335。 https://doi.org/10.1002/2016JA023475
Ern, M., Preusse, P., Alexander, M. J., & Warner, C. D. (2004). Absolute values of gravity wave momentum flux derived from satellite data. Journal of Geophysical Research, 109. https://doi.org/10.1029/2004JD004752 Ern, M., Preusse, P., Alexander, M. J., & Warner, C. D. (2004). 从卫星数据推导的重力波动量通量的绝对值. 地球物理研究杂志, 109. https://doi.org/10.1029/2004JD004752
Fedorova, A. A., Montmessin, F., Korablev, O., Luginin, M., Trokhimovskiy, A., Belyaev, D. A., et al. (2020). Stormy water on Mars: The distribution and saturation of atmospheric water during the dusty season. Science, 367, 297-300. https://doi.org/10.1126/science.aay9522 费多罗娃,A. A.,蒙特梅辛,F.,科拉布列夫,O.,卢金宁,M.,特罗希莫夫斯基,A.,别利亚耶夫,D. A.,等(2020)。火星上的暴风水:尘埃季节期间大气水的分布和饱和度。《科学》,367,297-300。 https://doi.org/10.1126/science.aay9522
Forget, F., Montmessin, F., Bertaux, J.-L., González-Galindo, F., Lebonnois, E., Quémerais, S., et al. (2009). Density and temperatures of the upper Martian atmosphere measured by stellar occultations with Mars Express SPICAM. Journal of Geophysical Research, 114(E01004). https://doi.org/10.1029/2008JE003086 Forget, F., Montmessin, F., Bertaux, J.-L., González-Galindo, F., Lebonnois, E., Quémerais, S., et al. (2009). 通过火星快车 SPICAM 的恒星掩蔽测量的火星上层大气的密度和温度. 地球物理研究杂志, 114(E01004). https://doi.org/10.1029/2008JE003086
Fritts, D. C., & Alexander, J. M. (2003). Gravity wave dynamics and effects in the middle atmosphere. Reviews of Geophysics, 41(1), 1003. https://doi.org/10.1029/2001rg000106 Fritts, D. C., & Alexander, J. M. (2003). 中层大气中的重力波动力学及其影响. 地球物理学评论, 41(1), 1003. https://doi.org/10.1029/2001rg000106
Fritts, D. C., Tsuda, T., Kato, S., Sato, T., & Fukao, S., (1988). Observational evidence of a saturated gravity wave spectrum in the troposphere and lower stratosphere. Journal of the Atmospheric Sciences, 45(12), 1741-1759. https://doi.org/10.1175/1520-0469(1988)045<1 741:OEOASG> 2.0.CO;2 Fritts, D. C., Tsuda, T., Kato, S., Sato, T., & Fukao, S., (1988). 对对流层和下平流层中饱和重力波谱的观测证据. 大气科学杂志, 45(12), 1741-1759. https://doi.org/10.1175/1520-0469(1988)045<1 741:OEOASG> 2.0.CO;2
Fritts, D. C., Wang, L., & Tolson, R. H. (2006). Mean and gravity wave structures and variability in the Mars upper atmosphere inferred from Mars Global Surveyor and Mars Odyssey aerobraking densities. Journal of Geophysical Research, 111(A12304). Fritts, D. C., Wang, L., & Tolson, R. H. (2006). 从火星全球探测器和火星奥德赛的气动制动密度推断的火星上层大气中的平均和重力波结构及其变异性。地球物理研究杂志, 111(A12304)。https://doi. org/10.1029/2006ja011897
Gamache, R. R., arese, M., & Renaud, C. L., (2016). A spectral line list for water isotopologues in the 1100-4100 cm-1 region for application to CO2-rich planetary atmospheres. Journal of Molecular Spectroscopy, 326, 144-150. https://doi.org/10.1016/j.jms.2015.09.001 Gamache, R. R., arese, M., & Renaud, C. L., (2016). 一份用于 CO2 丰富的行星大气的 1100-4100 cm-1 区域水同位素谱线列表。分子光谱学杂志, 326, 144-150. https://doi.org/10.1016/j.jms.2015.09.001
Gordon, I., Rothman, L., Hill, C., Kochanov, R., Tan, Y., Bernath, P., & Boudon, V. (2017). The HITRAN2016 molecular spectroscopic database. Journal of Quantitative Spectroscopy and Radiative Transfer, 203, 3-69. https://doi.org/10.1016/j.jqsrt.2017.06.038
Gröller, H., Montmessin, F., Yelle, R. V., Lefèvre, F., Forget, F., Schneider, N. M., et al. (2018). MAVEN/IUVS stellar occultation measurements of Mars atmospheric structure and composition. Journal of Geophysical Research, 123, 1449-1483-1483. Gröller, H., Montmessin, F., Yelle, R. V., Lefèvre, F., Forget, F., Schneider, N. M., et al. (2018). MAVEN/IUVS 星体掩蔽测量火星大气结构和成分。地球物理研究杂志, 123, 1449-1483-1483。https://doi. org/10.1029/2017JE005466
Heavens, N. G., Kass, D. M., Kleinböhl, A., & Schofield, J. T. (2020). A multiannual record of gravity wave activity in Mars’s lower atmosphere from on-planet observations by the Mars Climate Sounder. Icarus, 341, 113630. https://doi.org/10.1016/j.icarus.2020.113630 Heavens, N. G., Kass, D. M., Kleinböhl, A., & Schofield, J. T. (2020). 一项来自火星气候探测器的行星观测的火星下层大气中重力波活动的多年记录。Icarus, 341, 113630. https://doi.org/10.1016/j.icarus.2020.113630
Hinson, D. P., Simpson, R. A., Twicken, J. D., Tyler, G. L., & Flasar, F. M. (1999). Initial results from radio occultation measurements with Mars Global Surveyor. Journal of Geophysical Research, 104(E11), 26997-27012. https://doi.org/10.1029/1999JE001069
Irwin, P., Teanby, N., de Kok, R., Fletcher, L., Howett, C., Tsang, C., et al. (2008). The NEMESIS planetary atmosphere radiative transfer and retrieval tool. Journal of Quantitative Spectroscopy and Radiative Transfer, 109(6), 1136-1150. https://doi.org/10.1016/j.jqsrt.2007.11.006
Jesch, D., Medvedev, A. S., Castellini, F., Yiğit, E., & Hartogh, P. (2019). Density fluctuations in the lower thermosphere of mars retrieved from the exomars Trace Gas Orbiter (TGO) aerobraking. Atmosphere, 10(10). Jesch, D., Medvedev, A. S., Castellini, F., Yiğit, E., & Hartogh, P. (2019). 从 ExoMars 追踪气体轨道器(TGO)气动制动中获取的火星下热层的密度波动。大气,10(10)。https://doi.org/10.3390/atmos10100620
John, S. R., & Kumar, K. K., (2013). A discussion on the methods of extracting gravity wave perturbations from space-based measurements. Geophysical Research Letters, 40, 2406-2410. https://doi.org/10.1002/grl.50451
Keating, G. M., Bougher, S. W., Zurek, R. W., Tolson, R. H., Cancro, G. J., Noll, S. N., et al. (1998). The structure of the upper atmosphere of Mars: In situ accelerometer measurements from Mars global surveyor. Science, 279, 1672-1676. Keating, G. M., Bougher, S. W., Zurek, R. W., Tolson, R. H., Cancro, G. J., Noll, S. N., et al. (1998). 火星上层大气的结构:来自火星全球勘测者的原位加速度计测量。科学, 279, 1672-1676.https://doi.org/10.1126/ science.279.5357.1672
Korablev, O., Montmessin, F., Trokhimovskiy, A., Fedorova, A., Shakun, A., Grigoriev, A., et al. (2018). The Atmospheric Chemistry Suite (ACS) of three spectrometers for the ExoMars 2016 Trace Gas Orbiter. Space Science Reviews, 214. Korablev, O., Montmessin, F., Trokhimovskiy, A., Fedorova, A., Shakun, A., Grigoriev, A., et al. (2018). ExoMars 2016 追踪气体轨道器的三台光谱仪的大气化学套件(ACS)。空间科学评论,214。https://doi.org/10.1007/ s11214-017-0437-6
Leelavathi, V., Venkateswara Rao, N., & Rao, S. V. B. (2020). Interannual variability of atmospheric gravity waves in the martian thermosphere: Effects of the 2018 planet-encircling dust event. Journal of Geophysical Research: Plan, 125(12), e2020JE006649. Leelavathi, V., Venkateswara Rao, N., & Rao, S. V. B. (2020). 火星热层大气重力波的年际变异性:2018 年行星环绕尘埃事件的影响。地球物理研究杂志:行星, 125(12), e2020JE006649.https://doi. org/10.1029/2020JE006649
Marquardt, D. W. (1963). An algorithm for least-squares estimation of nonlinear parameters. Journal of the Society for Industrial and Applied Mathematics, 11(2), 431-441. Marquardt, D. W. (1963). 一种用于非线性参数最小二乘估计的算法。工业与应用数学学会杂志, 11(2), 431-441.https://doi.org/10.1137/0111030
Medvedev, A. S., & Klaassen, G. P. (2000). Parameterization of gravity wave momentum deposition based on nonlinear wave interactions: Basic formulation and sensitivity tests. JASTP, 62, 1015-1033. https://doi.org/10.1016/S1364-6826(00)00067-5 Medvedev, A. S., & Klaassen, G. P. (2000). 基于非线性波相互作用的重力波动量沉积参数化:基本公式和灵敏度测试。JASTP, 62, 1015-1033. https://doi.org/10.1016/S1364-6826(00)00067-5
Medvedev, A. S., Nakagawa, H., Mockel, C., Yiğit, E., Kuroda, T., Hartogh, P., et al. (2016). Comparison of the Martian thermospheric density and temperature from IUVS/MAVEN data and general circulation modeling. Geophysical Research Letters, 43(7), 3095-3104. 梅德韦杰夫,A. S.,中川,H.,莫克尔,C.,伊吉特,E.,黑田,T.,哈托赫,P.,等(2016)。基于 IUVS/MAVEN 数据和一般环流模型的火星热层密度和温度比较。《地球物理研究快报》,43(7),3095-3104。https://doi.org/10.1002/2016GL068388
Medvedev, A. S., & Yiğit, E. (2019). Gravity waves in planetary atmospheres: Their effects and parameterization in global circulation models. Atmosphere, 10(9), 531. Medvedev, A. S., & Yiğit, E. (2019). 行星大气中的重力波:它们的影响及在全球环流模型中的参数化。大气,10(9),531。https://doi.org/10.3390/atmos10090531
Medvedev, A. S., Yiğit, E., Hartogh, P., & Becker, E. (2011). Influence of gravity waves on the Martian atmosphere: General circulation modeling. Journal of Geophysical Research, 116(E10004). Medvedev, A. S., Yiğit, E., Hartogh, P., & Becker, E. (2011). 重力波对火星大气的影响:一般环流建模。地球物理研究杂志, 116(E10004)。https://doi.org/10.1029/2011JE003848
Millour, E., Forget, F., Spiga, A., Vals, M., Zakharov, V., Montabone, L., et al. (2018). The mars climate databaseScience workshop “from Mars express to ExoMars”. Retrieved from https://www.cosmos.esa.int/web/mars-science-workshop-2018
Nakagawa, H., Terada, N., Jain, S. K., Schneider, N. M., Montmessin, F., Yelle, R. V., et al. (2020). Vertical propagation of wave perturbations in the middle atmosphere on Mars by MAVEN/IUVS. Journal of Geophysical Research: Planets, 125, e2020JE006481. Nakagawa, H., Terada, N., Jain, S. K., Schneider, N. M., Montmessin, F., Yelle, R. V., et al. (2020). 火星中层大气中波动扰动的垂直传播,MAVEN/IUVS. 地球物理研究杂志:行星, 125, e2020JE006481.https://doi. org/10.1029/2020JE006481
Sen, P. K., (1968). Estimates of the regression coefficient based on Kendall’s Tau. Journal of the American Statistical Association, 63, 13791389. Sen, P. K., (1968). 基于肯德尔τ的回归系数估计。美国统计协会杂志, 63, 1379-1389.https://doi.org/10.2307/2285891
Siddle, A., Mueller-Wodarg, I., & Bruinsma, J.-C. M. (2020). Density structures in the Martian lower thermosphere as inferred by Trace Gas Orbiter accelerometer measurements. Icarus. https://doi.org/10.1016/j.icarus.2020.114109 Siddle, A., Mueller-Wodarg, I., & Bruinsma, J.-C. M. (2020). 通过追踪气体轨道器加速度计测量推断的火星下热气层的密度结构。Icarus. https://doi.org/10.1016/j.icarus.2020.114109
Spiga, A., Teitelbaum, H., & Zeitlin, V. (2008). Identification of the sources of inertia-gravity waves in the Andes Cordillera region. Annales Geophysicae, 26(9), 2551-2568. https://doi.org/10.5194/angeo-26-2551-2008 Spiga, A., Teitelbaum, H., & Zeitlin, V. (2008). 安第斯山脉地区惯性重力波源的识别。地球物理年鉴, 26(9), 2551-2568. https://doi.org/10.5194/angeo-26-2551-2008
Starichenko, E. D., (2021). Gravity wave activity in the Martian atmosphere at altitudes 20-160 km from ACS/TGO occultation measurements. https://doi.org/10.17632/gk6sr9fwzp. 1 Starichenko, E. D., (2021). 火星大气中 20-160 公里高度的引力波活动来自 ACS/TGO 掩星测量。 https://doi.org/10.17632/gk6sr9fwzp. 1
Terada, N., Leblanc, F., Nakagawa, H., Medvedev, A. S., Yiğit, E., Kuroda, T., et al. (2017). Global distribution and parameter dependences of gravity wave activity in the Martian upper thermosphere derived from MAVEN/NGIMS observations. Journal of Geophysical Research, 122(2), 2374-2397. https://doi.org/10.1002/2016ja023476 Terada, N., Leblanc, F., Nakagawa, H., Medvedev, A. S., Yiğit, E., Kuroda, T., et al. (2017). 从 MAVEN/NGIMS 观测中得出的火星上层热层重力波活动的全球分布和参数依赖性。地球物理研究杂志, 122(2), 2374-2397. https://doi.org/10.1002/2016ja023476
Theil, H. (1950). A rank-invariant method of linear and polynomial regression analysis. Proceedings of the Koninklijke Nederlandse Akademie van Wetenschappen, 53, 386-392. Theil, H. (1950). 一种秩不变的线性和多项式回归分析方法。荷兰皇家科学院院刊, 53, 386-392。
Tikhonov, A. N., & Arsenin, V. Y. (1977). Solutions of ill-posed problems. V. H. Winston.
Tolson, R., Bemis, E., Hough, S., Zaleski, K., Keating, G., Shidner, J., et al. (2008). Atmospheric modeling using accelerometer data during Mars Reconnaissance Orbiter aerobraking operations. Journal of Spacecraft and Rockets, 45. https://doi.org/10.2514/1.34301 Tolson, R., Bemis, E., Hough, S., Zaleski, K., Keating, G., Shidner, J., et al. (2008). 使用加速度计数据进行大气建模,研究火星侦察轨道器的气动制动操作。《航天器与火箭杂志》,45。https://doi.org/10.2514/1.34301
Tolson, R. H., Keating, G. M., George, B. E., Escalera, P. E., & Werner, M. R., (2005). Application of accelerometer data to Mars Odyssey aerobraking and atmospheric modeling. Journal of Spacecraft and Rockets, 45. https://doi.org/10.2514/1.15173 Tolson, R. H., Keating, G. M., George, B. E., Escalera, P. E., & Werner, M. R., (2005). 加速度计数据在火星奥德赛气动制动和大气建模中的应用. 航天器与火箭杂志, 45. https://doi.org/10.2514/1.15173
Vals, M., Spiga, A., Forget, F., Millour, E., Montabone, L., & Lott, F. (2019). Study of gravity waves distribution and propagation in the thermosphere of Mars based on MGS, ODY, MRO and MAVEN density measurements. Planetary and Space Science, 178, 104708. https:// doi.org/10.1016/j.pss.2019.104708 Vals, M., Spiga, A., Forget, F., Millour, E., Montabone, L., & Lott, F. (2019). 基于 MGS、ODY、MRO 和 MAVEN 密度测量的火星热层重力波分布和传播研究。行星与空间科学, 178, 104708. https:// doi.org/10.1016/j.pss.2019.104708
Whiteway, J., & Carswell, A. (1995). Lidar observations of gravity wave activity in the upper stratosphere over Toronto. Journal of Geophysical Research, 100(D7), 14113-14124. https://doi.org/10.1029/95jd00511 Whiteway, J., & Carswell, A. (1995). 多伦多上层平流层重力波活动的激光雷达观测. 地球物理研究杂志, 100(D7), 14113-14124. https://doi.org/10.1029/95jd00511
Withers, P. (2006). Mars Global Surveyor and Mars odyssey accelerometer observations of the Martian upper atmosphere during aerobraking. Geophysical Research Letters, 33. https://doi.org/10.1029/2005GL024447 Withers, P. (2006). 火星全球探测器和火星奥德赛加速度计对火星上层大气在气动制动期间的观测。地球物理研究快报, 33. https://doi.org/10.1029/2005GL024447
Wright, C. J. (2012). May). A one-year seasonal analysis of Martian gravity waves using MCS data. Icarus, 219(1), 274-282. https://doi.org/10.1016/j.icarus.2012.03.004 Wright, C. J. (2012). 五月). 一年季节性分析火星重力波的 MCS 数据. Icarus, 219(1), 274-282. https://doi.org/10.1016/j.icarus.2012.03.004
Yiğit, E., England, S. L., Liu, G., Medvedev, A. S., Mahaffy, P. R., Kuroda, T., & Jakosky, B. M. (2015). High-altitude gravity waves in the Martian thermosphere observed by MAVEN/NGIMS and modeled by a gravity wave scheme. Geophysical Research Letters, 42(21), 8993-9000. https://doi.org/10.1002/2015GL065307 Yiğit, E., England, S. L., Liu, G., Medvedev, A. S., Mahaffy, P. R., Kuroda, T., & Jakosky, B. M. (2015). 火星热层中的高空重力波通过 MAVEN/NGIMS 观测并通过重力波模型进行建模。地球物理研究快报, 42(21), 8993-9000. https://doi.org/10.1002/2015GL065307
Yiğit, E., & Medvedev, A. S. (2015). Internal wave coupling processes in Earth’s atmosphere. Advances in Space Research, 55(4), 983-1003. https://doi.org/10.1016/j.asr.2014.11.020 Yiğit, E., & Medvedev, A. S. (2015). 地球大气中的内部波耦合过程。空间研究进展, 55(4), 983-1003. https://doi.org/10.1016/j.asr.2014.11.020
Yiğit, E., & Medvedev, A. S. (2019). Obscure waves in planetary atmospheres. Physics Today, 6, 40-46. https://doi.org/10.1063/PT.3.4226 Yiğit, E., Medvedev, A. S., Benna, M., & Jakosky, B. M. (2021). Dust storm-enhanced gravity wave activity in the Martian thermosphere observed by MAVEN and implication for atmospheric escape. Geophysical Research Letters, 48. https://doi.org/10.1029/2020GL092095 Yiğit, E., Medvedev, A. S., & Hartogh, P. (2018). Influence of gravity waves on the climatology of high-altitude Martian carbon dioxide ice clouds. Annales Geophysicae, 36(6), 1631-1646. https://doi.org/10.5194/angeo-36-1631-2018 Yiğit, E., & Medvedev, A. S. (2019). 行星大气中的模糊波。今日物理,6,40-46。https://doi.org/10.1063/PT.3.4226 Yiğit, E., Medvedev, A. S., Benna, M., & Jakosky, B. M. (2021). 火星热层中尘暴增强的重力波活动由 MAVEN 观测到及其对大气逃逸的影响。地球物理研究快报,48。https://doi.org/10.1029/2020GL092095 Yiğit, E., Medvedev, A. S., & Hartogh, P. (2018). 重力波对高空火星二氧化碳冰云气候学的影响。地球物理年鉴,36(6),1631-1646。https://doi.org/10.5194/angeo-36-1631-2018