The paper proposes a multidimensional flux-form semi-Lagrangian (FFSL) transport scheme to extend one-dimensional (1D) forward-in-time, upstream-biased, flux-form transport schemes to higher dimensions. The FFSL scheme is designed to be conservative, upstream-biased, and monotonic, while preserving tracer correlations. It is implemented efficiently on spherical geometry and validated through idealized tests and realistic three-dimensional global transport simulations using wind data from data assimilation systems. The scheme's stability is analyzed using both von Neuman approaches and empirical methods on a 2D Cartesian plane. The FFSL scheme is shown to be accurate, stable, and efficient, with the ability to handle large time steps without stringent stability limitations, making it suitable for applications in atmospheric science and other geophysical fields.The paper proposes a multidimensional flux-form semi-Lagrangian (FFSL) transport scheme to extend one-dimensional (1D) forward-in-time, upstream-biased, flux-form transport schemes to higher dimensions. The FFSL scheme is designed to be conservative, upstream-biased, and monotonic, while preserving tracer correlations. It is implemented efficiently on spherical geometry and validated through idealized tests and realistic three-dimensional global transport simulations using wind data from data assimilation systems. The scheme's stability is analyzed using both von Neuman approaches and empirical methods on a 2D Cartesian plane. The FFSL scheme is shown to be accurate, stable, and efficient, with the ability to handle large time steps without stringent stability limitations, making it suitable for applications in atmospheric science and other geophysical fields.