NOx_stored(t)=NOx_stored(t???1)+FG ???NO x_rate???NOx efficency(t)?????t wherein FG_NOx_rate is a feedgas NOx generation rate from the engine determined based on engine operating conditions, such as engine speed, load, air-fuel ratio, mass airflow, etc.