Development of New Observation Operators for Assimilating GOES-R Geostationary Lightning Mapper Flash Extent Density Data Using GSI EnKF: Tests with Two Convective Events over the United States

In a prior study, GOES-R Geostationary Lightning Mapper (GLM) flash extent density (FED) data were assimilated using ensemble Kalman filter into a convection-allowing model for a mesoscale convective system (MCS) and a supercell storm. The FED observation operator based on a linear relation with column graupel mass was tuned by multiplying a factor to avoid large FED forecast bias. In this study, new observation operators are developed by fitting a third-order polynomial to GLM FED observations and the corresponding FED forecasts of graupel mass of the MCS and/or supercell cases. The new operators are used to assimilate the FED data for both cases, in three sets of experiments called MCSFit, SupercellFit, and CombinedFit, and their performances are compared with the prior results using the linear operator and with a reference simulation assimilating no FED data. The new nonlinear operators reduce the frequency biases (root-mean-square innovations) in the 0–4-h forecasts of the FED (radar reflectivity) relative to the results using the linear operator for both storm cases. The operator obtained by fitting data from the same case performs slightly better than fitting to data from the other case, while the operator obtained by fitting forecasts of both cases produce intermediate but still very similar results, and the latter is considered more general. In practice, a more general operator can be developed by fitting data from more cases.