The model evaluation of 2009 to 2012 simulations shows an overall good performance in terms of annual mean predictions for most meteorological variables such as temperature at 2-m (T2) and wind speed at 10-m (WSP10), radiative variables such as shortwave and longwave radiation, and cloud variables such as cloud fraction (CF) and cloud droplet number concentration. For example, the annual mean biases (MBs) for T2 and WSP10 are -0.7 °C to -0.1 °C and 0 to 1.1 m s-1, respectively. The annual mean normalized mean biases (NMBs) for CF are 0.6-0.8%. The model also shows a good performance for maximum 8-h ozone (O3) with NMBs of -22.8% to -6.0% and for fine particulate matter (PM2.5) with NMBs of 4.5-22.9%. However, moderate-to-large wet biases are found for precipitation with NMBs of 16.7-64.6%, which is mainly caused by to the limitations of the Grell 3D cumulus parameterization used. The cloud optical thickness (COT) is also largely underpredicted with annual mean NMBs of -51.2% to -45.6% due to poor performance of some cloud variables and limitations in COT calculations. In addition, large biases are found for PM2.5 composition such as sulfate and organic carbon and for tropospheric ozone residual. These large biases may be due to inaccurate meteorological predictions (e.g., precipitation), emissions (e.g., primary PM species), and boundary conditions for some species (e.g., O3). These results will provide valuable information to further enhance the skills of WRF/Chem-ROMS for the coupled atmosphere-ocean modeling through improving the model inputs, configurations, and representations of identified physical and chemical processes.