ABSTRACT An objective analysis dataset of the Venus atmosphere, ALERA‐V version 1.0, has been released. This objective analysis provides the best estimate of the state of the Venus atmosphere produced by combining observations and a forecast model weighted according to their respective levels of uncertainty. The accumulation of frequent observations by the Venus Climate Orbiter ‘Akatsuki’ has, for the first time, enabled an objective analysis of the Venus atmosphere. ALERA‐V stands for AFES‐LETKF experimental ensemble objective (re)analysis of the Venus atmosphere, and is generated by ALEDAS‐V, the AFES‐LETKF data assimilation system for the Venus atmosphere. ALEDAS‐V uses AFES‐Venus (Atmospheric General Circulation Model for the Earth Simulator for Venus) as a forecast model and the LETKF (Local Ensemble Transform Kalman Filter) for data assimilation. The observations of the zonal and meridional winds obtained by a cloud tracking technique from images taken by Akatsuki's UVI (Ultraviolet Imager) are assimilated to produce ALERA‐V version 1.0. The dataset consists of atmospheric variables defined on a grid of 128 longitude and 64 latitude points, and 60 vertical levels, with output every 6 Earth hours from September to December 2018, including the intensive observation period of November. ALERA‐V is expected to be useful for both scientific and engineering research, such as understanding the dynamical mechanisms of various atmospheric phenomena and planning for future satellite missions, with appropriate consideration of its quality and limitations.

The formation of the stable layer below about 2 × 106 Pa pressure level (about 20 km altitude) of the atmosphere of Venus detected by in situ observations is investigated by the use of a radiative-convective equilibrium model. We demonstrate that, assuming mixing ratio profiles of absorbers to be at the upper limits of the observed ranges for H2O and SO2 and the lower limit for CO, a stable layer forms as a radiative-convective equilibrium state, but its stability is lower than the observed one. Also, increasing the continuum absorption coefficient of CO2 and/or H2O, which are not well constrained observationally or experimentally, results in the formation of a stable layer whose stability is comparable to the observed one. These results suggest a practical method to form the stable layer in the dynamical models of the Venus atmosphere. Further, these results indicate that the important targets of future observations and laboratory measurements are to obtain more precise profiles of the mixing ratios of H2O, CO, and SO2 in the Venus atmosphere, and to determine the continuum absorption coefficients of those.

A set of line-by-line and correlated k-distribution radiation models are developed aiming for applications in simulations and examinations of Venus and Mars-like planetary atmospheres. Our line-by-line model is validated by comparing the results with observations and those of previous studies under conditions of Venus, and present and possible early Mars. The radiation fields calculated by our line-by-line model agree well with observed profiles and are within the acceptable range from those presented in previous studies. The results obtained by our line-by-line model are then processed to generate a series of parameters for our correlated k-distribution model. It is confirmed that the radiation fields calculated with those sets of parameters by our correlated k-distribution model sufficiently agree with those by our line-by-line model for the atmospheres with a wide range of surface pressure. By the use of our correlated k-distribution model implemented with those sets of parameters, we evaluate the radiation field for Venus and calculate radiative-convective equilibrium profiles for Venus and Mars. The obtained vertical thermal structures for Venus are qualitatively consistent with observations, and the behaviors of surface pressure and surface temperature for Mars are similar to those reported by previous studies. Those results demonstrate that our models including the procedure for generating tables of radiation parameters are applicable to examine climates of CO2 dominant atmospheres in our solar and exoplanetary systems.

Dependence of the radiative-convective equilibrium structure of the lower atmosphere of Venus on the speci-fication of an atmospheric thermodynamic model is investigated. A series of thermodynamic models including ideal gases, van der Waals gases, and real gases are introduced by the use of the Helmholtz energy given by the EOS-CG mixture model (EOS-CG: Equation of State for Combustion Gases and Combustion Gas-like Mixtures). It is demonstrated that the radiative-convective equilibrium profile for the real gas differs significantly from that for the ideal gas with temperature-dependent specific heat by an increase of about 7 K in the surface temperature. This difference is caused by the fact that the adiabatic lapse rate evaluated with the thermodynamic model of real gas is larger than that of ideal gas, since the non-ideality of gas increases the thermal expansion coefficient, which overwhelms the increases in density and specific heat. It is confirmed that, in order to obtain better calcu-lations of atmospheric circulations including the lower atmosphere of Venus, the ideal gas with a constant spe-cific heat should be abandoned. The ideal gas with a temperature-dependent specific heat may not be enough. A promising method is to use the ideal gas but with the temperature-dependent specific heat such that its adiabatic lapse rate profile mimics that for the real gas.

<p>Observations of the Venus Climate Orbiter “Akatsuki” provide us with horizontal distributions of the horizontal winds derived from cloud tracking of the Ultraviolet Imager (UVI) and of temperature observed by the Longwave Infrared Camera (LIR). However, these observations are limited in altitude, local time (day or night side), and frequency. Then it is difficult to elucidate the general circulation of the Venus atmosphere, including various temporal and spatial scales, only from observations. In this study, we produced a Venus dataset (analysis) that has high temporal and spatial resolutions by assimilating horizontal winds derived by the Akatsuki observations. At the top of the cloud layer of Venus, there are planetary-scale atmospheric waves that are excited by the solar heating and move with the sun, called the thermal tides. In this presentation, we focused on thermal tides to verify the analysis.</p> <p>We use the Venus atmospheric data assimilation system “ALEDAS-V" (Sugimoto et al., 2017) [1] for assimilation and the Venus atmospheric general circulation model “AFES-Venus" (Sugimoto et al., 2014) [2] for ensemble forecasts. AFES-Venus is a full nonlinear dynamical GCM on the assumption of hydrostatic balance, designed for the Venus atmosphere. ALEDAS-V uses the Local Ensemble Transform Kalman Filter, and is the first data assimilation system for the Venus atmosphere. We assimilated the cloud top (~70km) zonal and meridional winds obtained by tracking morphology, using Akatsuki UVI data (Horinouchi et al., 2021) [3] from September 1st to December 31st, 2018. The assimilation data (analysis) from October 1st to November 30th, 2018, is analyzed, because the root-mean-square-deviations (RMSD) from FR (free run; the case without data assimilation) are stable.</p> <p>Figures (a) and (d) show the observed zonal and meridional winds, respectively. The zonal wind has a local minimum near 11 LT (local time) around the equator (Figure a). The meridional wind is the weakest at the equator and increases with latitude, and the amplitude is maximum around noon (Figure d) in the local time direction. Note that these winds obtained from observations exist only the dayside equatorward of 50° latitudes (Figure a and d).</p> <p>Figures (b) and (e) show the deviations from the zonal means of zonal and meridional winds at an altitude of 70 km in the FR, respectively. For zonal wind, diurnal (zonal wavenumber 1) and semidiurnal (zonal wavenumber 2) tides are dominant at latitudes poleward and equatorward of 30, respectively (Figure b). The zonal wind deviation has a local minimum at 14-15LT, which is ~ 2 hours behind the observation (Figures a and b). The meridional wind deviation is polar and equatorial on the dayside and nightside, respectively (Figure e), and this distribution is consistent with Akatsuki's observation (Figure d).</p> <p>Figures (c) and (f) show the zonal and meridional winds as a result of assimilation, respectively. The zonal wind in the equatorial region have a local minimum near 11 LT. The assimilation improved the semidiurnal tide closer to the observations (Figures a and c). The meridional wind is not so different from FR. This is probably because FR was originally very similar to observations (Figures d and f). These results are consistent with a previous study by Sugimoto et al. (2019) [4]. In addition, while the observed winds exist only on the dayside, the results of assimilation show that the horizontal winds field is modified significantly even on the nightside. It is suggested that spatially limited data assimilation can improve the general circulation of GCM.</p> <p>In the future work, we are planning to release the assimilation dataset as the “objective analysis data” of Venus for the first time in the world.</p> <p><img src="" alt="" width="743" height="482" /></p> <p>[1] Sugimoto, N., et al. Development of an ensemble Kalman filter data assimilation system for the Venusian atmosphere. Scientific Reports 7(1), 9321 (2017).</p> <p>[2] Sugimoto, N., et al. Baroclinic instability in the Venus atmosphere simulated by GCM. J. Geophys. Res. Planets 119, 1950–1968 (2014).</p> <p>[3] Horinouchi, T., et al. Venus Climate Orbiter Akatsuki Cloud Motion Vector Data Set v1.0, JAXA Data Archives and Transmission System (2021).</p> <p>[4] Sugimoto, N., et al. Impact of data assimilation on thermal tides in the case of Venus Express wind observation. Geophys. Res. Lett. 46, 4573–4580 (2019).</p>

The planetary missions including the Venus Climate Orbiter 'Akatsuki' provide new information on various atmospheric phenomena. Nevertheless, it is difficult to elucidate their three-dimensional structures globally and continuously only from observations because satellite observations are considerably limited in time and space. We constructed the first 'objective analysis' of Venus' atmosphere by assimilating cloud-top horizontal winds on the dayside from the equator to mid-latitudes, which is frequently obtained from Akatsuki's Ultraviolet Imager (UVI). The three-dimensional structures of thermal tides, found recently to play a crucial role in maintaining the super rotation, are greatly improved by the data assimilation. This result is confirmed by comparison with Akatsuki's temperature observations. The momentum transport caused by the thermal tides and other disturbances are also modified by the wind assimilation and agrees well with those estimated from the UVI observations. The assimilated dataset is reliable and will be open to the public along with the Akatsuki observations for further investigation of Venus' atmospheric phenomena.

Ishiwatari et al. (2007) (https://doi.org/10.1029/2006JD007368) investigated multiple equilibrium solutions of a gray atmosphere for various values of solar constant, utilizing two types of models, namely, a one‐dimensional energy balance model and a general circulation model (GCM) with simplified physical processes, both of which permit existence of the runaway greenhouse state. The study was retracted, since there was a bug in their GCM that affected quantitative aspects of the study. Here, we revise the study with re‐performing all of the GCM experiments using an appropriately corrected model. The results of re‐experiments show that the main features of the climate regime diagram drawn in the solar constant ‐ ice line latitude plane are mostly unchanged, except that the ice‐free equilibrium state, which existed in Ishiwatari et al. (2007) (https://doi.org/10.1029/2006JD007368) now disappears. It is confirmed that there are the partially ice‐covered state, the globally ice‐covered state, and the runaway greenhouse state. These three states coexist for intermediate values of solar constant. The existence of the large ice cap instability is also confirmed, although the critical latitude of the partially ice‐covered state extends equatorward. Also same as before, the small ice cap instability does not seem to appear.

The existence of liquid water within an oxidized environment on early Mars has been inferred by the Mn-rich rocks found during recent explorations on Mars. The oxidized atmosphere implied by the Mn-rich rocks would basically be comprised of CO2 and H2O without any reduced greenhouse gases such as H-2 and CH4. So far, however, it has been thought that early Mars could not have been warm enough to sustain water in liquid form without the presence of reduced greenhouse gases. Here, we propose that H2O2 could have been the gas responsible for warming the surface of the oxidized early Mars. Our one-dimensional atmospheric model shows that only 1 ppm of H2O2 is enough to warm the planetary surface because of its strong absorption at far-infrared wavelengths, in which the surface temperature could have reached over 273 K for a CO2 atmosphere with a pressure of 3 bar. A wet and oxidized atmosphere is expected to maintain sufficient quantities of H2O2 gas in its upper atmosphere due to its rapid photochemical production in slow condensation conditions. Our results demonstrate that a warm and wet environment could have been maintained on an oxidized early Mars, thereby suggesting that there may be connections between its ancient atmospheric redox state and possible aqueous environment.

Direct-imaging techniques of exoplanets have made significant progress recently and will eventually enable monitoring of photometric and spectroscopic signals of Earth-like habitable planets. The presence of clouds, however, would remain as one of the most uncertain components in deciphering such direct-imaged signals of planets. We attempt to examine how the planetary obliquity produces different cloud patterns by performing a series of general circulation model simulation runs using a set of parameters relevant for our Earth. Then we use the simulated photometric lightcurves to compute their frequency modulation that is due to the planetary spin–orbit coupling over an entire orbital period, and we attempt to see to what extent one can estimate the obliquity of an Earth twin. We find that it is possible to estimate the obliquity of an Earth twin within the uncertainty of several degrees with a dedicated 4 m space telescope at 10 pc away from the system if the stellar flux is completely blocked. While our conclusion is based on several idealized assumptions, a frequency modulation of a directly imaged Earth-like planet offers a unique methodology to determine its obliquity.