You can get an approximate answer to this request using the Process Analysis features in CMAQ (CMAQ/CMAQ_UG_ch09_process_analysis.md at main · USEPA/CMAQ · GitHub). It is able to give you the production and loss of individual compounds, or user-defined families for processes that you mention. These data are processed for the model grid and are available in concentration units, not total mass gained or lost. In order to get domain-wide totals, you would need to convert the concentrations or mixing ratios to total mass, and sum across the domain.
In the next release of CMAQ, there will be an option for outputing a ‘BUDGET’ file which will include domain-wide mass gains and losses for all processes for any species or family you wish. These values will not include the spatial detail of process analysis, but will have the same hourly temporal resolution.