Introduction {#Sec1}
============

Biological communities in agricultural crop fields that have received low or no chemical application are more diverse than those of conventionally farmed fields (e.g. Bengtsson et al. [@CR2]). In particular, rice paddies maintained by biodiversity-friendly farming techniques are expected to function as important substitute habitats for organisms in natural wetlands, which have declined due to the development of urban and agricultural fields (Fasola and Ruiz [@CR5]; Washitani [@CR31]).

The farming practices of biodiversity-friendly paddy fields include the winter flooding of paddies to enhance the habitat conditions for wintering waterfowl (Elphic and Oring [@CR4]; Tourenq et al. [@CR30]; Kurechi [@CR12]). In the paddies at the Ramsar wetland site that surrounds the natural Kabukuri-numa Wetlands, which is one of the most important wintering sites for geese in Japan, the winter-flooding regime combined with no agricultural chemical applications has become a growing trend. This farming method is attracting attention as a new biodiversity-friendly approach throughout Japan and East Asia (Kurechi [@CR12]; Washitani [@CR31]). In fact, to restore paddy biodiversity including wintering waterfowl, no agricultural chemicals have been applied to the winter-flooded (WF) paddies in this region since around 1998. Moreover, WF paddies are also expected to enhance ecosystem services for rice production, including the biological control of insect pests by generalist predators.

In this study, the diversity of spiders----ubiquitous predators in rice paddies----was compared between WF and conventional paddies through field measurements conducted in 31 WF paddies and 7 conventional ones. In addition, factors limiting the abundance of spider populations in WF paddies were also examined. Crop yield in biodiversity-friendly fields can be low compared with that in paddies that receive conventional chemical applications (de Ponti et al. [@CR3]). Thus, effective pest management techniques other than the use of insecticides are needed to sustain growth in biodiversity-friendly farming (Zehnder et al. [@CR36]). Identifying the limiting factors of generalist predators in the WF paddies can help with the construction of effective management plans for insect pest control, because our previous studies showed that the spiders predate insect pests (Kobayashi et al. [@CR11]; Takada et al. [@CR29]), and suppress the density and the consequent damage in organic paddy fields (Takada et al. [@CR28]). Spider population growth is known to be constrained by habitat complexity (Langellotto and Denno [@CR13]; Takada and Miyashita [@CR26], [@CR27]) and/or prey availability (Wise [@CR33]; Miyashita [@CR14]). In this study, we focused on the vegetation structure of rice and abundant weeds, as well as the biomass of abundant arthropod groups, such as chironomidae and other nematocera, brachycera, and hoppers including Cicadellidae and Delphacidae, and the pecky rice bug *Stenotus rubrovittatus*, as potential limiting factors of the spiders.

Methods {#Sec2}
=======

Study site {#Sec3}
----------

A field survey was conducted in the area surrounding Kabukuri-numa at Osaki City, Miyagi Prefecture, Northern Japan (38°37′N, 141°07′E), where many WF paddies (total area: about 0.2 km^2^) have been cultivated (Kurechi [@CR12]) since around 1998. The mean annual precipitation from 2002 to 2007 was approximately 1166 mm, and the mean temperature in August during the same period was 23.6°C (Japan Meteorological Agency [@CR9]). In this region, insecticides for *Stenotus rubrovittatus* (Hemiptera: Miridae), the most abundant and influential rice pest (Yoshioka et al. [@CR34]; Takada et al. [@CR28]; [@CR35]), are applied to most of the conventional paddy fields after heading of rice plants in early August, when spider diversity in paddies is at the highest level (Oyama et al. [@CR19]). Field measurements were conducted twice, just before and after application of the insecticides to conventional fields, with the expectation that spider diversity decreases in conventional paddies, while remaining unchanged in WF ones after the application of insecticide to conventional fields.

In this area, clusters of WF paddies were distributed among conventionally managed paddies. The WF fields have been maintained by using similar management practices that include the application of natural fertilizers, such as fertilizer made from fish parts, and exclude chemical fertilizers, herbicides, and insecticides. Before seedling transplanting in the middle of May, the natural fertilizers are applied to WF fields. Some of the fields that did not receive herbicide application were dominated by weeds, such as *Monochoria vaginalis* var. *plantaginea*. So farmers do the mechanical weeding once or twice in June*.* WF fields are flooded all year. Even in winter (not rice-growing season) WF fields are flooded for providing suitable habitats for waterfowl (Kurechi [@CR12]). A typical farming schedule in the conventional paddies is as follows. Fungicides (Tetrachloroisophthalonitrile or hydroxyisoxazole and Metalaxyl) are applied to young rice plants in seedling boxes, and herbicides (indanofan, bensulfuronmethyl, and clomeprop (5 kg ha^-1^) or Pyriminobac-methyl, bromobutide, bensulfuronmethyl, and pentoxazone (10 kg ha^-1^)) are applied directly to the field after seedling transplanting in early May. Insecticides (etophenprox (2-3 l ha^-1^)) for rice leaf beetle *Oulema oryzae* Kuwayama are applied to the fields and then water is drained from the field for about ten days to stabilize rice yield and increase soil hardness in early June. After heading of rice plants in early August, insecticides for pecky rice bugs (dinotefuran (30 kg ha^-1^)) are applied once or twice to the field. The insecticides are applied to many of the conventional rice fields by radio-controlled helicopters. Most conventional paddies are seasonally dried-up after harvesting.

Field survey {#Sec4}
------------

A total of 31 WF and 7 conventional paddies were selected in an area of approximately 20 km^2^. The mean (SD) area of the study fields was 2100.8 ± 438.2 m^2^ and 1780.3 ± 523.7 m^2^ for the WF and conventional fields, respectively. Rice seedlings were transplanted in May 2007 at a mean (SD) density of 21.20 ± 2.0 hills/m^2^ and 23.24 ± 2.0 hills/m^2^ in the WF and conventional fields, respectively. In August 2007, net samplings (20 sweeps with a 36-cm-diameter insect net) for spiders and their potential prey were performed at the center of each field just before and after using insecticides for pecky rice bugs in conventional fields (insecticide application: 6, 18--19 August, first survey: 1--4 August; Second survey: 15--19 August). The sweeping sampling is an excellent method to collect invertebrates inhabiting the uppermost vegetation layer in paddies easily. Invertebrate body mass was calculated as mass (mg) = 0.0305 × (body length (mm))^2.62^ (Rogers et al. [@CR22]) to measure the abundance of potential prey for spiders. As indices of habitat complexity, we measured rice height of six randomly chosen hills and counted the number of individuals of the most abundant weed (*M. vaginalis*) within a 1 × 3 m^2^ quadrant at the center of each field (e.g., Rypstra et al. [@CR23], Miyashita and Takada [@CR15]).

Statistical analysis {#Sec5}
--------------------

To evaluate the effects of the farming practices (WF vs. conventional) and the insecticides applied to conventional fields (before vs. after) on the numbers of individual spiders and spider species in paddies of both types, we applied generalized linear mixed effect models with Poisson error distribution because the dependent variables are count data. Fixed factors in the models were the main effects of the farming practice and the insecticide applications and their interaction, and the random factor was the effect of the paddy. The significance of each fixed effect was checked by using analysis of deviance and by checking for the reduction in residual deviance by reducing each effect by means of the Chi-squared test.

The factors limiting the density and diversity of spiders in the WF paddies were analyzed by using generalized linear mixed effect models with Poisson error distribution because the dependent variables are count data. We examined the limiting factors of spider density and species richness using only the data of the first survey to avoid redundant analyses, because spider density in the WF paddies tended to be higher in the first field survey than in the second (see Results) and the results obtained from preliminary analyses using the data of the second survey were similar to those of the first survey. This model included the number of individual spiders or spider species richness as the dependent variable, and two indices of habitat complexity (mean rice height and *M. vaginalis* density) and the biomass of five prey groups as fixed factors (Table [1](#Tab1){ref-type="table"}), with the paddy field as a random effect. All tolerance values between the fixed factors were greater than the critical value of 0.1 (Table [1](#Tab1){ref-type="table"}), indicating no significant collinearity between them (Quinn and Keough [@CR20]). Model selection for the generalized linear mixed effect models was performed using Akaike's Information Criterion (AIC). AIC was calculated for each possible combination of fixed factors, and the models with small AIC were chosen as optimal models. All statistical analysis was carried out using *R* for Windows 2.13.1 (*R* Development Core Team [@CR21]).Table 1**Independent variables used in the generalized linear mixed effect models for identifying limiting factors of spider density and species richness in winter-flooded (WF) paddies**Independent variablesMeanRangeTolerance valuesHabitat complexityMean rice height (cm)67.3459.67-77.830.903*Monochoria vaginalis* (no.individuals/3 m^2^)67.420-3930.794Prey availability (mg/20 sweeps with an insect net)Chironomid130.6814.23-975.860.837Other nematocera32.076.58-100.540.913*S. rubrovittatus*7.080-65.370.746Brachycera5.230.38-29.900.831Hopper3.980.79-10.170.823

Results {#Sec6}
=======

A total of 234 spiders in 15 species or families were observed. Dominant spider species were *Tetragnatha caudicula*, *T. praedonia,* and *T. extensa* for web-building spiders and *Pachygnatha quadrimaculata*, *P. clercki*, and Lycosidae spp. for cursorial spiders (Table [2](#Tab2){ref-type="table"}). Juveniles were identified only to family.Table 2**Mean density of each spider species or family per field for WF and conventional paddies**Species/familyBefore insecticide applicationAfter insecticide applicationWFConventionalWFConventional*Tetragnatha caudicula***2.7740.1430.4060***T. praedonia***0.7740.1430.4690***T. extensa***0.1940.1430.04060***T. maxillosa***0.48400.0940***T. squamata***0.09700.0310***Pachygnatha quadrimaculata***0.51600.4060***P. clercki***0.16100.3440Lycosidae0.06500.2500***Neoscona adianta***0.065000***Argiope bruennichii***0.0320.28600***Larinioides cornutus***000.0630Linyphiidae0.03200.2810Thomiside0.0320.1430.0310Theridiidae000.0310Unknown spiders01.7810.143**

The main effect of farming practice (WF vs. conventional) on spider density (χ^2^~1~ = 23.369, *P* \< 0.001) and species richness (χ^2^~1~ = 17.817, *P* \< 0.001) was highly significant, while the interaction effect of farming practice and insecticide application was significant on neither spider density nor species richness (density: χ^2^~1~ = 2.792, *P* = 0.095; species richness: χ^2^~1~ = 3.204, *P* = 0.073). This result indicates that spider density and species richness were higher in WF paddies than in conventional ones both before and after insecticide applications to conventional paddies (Figure [1](#Fig1){ref-type="fig"}). Since the main effect of insecticide application (before vs. after) was significant only for spider density (density: χ^2^~1~ = 5.920, *P* = 0.015; species richness: χ^2^~1~ = 1.417, *P* = 0.234), insecticide application to conventional paddies was suggested to reduce spider density in both the WF and conventional paddies (Figure [1](#Fig1){ref-type="fig"}). Moreover, the variations in spider density and species richness in the WF paddies seemed to increase after applying insecticide: the interquartile range----the difference between upper (75% quartile) and lower (25% quartile) hinges----of spider density was 5.5 before insecticide application, but 6 after (Figure [1](#Fig1){ref-type="fig"}a), and the interquartile range of spider species richness was 2 before insecticide application, but 4 after insecticide application (Figure [1](#Fig1){ref-type="fig"}b).Figure 1**Spider diversity in WF and conventional paddy fields.** Boxplots for the spider density **(a)** and species richness **(b)** in the WF and conventional paddy fields before (filled bars) and after (open bars) applying insecticide to conventional paddies in early August.

The best model for explaining the spider density and species richness in WF paddies included the biomass of chironomids and other nematocera (Table [3](#Tab3){ref-type="table"}) but did not include any factors of habitat complexity, such as mean rice height or *M. vaginalis* density. These two groups of prey increased both spider density and species richness (Figure [2](#Fig2){ref-type="fig"}).Table 3**Estimated coefficients(SE) and information theoretic statistics of the top five and null (intercept only) models explaining spider density and species richness in WF paddies**RankInterceptRice height***M. vaginalis***ChironomidOther nematocera***S. rubrovittatus***BrachyceraHopperAICΔ AIC***w*** ~***i***~Spider density1-0.941 (0.457)\--0.473 (0.092)0.014 (0.004)\-\--52.1000.1192-0.717 (0.471)\--0.460 (0.090)0.013 (0.004)\--0.028 (0.022)-52.450.340.10030.844 (1.736)-0.027 (0.26)-0.487 (0.93)0.013 (0.004)\-\--52.970.870.0774-0.755 (0.543)\--0.031 (0.53)0.445 (0.0102)0.014 (0.004)\-\--53.751.650.0525-0.473 (0.091)0.013 (0.004)\--0.025 (0.023)-53.801.700.051Null1.490 (0.164)\-\-\-\-\--73.4421.330Spider species richness1-0.535 (0.481)\--0.278 (0.098)0.008 (0.005)\-\--21.1300.0742-0.286 (0.517)\--0.262 (0.097)0.009 (0.005)\-\--0.053 (0.044)21.630.500.0573-0325 (0.451)\--0.294 (0.096)\-\-\--21.870.740.0514-0372 (0.515)\--0.268 (0.099)0.007 (0.005)\--0.020 (0.25)-22.401.270.0395-0.138 (0.481)\--0.278 (0.097)\-\--0.025 (0.025)-22.731.600.033Null0.960 (0.111)\-\-\-\-\-\--28.687.550.002Ranking of the sub models is based on Akaike's information criterion.Figure 2**Spider diversity and their potential prey in WF paddy fields.** Relationships between total biomass of important prey (chironomids and other nematocera) and spider density **(a)** and species richness **(b)** in WF paddies before applying insecticide to conventional fields. Poisson regression lines are also shown.

Discussion {#Sec7}
==========

We revealed that spider density and species richness were higher in the WF paddies than in the conventional ones both before and after insecticide application to conventional paddies for pecky rice bugs in early August. Moreover, our results suggest that the insecticide application decreased spider density in not only conventional paddies but also WF ones. The spider diversity in the WF paddies was significantly affected by the availability of two prey groups, chironomids and other nematocera.

The significant main effect of farming practice (WF vs. conventional) on spider density and species richness appears to indicate that these values were already lower in conventional fields than in WF ones before the insecticide application to conventional paddies (Table [2](#Tab2){ref-type="table"}, Figure [1](#Fig1){ref-type="fig"}). At least four possible explanations for the farming practice effect can be identified. First, fungicides applied to young rice plants in seedling boxes, herbicides and insecticides early in the rice-growing season may be responsible for the reduced spider diversity in conventional fields during the early-growing season. Such agricultural chemicals are known to affect invertebrates (Oyama and Kidokoro [@CR18]; Wilson et al. [@CR32]; Amano et al. [@CR1]) and were probably responsible for decreasing the number of spider prey in the paddies. The second mechanism for the farming practice effect on the spiders is due to the drainage of conventional paddies in early June. Drainage dries and hardens the soil in conventional paddies (see Study site) and likely decreases the number of aquatic invertebrate, including juveniles of dipterous insects, which are important prey for spiders (Ishijima et al. [@CR8]; Tahir and Butt [@CR25]). Third, organic fertilizer applied to the WF fields before seedling transplanting might increase spider diversity due to increasing the availability of important prey for spiders. Previous studies have reported that detrital subsidy, such as organic fertilizer, enhances spider populations indirectly through increasing detritivores, such as chironomids and collembolas (Murata [@CR16]; Settle et al. [@CR24]; Halaj and Wise [@CR6]). Finally, tubificid worms as well as waterfowl fertilize WF paddy soil with their feces (Kurechi [@CR12]), which may increase spider diversity due to increasing the availability of important prey for spiders. Because WF fields are flooded all year, tubificid worms are known to active even in winter, fertilizing the soil in WF fields. On the other hand, conventional fields are dried-up after harvesting. Clarification of these causalities requires manipulative experiments controlling each component with many replications.

The main effect of insecticide application (before vs. after) on spider density was also significant, suggesting the insecticides applied to conventional fields may decrease spider density in not only conventional but also WF paddies, despite not applying insecticides to WF paddies directly. One possible explanation is as follows. WF and conventional paddy fields are spatially intermingled in the study area, so insecticides for pecky rice bugs sprayed onto the conventional rice fields by radio-controlled helicopters may easily diffuse into adjacent WF paddies and consequently affect spider density. Increased variation in the spider density and species richness in WF paddies after insecticide application (Figure [1](#Fig1){ref-type="fig"}) appears to support this possibility, because the intensity of the negative effect of the insecticides on spiders inhabiting the WF paddies may vary depending on the distance between the WF and conventional paddies. If so, increasing the spatial gathering of WF paddies will lead to an increase in the number of spiders in each WF paddy due to avoiding the diffusion of insecticides sprayed by radio-controlled helicopters. Enhancement of spider populations by strategic spatial placement of WF paddies may also suppress damage by pecky rice bugs in WF paddies. Spider species with high abundance in the WF paddies were also shown to predate pecky rice bugs through previous DNA-based gut-content analysis (Kobayashi et al. [@CR11]; Takada et al. [@CR29]). Moreover, *Tetragnatha* spp. spiders were revealed to suppress bug density and the consequent damage in organic paddies fields (Takada et al. [@CR28]).

Spider diversity in WF paddies increased with the increasing abundance of chironomids and other nematocera. This was revealed by statistical analysis of the limiting factors of spider diversity (Table [3](#Tab3){ref-type="table"}). Because we used a comparative approach, the casual relationship between spiders and these insects was no proven in this study. However, previous studies have suggested that these dipterous insects are frequently consumed by the spiders in paddies (Ishijima et al. [@CR8]; Tahir and Butt [@CR25]). Taken together, the dipterous insects most likely limit spider diversity in WF paddy fields. Our study is one of the few that provides evidence of the correlation between the abundance of important prey and increased spider diversity in paddies, with many replications.

Spiders and their potential prey were collected only by one method (i.e. sweeping sampling) in this study. Sweeping sampling is an excellent method to collect invertebrates inhabiting the uppermost vegetation layer in paddies easily. But it is necessary to use other methods (e.g., direct counting), to estimate spiders inhabiting the bottom layer, such as Lycosidae.

Many studies have shown that flooding rice paddies can provide suitable habitats for waterfowl (Elphick and Oring [@CR4]; Tourenq et al. [@CR30]; Kurechi [@CR12]). The present study suggests that in WF paddies the diversity of generalist predators is higher than in conventional paddies during the rice-growing season. Moreover, a combination of management at both the landscape and field level seems to be more effective for increasing spider diversity in WF paddies. Spiders are not only important natural enemies of various insect pests (Kiritani et al. [@CR10]; Oraze et al. [@CR17]; Takada et al. [@CR28]), but also important prey for larger predators, such as frogs (e.g., Hirai and Matsui [@CR7]) and larger spiders (e.g., Ishijima et al. [@CR8]) in paddies. Taken together, the enhancement of spider diversity in WF paddy fields might achieve both biological control and biodiversity conservation in rice paddy ecosystems.

**Competing interests**

The authors declare that they have no competing interests.

**Authors' contributions**

MBT conducted field surveys and wrote the paper. ST, SI, and TM collected field data. IW designed field surveys and wrote the paper. All authors discussed the results and commented on the paper. All authors read and approved the final manuscript.

This study was made possible by the help of the farmers of the paddy fields in Osaki City and by farmers at the Osaki City Hall. We would also like to thank J. Oyama, A. Yoshioka, and Y. G. Baba for providing helpful comments and A. Tanikawa and Y. G. Baba for identification of the spiders. This study was supported by a Grant-in-Aid for Young Scientists (B) (23780049) and the 21st Century COE Program at the University of Tokyo "Biodiversity and Ecosystem Restoration Research Project".
