Masaki Ishiwatari
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Masaki Ishiwatari has a total of 25 co-authors. Here are the top co-authors by number of shared publications:
- George L. Hashimoto
Department of Earth Sciences Okayama University Okayama Japan
4 shared publications - Yoshi-Yuki Hayashi
Kobe University
4 shared publications - Yoshiyuki O. Takahashi
Graduate School of Science Kobe University Kobe Japan
3 shared publications - Hiroki Kashimura
Kobe University
2 shared publications - Kensuke Nakajima
Faculty of Sciences, Kyushu University, Fukuoka, Japan
2 shared publications - Kiyoshi Kuramoto
Hokkaido University
2 shared publications
Latest Publications
- Open Access
- 12 May 2026
- Published by Wiley in Geoscience Data Journal
- Vol. 13 (3) , e70080
ALERA ‐V Version 1.0: An Objective Analysis Dataset of the Venus Atmosphere
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.
- 12 May 2026
- Published by Wiley in Geoscience Data Journal
- Vol. 13 (3) , e70080
Open AccessALERA ‐V Version 1.0: An Objective Analysis Dataset of the Venus Atmosphere
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.
- Open Access
- 23 September 2022
- Published by Copernicus GmbH
Thermal tides reproduced in the assimilation results of horizontal winds obtained from Akatsuki UVI observations
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. 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. 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). 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). 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. 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. [1] Sugimoto, N., et al. Development of an ensemble Kalman filter data assimilation system for the Venusian atmosphere. Scientific Reports 7(1), 9321 (2017). [2] Sugimoto, N., et al. Baroclinic instability in the Venus atmosphere simulated by GCM. J. Geophys. Res. Planets 119, 1950–1968 (2014). [3] Horinouchi, T., et al. Venus Climate Orbiter Akatsuki Cloud Motion Vector Data Set v1.0, JAXA Data Archives and Transmission System (2021). [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).
- 23 September 2022
- Published by Copernicus GmbH
Open AccessThermal tides reproduced in the assimilation results of horizontal winds obtained from Akatsuki UVI observations
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. 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. 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). 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). 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. 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. [1] Sugimoto, N., et al. Development of an ensemble Kalman filter data assimilation system for the Venusian atmosphere. Scientific Reports 7(1), 9321 (2017). [2] Sugimoto, N., et al. Baroclinic instability in the Venus atmosphere simulated by GCM. J. Geophys. Res. Planets 119, 1950–1968 (2014). [3] Horinouchi, T., et al. Venus Climate Orbiter Akatsuki Cloud Motion Vector Data Set v1.0, JAXA Data Archives and Transmission System (2021). [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).
- Open Access
- 7 June 2021
- Published by American Geophysical Union (AGU) in Journal of Geophysical Research: Atmospheres
- Vol. 126 (11)
Revision of “Dependence of Climate States of Gray Atmosphere on Solar Constant: From the Runaway Greenhouse to the Snowball States” by Ishiwatari et al. (2007)
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.
- 7 June 2021
- Published by American Geophysical Union (AGU) in Journal of Geophysical Research: Atmospheres
- Vol. 126 (11)
Open AccessRevision of “Dependence of Climate States of Gray Atmosphere on Solar Constant: From the Runaway Greenhouse to the Snowball States” by Ishiwatari et al. (2007)
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.
- 28 July 2020
- Published by American Astronomical Society in The Astrophysical Journal
- Vol. 898 (2) , 95
Obliquity of an Earth-like Planet from Frequency Modulation of Its Direct-imaged Lightcurve: Mock Analysis from General Circulation Model Simulation
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.
- 28 July 2020
- Published by American Astronomical Society in The Astrophysical Journal
- Vol. 898 (2) , 95
Obliquity of an Earth-like Planet from Frequency Modulation of Its Direct-imaged Lightcurve: Mock Analysis from General Circulation Model Simulation
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.
- 1 April 2020
- Published by American Astronomical Society in The Astrophysical Journal
- Vol. 893 (2) , 168
H2O2-induced Greenhouse Warming on Oxidized Early Mars
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 H2 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.
- 1 April 2020
- Published by American Astronomical Society in The Astrophysical Journal
- Vol. 893 (2) , 168
H2O2-induced Greenhouse Warming on Oxidized Early Mars
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 H2 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.