The vertical distribution of biomass burning aerosol (BBA) is important in regulating their impacts on weather and climate. The plume‐rise process affects the injection height of BBA and interacts with the air parcel lifting and cloud processes. However, these processes are not represented in most global climate models. In this study, we replaced the fixed vertical profiles of monthly BBA emissions in the Department of Energy's Energy Exascale Earth System Model version 1 (E3SMv1) with an interactive fire plume‐rise model. The vertical distribution of BBA emissions was calculated as a function of ambient thermodynamic conditions from the host E3SMv1, with distributions of fire sizes and sensible heat fluxes derived from the observations. The maximum fire radiative power (FRP) technique was used to determine the fire size. Scaling‐FRP technique is used to calculate the wildfire heat release. Daily BBA emission, superimposed with a fire diurnal cycle retrieved from the satellite observation, was included in model simulations. The model shows improved agreement with satellite retrievals and in situ measurement during the National Oceanic and Atmospheric Administration Wildfire Experiment for Cloud chemistry, Aerosol absorption, and Nitrogen campaign. The model‐observation comparison demonstrates the importance of the plume‐rise model and fire diurnal cycle assumption in determining the BBA fields. We also find that E3SMv1 with new features produces a larger carbonaceous aerosol burden, leading to 0.13 W m−2 warming at the top of atmosphere compared to the default E3SMv1. This highlights the importance of accurately representing the BBA injection height and indicates a no‐linear nature in the BBA‐induced radiative effect.