In this presentation, the simulations of the Weather Research and Forecast (WRF) and Community Multiscale Air Quality (CMAQ) Models applied to the New York City (NYC) area are assessed with the aid of vertical profiling and column integrated remote sensing measurements. First, we find that when turbulent mixing processes are dominant, the WRF derived planetary boundary layer (PBL) height exhibits a strong linear correlation (R>0.85) with lidar derived PBL height. In these comparisons, we estimate the PBL height from the lidar measurements using a Wavelet Covariance Transform (WCT) approach which is modified to better isolate the convective layer from the residual layer (RL). Furthermore, the WRF-Lidar PBL height comparisons are made using different PBL parameterization schemes including the Asymmetric Convective Model-version2 (ACM2) and the Modified Blackadar (BLK) scheme which are both runs using hindcast data, and the Mellor-Yamada-Janjic (MYJ) scheme run in forecast mode. Our findings show that the correlations for these runs were high (>0.8), but the hindcast runs exhibit smaller overall dispersion (≈0.1) than the forecast runs. In addition, we apply continuous 24hours/7days vertical ceilometer measurements to assess WRF-CMAQ model forecasts of surface PM2.5 (particulate matter has aerodynamic diameter <2.5µm). Particularly, strong overestimations in the surface PM2.5 mass that are observed in the summer prior to sunrise are shown to be strongly connected to underestimations of the PBL heights and less to enhanced emissions. This interpretation is consistent with observations that TEOM (Tapered Element Oscillating MicroBalance) PM2.5 measurements are better correlated to path integrated CMAQ PM2.5 than the near surface measurements during these periods.