“…The values of factors (6)-(8) were obtained using the GIS software ''ArcGIS'' (ESRI Japan, Tokyo). With regards to factor (6), the spatial extent at which the variables should be extracted was not clear, so we generated buffers of given radii (200, 400, 600, 800, 1,000, 1,200, 1,400 and 1,600 m) around each site and then calculated the proportion of paddy field area within the buffered area (e.g., Miyashita et al 2008). With regards to factor (7), for the same reason, the slope variables were obtained from the average slopes in 200, 400,…,1,600 m rasterized grid squares using the GIS software.…”