Long-Term E ﬀ ects of Fertilizers with Regional Climate Variability on Yield Trends of Sweet Corn

: Climate change a ﬀ ects global crop production year after year. Exploring the impact of di ﬀ erent fertilization methods on crop yield stability has become an extremely important topic in sustainable agriculture. The objective of this study is to explore the e ﬀ ects of various fertilization regimes with climate variability on yield stability for sweet corn production in southern Taiwan. Three fertilization treatments composed of chemical fertilizer only (CF), integrated fertilizer (half organic / half chemical fertilizer) (IF), and organic fertilizer only (OF) were implemented from 2009 to 2018 based on the well-maintained soils since 1988. While the same amounts of these fertilizers were applied during the period, we found that di ﬀ erent fertilization changed the marketable yields of fresh fruit (ear), which slightly increased for organic fertilizer, but substantially decreased for both chemical ( p = 0.0001) and integrated ( p = 0.0061) fertilizer. Thus, based on these 10 years of observation, yields among fertilization treatments were analyzed with weather and soil parameters to determine the possible factors involved. Both multiple linear regression equation ( p < 0.0001, adj. R 2 > 0.57) and regression tree analysis illustrated signiﬁcantly negative correlations between average ear weight and relative humidity under the chemical fertilizer treatment. In this study, we show for the ﬁrst time that chemical fertilizer had the lowest yield resilience in response to regional relative humidity change compared to organic and integrated fertilizers. Our results also indicate that speciﬁc soil microbes have the potential to help sweet corn face environmental vulnerability in subtropical regions. C.-N.C.; preparation, P.-F.H.;


Introduction
Yield stability is one of the critical indicators for agricultural sustainability. Environmental complexity includes weather, soil, fertilization, and tillage methods diversely influence stable crop production [1]. Since the 1990s, the observed trend of global climate warming has become a great challenge to crop production in many locations around the world [2]. Soil quality has also suffered from degradation and chemical abuse in many places throughout the world [3]. Until now, these abiotic stresses continue to impact crop production, posing a threat to production systems.
In Taiwan (geographical position in Southeast Asia: latitude, 23 • 33 N; longitude, 121 • 01 E; koppen climate type), climate warming accelerated conspicuously at a rate of 0.29 • C per decade, and collected from 1990 until the present at an agricultural meteorological station: 72 V14 (22.86 • N, 120.51 • E) (Central Meteorological Administration, Kaohsiung, Taiwan), which is 24 m above sea level. Temperature records were divided into average temperature (Tave), maximum temperature (Tmax), minimum temperature (Tmin). The linear distance from the meteorological station to the experimental field is less than 1 km. We transformed the daily weather data as average data for fitting the spring cropping season growth period of sweet corn (February 1-May 31) for Pearson correlation and regression analysis during the period of 2009-2018. The range between daily Tmax and Tmin was defined as Tran.

Sweet Corn Experimental Field
The original experimental field soil was hyperthermic, udic, haplaquept, mixed, and calcareous, with a silty loam texture measured in 1988. The original main plots field was designed for two kinds of crop rotation systems. Different crops were planted in rotation systems annually, but the same variety of sweet corn was planted only in an upland rotation system in every spring cropping season of 2009-2018. Commercial sweet corn seed (Zea mays scharata var. wha-chen), a non-genetically modified organism, was purchased from Known-You Seed Co., Ltd. in Taiwan. The planting space of each sweet corn strip is 80 cm wide and 30 cm long. The plot size of each fertilization treatment was around 0.1 hectare (CF: 880 square meters; IF: 881 square meters; OF: 864 square meters). From 2001 to present, no chemical pesticides or herbicides have been applied to any of the three fertilization treatment plots. In general, the corn was planted around February 10, and first and second fresh ear were harvested at approximately 75 days and 84 days (about April 25 and May 5), respectively, after seeding. All crop residues were plowed back into the soil after harvest, regardless of fertilization treatments.

Fertilization
Each main rotation system contained three fertilization treatments: application of organic fertilizer only (OF), chemical fertilizer only (CF), and a combined application of half organic fertilizer and half chemical fertilizer, based on the same amount of nitrogen (N) as in OF and CF, called integrated fertilizer (IF). The amounts of chemical fertilizers and composted manure (or commercial organic fertilizer) applied to the sweet corn were described in a previous paper [22]. The application rate of chemical fertilizer for sweet corn in the chemical fertilizer treatment is shown in Table 1. Half of the N and K, and all of the P in the chemical fertilizer treatment, were banded beneath the surface as basal fertilizer and applied one day before seeding. Table 1. Amount of chemical fertilizers (N, P, and K added as urea, superphosphate, and potassium oxide, respectively) and composted manure applied to sweet corn production.
The remaining half of the N and K was equally split and used as a top dressing on the 30th day after seeding. The application rate of organic fertilizer was based on the assumption that 50% of the N in the organic fertilizer would become available to plants during the cropping season. Therefore, vegetable, farmyard, and bone dust mixture compost as organic fertilizer had been applied from 2005 to present. Each corresponding treatment is kept constant for every spring cropping season. Before this experiment set a record, each fertilizer treatment had been applied over 20 years in the same area as uniform as possible to obtain the homogeneous soil properties. Three continuous cropping schedules each year also provided reliable support for farmland management. As a result, we believe it reasonable to assume any significant different yield trends observed between plots to have been caused by the different fertilizer treatments. For instance, this point of view was also carried out by Yang et al., who conducted a series of analysis for such long-term data with sampling replicates [23].

Soil Analysis
Soil samples: Sampling was done in February, after the winter crop was harvested. The soil was sampled from the superficial layer (0-15 cm) and passed through 10 mesh sieves to remove stones and plant debris. Each plot was divided into four subplots for sampling among fertilizer treatments for soil properties each year. Soil chemistry: The methods for measuring soil chemical properties were in accordance with previous reports [24,25]. Ionic form of K, Ca, Mg, Fe, Mn, Cu, Zn, Na was detected by inductively coupled plasma optical emission spectrometer (ICP-OES) (Agilent Technologies, Inc., Santa Clara, USA). Soil microbial community: To investigate the relation between yield and microbes, the four sampling average ear weights in 2014 were further correlated with the previous four sampling microbial data from metagenomic sequences in the same year [15].

Yield Sampling and Analysis
Each replicate consisted of one strip (27 square meters, approximately 112-115 plants). All fresh fruit were harvested from each strip and every fruit was weighed. Above 200 g per fruit was characterized as qualified fruit, otherwise it was defined as unqualified. Unqualified fruit was ultimately grouped by insect pest or disease. Total N from the harvest sweet corn leaves were determined by combustion (PerkinElmer 2400 II elemental analyzer, PerkinElmer, Shelton, CT, USA).

Statistical Analysis
The SAS-EG software package (v7.1) (SAS Institute Inc., Cary, USA) was used to perform all statistical analyses. To minimize the heterogeneity of soil and climate change over space and time, 10 years of continuously observed data for estimating the trends and correlations of these parameters were considered.
Statistical analysis was performed by ANOVA, followed by multiple range test to determine the differences of sweet corn yields and soil properties among treatments. Principal component analysis was performed based on the soil properties (12 parameters-pH, OM, P, K, Ca, Mg, Fe, Mn, Cu, Zn, Na, EC) collected in 2009 and 2018 to confirm whether the soils were distinguishable and representative. Least-squares linear regressions of climate parameters, soil properties, and sweet corn yield against time series were calculated to test the hypothesis that, throughout the experimental period, these variables did not change. The relationships among yield, climate parameters, and soil properties were evaluated with simple Pearson correlation analyses. To predict sweet corn yield by extension, multiple linear regression analysis was used with the backward elimination procedure to identify the subset of climate parameters and soil properties.
To further investigate the effects of the climate parameters and soil properties on sweet corn yield under CF, a regression tree analysis was performed.

Weather Trends
Daily climate parameters were collected and divided into three cropping seasons: spring (1 February-31 May), summer (1 June-30 September), and autumn (1 October-31 January), to match the growth of sweet corn in the spring cropping season. The trends of Tave (R square = 0.0991), Tmin (R square = 0.1093), and Tmax (R square = 0.0978) did not change significantly. Relative humidity and sunshine hours ( Figure 1A were uneven distribution patterns of precipitation among these periods, and especially heavy rainfall occurred in 2013 (7.04 mm per day in 2013 compared to 2.38 mm per day on average in 2009-2018) ( Figure 1B). The linear trend equation for sunshine hours had a marked rising trend of 0.1175 h per year, while the one for relative humidity had a moderately increasing trend of 0.5778% per year. Although the relationship between sunshine hours and relative humidity was not confirmed in this study, decreasing wind speed was observed in Taiwan [4], which might be the cause of the increase in relative humidity.
Sustainability 2020, 12, x FOR PEER REVIEW 5 of 18 ( Figure 1B). The linear trend equation for sunshine hours had a marked rising trend of 0.1175 h per year, while the one for relative humidity had a moderately increasing trend of 0.5778% per year. Although the relationship between sunshine hours and relative humidity was not confirmed in this study, decreasing wind speed was observed in Taiwan [4], which might be the cause of the increase in relative humidity.

Soil Properties
Most of the soil chemical properties show stepwise feature from highest to lowest in the following order: OF, IF, and CF, except elements of Fe, Cu, and EC. Meanwhile, these elements had significant differences among different fertilizer treatments from 2009 to 2018.

Soil Properties
Most of the soil chemical properties show stepwise feature from highest to lowest in the following order: OF, IF, and CF, except elements of Fe, Cu, and EC. Meanwhile, these elements had significant differences among different fertilizer treatments from 2009 to 2018.
Treatment had a significantly higher pH value on average and higher accumulation of organic matter, P, K, Ca, Mg, Fe, Mn, Cu, Zn, Na, and EC in OF than in CF soil. Compared with the original soil in 1988, the OF soil was more enriched in organic matter (+61%), P (+37%), K (+196%), Ca (+313%), and Mg (+313%), but CF soil was depleted in organic matter (−24%) and P (−31%) in 2018 (Table 2). In Figure 2, there were clear trends of decreasing pH and EC in CF, IF, and OF soils. Although natural soil is proverbially believed to be formatted by acid rain, the real causes of the decrease of soil pH values in this experiment are uncertain. Acid rain (pH < 5.0) was frequently (21% of observation days) recorded in Kaohsiung in 2018 from the EPA (Environmental Protection Administration), Taiwan, so it is a possible cause of increasing soil acidity. Even though application of chemical fertilizer has obvious effects on soil acidification [19], it is notable that the acidification rate of the organic fertilizer treatment was the least steep among the three; OF was −0.0529 pH units per year, compared to −0.1009 pH units per year for CF ( Figure 2A). That is, the CF acidification rate was almost twice as severe as the OF rate. In addition, EC is an integral factor in soil fertility for agricultural crops. Studies' reports indicate that extractable amounts of Na and Mg obviously affect the EC [26], as does as the soil texture [27]. The decreasing EC trend may be the cause of increasing relative humidity or the integrated intention of acidification of these three soil treatments. The negative correlation between EC and relative humidity is significant in all three fertilization treatment soils (Pearson correlation coefficient: −0.47743~−0.60472, p ≤ 0.0018), which means the higher the relative humidity, the lower the EC. This point of view is somewhat of a concern to Krug, who emphasizes that, in acidification, acid rain is responsible for leaching of nutrients, release of cation [28].
Sustainability 2020, 12, x FOR PEER REVIEW 1 of 18 In Figure 2, there were clear trends of decreasing pH and EC in CF, IF, and OF soils. Although natural soil is proverbially believed to be formatted by acid rain, the real causes of the decrease of soil pH values in this experiment are uncertain. Acid rain (pH < 5.0) was frequently (21% of observation days) recorded in Kaohsiung in 2018 from the EPA (Environmental Protection Administration), Taiwan, so it is a possible cause of increasing soil acidity. Even though application of chemical fertilizer has obvious effects on soil acidification [19], it is notable that the acidification rate of the organic fertilizer treatment was the least steep among the three; OF was −0.0529 pH units per year, compared to −0.1009 pH units per year for CF ( Figure 2A). That is, the CF acidification rate was almost twice as severe as the OF rate. In addition, EC is an integral factor in soil fertility for agricultural crops. Studies' reports indicate that extractable amounts of Na and Mg obviously affect the EC [26], as does as the soil texture [27]. The decreasing EC trend may be the cause of increasing relative humidity or the integrated intention of acidification of these three soil treatments. The negative correlation between EC and relative humidity is significant in all three fertilization treatment soils (Pearson correlation coefficient: −0.47743~−0.60472, p ≤ 0.0018), which means the higher the relative humidity, the lower the EC. This point of view is somewhat of a concern to Krug, who emphasizes that, in acidification, acid rain is responsible for leaching of nutrients, release of cation [28]. At the beginning and end of the study in 2009 and 2018, respectively, the chemical properties of the soil were keeping it steady and distinguishable from the three fertilizers' soil via principal component analysis (Figure 3), that is, the soil samples and sweet corn yields from different fertilizer treatments belonging to CF, IF, and OF in their representatives. At the beginning and end of the study in 2009 and 2018, respectively, the chemical properties of the soil were keeping it steady and distinguishable from the three fertilizers' soil via principal component analysis (Figure 3), that is, the soil samples and sweet corn yields from different fertilizer treatments belonging to CF, IF, and OF in their representatives. Sustainability 2020, 12, x FOR PEER REVIEW 2 of 18

Sweet Corn Yield
There are five survey items of sweet corn yield. The marketable yield of sweet corn was composed of qualified and unqualified ear yield, which can be further separated as first ear and second ear. Although these yields were not significantly different during the period of 2009-2018, the sources of variation belong to "Years," and "Treatments × Years" had significantly different yields among treatments (Table 3, Table 4).  Years (Y) *** *** *** *** *** T × Y *** *** *** *** ***
Consequently, the yield among treatments with each year should be separately compared. Therefore, the results of yield survey showed that the year 2013 was a turning point. Fresh marketable yield and components of sweet corn yield and weight presented significant differences in CF treatment in the 2009-2013 period, while OF treatment presented a statistical significance in 2013-2017 (Table 5). Some anomalies and events, such as typhoons in August-December, have arisen as an inferred result of influences on maize growth in February-May of the next year of the 2013-2017 period. Table 5. Analysis of variance between CF and OF for fresh marketable yield and components of sweet corn yield and weight # . A further comparison of yield trends on scatter plot among treatments was displayed in Figure 4. The marketable yield trends of three fertilization treatments had different distribution patterns; the trend of marketable yield was significantly decreased under CF (p = 0.0001) and IF (p = 0.0061), but slightly increased under OF. significance in 2013-2017 (Table 5). Some anomalies and events, such as typhoons in August-December, have arisen as an inferred result of influences on maize growth in February-May of the next year of the 2013-2017 period. Table 5. Analysis of variance between CF and OF for fresh marketable yield and components of sweet corn yield and weight # . A further comparison of yield trends on scatter plot among treatments was displayed in Figure  4. The marketable yield trends of three fertilization treatments had different distribution patterns; the trend of marketable yield was significantly decreased under CF (p = 0.0001) and IF (p = 0.0061), but slightly increased under OF.  The yield was further evaluated by dividing marketable yield into qualified yield of first ear and second ear ( Figure 5A). Meanwhile, the distribution pattern of average ear weight of first ear and second ear is displayed in Figure 5B.

Sources of Variation
The yield was further evaluated by dividing marketable yield into qualified yield of first ear and second ear ( Figure 5A). Meanwhile, the distribution pattern of average ear weight of first ear and second ear is displayed in Figure 5B. The qualified first ear yield trends slightly declined for OF (p = 0.3884), but significantly decreased for IF (R square = 0.3094, p = 0.0002) and CF (R square = 0.4978, p < 0.0001). Comparatively, the second ear qualified yield trends notably rose in OF (R square = 0.2051, p = 0.0033), but significantly declined in CF (R square = 0.2284, p = 0.0018), and did not have clear trends in IF (R square = 0.0113, p = 0.5139). In addition, to comprehend the overview of each ear weight, we The qualified first ear yield trends slightly declined for OF (p = 0.3884), but significantly decreased for IF (R square = 0.3094, p = 0.0002) and CF (R square = 0.4978, p < 0.0001). Comparatively, the second ear qualified yield trends notably rose in OF (R square = 0.2051, p = 0.0033), but significantly declined in CF (R square = 0.2284, p = 0.0018), and did not have clear trends in IF (R square = 0.0113, p = 0.5139). In addition, to comprehend the overview of each ear weight, we calculated the qualified yield into average ear weight for the first and second ears. Although, the trends of average ear weight from first ear and second ear had differentiable decreasing patterns among three different fertilization treatments, which was significantly decreasing in CF and IF (p < 0.001), whereas not in OF (p = 0.4884). That is, the distribution patterns of average first and second ear weight of OF were more stable than CF and IF. Therefore, we conclude that the decrease in marketable yield was caused by the drop in first and second ear weights for CF and IF, but the yield increase for OF mainly resulted from the second ear.
The results indicate that the application of chemical fertilizer (CF) was unable to maintain the sweet corn yield in the period of 2009-2018, especially for second ear yield (CF average first ear weight trend: y = −3.4461x + 321.08; average second ear weight trend: y = −7.1037x + 251.24). These results also imply that the three fertilization treatments were not adequate for the second fruit growth on CF and IF under this circumstance, but organic fertilizer may enhance soil productivity, over and above nutrient content effects, when large inputs are applied over many years [29].

Correlation Analysis among Climate Variables, Soil Properties, and Crop Yield
To summarize, the crop yield distribution patterns were different among the three treatments, especially for second ear yield of OF. Analysis focused on the decreasing yield of CF and increasing yield of OF second ear. Unexpectedly, conspicuous correlations were found between average ear weight and both EC and relative humidity for the CF and IF treatments, but not for OF, via Pearson correlation analysis (Table 6).  To confirm whether the yield remarkably correlated with these parameters, multiple regression analysis was performed (Table 7). Different average ear weight regression equations were constructed for the three treatments. In the equation for average ear weight from CF, the EC and relative humidity were significant predictors of average first ear weight, while relative humidity and Tmin were predictors of average second ear weight. The ear weight equation model for CF had p < 0.0001 and adjusted R square > 0.5722. In summary, more than 57.2% of sweet corn ear weight could be predicted under the CF treatment. This result confirms that the yield was broadly affected by relative humidity, particularly under the CF treatment. Moreover, the average first ear weight of IF was also influenced by relative humidity, but there was no extraordinary predictor in the yield (weight) regression equation under the OF treatment (adjusted R square < 0.4265). Consequently, we conclude that the increasing relative humidity had negative impacts on the sweet corn yield of CF first and second ear weight, as the same amount of chemical fertilizer was used annually (Table 1).   The Family Chitinophagaceae (Nos. 1, 2, 6, 10, 12) (facultative aerobic bacteria, some active even at very low pH) co-occurred in negative correlation with ear weight as an obvious indicator, whereas the Classes of Acidobacteria (Nos. 16, 20, 21, 26) (aerobic bacteria, low pH) and Betaproteobacteria (Nos. 25,29) (in groundwater ecosystems) were the two differentiable dominant microbes in positive correlation with ear weight. This result suggests that microbial groups influence the decrease in ear weight in CF because microorganisms are affected by chemical components at high RH and soil moisture. Betaproteobacteria, one of the Proteobacteria subgroups, was approved as abundant specific microbes on organic fertilizer treatment soil compared to chemical fertilizer soil [15]. A previous study showed that a 54% yield increase of maize was conferred by microorganisms in no-till strip field [30]. Above all, the Pearson correlation and multiple regression analyses both illustrated the same conclusion, that relative humidity was the key factor determining the yield under the CF and IF treatments, especially under the CF treatment.
Furthermore, to confirm the effect of these factors on the yield, regression tree analysis was performed. Interestingly, the outcome also shows the major importance of relative humidity: 70% of the yield is produced under the relative humidity above 69% and the average yield is far lower than the average yield under relative humidity below 69% (30% of the yield). This is the case for both average weights of the first ear and second ear, under the CF treatment ( Figure 6). Therefore, these results are similar to our previous ones, i.e., increase of relative humidity had a negative impact on the average ear weight of CF. Regression tree analysis also shows that, with the relative humidity below 69% and the soil Cu below 4.6 mg kg −1 , the average yield for the first ear reaches high average weight of 331 g ( Figure 6A). With the relative humidity below 69% and the soil Na above 29 mg kg −1 , the average yield for the second ear reaches high average weight of 266 g ( Figure  6B). This indicates that there are interaction effects of the climate parameters and soil properties with the average yields under CF. The maize yield is more sensitive to climate change under long-term chemical fertilizer treatment, because Southeast Asia and Taiwan have a tropical/subtropical monsoon climate with frequent typhoons and storms. CF and IF treatments gave higher yield values in 2009-2012 with climate anomalies (2013 with heavy precipitations), but the OF treatment was more sustainable in the 2009-2018 period. Through the regression tree analysis, we can find the most important variables of the climate parameters and soil properties which will promote or reduce the sweet corn average yield. Furthermore, the critical interactions of the climate parameters and soil properties to sweet corn average yield are also identified.
As the environment changes, various factors are widely known to affect the maize yield. One study determined that warming temperatures were positively correlated with maize yield at Zhengzhou [9]. Recently, the CERES-Maize model was used to determine optimum planting dates [31]. Lana et al. compared the yields of four maize varieties to choose the best one for coping with the impacts of climate change. They found that change simulation in increasing precipitation was beneficial to the MPA01 cultivar with higher yield stability [32]. This result implied that precipitation might be the limiting factor affecting yield stability. Interestingly, relative humidity had a close relationship with precipitation (Pearson correlation coefficient: 0.40917) in this study. We also noticed the lowest yield of sweet corn in CF compared with the highest yield in OF (on average, CF: 10.9 t ha −1 ; IF: 12.9 t ha −1 ; OF: 13.5 t ha −1 ), as the heavy precipitation occurred in 2013 ( Figure 1B). Nevertheless, the effects of relative humidity on the yield of maize or sweet corn are rarely reported. One study indicated that high relative humidity facilitated reduced iron deficiency chlorosis symptoms to obtain higher plant dry weight and increasing plant height, when comparing lower (60%) and high (90%) relative humidity effects on iron deficiency chlorosis in soybean [33]. Meanwhile, if adequate nutrition is supplied, high relative humidity also increases yield on peanut [34]. Thus, increase of relative humidity seems beneficial to crops when adequate nutrition is applied. However, we observed decreasing yield for CF due to the decline of average weight of sweet corn, especially the second ear weight under long-term application of fertilizer. This might have resulted from the lack of nutrition supplied under the CF treatment or was due to different varieties of crop.
Alternatively, vegetable crops require an adequate and continuous supply of N for proper growth and maximum high-quality yields. As the crop grows, N is the most important of all the essential nutrients for plants. In this study, the accumulated total N in leaves during the harvest stage Therefore, these results are similar to our previous ones, i.e., increase of relative humidity had a negative impact on the average ear weight of CF. Regression tree analysis also shows that, with the relative humidity below 69% and the soil Cu below 4.6 mg kg −1 , the average yield for the first ear reaches high average weight of 331 g ( Figure 6A). With the relative humidity below 69% and the soil Na above 29 mg kg −1 , the average yield for the second ear reaches high average weight of 266 g ( Figure 6B). This indicates that there are interaction effects of the climate parameters and soil properties with the average yields under CF. The maize yield is more sensitive to climate change under long-term chemical fertilizer treatment, because Southeast Asia and Taiwan have a tropical/subtropical monsoon climate with frequent typhoons and storms. CF and IF treatments gave higher yield values in 2009-2012 with climate anomalies (2013 with heavy precipitations), but the OF treatment was more sustainable in the 2009-2018 period. Through the regression tree analysis, we can find the most important variables of the climate parameters and soil properties which will promote or reduce the sweet corn average yield. Furthermore, the critical interactions of the climate parameters and soil properties to sweet corn average yield are also identified.
As the environment changes, various factors are widely known to affect the maize yield. One study determined that warming temperatures were positively correlated with maize yield at Zhengzhou [9]. Recently, the CERES-Maize model was used to determine optimum planting dates [31]. Lana et al. compared the yields of four maize varieties to choose the best one for coping with the impacts of climate change. They found that change simulation in increasing precipitation was beneficial to the MPA01 cultivar with higher yield stability [32]. This result implied that precipitation might be the limiting factor affecting yield stability. Interestingly, relative humidity had a close relationship with precipitation (Pearson correlation coefficient: 0.40917) in this study. We also noticed the lowest yield of sweet corn in CF compared with the highest yield in OF (on average, CF: 10.9 t ha −1 ; IF: 12.9 t ha −1 ; OF: 13.5 t ha −1 ), as the heavy precipitation occurred in 2013 ( Figure 1B). Nevertheless, the effects of relative humidity on the yield of maize or sweet corn are rarely reported. One study indicated that high relative humidity facilitated reduced iron deficiency chlorosis symptoms to obtain higher plant dry weight and increasing plant height, when comparing lower (60%) and high (90%) relative humidity effects on iron deficiency chlorosis in soybean [33]. Meanwhile, if adequate nutrition is supplied, high relative humidity also increases yield on peanut [34]. Thus, increase of relative humidity seems beneficial to crops when adequate nutrition is applied. However, we observed decreasing yield for CF due to the decline of average weight of sweet corn, especially the second ear weight under long-term application of fertilizer. This might have resulted from the lack of nutrition supplied under the CF treatment or was due to different varieties of crop.
Alternatively, vegetable crops require an adequate and continuous supply of N for proper growth and maximum high-quality yields. As the crop grows, N is the most important of all the essential nutrients for plants. In this study, the accumulated total N in leaves during the harvest stage was measured as 3.03, 2.81, and 2.56 ppm on average for CF, IF, and OF, respectively. The leaf N accumulation was significantly different (Table 9) and negatively correlated with the average ear weight (Table 6). Namely, the higher the accumulated N in leaves, the lower the ear weight. This result corresponds to a previous study, which showed that higher relative humidity produces heavier leaf fresh weight and larger leaf area [34]. Another study showed that excising the uppermost two leaves promoted N remobilization from vegetative organs to kernels on maize [35].
Here, we showed that although the trend of N content in the leaves was rising under both CF and OF, N accumulation was higher in CF than in OF on average ( Table 9). The immediate N supply might not be enough or as available to sweet corn in OF as in CF. In addition, due to the slow release of N components from organic matter, a decline in crop yield was observed in the initial years compared with chemical fertilizer [36], but N affects the uptake of other nutrients, which are constituents of organic compounds [37]. As a result, we suggest that after the long-term application of the organic fertilizer, the slow release of N and other essential elements may provide a lasting nutrition supply to the sweet corn under increasing relative humidity, compared to an instant supply using CF. These findings further support the ideas that organic farming may result in greater spatial stability of soil biotic and abiotic properties [38]. Indeed, although the yield of sweet corn from OF was not notably higher than that from CF, the yield trend had higher resilience under OF than CF under increasing relative humidity. Therefore, organic fertilization is certainly necessary in sustainable agriculture.

Conclusions
Organic farming is generally considered more sustainable but less productive than conventional farming. Here, we observed that the yields of sweet corn under organic fertilizer had the potential to surpass those under chemical fertilizer treatment, as the same amount of fertilizer was applied during the period of 2009-2018. From these long-term observations and analyses of the correlations among the yield of sweet corn and environmental factors, this is the first study in which increase of relative humidity is verified as the key factor that impacts the production of sweet corn under chemical fertilizer-treated soil. These consequences suggest that increasing relative humidity may reduce the soil EC, which is a cause of soil degradation, lowering the soil fertility results in decreasing yield under limited nutrition conditions via correlation and regression analysis. These results also imply that higher relative humidity without sufficient nutrition would result in reduced yield. In addition, as the trend of increasing relative humidity continues in company with uneven precipitation, the yield stability of sweet corn with chemical fertilizer was obviously affected compared with organic fertilizer. However, long-term application of organic fertilizer could accumulate more essential chemicals and harbor more specific microbes under soil, which may play considerable roles in helping sweet corn mitigate climate variability.