Sub-grid variability (SGV) in atmospheric trace gases within satellite pixels is a key issue in satellite design and interpretation and validation of retrieval products. However, characterizing this variability is challenging due to the lack of independent high-resolution measurements. Here we use tropospheric NO.sub.2 vertical column (VC) measurements from the Geostationary Trace gas and Aerosol Sensor Optimization (GeoTASO) airborne instrument with a spatial resolution of about 250 mx250 m to quantify the normalized SGV (i.e., the standard deviation of the sub-grid GeoTASO values within the sampled satellite pixel divided by the mean of the sub-grid GeoTASO values within the same satellite pixel) for different hypothetical satellite pixel sizes over urban regions. We use the GeoTASO measurements over the Seoul Metropolitan Area (SMA) and Busan region of South Korea during the 2016 KORUS-AQ field campaign and over the Los Angeles Basin, USA, during the 2017 Student Airborne Research Program (SARP) field campaign. We find that the normalized SGV of NO.sub.2 VC increases with increasing satellite pixel sizes (from â¼10 % for 0.5 kmx0.5 km pixel size to â¼35 % for 25 kmx25 km pixel size), and this relationship holds for the three study regions, which are also within the domains of upcoming geostationary satellite air quality missions. We also quantify the temporal variability in the retrieved NO.sub.2 VC within the same hypothetical satellite pixels (represented by the difference of retrieved values at two or more different times in a day). For a given satellite pixel size, the temporal variability within the same satellite pixels increases with the sampling time difference over the SMA. For a given small (e.g., less than or equal to4 h) sampling time difference within the same satellite pixels, the temporal variability in the retrieved NO.sub.2 VC increases with the increasing spatial resolution over the SMA, Busan region, and the Los Angeles Basin.