Quantitative Analysis of the Vertical Interactions between Dust, Zonal Wind, and Migrating Diurnal Tide on Mars and the Role of Gravity Waves

: In the atmospheric system of Mars, vertical interactions are crucial, yet quantitative studies addressing this issue remain scarce. Based on simulations using the Mars PCM-LMDZ, we pre-sent the first frequency-domain quantitative analysis of the vertical interactions among Martian atmospheric dust, zonal circulation, and the migrating diurnal tide (DW1), employing Partial Directed Coherence (PDC) techniques to quantify the strength of associations between different variables. Our findings reveal a chain of influence where sub-seasonal-scale dust signals in the troposphere, through the Doppler effect of middle atmospheric zonal winds, transmit modulated energy to the DW1 in the upper mesosphere, thereby facilitating interlayer atmospheric interactions. The radiative heating from dust activities enhances the residual mean meridional circulation, which, under the influence of the Coriolis force, further accelerates the westerlies. Although gravity wave activity also contributes to the acceleration of the westerlies, its forcing generally remains below 5 m/s, which is relatively weak compared to the impact of intense dust activities in the warm scenario experiments (approximately 20 m/s). Overall, this study quantifies the interactions among atmospheric layers by means of PDC technology and analytically demonstrates how dust energy is transferred to mesospheric tides by shaping the zonal winds in between.


Introduction
Substantial qualitative research has been conducted on the interlayer interactions in the Martian atmosphere, but quantitative studies are less common.Several observationbased studies suggest that atmospheric layer interactions on Mars are at least as significant as those on Earth, if not stronger [1][2][3].The importance of atmospheric layer interactions on Mars is attributed to the more extreme weather conditions, such as a lower atmospheric pressure, a rougher terrain, a smaller planetary radius, stronger zonal winds, and more persistent diurnal slope flows [4].Martian mid-atmospheric interactions with the lower atmosphere are facilitated through waves (e.g., gravity waves, planetary waves, and tides), dust, circulation, and atmospheric seasonal expansion/contraction related to the Martian carbon cycle.Vertical interactions on Mars lead to the redistribution of chemical species (such as water and ozone), thermal structures (such as polar warming), and mesoscale clouds (such as water and carbon dioxide).While these conclusions have a qualitative physical basis, a quantitative frequency-domain analysis is necessary to enhance our understanding of the strength and periodicity of interlayer interactions in the Martian atmosphere.
The lack of quantitative frequency-domain analyses in studies of vertical interactions has left the frequency-dependent interactions between dust, migrating diurnal tide (DW1), and zonal circulation inadequately explored.The meridional circulation in the Martian atmosphere typically receives more attention than the zonal circulation, primarily due to the more observable dynamical phenomena, including the global-scale Hadley circulation [1,3,[5][6][7] and its thermodynamic effects in the thermosphere through polar warmings [8][9][10].
Zonal winds have also been found to be critical in the vertical interactions of the Martian atmosphere, possibly playing a mediating role in dynamical processes.As a major momentum-containing component of atmospheric circulation [5], zonal winds are crucial for distributing heat and momentum.They modulate DW1 through the Doppler frequency shift effect [2,11]; affect the outbreak, propagation, and development of dust storms [7,12]; and absorb momentum and energy from gravity wave drag [13,14], thereby filtering gravity waves to specific frequencies [15].Gravity waves, dust, and DW1, in turn, influence zonal winds.The convergence of tidal heat and momentum fluxes exerts a strong forcing on zonal winds [16][17][18][19].These forcings on Mars are stronger than on Earth, especially during the Martian dust storm season, when the zonal-mean convergence of tidal fluxes is enhanced globally by dust radiative forcing [20].Furthermore, the damping effect of gravity wave drag is closely related to zonal winds [21,22], with Gilli finding particularly strong gravity wave drag where there is a large difference between the mean flow speed and the given phase speed (the speed of the flow being decelerated/accelerated), such as the upper part of the zonal winter jet stream [22].Miyamoto suggests that zonal winds are significantly influenced by planetary-scale dust events, transmitting their effects from the lower to the middle atmosphere [23].
These vertical interactions are typically nonlinear and challenging to discern in the time domain, where significant phenomena may not be readily apparent.In contrast, frequency-domain analysis excels at isolating unrelated signals, thereby enhancing the visibility of the interactions of interest.This analytical shift allows for a clearer observation of the dynamics at play, providing a more nuanced understanding of the complex interactions between atmospheric layers.However, frequency-domain analyses of dust, DW1, and zonal circulation have not yet been conducted.
Consequently, while traditional views have qualitatively demonstrated the vertical interactions, we will quantitatively address the frequency domain, conditions, and strength of interactions between signals of dust and DW1 and the specific role of zonal winds in this context.Based on the experiments conducted using the Martian planetary climate model of the Laboratoire de météorologie dynamique Zoom (Mars PCM-LMDZ) [24], our Partial Directed Coherence (PDC) analysis revealed sub-seasonal-scale dust signals caused by variations in dust activity.Here, the sub-seasonal scale indicates signals with periods shorter than a Martian season [24], which lasts approximately 193.47 sols in spring, 178.52 sols in summer, 142.77 sols in fall, and 154.24 sols in winter.Through zonal winds as an intermediary, sub-seasonal dust signals are transmitted to the DW1.When the background dust content is high enough, the dust signal strengthens to be equally important to the semiannual oscillation's signal.Gravity wave drag can reduce zonal wind signals in the sub-seasonal range, countering their modulation by dust, but it is not significant on semiannual and longer scales.In summary, in the vertical interactions between dust and DW1, zonal winds are a crucial link in signal transmission, while gravity waves counteract the modulation of zonal winds by dust by decelerating the winds.
Addressing the current shortfall in research on vertical interactions within the Martian atmospheric structure, we conduct targeted experiments using the Mars PCM-LMDZ to investigate the interactions between dust, DW1, and zonal circulation.Through the application of PDC analysis, this study explores sub-seasonal-scale dust signals induced by variations in dust activity.A focal point of this research is the transmission of dust signals to DW1 through zonal winds as the intermediary medium, alongside an examination of the effects of gravity waves on this interaction.We also discuss the roles that zonal winds and gravity waves play in these vertical interactions.

Experimental Design
In this study, we conducted six Martian atmospheric simulation experiments spanning 10 Martian years each, utilizing the Mars PCM-LMDZ model to simulate the interactions between atmospheric layers.These simulations were based on three initial conditions (cold, warm, and climatology states) derived from the widely used Martian Climate Database [25] (MCD).By toggling the presence of gravity waves (with/without gravity wave influences) through the control of the gravity wave module, we aimed to simulate the interactions of signals between atmospheric layers.
Secondly, drawing on the simulation outcomes, we employed generalized Partial Directed Coherence [26,27] (gPDC) to identify correlations between the dust content within the troposphere (<~40 km) and zonal winds in the upper troposphere and lower mesosphere (UTLM, 20~50 km), as well as DW1 within the middle atmosphere (50~100 km).This exploration seeks to characterize the representation of layer interactions within the frequency domain by comparing these dynamics against inherent atmospheric phenomena.Considering the variable DW1 as  1 (), dust as  2 (), and zonal winds as  3 (), we used Equation (1) with the input of () = [ 1 (),  2 (),  3 ()]  , while   () and   () are adjusted accordingly (detailed in Section 2.3).All data were vertically weighted by mass, with a meridional weighting factor of  , where  is the latitude before conducting the PDC analysis.Then, utilizing a vector autoregressive (VAR) [3] process of gPDC, we aimed to elucidate the interconnections among Martian atmospheric processes (dust, zonal winds, and DW1).Herein, gPDC calculations were conducted using the AsympPDC version 3.0 package in the MATLAB toolbox.

Numerical Simulation Experiments
In our study, we performed six global simulations using the Mars PCM-LMDZ model, encompassing a grid with 64 × 48 × 29 points across two horizontal and one vertical dimension (with the model top at approximately 102.5 km above the surface).The resolution was as follows: spatial resolution at 5.625° longitude by 3.75° latitude, with 29 hybrid sigma-pressure layers; the temporal resolution was set at every 1.5 Martian hours.The core version of the model was LMDZ5, which is more stable, reliable (verified through years of usage), and operationally efficient than the newest LMDZ6, and it was enough for our experiments.
The initial conditions for the six simulations were based on three dust scenariosclimatology, cold, and warm dust backgrounds-as modeled for the seasonal evolution and spatial variability of dust column opacity in the MCD [25].The standard 'climatology' dust scenario synthesizes observations from Martian years without global dust storms (i.e., MY 24,26,27,29,30,and 31), creating a composite scenario from these observationally average years.In the 'cold' scenario, dust opacity at each location is set to 50% of the minimum values observed during MY 24-31, representing an exceptionally clear atmosphere.The 'warm' scenario corresponds to atmospheric layers with dust but without global dust storms, where opacity due to sunlight at each location and throughout the year is set to the maximum values observed during MY 24-31.Dust observations for these Martian years in MCD v6.1 were compiled by Montabone [28].
The terrain used to generate orographic gravity waves in the PCM was smoothed to the model's horizontal resolution.Sub-grid-scale mountain effects were parameterized along with low-level drag schemes [24,29] and gravity wave drag schemes [30,31].These schemes assume that all observed stress reductions above mountaintops are due to wave breaking, although, at times, partially trapped non-orographic gravity waves also explain the stress reduction.Furthermore, it is assumed that the area-averaged momentum flux is comparable to the line average measured above the central cross-section.Thus, the scheme's simulation of momentum flux aligns qualitatively with observations.The gravity wave drag scheme employs a single gravity wave hypothesis, with gravity waves propagating vertically from the model surface in the vertical plane.We adopted a stochastic parameterization scheme for non-orographic gravity waves on the Mars PCM-LMDZ as described by Liu [32], which is an improvement upon the original scheme described by Gilli [22].This scheme assumes that the sources are located above typical convective cells (~250 Pa) and operate continuously.
Overall, we conducted six experiments ('cold', 'warm', and 'climatic state' scenarios with and without gravity waves, Experiments 1 to 6) using the Mars PCM-LMDZ, as outlined in Table 1.This detailed configuration enables a comprehensive analysis of Martian atmospheric dynamics by contrasting conditions across a range of dust scenarios and gravity wave influences, providing insights into the complex interactions within the Martian atmosphere.

Partially Directed Coherence Analysis
Mars's extreme weather, such as lower atmospheric pressure and stronger zonal winds, suggests more robust vertical interactions than Earth.These interactions, involving waves, dust, and seasonal atmospheric changes, affect the distribution of chemical species and thermal structures.Quantitative studies on interlayer interactions in the Martian atmosphere are scarce despite qualitative research indicating their significance.To address this gap, our study employed PDC analysis on data from the Mars PCM-LMDZ model.PDC is well-suited for this task as it effectively isolates and quantifies directional influences between multiple time series data in the frequency domain.This method allows us to discern complex, nonlinear interactions that are difficult to detect in the time domain.By employing PDC analysis, we quantitatively elucidated the frequency, conditions, and strength of interactions between dust, DW1, and zonal winds, enhancing our understanding of Martian atmospheric dynamics.
PDC is a multivariate time series technique derived from the factorization of partial coherence [33].Research indicates that PDC captures certain types of nonlinear interactions [34,35].It has been successfully applied in highly nonlinear time series data in fields such as neuroscience [36,37] and meteorology [38].The definition of causality related to the predictability of two or more different processes is currently known as Granger causality [34].Granger causality can determine the direction and spectrum of interactions.For instance, interactions between the stratosphere and the tropical troposphere on Earth have been studied based on Granger causality [38].PDC, derived from the factorization of partial coherence [33], closely embodies the concept of Granger causality in the frequency domain [34].When analyzing interactions and identifying correlations between different spectral characteristic variables (e.g., dust and wind on Mars) across atmospheric layers, gPDC is an effective tool [27].PDC's multivariate analysis evaluates each possible direction of interaction between the mesosphere and the troposphere, revealing influences received or transmitted by the atmosphere, thereby even identifying feedback effects.Notably, PDC shares duality with Granger causality, where a zero   () corresponds to the absence of a connection, similar to the conditions of Granger causality.
The following definitions were used for the three time series variables in this study, i.e., DW1  1 (), dust  2 (), and zonal winds  3 ().Here, we consider a vector-valued signal using a p-order vector autoregressive model, VAR(p), defined by the following equation: where   () are the VAR(p) coefficients representing the k-order lagged effect of   () on   (), with t denoting the time series variable.The random component   () is a zeromean novelty process with a covariance matrix  = [  ]: If   () ≠ 0 and is statistically significant, then   () is a Granger cause of   ().This causality implies that   () indeed helps predict   (), and vice versa.The VAR(p) model was used to formulate and test the null hypothesis.
0 :   () = 0,  = 1, . .,   1 :  ∈ {1, . ., },   () ≠ 0 If   () in the past helps predict   () in the future, then   () is considered a Granger cause of   ().Again, we consider a binary time series with a p-order vector autoregressive process VAR(p), as defined in Equation ( 1), where coefficients   () describe the lagged effects of   () on   ().The frequency-domain representation of A(k) can be denoted as  ̄(), with where   is the Kronecker delta,   = √−1, p is the model order, and  is the Fourier frequency (Hz).The gPDC information flow of   () to   () at a specific frequency  is denoted as   (), indicating where  is   −1/2 , the columns of  ̄() are denoted as   (),    () represents the Hermitian transpose, and  is (  ⊙ ) −1 for PDC, where   is a  ×  identity matrix (with  = 3 being the variable number of signals), and  is the covariance matrix of random components.The 95% confidence interval is computed based on the asymptotic normal distribution of |  ()| 2 : where   is the number of available observations, and the frequency-related  2 =        .When the covariance matrix  is given as being a priori known from the previous [  ], it is known that and For in which   and   are  2 ×  2 matrices.The vector   is given by where Here,   is an N×N identity matrix,   has coefficients equal to 1 only at ( , ) = ((  − 1)) + , (  − 1) + ) , and   has non-zero coefficients only at (, ) = (  − 1) + 1 ≤  =  ≤ .The operator  denotes the matrix-to-column stacking operator, and () had been defined in Equation ( 1).
The significance level is calculated through the convergence of the following distribution: where  1 and  2 are weights that can be estimated from the data, and  1 and  2 are two independent  1 2 distributed random variables, presented under the null hypothesis  0 .PDC analysis, while sometimes underperforming in the context of linear interactions [35], has been shown to effectively reflect the characteristics of nonlinear interactions despite its various limitations with long time intervals or non-Gaussian processes.PDC operates similarly to the concept of the Mutual Information Rate (MIR) in Gaussian processes, elucidating all information flows between two or more signals, thereby quantitatively reflecting the mutual dependencies between the signals.Conversely, in non-Gaussian scenarios, a boundary exists between the estimated MIR by PDC and the actual MIR, within which the differences are deemed acceptable, allowing PDC results to reflect real interactions to a certain extent [39].Moreover, PDC has been successfully applied in the analysis of highly nonlinear time series data in fields such as neuroscience, demonstrating its efficacy in revealing the intricate dynamics of complex systems [37].Thus, PDC is also suitable for analyzing nonlinear problems like the vertical interactions in the Martian atmosphere, offering a quantitative depiction of how different layers within the Martian atmosphere interact and influence each other.

Numerical Simulation Experiment Results
Before proceeding with the PDC analysis, we first examined the general distribution of dust, DW1, and zonal circulation.As outlined in the experimental design, our primary focus was on the distribution of tropospheric dust, the zonal winds in the UTLM, and DW1 in the mesosphere.
For DW1, the simulation results within both solar longitude (Ls) ranges are broadly consistent with previous studies.In the 'climatology' scenario, for Ls = 180°-210° (Figure 1i) and Ls = 270°-300° (Figure 2i), the morphology, structure, and amplitude distribution of DW1 closely resemble those depicted in Figure 2c of Forbes et al. [2], Figure 5 of Lee et al. [40], and Figure 4 of Scott et al. [41].It is illustrated that during the southern spring, the equatorial amplitude maxima of DW1 are located around 70 km, with smaller non-equatorial maxima at high latitudes at the same altitude (Figure 1).The seasonal variation in the height of the amplitude maxima of DW1, as observed in the MCS data (Wu et al. [42]), is reproduced by the model, with the springtime maxima in the Southern Hemisphere (Figure 2i) being lower than the summertime maxima (Figure 1i).Additionally, the modeled distribution of dust aligns well with observations.At Ls = 180°-210° and 270°-300°, dust development in the Southern Hemisphere reaches higher altitudes compared to the Northern Hemisphere (Figures 1g  and 2g), the latter experiencing more vigorous growth due to stronger solar radiation.Our findings also indicate that under a 'warm' dust background (Figure 1m,p), regardless of the presence of gravity waves, the average dust mass mixing ratio can reach nearly 5 × 10 −6 kg/kg at 40 km altitude.The dust reaching the upper levels of the tropical troposphere imposes significant thermodynamic forcing on the circulation and tidal patterns, as supported by [43,44].

Vertical Interactions among Dust, Zonal Winds, and DW1
In the 'warm' scenario with a high background dust content, the energy of zonal winds is modulated by dust due to the impact of temperature changes caused by dust absorption of solar radiation.The resulting zonal winds allow DW1 of specific frequency ranges to propagate upwards, filtering out other tides in the process.During this filtering, energy from the zonal winds is transferred to the tides, a portion of which originally stems from dust.Filtering occurs regardless of dust presence, but dust alters the effectiveness of filtering within certain frequency ranges, thereby facilitating interactions across different atmospheric layers.To study the vertical interactions between dust, zonal winds, and DW1, we employ the PDC technique for frequency domain analysis.
In the Martian atmosphere, the primary sources of thermal tide generation are dust and water ice clouds, with dust being the more significant source due to the sparse and cold nature of the Martian atmosphere that limits the formation of water ice clouds [45].Moreover, the forcing effect of dust on zonal winds and diurnal migrating tidal wave DW1 is stronger than their mutual opposing effects.Our study thus focuses on the interactions among dust, zonal winds, and tides beyond seasonal factors within the same frequency range.We identify energy transfer from dust to zonal winds and tides, with sub-seasonal scales (periods less than 200 Martian days) of dust energy successively transferred to zonal winds and DW1.To corroborate the specific relationship with dust, controlled variable experiments were conducted.An increase in the background dust content led to a simultaneous increase in the Power Spectral Density (PSD) of dust, zonal winds, and DW1 at sub-seasonal scales (Figure 3).The increase is more pronounced without disturbances (Figure 4), as gravitational wave drag partially offsets dust energy, reducing its impact on zonal winds.Besides sub-seasonal dust energy, there are also annual and semi-annual interactions and energy peaks, primarily driven by Mars's orbital and climatic system rhythms, which are not the focus of this study.
Figures 3 and 4 illustrate the distribution of PDC values and PSD in the frequency domain for experiments with and without gravity waves, respectively.The middle layer's zonal winds act as an intermediate bridge in the energy transfer chain, conveying energy upwards while retaining little themselves.In high dust background scenarios, multiple energy peaks at sub-seasonal scales were observed for dust and DW1, with zonal winds showing few such peaks.Despite absorbing minimal energy at sub-seasonal scales, zonal winds significantly affect DW1 at these frequencies, acting like an energy transfer bridge.The convective layer's dust energy cannot directly reach the middle layer's DW1 energy; an intermediate energy transfer bridge is required to convey energy upwards, with zonal winds serving as that bridge.Without gravity wave interference, the PSD consistency between dust and DW1 at sub-seasonal frequencies improves.According to the theory of the Doppler filtering effect, zonal winds modulate the wavelength of thermal tidal waves, facilitating their dissipation when the propagating directions coincide and aiding their upward propagation when the directions oppose each other.Thus, mid-layer zonal winds convey sub-seasonal-scale dust storm signals from the lower atmosphere (below 40 km) to mid-layer tides (50-100 km) via Doppler shifting the upward-propagating DW1's frequency.With a higher background dust content, enhanced dust radiative heating transfers more sub-seasonal-scale dust signal energy upwards, making this impact chain more pronounced.This is observed in the PDC value, which indicates the consistent effects of dust on zonal winds and of zonal winds on DW1 at sub-seasonal scales.

Mechanism Analysis
In this study, we focus on the Southern Hemisphere's spring season (Ls = 180°-270°) to elucidate the entire chain of impacts, as the effects of chosen dust scenarios in our experiments are most pronounced during the dust seasons.
The zonal wind's filtering effect on DW1 can be observed by comparing the activity of lower-level DW1 (Figure 5c) against upper-level DW1 (Figure 5d).The unfiltered lowerlevel DW1 is predominantly driven by dust radiative heating wave sources, as evidenced by the significant rise in the average DW1 amplitude with an increasing background dust content ('Cold': 2.3 K, 'Climatology': 4.1 K, 'Warm': 6.1 K) (Figure 5c).However, DW1 reaching the upper levels is primarily influenced by mid-level zonal winds, rather than solely by lower-level dust heating.Prior to Ls = 180°-270°, the 'Warm' scenario exhibits the strongest westerlies, thereby shifting the westward-propagating DW1 to the highest frequencies (Figure 3l), aligning closely with the sub-seasonal-scale dust's intrinsic frequency (Figure 3o).Similarly, post-Ls = 210°, scenarios with the highest zonal wind speeds consistently correspond with the greatest DW1 magnitudes, as evidenced by the alignment of the maximum bars (ordinate axis labels) in Figure 5b,d.The Doppler shift effect is encapsulated by Equation ( 13): with σ denoting the Doppler frequency shift, Ω denoting the rotational rate of Mars,  denoting the zonal wind speed, θ denoting the co-latitude, and  0 denoting the equatorial wave phase speed.A simplified version of the dispersion relation for tidal waves in a rotating atmosphere [46] is where  is the horizontal wavenumber,  is the Coriolis parameter, and  is the buoyancy frequency.The vertical wavenumber  is related to the vertical wavelength   by   = 2/, so Equation ( 14) means that   is positively related to frequency .Given that Ω is constant and sin(θ) is always positive, if the zonal wind  and wave phase speed  0 are in opposite directions, DW1 will experience a shift to higher frequencies and longer vertical wavelengths upon traversing the UTLM region, ultimately resulting in greater amplitudes by the time it reaches the upper stratosphere.The first link in the impact chain is the influence of dust on zonal winds.As depicted in Figure 5, lower-level dust can only directly affect the lower-level DW1 and mid-level zonal winds, not directly impacting DW1 in the middle to upper mesosphere.Hence, the upward propagation of dust signals must first modulate the UTLM's zonal winds.Dust influences by enhancing the Hadley circulation, altering vorticity advection and thus modulating zonal winds, and providing an eastward torque to the equatorial mean zonal winds at the solstice [12,47], as shown in Equation (15).This transformation in the form of the Transformed Eulerian Mean [48] (TEM) equation estimates the impact of meridional (̅ * ), vertical circulation ( ̅ * ), and gravity waves on the acceleration of the mean zonal winds: The terms −̅ *  ̅   −1 and − ̅ *  ̅  represent the acceleration of westerlies due to meridional and vertical circulations [49,50].The impact of meridional circulation is at least two orders of magnitude greater than that of vertical circulation, as the average ̅ * in Figure 6 is almost two orders of magnitude larger than  ̅ * in Figure 7.At the onset of the Northern Hemisphere's autumn, westerlies dominate (Figure 5b); however, given that  ̅  is positive in the Northern Hemisphere's autumn [23] and  −1 is a constant, Equation (15) suggests that a larger near-equatorial ̅ * would result in stronger easterlies, evident from the concurrent rise in ̅ * in Figure 6a and easterly wind speeds in Figure 5b, correlating with increasing background dust levels.Furthermore, ̅ * shows a positive correlation with the background dust content, consistent with the enhancements in meridional circulation during periods of intense dust activity found in previous studies [50].Thus, dust's modulation of zonal winds primarily originates from dust radiative heating-induced changes in meridional winds.The influence of gravity waves on the acceleration of mean zonal winds, represented by (ρ 0 ) −1 ∇ ⋅ , can be estimated by directly controlling the gravity wave parameterization scheme in the model.Gravity waves facilitate easterly acceleration in the 'Climatological' and 'Warm' dust scenarios, while they decelerate easterlies in the 'Cold' scenario, though the forcing of gravity waves on the zonal wind field typically remains below 5 m/s, significantly lower than the forcing from dust on the zonal wind field (approximately 20 m/s), as shown in Figure 8. Notably, while the forcing from gravity waves plays a less significant role in shaping the zonal wind field, it markedly alters the ̅ * and  ̅ * in the 'Climatological' scenario, as depicted in Figures 6c and 7c  In summary, sub-seasonal-scale dust energy in the low-latitude regions, through thermodynamic effects, enhances the residual mean meridional circulation (̅ * ), strengthens the Hadley circulation, and alters vorticity advection.This provides an eastward torque to the equatorial mean zonal winds, augmenting easterlies and leading to a greater dissipation of DW1 in the UTLM region through the Doppler shift effect.Thus, although an increase in the background dust content boosts the excitation of lower-level thermal tidal fluctuations, the thermal tides that reach the middle to upper mesosphere are reduced.In this process, DW1 reaching upper levels undergoes more extensive modulation by zonal winds, with the modulating energy initially stemming from dust energy in the lower-to-middle troposphere.This outlines the complete chain of impacts among dust, zonal winds, and tides.

Discussion
This study quantitatively analyzes the vertical interactions between dust, zonal winds, and the migrating diurnal tide DW1 in the Martian atmosphere using a frequencydomain approach based on PCM-LMDZ simulations and Partial Directed Coherence (PDC) techniques.Our findings reveal an important chain of influence, where sub-seasonal-scale dust signals originating from the troposphere are transmitted upward through the modulation of middle atmospheric zonal winds, facilitating the transfer of energy to DW1 in the upper mesosphere.This mechanism highlights the crucial role of zonal winds in linking different atmospheric layers on Mars.
While previous studies have qualitatively demonstrated such vertical coupling, our work provides quantitative insights into the frequency dependence, conditions, and relative strengths of these interactions.Notably, the influence of dust on DW1 becomes increasingly pronounced under high dust-loading scenarios, rivaling the prominent semiannual oscillation signal.Gravity waves, on the other hand, exert a relatively weak forcing on zonal winds (below 5 m/s), partially counteracting the modulating effects of dust (approximately 20 m/s) at sub-seasonal scales.This distinction underlines dust significance and the varying scales of impact that different atmospheric phenomena have on Martian atmospheric dynamics.
Despite these significant findings, this study has several limitations that should be addressed in future research.First, our simulations rely on parameterized representations of dust activities [28] and gravity wave processes [32], which may not fully capture the complexity of these phenomena.Incorporating more observational constraints and advanced parameterization schemes could improve the model's fidelity.Second, our analysis primarily focuses on the interactions between dust, zonal winds, and DW1, while other secondary potential mechanisms, such as the role of planetary waves, are not considered.A more comprehensive investigation of various coupling pathways is warranted.Furthermore, our study primarily examines the general behavior and averages over latitudes.However, it is crucial to recognize that the Martian atmosphere exhibits longitudinal variations due to topographical characteristics, which may also influence the strength and patterns of vertical interactions.Although topography is not the focus of this study, future studies could explore these spatial and temporal differences in detail, potentially revealing new insights.
Looking ahead, several promising research directions emerge.Combining our modeling efforts with enhanced observational data from ongoing and future Mars missions could lead to improved parameterizations and validation of the simulated processes.Additionally, exploring the impacts of other atmospheric coupling mechanisms, such as planetary waves, on zonal winds, tides, and dust dynamics would provide a more holistic understanding of the Martian atmospheric system.Moreover, investigating the scale dependence of atmospheric layer coupling could yield valuable insights.Our study primarily focused on sub-seasonal scales, but interactions at longer timescales, such as interannual or decadal variations, may exhibit distinct characteristics and influence large-scale climate patterns on Mars.Conversely, examining the coupling mechanisms at smaller spatial and temporal scales could unravel the intricate dynamics governing mesoscale phenomena like dust devils and slope winds.However, scaling up requires support from longer-duration observational data, while scaling down demands higher-resolution data.
In conclusion, this study represents a significant step forward in quantifying the vertical interactions within the Martian atmosphere.By employing advanced frequency-domain techniques, we elucidated the crucial role of zonal winds in facilitating energy transfer from dust activities to the migrating diurnal tide.However, our findings also underscore the need for further research to address remaining uncertainties, incorporate additional coupling mechanisms, and explore the scale dependence of these interactions.Such efforts will not only deepen our understanding of the Martian climate system but also contribute to the broader field of comparative planetary atmospheric dynamics.

Conclusions
This work addresses the existing shortcomings in quantitative studies on vertical interactions within the Martian atmospheric structure by conducting targeted experiments with the Mars PCM-LMDZ.Through PDC analysis, it explores sub-seasonal-scale dust signals induced by variations in dust activity.It focuses on the process by which dust signals are transmitted to DW1 through zonal winds as the intermediary medium and examines the effects of gravity waves on this interaction.The roles of zonal winds and gravity waves in these vertical interactions are also discussed.
Our findings emphasize the importance of employing frequency-domain analysis to quantify vertical interactions in studies of Martian atmospheric dynamics.Through detailed analysis of the interactions among dust, DW1, and zonal circulation, this study not only advances our understanding of the dynamic interplay between different atmospheric layers on Mars but also lays a novel methodological foundation for future Martian climate modeling.Moreover, we showcase the potential of PDC techniques in analyzing complex atmospheric systems' nonlinear interactions, offering a valuable tool for similar studies in other planetary atmospheres.

Figure 1 .
Figure 1.Vertical profiles of dust, zonal winds, and DW1 amplitude during the first Martian year at Ls = 180°-210° across Experiments 1 to 6.The profiles include the average dust mass mixing ratio (10 −6 kg/kg, left), zonal winds (m/s, middle), and amplitude of the diurnal migrating tide DW1 (K, right), with (a-c) showing results from EXP1, (d-f) from EXP2, (g-i) from EXP3, (j-l) from EXP4, (m-o) from EXP5, and (p-r) from EXP6.The red dashed boxes in the images highlight the data ranges selected for the PDC analysis: dust at 5-40 km within ±60° latitude, zonal winds at 20-50 km within ±30° latitude, and DW1 at 50-97.5 km within ±30° latitude.Temperatures below 50 km (contour lines) are indicated by the shaded area of dust mass mixing ratio in the left column, while zonal winds and tidal effects below 100 km are shown in the middle and right columns, respectively.Ls = 180° corresponds to the 372nd Martian sol, and Ls = 210° corresponds to the 422nd Martian sol.

Figure 2 .
Figure 2. Similar to Figure 1 but at Ls = 270°-300° during the first Martian year in Experiments 1 to 6. Ls = 270° corresponds to the 515th Martian sol, and Ls = 300° to the 562nd Martian sol.

Figure 3 .
Figure 3.The distribution of PDC values and PSD in the frequency domain in Experiments 1 (left), 3 (center), and 5 (right) involving gravity waves.The sequence from top to bottom (a-c) shows the PDC values from dust to zonal winds (u); (d-f) illustrate the PDC values from zonal winds to the diurnal migrating tidal wave DW1; (g-i) depict the PDC values from dust to DW1; (j-l) display DW1's PSD; (m-o) showcase the PSD of dust; (p-r) present the PSD of zonal winds.The data represented in these figures correspond to the averages within the areas of focus described in the experimental chapters.In the top three rows, red (green) lines indicate PDC values with (without) Granger statistical significance, with the shaded gray area representing the 99% confidence interval and the black curved dashed line marking the confidence lower limit.The bottom three rows feature blue lines for the PSD.Dashed vertical red and black lines denote semiannual and annual cycles, respectively.

Figure 7 .
Figure 7.These numerical simulation results are analogous to Figure 6 but focus on the residual mean vertical wind ( ̅ * ).

Figure 8 .
Figure 8. Numerical simulation results within the UTLM across 30°S-30°N for 'cold', 'climatological', and 'warm' scenarios, with/without gravity wave impacts.The legend indicates the change in zonal wind speeds from 180 Ls to 270 Ls.