Nitrogen oxides (NOx = NO + NO2 ) play a crucial role in the formation of ozone and secondary inorganic and organic aerosols, thus affecting human health, global radiation budget, and climate. The diurnal and spatial variations in NO2 are functions of emissions, advection, deposition, vertical mixing, and chemistry. Their observations, therefore, provide useful constraints in our understanding of these factors. We employ a Regional chEmical and trAnsport model (REAM) to analyze the observed temporal (diurnal cycles) and spatial distributions of NO2 concentrations and tropospheric vertical column densities (TVCDs) using aircraft in situ measurements and surface EPA Air Quality System (AQS) observations as well as the measurements of TVCDs by satellite instruments (OMI: the Ozone Monitoring Instrument; GOME-2A: Global Ozone Monitoring Experiment – 2A), ground-based Pandora, and the Airborne Compact Atmospheric Mapper (ACAM) instrument in July 2011 during the DISCOVER-AQ campaign over the Baltimore– Washington region. The model simulations at 36 and 4 km resolutions are in reasonably good agreement with the regional mean temporospatial NO2 observations in the daytime. However, we find significant overestimations (underestimations) of model-simulated NO2 (O3 ) surface concentrations during nighttime, which can be mitigated by enhancing nocturnal vertical mixing in the model. Another discrepancy is that Pandora-measured NO2 TVCDs show much less variation in the late afternoon than simulated in the model. The higher-resolution 4 km simulations tend to show larger bi-