Capabilities are developed to assimilate GOES-R GLM FED data within the operational GSI framework, using ensemble Kalman filter (EnKF), and ensemble-3DVar (En3DVar) hybrid methods. FED observation operators based on graupel mass (FEDM) or graupel volume (FEDV) are used in the EnKF and while FEDM is used in the variational framework. The capabilities are tested with a mesoscale convective system (MCS) case and a case with multiple supercells. When tested with the MCS case with EnKF, FEDM and FEDV observation operators are found to perform similarly, and the assimilation of FED data shows clear positive impact. FED DA is effective in identifying and properly initializing regions of intense convection and helps improve the forecasts up to 3 hours. Direct adjustment to graupel as well as to other model states (such as vertical velocity) through flow-dependent cross-covariances work together to produce more consistent analyses then improved forecasts. 3DVar, EnKF and En3DVar are applied to the supercell case, and FED DA again leads to positive impacts on analyses and forecasts of the storms. EnKF and hybrid En3DVar are found to perform similarly while 3DVar performs quite a bit worse. Overall, the 2D FED data show great promise for initializing convective storms via advanced DA methods, especially in the absence of 3D volumetric radar data.