In suburban regions, vacant lots potentially offer significant opportunities for biodiversity conservation. Recently, in Japan, due to an economic recession, some previously developed lands have become vacant. Little is known, however, about the legacy of earlier earthmoving, which involves topsoil removal and ground levelling before residential construction, on plant community composition in such vacant lots. To understand (dis)assembly processes in vacant lots, we studied 24 grasslands in a suburban region in Japan: 12 grasslands that had experienced earthmoving and 12 that had not. We surveyed plant community composition and species richness, and clarified compositional turnover (replacement of species) and nestedness (non-random species loss) by distance-based β-diversities, which were summarized by PCoA analysis. We used piecewise structural equation modeling to examine the effects of soil properties, mowing frequency, past and present habitat connectivities on compositional changes. As a result, past earthmoving, mowing frequency, soil properties, and past habitat connectivity were found to be the drivers of compositional turnover. In particular, we found legacy effects of earthmoving: earthmoving promoted turnover from native grassland species to weeds in arable lands or roadside by altering soil properties. Mowing frequency also promoted the same turnover, implying that extensive rather than intensive mowing can modify the negative legacy effects and maintain grassland species. Decrease in present habitat connectivity marginally enhanced non-random loss of native grassland species (nestedness). Present habitat connectivity had a positive effect on species richness, highlighting the important roles of contemporary dispersal. Our study demonstrates that community assembly is a result of multiple processes differing in spatial and temporal scales. We suggest that extensive mowing at local scale, as well as giving a high conservation priority to grasslands with high habitat connectivity at regional scale, are the promising actions to maintain endangered native grassland species in suburban landscapes with negative legacy effects of earthmoving.

Species occurrence data (vegetation_data.csv) was first analyzed by PCoA, and the resultant ordination axes were used as response variables in piecewise SEM analysis. Principal components of soil variables (soil_properties.csv), past and present habitat connectivities (site_area_and_connectivity.csv), and land-use history (landuse_and_mowing.csv) were all used as explanatory variables in the SEM analysis.