Intercropping of Rice and Water Mimosa ( Neptunia oleracea Lour.): A Novel Model to Control Pests and Diseases and Improve Yield and Grain Quality while Reducing N Fertilizer Application

: Cereal/legume intercropping is an effective agricultural practice for pest and disease control and crop production. However, global research on rice and aquatic legume intercropping is relatively rare. A ﬁeld experiment during two seasons (2018 late season and 2019 early season) was conducted to explore the effects of rice and water mimosa intercropping on rice canopy microclimate, pest and disease, yield, grain quality, and economic income. Two cultivation patterns including rice/water mimosa intercropping and rice monocropping were employed, and three nitrogen (N) fertilizer application levels, including zero N (ZN, 0 kg ha − 1 N), reduced N (RN, 140 kg ha − 1 N), and conventional N (CN, 180 kg ha − 1 N) levels, were applied for the above two cultivation patterns. The results showed that rice/water mimosa intercropping formed a canopy microclimate of rice with higher temperature and lower relative humidity and dew point temperature. In addition, there was a signiﬁcant reduction in the occurrences of rice leaf blast by 15.05%~35.49%, leaf folders by 25.32%~43.40%, and sheath blight by 16.35%~41.91% in the intercropping treatments. Moreover, rice/water mimosa intercropping increased rice per unit yield by 43.00%~53.10% in the late season of 2018 and 21.40%~26.18% in the early season of 2019. Furthermore, rice grain quality was totally improved, among which brown and head rice rates increased but rice chalky rate and chalkiness degree decreased in the intercropping system. We suggest that combining rice/water mimosa intercropping and N fertilizer reduction can be used as an environmentally friendly eco-farming technique because it can decrease N fertilizer application by approximately 40 kg · ha − 1 . This combination would not only mitigate nonpoint source pollution but also obtain advantages for controlling rice pests and diseases that would alleviate pesticide usage and improve rice yield and grain quality, which can be extended for green rice production to increase income for producers.


Introduction
Globally, rice (Oryza sativa L.) is one of the oldest and most important staple foods and is consumed by approximately 50% of the world's population and 60% of China's population [1][2][3][4][5]. In addition, modern, intensive agriculture has resulted in a series of grain quality, and economic income in paddy fields. In addition, previous studies on the microclimate in rice/water mimosa intercropping systems are relatively rare. Therefore, the present study was conducted to explore whether the microclimate of the rice canopy would be changed by rice/water mimosa intercropping and whether this intercropping can control pests and diseases and increase rice yield and grain quality while reducing N fertilizer application. Our hypotheses were as follows: (1) Rice/water mimosa intercropping can modify the rice canopy microclimate; and (2) rice/water mimosa intercropping can control diseases and pests and then increase yield and grain quality while reducing N fertilizer application.

Study Site and Materials
A field experiment was performed at Zengcheng Teaching and Research Farm (23 • 14 N, 113 • 38 E), South China Agricultural University, Guangzhou, Guangdong Province, China. The area has a subtropical monsoon climate with a warm winter and hot summer ( Figure S1). In addition, rice (Oryza sativa L., Huanghuazhan), a typical local crop that is widely cultivated in southern China for its good taste, grain quality, strong lodging resistance, and steady productivity, and water mimosa (Neptunia oleracea Lour.), an aquatic leguminous plant, were selected for the experiment. The basic physical and chemical properties of the soil were as follows: sandy loam, pH 4.88, containing 15.80 g kg −1 organic matter, 2.27 g kg −1 total N, 0.51 g kg −1 total phosphorous, 10.91 g kg −1 total potassium, 9.06 mg kg −1 ammonium nitrogen (N), 4.69 mg kg −1 nitrate N, 43.82 mg kg −1 available phosphorus, and 47.56 mg kg −1 available potassium.

Cultivating Experimental Design
A field experiment was conducted during 2018 and 2019 consisting of two rice growing seasons. The 2018 late growing season was from August to November, and the 2019 early growing season was from April to July. Six treatments with four replicates were applied to the experiment: rice monocropping (0, 140, and 180 kg·ha −1 N) and rice/water mimosa intercropping (0, 140, and 180 kg·ha −1 N). A completely randomized design was used in the field experiment. Each plot area was 35 m 2 (5 m × 7 m) and irrigated and drained independently. The seedlings' density of rice was 250,000 holes/ha and 166,667 holes/ha in monocropping and intercropping treatments. The density of water mimosa was 83,333 plants/ha in intercropping treatments. The cultivation standard per strip in the intercropping treatments was performed every four rows of rice and then connected to three rows of water mimosa. Within the rice monocropping treatments, rice row spacing was 0.2 m. In intercropping treatments, the intrarow spacing of rice and water mimosa was 0.2 m and 0.15 m, respectively, and the row distance between rice and water mimosa was 0.25 m (Figure 1). The area ratio of rice and water mimosa was 2:1 in the intercropping treatment. There were 2/3 and 1/3 areas of each plot occupied by rice and water mimosa in intercropping treatments, respectively. Urea (CO(NH 2 ) 2 ) was applied as N fertilizer before transplanting and at the rice tillering, heading, and filling stages with proportions of 40%, 20%, 30%, and 10%, respectively ( Figure 1). In addition, the application method of phosphate and potassium fertilizer was the same in each treatment. Calcium superphosphate (P 2 O 5 12%), as the phosphate fertilizer, was applied only as the base fertilizer at 45 kg ha −1 . Potassium chloride (KCl 60%) was used as potassium fertilizer at 135 kg ha −1 ; one half was applied as the base, and the other half was applied as heading stage fertilizer. Seeds were soaked in water for 24 h at room temperature and then germinated under moisture conditions for seedling preparation. Germinated seeds were sown on 25 July (late growing season of 2018) and on 10 March (early growing season of 2019) for rising nurseries. In addition, 15 kg ha −1 urea fertilizer was applied during the rice seedling stage. Water mimosa (length of 0.3 m) and rice seedlings (three-leaf-and-one-leaflet stage) were then concurrently transplanted to the paddy field. Fertilization management methods were the same for the two seasons. The drainage and irrigation managements were to keep the water layer at 6~8 cm throughout the rice-growing season in the field because water mimosa grows better in water flooded condition, but irrigation was stopped 1 week before the rice harvest. We did not apply any pesticides, herbicides, or weed control practices for any treatments.

Canopy Microclimate
The canopy microclimate of rice in the monocropping and intercropping treatments under the reduced N fertilizer level was measured by a small, portable climate instrument UA-002064 (Onset, Bourne, MA, USA). On the 50th, 55th, 60th, 65th, 70th, and 75th days after transplanting, daily changes of dew point temperature, air temperature, and relative humidity (RH) in the rice canopy were measured for the four plots in each treatment. A small, portable climate instrument was placed among the rows of rice. The instrument was tied on a pole at the height of the rice canopy, and its height was changed according to the height of the rice canopy. The data were recorded every 20 min, and then the dew point temperature, air temperature, and RH value were calculated to obtain an average value per hour.

Pest and Disease Survey
Rice leaf blast, leaf folders, and sheath blight were investigated on the 50th, 55th, 60th, 65th, 70th, and 75th days after transplanting in each season. Three holes of rice in different rows were selected randomly as the initiation points, and then five holes were extended in each row when using a parallel jumping sampling method. The calculation methods of rice leaf blast, leaf folders, and sheath blight were as follows.
The damaged leaf number caused by rice leaf blast and leaf folders and the total leaf number of the 15-hole rice were counted. The incidences of rice leaf blast and leaf folders were calculated by using the following Equation (1) [1]: The incidence of rice leaf blast (leaf folders) (%) = damaged leaf number/total leaf number × 100% (1) Rice sheath blight was surveyed, according to Wu et al. (2012) [26], and the disease categories of rice sheath blight were as follows: The numbers of infected stems in different disease categories and total stems were counted. The disease index was calculated according to Equation (2)  Rice and water mimosa yields were measured from four plots with an area of 1 m × 1 m for each replicate, and each sample contained insiders and outsiders of the plants. In addition, rice was threshed manually and sun-dried (adjusted to 14% moisture content) to obtain the grain yield, and the water mimosa yield was determined by harvesting the top tender part. Per unit yield indicated the rice yield which was calculated in the 3 × 1 m 2 plots, and the actual yield of rice and water mimosa represented the yield considering the area ratio of rice and water mimosa in the intercropping treatment. Yield components were measured according to   [27] and determined from five hills for each replicate. Effective panicle numbers per hill were counted and calculated per unit effective panicle numbers and actual panicle numbers. Then, the panicles were threshed manually and divided into filled grains and unfilled grains. Five subsamples of filled and unfilled grains were used to estimate the grain numbers per panicle and seed-setting rate. The 1000-grain weight was also calculated from the sampled grains.
Rice grain quality was measured according to   [27] and determined from five samples for each replicate. A 100-g sample of rice grain was passed through a dehusker for polishing and then divided into broken and unbroken grains. The brown, milled, and head rice rates were calculated as the percentages of the total (100 g) rice grains. Amylose content, soluble protein content, and alkali value were measured using an Infratec-1241 grain analyzer (FOSS-TECATOR, Hilleroed, Denmark). The chalky rice rate, chalkiness degree, and length/width were scanned with Plant Mirror Image Analysis (MICROTEK, Hsinchu, Taiwan, China), and then the resulting images were processed with SC-E software (HangzhouWanshen Detection Technology Co., Ltd., Hangzhou, Zhejiang, China).
Net income for each treatment was estimated using Equation (3) [28]: Here, the total cost of cultivation comprised the costs of inputs and labors. The costs of inputs (seeds and fertilizers) were based on the local market prices. Water mimosa is a perennial herb that can be reproduced in the field so that its plant prices are free [20,21]. The costs of labors were calculated by the cultural activities (land preparation, seedling, transplanting, applying fertilizer, and harvesting) and paid at the rate of 120 yuan (Chinese yuan) person-day −1 of 8 h. Gross income was calculated as the total value of economic yield (water mimosa and rice) per treatment. The market prices of rice and water mimosa were 4 yuan kg −1 and 7 yuan kg −1 , respectively.

Data Analysis
Data are expressed as the mean value ± standard error and subjected to two-way ANOVA in the two cultivation patterns with three different N fertilizer levels (p < 0.05). In addition, treatment differences in the same cropping system were statistically assessed by using Duncan's method of one-way ANOVA for multiple comparisons when the data met the normality and homoscedasticity hypotheses; otherwise, the data were evaluated through a Games-Howell method (p < 0.05). The differences between intercropping and monocropping were statistically evaluated by independent T-tests (p < 0.05). All statistical analyses were performed using SPSS 17.0. The data were tested by the Pearson correlation analysis using the "corrplot" package in R [29].

Rice Canopy Microclimate
On the 50th, 55th, 60th, and 75th days in the late season of 2018 and on the 55th, 65th, and 75th days in the early season of 2019, the air temperatures were significantly higher under intercropping treatments than under the monocropping treatments, by 0.07~2.30 • C ( Figure S2a,b). In contrast, in the late season of 2018, in the middle of the day, the dew point temperature and RH in the intercropping treatments were generally significantly lower, by 0.07~1.77 • C and 0.88~11.85% compared with those in the rice monocropping treatments, respectively ( Figure S2c,e).

Rice Pest and Disease
During the whole investigation, intercropping treatments reduced the occurrence of rice pests and diseases. For instance, the incidences of rice leaf blast in intercropping treatments were generally significantly lower than those in monocropping treatments except on the 50th day after transplanting in the late season of 2018 (Figure 2a,b, Table S1). Intercropping treatments also significantly reduced the occurrence of rice leaf folders compared with the monocropping treatments on the 60th~75th days after transplanting in both seasons (Figure 2c,d, Table S3). Moreover, there was a significant reduction in the incidence of rice sheath blight in the intercropping treatment on the 55th and 70th days after transplanting in the 2018 late season and on the 60th, 65th, and 75th days after transplanting in the 2019 early season (Figure 2e,f, Table S5). Additionally, the incidences of rice leaf blast, leaf folders, or sheath blight in 2018 late season generally were significantly higher than those in 2019 early season after the 70th day of transplanting (Tables S2, S4 and S6).
More importantly, the incidences of rice leaf blast, leaf folders, and sheath blight in the intercropping with reduced N treatments were significantly lower than those in the monocropping with zero, reduced, and conventional N treatments (Tables S1, S3 and S5). No significant difference was found between the N fertilizer treatments in pests and diseases.

Rice Yield and Yield Components
There was a significant interactive effect between cultivation pattern and N fertilizer on rice yield and per unit effective panicle number (Table S7). The per unit yield, actual yield, and per unit effective and actual panicle number of rice and the yield of water mimosa in the CN and RN treatments were significantly higher than those in the ZN treatments. In addition, intercropping treatments significantly increased the rice's per unit yield and per unit effective panicle numbers by 43.00%~53.10% and 22.15%~41.41% in the late season of 2018 and 21.40%~26.18% and 15.72%~41.54% in the early season of 2019, respectively. However, no significant difference was found in grain numbers per panicle, seed-setting rate, or 1000-grain weight (Table 1). Additionally, rice yield and yield components in 2019 early season were generally significantly higher than those in 2018 late season except the grains per panicle (Table S8).
Furthermore, compared with rice monocropping with conventional N application treatments, intercropping with reduced N treatments significantly increased the rice's per unit yield by 35.64% in the 2018 late season and increased the per unit effective panicle numbers by 17.01% in the 2019 early season (Table 1).

Rice Grain Quality
The interactions between the two factors (cultivation pattern and N fertilizer) had significant effect on rice chalky rice rate and chalkiness degree (Table S9). Compared with the ZN treatments, the RN and CN treatments significantly increased the brown rice rate, protein content, length/width, and alkali value of rice in the early season of 2019 ( Table 2). As our prediction, brown and head rice rates in intercropping treatments were slightly higher than those in monocropping treatments in the late season of 2018. There was a significant reduction in chalky rice rate and chalkiness degree in intercropping treatments compared with monocropping treatments during the two seasons (Table 2). No significant difference was found in the milled rice rate or amylose content (Table 2). Additionally, rice grain quality in 2019 early season was generally significantly higher than that in 2018 late season except the chalky rice rate, chalkiness degree, and length/width (Table S10).
The chalky rice rate in the intercropping with reduced N treatment was significantly lower than that in the monocropping with the conventional N treatment ( Table 2). In addition, the length/width and alkali value in the intercropping with zero N treatments were significantly lower than those in the monocropping with the conventional N treatment in the early season of 2019 (Table 2).

Correlation Analysis
In the late season of 2018, the incidences of rice blast, leaf folders, and sheath blight generally had significant negative correlations with the air temperature, rice per unit yield, per unit effective panicle number, and head rice rate, but had significant positive correlations with the dew point temperature, RH, chalky rice rate, and chalkiness degree ( Figure S3a). Although the incidences of pathogens did not have obvious correlations with microclimatic parameters, the incidences of pathogens had negative correlations with rice per unit yield, effective number, and head rice rate and had significant positive correlations with chalky rice rate and chalkiness degree in the early season of 2019 ( Figure S3b).

Economic Analysis
In the two seasons, intercropping with reduced N treatments had the highest economic income and accounted for 18,234 and 18,830 yuan ha −1 in the 2018 late season and 2019 early season, respectively (Table 3). Thus, combining rice/water mimosa intercropping and N reduction could obtain the maximum benefit for farmers and producers.
The prices of rice and water mimosa are 4 and 7 yuan kg −1 , respectively. ZN, RN and CN indicate the 0, 140 and 180 kg ha −1 N fertilizer levels, respectively. "Mono"and "Inter" denote the mono-cropping and intercropping treatments. All the data presented are the means of four replicates ± standard errors. Different capital letters indicate significant differences between cultivation patterns under the same N fertilizer level (p < 0.05). Different lowercase letters indicate significant differences between N fertilizer levels in the same cultivation pattern (p < 0.05). Actual yield and effective panicle denote the indexes according to the area ratio of the two crops in the intercropping treatment. ZN, RN and CN indicate the 0, 140 and 180 kg ha −1 N fertilizer levels, respectively. "Mono"and "Inter" denote the monocrop-ping and intercropping treatments. The "*" represents a significant difference between mono-cropping with the conventional N treatments and intercropping with zero and reduced N treat-ments (p < 0.05). All the data presented are the means of four replicates ± standard errors. Different capital letters indicate significant differences between cultivation patterns under the same N fertilizer level (p < 0.05). Different lowercase letters indicate significant differences between N fertilizer levels in the same cultivation pattern (p < 0.05). ZN, RN and CN indicate the 0, 140 and 180 kg ha −1 N fertilizer levels, respectively. "Mono"and "Inter" denote the monocropping and intercropping treatments. The "*" represents a significant difference between monocropping with the conventional N treatments and intercropping with zero and reduced N treatments (p < 0.05).

Canopy Microclimate Modifying
Intercropping is a simple and effective way to modify the microclimate of rice canopies [30][31][32]. In our two-season field experiments, intercropping treatments increased the canopy air temperature compared with the rice monocropping treatments ( Figure S2). These results were similar to those of a previous study that showed that canopy temperature increased in maize/cowpea intercropping systems relative to sole maize cropping [33,34]. Additionally, rice/water mimosa intercropping reduced the dew point temperature and relative humidity of rice canopy ( Figure S2). Similarly, many studies have demonstrated that intercropping reduces the canopy relative humidity and dew point temperature [30,31]. In the rice/water mimosa intercropping system, water mimosa prostrated growth, rice erected growth, and the height of rice was much higher than that of water mimosa [21]. It is known that higher plants have stronger competition for light in intercropping systems, and high irradiation is usually combined with high temperature and low humidity [14,32]. In addition, the height difference between the two crops exhibited better air circulation. Moreover, this intercropping system had a higher air temperature and lower relative humidity and dew point temperature and, thus, exhibited a special canopy microclimate of rice as a result (Figures 3 and S2).

Pest and Disease Control
It is believed that intercropping usually reduces the occurrence of harmful organisms due to the high biodiversity and stability of the ecosystem [1]. In the present study, rice/water mimosa intercropping reduced the incidences of rice leaf blast, leaf folders, and sheath blight (Figure 2, Tables S1, S3 and S5). A similar study on rice and water spinach intercropping in paddy fields also showed that pests and diseases were substantially lower in intercropping systems [4]. In addition, our study also showed that the occurrence of pathogens had a positive correlation with dew point temperature and relative humidity and a negative correlation with air temperature in the late season of 2018 ( Figure S3). High air temperature inhibited rice leaf roller egg hatching and reduced the longevity of this insect [1]. Hence, the high air temperature in this intercropping system might contribute to reducing the breeding of rice folders. Universally, dew is formed at night, and relative humidity can control fungal spores to germinate the conidial sporulation and dispersal process of pathogens [31,33]. Thus, the high temperature and low dew point and humidity probably had some functions to control pathogens in the present study. According to the characteristics of the rice pest and disease, we speculated and summarized that the causes of rice/water mimosa intercropping reduced the incidences of pathogens would be as follows. (1) The strip distribution in the intercropping system established a physical barrier that blocked the horizontal spread of the pathogens [1].
(2) Microclimatic factors such as high temperature, low relative humidity and dew point temperature, stronger solar radiation, and better air circulation in the intercropping system had some functions in resisting pests' and diseases' activity, such as insect migration and reproduction [4,34]. (3) Water mimosa was neither the food nor the host of the pathogens; hence, the density of pathogens was relatively diluted. Thus, the incidences of pests and diseases reduced in the rice/water mimosa intercropping system ( Figure 3).
Moreover, previous studies revealed that N fertilizer usually increases the disease susceptibility of rice [35][36][37]. However, in our study, there were no significant differences between N fertilizer levels in the incidences of rice leaf blast, sheath blight, and leaf folders. The reason for this phenomenon was probably that the application of N fertilizer was insufficient to increase the disease susceptibility of rice in the present study because 180 kg ha −1 N was the optimal N fertilizer for rice growth in paddy fields [1,4,[35][36][37]. We also found that the pests and diseases in the monocropping with conventional N treatments were significantly higher than those in the intercropping with zero, reduced, and conventional N treatments (Figure 2, Tables S1, S3 and S5). Therefore, the reduction in rice pest and disease rates in rice/water mimosa intercropping could mainly be attributed to the border effect.

Yield and Yield Components Advantages
Environmental conditions and management methods play a crucial role in yield and yield components [2]. The increase in rice per unit yield in the early season of 2019 might be attributed to the higher rainfall compared with that in the late season of 2018 ( Figure S1 and Table S8). This is consistent with a previous study that showed that sufficient irrigation could increase crop yield [38]. In the present study, we observed that the interaction between cultivation pattern and N fertilizer significantly affected the rice per unit yield and effective panicle numbers in the two seasons (Table S7). These results are in line with previous research that showed that the interactive effect between the intercropping pattern and N fertilization significantly affected plant growth in a groundnut (Arachis hypogaea Linn.) and sesame (Sesamum indicum Linn.) intercropping system [39]. In the present research, compared with zero N treatment, the reduced and conventional N treatments substantially increased the yield of rice and water mimosa ( Table 1). As a comparison, previous research found that crop yield increased with N fertilizer levels [40,41]. It is acknowledged that higher grain yield was obtained due to higher panicle numbers [40]. Similar to rice yield, reduced and conventional N treatments also had higher rice effective panicle numbers than zero N treatments (Table 1). Ju et al. (2019) [24] also reported that N fertilizer application increased the panicle numbers of rice. Consequently, both reduced and conventional N treatments were beneficial to rice yield and effective panicle numbers.
Intercropping generally increased the crop yield. In the present study, intercropping treatments increased rice per unit yield by 43.00%~53.10% in the late season of 2018 and 21.40%~26.18% in the early season of 2019. Meanwhile, per unit effective panicle numbers were higher in rice/water mimosa intercropping than in rice monocropping (Table 1). Our results were consistent with previous studies. Ning et al. (2017) [1] and Liang et al. (2016) [4] both found that rice/water spinach intercropping increased the rice yield and effective panicle numbers. In addition, intercropping treatments significantly increased the yield, while there was no significant difference in grain components between cultivation patterns during most of the time ( Table 1). The reason probably was that rice had higher per unit effective panicle numbers among the intercropping treatments, and the effective panicle number was one of the yield components. The increase in rice per unit yield in intercropping treatments could also be attributed to the promoting effect of water mimosa N fixation [23,42]. Moreover, the rice per unit yield and effective panicle numbers in the intercropping with reduced N treatments were higher than those in the monocropping with conventional N treatments (Table 1). These findings suggested that rice/water mimosa intercropping combined with reduced N fertilizer application would be a good practice choice to increase rice yield in paddy fields.

Grain Quality and Economic Incomes Improvement
Grain quality comprises grain appearance, milling, cooking, eating, and nutritional qualities and depends on genetic, environmental, and crop management factors [27,43]. Our results indicated that N fertilizer improved rice grain quality. For instance, the brown rice rate, protein content, length/width, and alkali value in the reduced and conventional N treatments were higher than those in the zero N treatments (Table 2). These results are consistent with previous studies which found that crop quality indexes improved with N fertilizer rates [25,44,45].
Moreover, intercropping had a positive and significant effect on brown and head rice rates (Table 2). An increase in crop quality in intercropping systems was also reported by Yusef et al. (2014) [37]. Low chalkiness is often associated with more translucent rice grain and represents a higher value and price of rice in markets [46,47]. In the present study, rice/water mimosa intercropping resulted in a significant reduction in the chalky rice rate and chalkiness degree ( Table 2). The lower chalky rate and chalkiness degree of rice probably bring more profit to farmers and producers. In the rice/water mimosa intercropping system, the reason that led to better grain quality could be attributed to the following aspects: (1) Organic fertilizer and biofertilizer supplied by water mimosa could improve grain quality by improving photosynthesis and nutrient uptake by cereals, and these nutrients were eventually transported to the seed and contributed to improving grain quality ( Figure 3) [48]. (2) Rice grain quality was associated with microclimatic parameters ( Figure S3). Similarly, intercropping had higher temperature and good ventilation and light transmission, which could lead to less inhomogeneous substances, such as white opacity and chalkiness in grains [49]. (3) Our study found that the incidences of pathogens had negative correlations with the head rice rate and had positive correlations with the chalky rice rate and chalkiness degree ( Figure S3). These results are consistent with a previous study which showed that lower pathogens in rice resulted in better grain quality ( Figure 3) [50,51].
Notably, compared with monocropping with conventional N treatments, intercropping with zero N treatments resulted in lower length/width and alkali values. However, intercropping with reduced N treatment decreased the chalky rice rate relative to monocropping with conventional N treatments (Table 2). Therefore, we could infer that intercropping with reduced N treatment is an optimal choice for rice grain quality.
For economic analysis, we found that the intercropping system with reduced N treatments had the highest economic income (Table 3). Higher net income has also been reported under potato/legume and sorghum/legume intercropping systems [28,52]. Rice/water mimosa intercropping with reduced N treatments achieved the maximum economic income which would be beneficial for farmers and producers. However, there were only two seasons of performance assessment in this study, and further long-term experiments are needed in the future.

Conclusions
Combining rice/water mimosa intercropping and N fertilizer reduction could not only decrease N fertilizer application by approximately 40 kg ha −1 , which would mitigate nonpoint source pollution, but also form a canopy microclimate of rice with higher temperature, lower relative humidity and dew point temperature, and better control of rice pests and diseases to a certain extent, which would alleviate pesticide usage, finally increasing rice productivity, grain quality, and economic income, and that can be extended for green rice production to increase income for farmers and producers. In a word, our study provided an environmentally friendly eco-farming technique in rice production.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/agriculture12010013/s1, Figure S1: Rainfall (RF) and average temperature (AT) per month during the two seasons in Zengcheng, Guangzhou, Guangdong, China. Figure S2: Air temperature, relative humidity (RH) and dew point temperature of rice canopy in the monocropping and intercropping under the reduced N fertilizer application level. Figure S3: Pearson correlations among daily air temperature (AT), dew point temperature (DPT) and relative humidity (RH) of rice canopy at 50th7 5th after transplanting, rice leaf blast (RLB), rice leaf folders (RLF) and rice sheath blight (RSB) at 75th day after transplanting, rice yield (RY), per unit effective panicle number (EPN), head rice rate (HRR), chalky rice rate (CRR) and chalkiness degree (CD) at maturity stage. Table S1: Incidence (%) of rice leaf blast in monocropping and intercropping treatments with the three different N fertilizer application levels. Table S2: Variance in incidence of rice leaf blast between the 2018 late season and 2019 early season. Table S3: Incidence (%) of rice leaf folders in monocropping and intercropping treatments with the three different N fertilizer application levels. Table S4: Variance in incidence of rice leaf folders between the 2018 late season and 2019 early season. Table S5: Disease index (%) of rice sheath blight in monocropping and intercropping treatments with the three different N fertilizer application levels. Table S6: Variance in incidence of rice sheath blight between the 2018 late season and 2019 early season. Table S7: Interaction effects between the two factors (cultivation pattern and N fertilizer) on rice yield and yield components by two-way ANOVA. Table S8: Variance in rice yield and yield components between the 2018 late season and 2019 early season. Table S9: Interaction effects between the two factors (cultivation pattern and N fertilizer) on rice grain quality determined by two-way ANOVA. Table S10: Variance in rice grain quality between the 2018 late season and 2019 early season.