Environment and climate change problems are very complicated, and their research and operational prediction heavily depend on powerful computer techniques. Even with today most powerful supercomputer, the climate and environmental models are still limited to very coarse resolution. In this paper, we report our recent effort in parallelizing our very-high-resolution numerical model systems. First, the mathematical equations, algorithms, and numerical schemes are designed and analyzed; then domain decomposition, data decomposition, and functional decomposition schemes are tested in our implementations on clusters of HP workstations and/or DEC Alpha stations with PVM; finally, the performance and load balance are analyzed. Shelterbelts cause significantly inhomogeneous computation distribution on the domain, therefore, common and easiest domain decomposition does not work well on our problem. Special care must be taken to treat computations around shelterbelts. With carefull design of algorithms, we found that cheap and still powerful workstations or PCs make it possible to run these models in clusters of workstations or PCs.