This paper describes the first official release (v1.0) of RTTOV-gb. RTTOV-gb is a FORTRAN 90 code developed by adapting the atmospheric radiative transfer code RTTOV, focused on satellite-observing geometry, to the ground-based observing geometry. RTTOV-gb is designed to simulate ground-based upward-looking microwave radiometer (MWR) observations of atmospheric downwelling natural radiation in the frequency range from 22 to 150 GHz. Given an atmospheric profile of temperature, water vapor, and, optionally, cloud liquid water content, and together with a viewing geometry, RTTOV-gb computes downwelling radiances and brightness temperatures leaving the bottom of the atmosphere in each of the channels of the sensor being simulated. In addition, it provides the sensitivity of observations to the atmospheric thermodynamical state, i.e., the Jacobians. Therefore, RTTOV-gb represents the forward model needed to assimilate ground-based MWR data into numerical weather prediction models, which is currently pursued internationally by several weather services. RTTOV-gb is fully described in a previous paper (De Angelis et al., 2016), while several updates are described here. In particular, two new MWR types and a new parameterization for the atmospheric absorption model have been introduced since the first paper. In addition, estimates of the uncertainty associated with the absorption model and with the fast parameterization are given here. Brightness temperatures (T_B) computed with RTTOV-gb v1.0 from radiosonde profiles have been compared with ground-based MWR observations in six channels (23.8, 31.4, 72.5, 82.5, 90.0, and 150.0 GHz). The comparison shows statistics within the expected accuracy. RTTOV-gb is now available to licensed users free of charge from the Numerical Weather Prediction Satellite Application Facility (NWP SAF) website, after registration. Coefficients for four MWR instrument types and two absorption model parameterizations are also freely available from the RTTOV-gb support website.

Total precipitable water (TPW), a column of water vapor content in the atmosphere, provides information on the spatial distribution of moisture. The high-resolution TPW, together with atmospheric stability indices such as convective available potential energy (CAPE), is an effective indicator of severe weather phenomena in the pre-convective atmospheric condition. With the advent of high performing imaging instrument onboard geostationary satellites such as Advanced Himawari Imager (AHI) onboard Himawari-8 of Japan and Advanced Meteorological Imager (AMI) onboard GeoKompsat-2A of Korea, it is expected that unprecedented spatiotemporal resolution data (e.g., AMI plans to provide 2 km resolution data at every 2 min over the northeast part of East Asia) will be provided. To derive TPW from such high-resolution data in a timely fashion, an efficient algorithm is highly required. Here, machine learning approaches—random forest (RF), extreme gradient boosting (XGB), and deep neural network (DNN)—are assessed for the TPW retrieved from AHI over the clear sky in Northeast Asia area. For the training dataset, the nine infrared brightness temperatures (BT) of AHI (BT8 to 16 centered at 6.2, 6.9, 7.3, 8.6, 9.6, 10.4, 11.2, 12.4, and 13.3 $\mathsf{\mu }\mathrm{m}$ , respectively), six dual channel differences and observation conditions such as time, latitude, longitude, and satellite zenith angle for two years (September 2016 to August 2018) are used. The corresponding TPW is prepared by integrating the water vapor profiles from InterimEuropean Centre for Medium-Range Weather Forecasts Re-Analysis data (ERA-Interim). The algorithm performances are assessed using the ERA-Interim and radiosonde observations (RAOB) as the reference data. The results show that the DNN model performs better than RF and XGB with a correlation coefficient of 0.96, a mean bias of 0.90 mm, and a root mean square error (RMSE) of 4.65 mm when compared to the ERA-Interim. Similarly, DNN results in a correlation coefficient of 0.95, a mean bias of 1.25 mm, and an RMSE of 5.03 mm when compared to RAOB. Contributing variables to retrieve the TPW in each model and the spatial and temporal analysis of the retrieved TPW are carefully examined and discussed.

As convective clouds in Northeast Asia are accompanied by various hazards related with heavy rainfall and thunderstorms, it is very important to detect convective initiation (CI) in the region in order to mitigate damage by such hazards. In this study, a novel approach for CI detection using images from Meteorological Imager (MI), a payload of the Communication, Ocean, and Meteorological Satellite (COMS), was developed by improving the criteria of the interest fields of Rapidly Developing Cumulus Areas (RDCA) derivation algorithm, an official CI detection algorithm for Multi-functional Transport SATellite-2 (MTSAT-2), based on three machine learning approaches—decision trees (DT), random forest (RF), and support vector machines (SVM). CI was defined as clouds within a 16 × 16 km window with the first detection of lightning occurrence at the center. A total of nine interest fields derived from visible, water vapor, and two thermal infrared images of MI obtained 15–75 min before the lightning occurrence were used as input variables for CI detection. RF produced slightly higher performance (probability of detection (POD) of 75.5% and false alarm rate (FAR) of 46.2%) than DT (POD of 70.7% and FAR of 46.6%) for detection of CI caused by migrating frontal cyclones and unstable atmosphere. SVM resulted in relatively poor performance with very high FAR ~83.3%. The averaged lead times of CI detection based on the DT and RF models were 36.8 and 37.7 min, respectively. This implies that CI over Northeast Asia can be forecasted ~30–45 min in advance using COMS MI data.

Hyperspectral imagery (HSI) provides substantial information on optical features of water bodies that is usually applicable to water quality monitoring. However, it generates considerable uncertainties in assessments of spatial and temporal variation in water quality. Thus, this study explored the influence of different optical methods on the spatial distribution and concentration of phycocyanin (PC), chlorophyll-a (Chl-a), and total suspended solids (TSSs) and evaluated the dependence of algal distribution on flow velocity. Four ground-based and airborne monitoring campaigns were conducted to measure water surface reflectance. The actual concentrations of PC, Chl-a, and TSSs were also determined, while four bio-optical algorithms were calibrated to estimate the PC and Chl-a concentrations. Artificial neural network atmospheric correction achieved Nash-Sutcliffe Efficiency (NSE) values of 0.80 and 0.76 for the training and validation steps, respectively. Moderate resolution atmospheric transmission 6 (MODTRAN 6) showed an NSE value >0.8; whereas, atmospheric and topographic correction 4 (ATCOR 4) yielded a negative NSE value. The MODTRAN 6 correction led to the highest R² values and lowest root mean square error values for all algorithms in terms of PC and Chl-a. The PC:Chl-a distribution generated using HSI proved to be negatively dependent on flow velocity (p-value = 0.003) and successfully indicated cyanobacteria risk regions in the study area.

The Geostationary Environment Monitoring Spectrometer (GEMS) onboard the Geostationary Korean Multi-Purpose Satellite 2B was successfully launched in February 2020. GEMS is a hyperspectral spectrometer measuring solar irradiance and Earth radiance in the wavelength range of 300 to 500 nm. This paper introduces the spectral calibration algorithm for GEMS, which uses a nonlinear least-squares approach. Sensitivity tests for a series of unknown algorithm parameters such as spectral range for fitting, spectral response function (SRF), and reference spectrum were conducted using the synthetic GEMS spectrum prepared with the ground-measured GEMS SRF. The test results show that the required accuracy of 0.002 nm is achievable provided the SRF and the high-resolution reference spectrum are properly prepared. Such a satisfactory performance is possible mainly due to the inclusion of additional fitting parameters of spectral scales (shift, squeeze, and high order shifts) and SRF (width, shape and asymmetry). For the application to the actual GEMS data, in-orbit SRF is to be monitored using an analytic SRF function and the measured GEMS solar irradiance, while a reference spectrum is going to be selected during the instrument in-orbit test. The calibrated GEMS data is expected to be released by the end of 2020.

The Geostationary Environment Monitoring Spectrometer (GEMS) on board GEO-KOMPSAT-2B was launched in February 2020 and has been monitoring atmospheric chemical compositions over Asia. We present the first evaluation of the operational GEMS formaldehyde (HCHO) vertical column densities (VCDs) during and after the in-orbit test (IOT) period (August–October 2020) by comparing them with the products from the TROPOspheric Monitoring Instrument (TROPOMI) and Fourier-transform infrared (FTIR) and multi-axis differential optical absorption spectroscopy (MAX-DOAS) instruments. During the IOT, the GEMS HCHO VCDs reproduced the observed spatial pattern of TROPOMI VCDs over the entire domain (r= 0.62) with high biases (10 %–16 %). We found that the agreement between GEMS and TROPOMI was substantially higher in Northeast Asia (r= 0.90), encompassing the Korean Peninsula and east China. GEMS HCHO VCDs captured the seasonal variation in HCHO, primarily driven by biogenic emissions and photochemical activities, but showed larger variations than those of TROPOMI over coastal regions (Kuala Lumpur, Singapore, Shanghai, and Busan). In addition, GEMS HCHO VCDs showed consistent hourly variations with MAX-DOAS (r= 0.77) and FTIR (r= 0.86) but were 30–40 % lower than ground-based observations. Different vertical sensitivities of GEMS and ground-based instruments caused these biases. Utilizing the averaging kernel smoothing method reduces the low biases by approximately 10 % to 15 % (normalized mean bias (NMB): −47.4 % to −31.5 % and −38.6 % to −26.7 % for MAX-DOAS and FTIR, respectively). The remaining discrepancies are due to multiple factors, including spatial collocation and different instrumental sensitivities, requiring further investigation using inter-comparable datasets.