Understanding the local transport of microplastics (MPs) emitted from the urban environment, such as those from vehicle tire wearing in streets and highways, is a necessary first step for quantifying their global transport cycle. By approximating microplastics as heavy particles, we conduct numerical simulations using large-eddy simulations (LESs) to understand how spatially organized sources and complex urban surface roughness affect their transport. Three sets of cases are considered, namely: (i) spatially uniform source and (ii) spatially organized source, and (iii) spatially organized source with explicitly resolved roughness elements, respectively. Results suggest that for a spatially organized source without buildings, source heterogeneity in streamwise direction only influences the vertical concentration profile up to z/Lz=0.11\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$z/L_z=0.11$$\end{document}. In contrast, that in spanwise direction influences the profile till z/Lz=0.58\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$z/L_z=0.58$$\end{document}, where z and Lz\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$L_z$$\end{document} are the vertical coordinate and the domain height, respectively. Simulations with buildings reveal that the buildings impede the transport of particles and particles accumulate on the leeward side of the buildings, which are characterized by wake turbulence and relatively quiescent flow. Within the canopy sub-layer, the gravitational settling effect of the particles becomes more significant because of the reduced wind. Because of that, the escape fraction of particles is smaller than in cases with no buildings and it decreases with increasing building height h and increasing building plan area fraction λp\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\lambda _p$$\end{document}. Finally, by finding suitable scalar displacement height ds\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$d_s$$\end{document} and scalar roughness length zos\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$z_{os}$$\end{document}, we find that similar to a passive scalar, an inertial sub-layer (ISL) still exists for heavy particles. This study highlights that for spatially organized particle sources, the momentum sinks due to urban roughness and the gravitational settling jointly affect the transport of heavy particles, which implies that the surface heterogeneity effect can be substantial in quantifying the atmospheric transport of microplastics of urban origins.
