Implementation of Bessel's method for solar eclipses prediction in the WRF-ARW model

Solar eclipses are predictable astronomical events that abruptly reduce the incoming solar radiation into the Earth's atmosphere, which frequently results in non-negligible changes in meteorological fields. The meteorological impacts of these events have been analyzed in many studies since the late 1960s. The recent growth in the solar energy industry has greatly increased the interest in providing more detail in the modeling of solar radiation variations in numerical weather prediction (NWP) models for the use in solar resource assessment and forecasting applications. The significant impact of the recent partial and total solar eclipses that occurred in the USA (23 October 2014) and Europe (20 March 2015) on solar power generation have provided additional motivation and interest for including these astronomical events in the current solar parameterizations.<br><br>Although some studies added solar eclipse episodes within NWP codes in the 1990s and 2000s, they used eclipse parameterizations designed for a particular case study. In contrast to these earlier implementations, this paper documents a new package for the Weather Research and Forecasting–Advanced Research WRF (WRF-ARW) model that can simulate any partial, total or hybrid solar eclipse for the period 1950 to 2050 and is also extensible to a longer period. The algorithm analytically computes the trajectory of the Moon's shadow and the degree of obscuration of the solar disk at each grid point of the domain based on Bessel's method and the <i>Five Millennium Catalog of Solar Eclipses</i> provided by NASA, with a negligible computational time. Then, the incoming radiation is modified accordingly at each grid point of the domain.<br><br>This contribution is divided in three parts. First, the implementation of Bessel's method is validated for solar eclipses in the period 1950–2050, by comparing the shadow trajectory with values provided by NASA. Latitude and longitude are determined with a bias lower than 5 × 10⁻³ degrees (i.e., ∼ 550 m at the Equator) and are slightly overestimated and underestimated, respectively. The second part includes a validation of the simulated global horizontal irradiance (GHI) for four total solar eclipses with measurements from the Baseline Surface Radiation Network (BSRN). The results show an improvement in mean absolute error (MAE) from 77 to 90 % under cloudless skies. Lower agreement between modeled and measured GHI is observed under cloudy conditions because the effect of clouds is not included in the simulations for a better analysis of the eclipse outcomes. Finally, an introductory discussion of eclipse-induced perturbations in the surface meteorological fields (e.g., temperature, wind speed) is provided by comparing the WRF–eclipse outcomes with control simulations