Transparent and Black Film Mulching Improve Photosynthesis and Yield of Summer Maize in North China Plain

Abstract

In order to clarify the influences of drip irrigation under different mulch materials on crop yield, field experiments were carried out in the North China Plain for two seasons in 2020 and 2021. The changes in field microenvironment, photosynthetic capacity, leaf biological factors, and maize growth indexes were analyzed under drip irrigation with transparent film (W), black film (B), and straw mulching (S), with a nonmulching field as control (CK). The results showed that compared with CK, the yield of W and B increased by 7.2–9.9% and 7.1–12.4%, and the yield of S did not change significantly. The increase in yield was related to the improvement of the field microenvironment and photosynthetic capacity and higher LAI. Compared with CK, the soil water content 0–40 cm below the soil surface of W, B, and S increased by 13.6%, 9.1%, and 4.6%, respectively, and the 5 cm effective accumulated soil temperature of W and B increased by 7.9–10.2% and 4.1–4.7%, respectively. The maximum carboxylation rate (V_max) of W, B, and S at the jointing stage was significantly increased by 3.5–17.3%, 12.7–17.6%, and 10.1–12.7% compared with CK. There was a significant linear correlation between V_max and N_mass, and the correlation was affected by mulching treatments. At the jointing stage, compared with the CK, the LAI of W and B significantly increased by 8.6–66.5% and 7.2–56.0%, but there was no significant difference between S and CK. In conclusion, the increase in yield of W and B resulted from the combined effect of increasing LAI, V_max, and soil water content and temperature.

1. Introduction

By 2050, food production needs to double to meet the growing demand for food [1]. Climate change intensifies soil drought and water shortage, which pose a great threat to crop yield [2,3]. In this context, the sustainable way to achieve food security is to improve crop yield and water use efficiency. The combination of agronomic measures, such as film or straw mulching, and efficient irrigation systems, such as drip irrigation, is conducive to reducing soil evaporation, improving water use efficiency, and achieving the purpose of saving water and increasing production [4,5,6].

The research results on crop growth and water use efficiency under different mulching and irrigation methods have been different. Transparent film, black film, and straw mulching are commonly applied in croplands to improve soil water conditions. Although each single one of the mulching methods has been reported applied in croplands, comparisons of transparent film, black film, and straw mulching effects on crops in one field have been relatively rare. Results from different studies may greatly vary because of the varied soil and air conditions of different research areas. First, the effects of transparent film and black film on crop growth and yield vary with the region and application conditions. Some researchers found that both transparent film and black film could improve yield, but the high temperature of transparent film would lead to tassel abortion of maize, which made the yield improvement effect of the transparent film not as good as that of the black film [7,8]. However, Xu, et al. [9] found that plants with black film grew higher and thinner than those with transparent film, resulting in lodging occurring at the booting stage under black film mulching with drip irrigation of spring maize in North China, which made the effect on yield improvement of the black film less apparent than that of the transparent film. Second, the impact of straw mulching on crop yield also varies with crop types and research areas. Some studies have shown that straw mulching increased soil moisture, reduced water consumption, and improved crop yield and water use efficiency [10]. In other studies, crop yield decreased because of the low soil temperature caused by straw mulching [11,12]. Therefore, the effects of black film, transparent film, and straw mulching on crop yield are inconsistent, and the influence mechanism remains unclear.

Soil water and temperature conditions influence the water status and physiological activity of crops [13]. After surface mulching, the changes in soil hydrothermal conditions cause changes in plant growth and leaf physiological factors, which affect the yield. The photosynthetic rate is positively correlated with temperature within an appropriate temperature range. When the temperature exceeds a certain threshold, the activities of key metabolic enzymes in leaves decrease, causing a decrease in the photosynthetic rate [14,15]. In addition, better physiological characteristics of maize leaves usually correspond to higher grain yield [16]. The improvement of photosynthetic parameters is crucial for improving crop yield and water use efficiency [17,18]. Leaf photosynthetic capacity is closely related to leaf biological factors. Photosynthetic capacity per leaf area is linearly correlated with leaf nitrogen (N) content per leaf area (N_area) and leaf N content per leaf mass (N_mass) [19,20,21]. Leaf N content determines the photosynthetic capacity of leaves, and research results can provide a reference for yield estimation [19,22]. Under drought stress, photosynthesis is restricted by stomata and the biochemical activities of photosynthetic enzymes [23,24]. The biochemical limitation of C3 plants is mainly reflected in the reduction in the maximum carboxylation rate of Rubisco enzyme (V_cmax) [23], while that of C4 plants is reflected in the maximum carboxylation rate of PEP carboxylase (V_pmax) and other remaining enzymes (including Rubisco enzyme, V_max) [25]. It was found that V_cmax of rice was significantly correlated with N_mass and chlorophyll content expressed with SPAD readings, but not with leaf mass per area (LMA) [26]. There was a significant linear relationship between N_mass and V_max of maize leaves [27]. However, it is still unclear whether field surface mulching affects leaf N content and photosynthetic capacity. Studies with physiological parameters of plants usually give descriptive results, but studies with quantified mathematical relationships are fewer. The leaf photosynthetic characteristic parameters of crops grown under different filed mulching should be quantified to uncover the internal causes of yield differences caused by mulching and irrigation modes. This has important scientific significance and practical value for understanding the photosynthetic physiological response of maize under different soil hydrothermal conditions.

In this study, summer maize in the North China Plain was selected, and three mulching methods were set up. Field measurements were carried out for two consecutive years, from 2020 to 2021. The field soil hydrothermal parameters, crop growth, yield, and leaf gas exchange parameters were determined to achieve the following three goals: (1) to examine the temporal and spatial variations in the field microenvironment, such as soil moisture, soil temperature, and leaf gas exchange parameters, under different treatments; (2) to determine the effects of different mulching treatments on crop growth, components of yield, and yield; (3) to explain the internal mechanisms of different mulching treatments affecting leaf photosynthetic capacity via the changed crop growth environment and consequently influencing crop growth and yield. We hypothesized that straw mulching, black film, and transparent film mulching affected plant growth and yield from mildly to severely because of the differences in their improvement effects on soil moisture and temperature.

2. Materials and Methods

2.1. Research Site and Experimental Design

The experiment was carried out in the field of the Agricultural Water Conservation Irrigation Experiment Station in the Daxing District of Beijing, China (39°39′ N, 116°15′ E). The local climate belongs to a typical semiarid continental monsoon climate, with an average annual rainfall of 540 mm. The soil texture of 0–00 cm soil layer in the experimental plot was loam with average field capacity and soil bulk density of 36.58% and 1.41 g/cm³, respectively. The field experiment was carried out from June to October every year from 2020 to 2021. The summer maize was sewn by machine, with wide row spacing of 80 cm, narrow row spacing of 40 cm, plant spacing of 20 cm, and planting density of 83,300 plants/ha. In 2020, the sowing and harvesting times were 30 June and 17 October, while those for 2021 were 26 June and 7 October.

The field was surface drip irrigated and treated with transparent film (W), black film (B), and straw (S) mulching compared with no mulching (CK) under drip irrigation. Three replicates of each treatment were set, with each area being 10.8 m × 4 m and arranged randomly. The film was polyethylene transparent film and black film with a thickness of 0.01 mm. All the straws of previous winter wheat were harvested and removed from the field. Part of the straw was crushed and covered over the driplines in the S plots. The coverage amount was 6000 kg ha⁻¹. Straw mulching was covered evenly across the field, while 80 cm-width films was covered two lines of maize and the narrow rows, leaving the wide row bared. The transparent film, black film, and straw were covered on the drip irrigation belt. The driplines were laid in the middle of narrow rows with a ’one dripline, two rows‘ layout mode.

Rain-fed and supplemental drip irrigations were applied as needed to remove the water stress of the crops. In 2020, because of the rainfall after sowing, the seedling water was not irrigated. It was irrigated three times, 5 mm every time, during topdressing. In 2021, there were three times of irrigation. The first time was 50 mm of seedling water, and the last two times were 5 mm each time of topdressing irrigation. The total rainfall in 2020 and 2021 was 309.63 mm and 341.88 mm, respectively. Rainstorms are usually heavy and frequently happen in the summer maize sowing season in the research area, leading to fertilizer leaching. In such cases, no base fertilizer but only topdressing was used. During the growth period of maize, the application rates of nitrogen, phosphorus (P₂O₅), and potassium (K₂O) were 240, 135, and 135 kg/ha, respectively. The specific water and fertilizer schedules are shown in

Table 1. Irrigation and fertilization schedules for 2020 and 2021.

Years	Times	Growth Stages	Date	Irrigation (mm)	Nutrients (kg/ha)
					N	P2O5	K2O
2020	1	Jointing stage	8/1	5	156	108	135
	2	Tasseling period	8/16	5	48	27
	3	Filling stage	9/13	5	36
	Total			15	240	135	135
2021	1	Seedling stage	6/27	50
	2	Jointing stage	8/2	5	156	108	135
	3	Tasseling period	8/22	5	84	27
	Total			60	240	135	135

.

2.2. Soil Water Content and Soil Temperature

The soil water content and temperature were monitored by EM50 soil moisture and temperature sensors and data collectors (Decagon devices Inc., Pullman, WA, USA) and collected every 30 min. Three depths of 5 cm, 15 cm, and 30 cm below the soil surface were selected at which to install the 5TM sensors for the two years. The sensor measurements were carried out two repetitions in each treatment. Data of one system were collected by a datalogger automatically as previously described, and those of the other system were collected manually with a precheck device (Decagon devices Inc., Pullman, WA, USA) at 18:00 every two days. The soil water content of the profile was also measured by the oven-drying method. The soil sampling depths were 5 cm, 15 cm, and 30 cm. The soil sampling times were during the growth stage to cover the different soil water content ranges. The times were 19 July, 18 August, 4 September, and 9 October in 2020 and 7 July, 11 September, and 8 October in 2021. By fitting the soil water contents from the oven-drying method and those from the sensors on the same day, the soil water content of the probe was corrected.

The effective accumulated temperature of soil is referred to the method of Zhang et al. [28]. The lower limit temperature (biological zero) for the growth and development of summer maize in North China is 10 °C, and the effective accumulated temperature is the sum of the daily average ground temperature exceeding biological zero during the growth period of maize.

2.3. Air Temperature and Relative Humidity

The air temperature and relative humidity between rows were monitored by automatic-datalogging sensor (Hobo Pro V2, ONSET, Bourne, MA, USA), which was installed 20 cm away from the ground surface for all treatments. The data acquisition interval was once every 30 min.

2.4. Photosynthetic Physiological Parameters

At the jointing, filling, and maturity stages, three sun leaves were randomly selected from each plot. The response of the net rate of photosynthetic CO₂ assimilation (A, μmol m⁻² s⁻¹) versus the intercellular CO₂ concentrations (A-Ci) was measured by a portable photosynthetic system (Li- 6800, LiCor Inc., Lincoln, NE, USA). The temperature and humidity of the leaf chamber were set to 26 °C and 60%, and the light intensity was set to 2100 μmol m⁻²s⁻¹. Measurements were taken at CO₂ concentrations of 400, 300, 200, 150, 100, 80, 40, 400, 400, 600, 800, 1000, and 1200 μmol mol⁻¹. According to Von Caemmerer [25], the A to Ci relationship at a Ci < 50 μmol mol⁻¹ was used to solve the maximum carboxylation capacity of phosphoenolpyruvate carboxylase (V_pmax, μmol m⁻² s⁻¹), and the CO₂ saturated photosynthetic capacity (V_max, μmol m⁻² s⁻¹) was the horizontal asymptote of the A/Ci curve. Next, the relative value of chlorophyll was measured by a SPAD instrument, and the measurement position was consistent with that of photosynthesis. Five points were selected in the middle of the leaf for reading. Then, avoiding the main vein, ten holes were punched in the leaf with a 20 mm diameter punch, and the leaves were weighed, and the leaf mass per area (LMA) measured, after drying. The main veins, tip, and base of the remaining leaves were removed from the middle part of the leaves. After drying and grinding, the leaves were digested with sulfuric acid, and the N content per leaf mass (N_mass) was determined with an automatic Kjeldahl apparatus (Kjeltec 8400; FOSS, Denmark).

2.5. Leaf Area Index and Yield

Leaf area index (LAI) was measured three to four times during the growing season. The length and width in the widest part of each leaf were measured using a tape with a precision of 1 cm, while the leaf area was obtained by multiplying the empirical coefficient of 0.74. LAI was calculated as the sampled leaf area divided by the product of the plant spacing (20 cm) and sampled length (60 cm). At the harvest, four subplots with areas of 3 m × 1.2 m were selected from each plot, and the ear length (cm), ear diameter (mm), kernels per ear (ear⁻¹), and the hundred grain mass (g) were determined and converted to yield of the plot (Y, t/ha).

2.6. Statistical Analysis

The soil water content and soil temperature among CK, W, B, and S treatments were compared via a paired sample T test in SPSS (v13.0, IBM SPSS Statistics, Chicago, IL, USA). The mean differences in V_max, V_pmax, LMA, N_mass, SPAD, and yield among different treatments were determined by one-way ANOVA. The data of V_max, V_pmax, LMA, N_mass, and SPAD in two years were analyzed using Origin 2021 and tested by the general linear model. If there was a significant difference in the slope and intercept of the regression line among V_max, V_pmax, and N_mass, four lines were fitted: CK, W, B, and S; otherwise, a single line was fitted for all treatments. All figures were drawn using SigmaPlot (v13.0, Systat Software Inc., San Jose, CA, USA).

3. Results

3.1. Field Microenvironment

The average value of 0–40 cm soil water content (θ) in the growing season of the three mulching treatments in 2020 was 0.23–0.25 cm³ cm⁻³ (Figure 1). Compared with CK, the θ of W, B, and S increased significantly by 13.6%, 9.1%, and 4.6% (p < 0.001), respectively. In 2021, the θ varied from 0.28 to 0.30 cm³ cm⁻³ among all the treatments. Compared with CK, θ of W, B, and S increased significantly by 7.4%, 11.1%, and 3.7% (p < 0.001), respectively. In 2020, there were seven days with θ below 60% field capacity (60% FC) in both CK and S treatments during the whole growing season. B and W treatments reduced the days with θ below 60% FC to four and zero days, respectively.

After the emergence of seedlings in 2020, the averaged soil temperatures at the depth of 5 cm below ground (T_s-avg-5) of CK, W, B, and S were 22.5 °C, 23.5 °C, 23.0 °C, and 22.3 °C, respectively (Figure 1). The 5 cm soil effective accumulated temperature (GDD₅) of W and B were 7.9% and 4.1% higher than that of CK, respectively, and that of S was 1.5% lower than that of CK. The results of soil temperature in 2021 were similar to those in 2020. The GDD₅ of W and B increased by 10.2% and 4.7%, respectively, and S decreased by 1.0%.

Figure 2 shows the seasonal dynamics of the differences in the daily maximum soil temperature at 5 cm below ground (T_s-max-5) between the three mulching treatments and CK. During the growth period, even with several exceptions in 2020, during which the T_s-max-5 difference between W and CK was negative because of the rainy weather or irrigation, W significantly increased the T_s-max-5 (p < 0.001), and S significantly reduced the T_s-max-5 (p < 0.001), for two years. The T_s-max-5 of B increased significantly in 2020 (p < 0.001), but it did not increase significantly in 2021 (Figure 2). W and B significantly increased T_s-avg-5 and minimum value of soil temperature at 5 cm below ground (T_s-min-5, p < 0.001). The T_s-avg-5 of S increased significantly in 2020 (p < 0.001), but it did not increase significantly in 2021, and S significantly increased the T_s-min-5 (p < 0.001) for two years (Figure S1).

The daily changes in air temperature and relative humidity at 20 cm on the surface showed the opposite trend. There was no significant difference in air temperature and relative humidity among all treatments (Figure S2).

3.2. Leaf Area Index

The effect of mulching on LAI occurred mainly before the jointing stage (Figure 3). Compared with CK, the LAI of W and B increased significantly by 8.6–40.5% and 7.2–38.7% (p < 0.001), respectively, during the seedling and jointing stages in 2020, but there was no significant difference between S and CK during the same stages (p > 0.1). There was no significant difference in the LAI among W, B, S, and CK after the jointing stage (p > 0.1). The changes in 2021 were similar to those in 2020. During the seedling and jointing stages, the LAI of W and B increased significantly by 44.3–66.5% and 51.9–56.0% (p < 0.001), respectively, but there was no significant difference among treatments (p > 0.1) after the jointing stage.

3.3. Yield and Water Use Efficiency

Compared with CK, the kernels per ear of W and B treatments increased by 4.2–12.5% and 5.0–8.5%; the hundred kernel mass of W and B treatment increased by 2.3–2.8% and 1.9–3.9%; and the yield of W and B increased significantly by 7.2–9.9% and 7.1–12.4%, respectively. Compared with W and B, the hundred kernel mass of S decreased by 5.1–8.4% and 6.8–8.3%, respectively. No significant difference was observed between the yield of S and CK (

Table 2. Yield and its components among treatments for the two years. Notes: Data are given as means with standard deviations in brackets.

Years	Treatments	Ear Length (cm)	Ear Diameter (mm)	Kernels/ear	Hundred Kernel Mass (g)	Yield (t/ha)
2020	CK	19.43	54.29a	481.29b	33.32ab	12.58b
	W	19.25	54.38a	501.63a	34.24a	13.48a
	B	19.25	54.13a	505.50a	33.95a	13.47a
	S	19.25	53.5b	489.38ab	32.39b	12.44b
	p value	0.903	0.064	0.073	0.052	0.027
2021	CK	17.63	49.63	410.75b	27.35bc	10.37b
	W	18.00	49.5	462.13a	27.97ab	11.40a
	B	17.75	49.75	445.50a	28.43a	11.66a
	S	17.50	49.38	409.50b	26.62c	10.16b
	p value	0.455	0.905	<0.001	0.007	0.001

Notes: Different lower-case letters (a, b, and c) in the table indicate statistically significant differences at 0.05 < p < 0.1.

).

Principal component analysis of yield components in 2020 and 2021 showed that the first two principal components in 2020 explained 80% of the total variation (Figure 4). PC1 explained 58.9% of the total variation and positively correlated with yield, kernels per ear, and ear diameter. In addition, PC2 explained 21.1% of the total variation. PC2 had a positive correlation with hundred kernel mass and a negative correlation with kernels per ear. In 2021, the first two principal components explained 79.9% of the total variation. PC1 explained 63.7% of the total variation and showed a positive correlation with yield, kernels per ear, ear diameter, and hundred kernel mass. PC2 explained 16.2% of the total variation and was positively correlated with ear length. The results showed that kernels per ear and ear diameter had important effects on maize yield.

3.4. Photosynthetic Capacity

In 2020 and 2021, the V_max was highest at the jointing stage and then gradually decreased (

Table 3. Changes in photosynthetic capacity among treatments.

Years	Parameters	Treatments	Jointing Stage	Filling Stage	Mature Stage
2020	Vmax	CK	42.7 ± 2.4a	40.9 ± 4.7	30.2 ± 5.7
	(µmol m⁻2 s⁻¹)	W	44.2 ± 4.7ab	42.3 ± 3.6	33.5 ± 3.1
		B	50.2 ± 1.8c	39.4 ± 3.5	36.6 ± 4.4
		S	47 ± 3.9bc	40.5 ± 4.5	35.1 ± 5.3
		p value	0.005	0.723	0.186
	Vpmax	CK	123.6 ± 19.6	115.4 ± 16.7	104.3 ± 16.1
	(µmol m⁻2 s⁻¹)	W	120.3 ± 25.4	121.7 ± 20.9	96.2 ± 18.6
		B	142.8 ± 15.4	120.8 ± 18.2	103.8 ± 14.4
		S	134.1 ± 24.1	116.2 ± 8.5	116.3 ± 11.2
		p value	0.281	0.905	0.328
2021	Vmax	CK	41.7 ± 6.3a	39.4 ± 1.8	27.9 ± 1.8a
	(µmol m⁻2 s⁻¹)	W	48.9 ± 3.0b	36.4 ± 4.7	31.5 ± 3.3bc
		B	47.0 ± 0.9b	38.6 ± 3.6	29.6 ± 3.0ab
		S	47.0 ± 4.1b	38.0 ± 5.5	33.1 ± 4.3c
		p value	0.066	0.685	0.033
	Vpmax	CK	173.1 ± 35.3	126.3 ± 3.4	98.4 ± 8.7
	(µmol m⁻2 s⁻¹)	W	222.9 ± 45.5	122.5 ± 20.4	92.0 ± 15.7
		B	191.5 ± 41.3	113.4 ± 15.0	94.7 ± 22.7
		S	195.3 ± 35.7	111.9 ± 19.3	109.0 ± 25.2
		p value	0.296	0.429	0.386

Notes: Different lower-case letters (a, b, and c) in the table indicate statistically significant differences at 0.05 < p < 0.1.

). At the jointing stage in 2020, the V_max of W, B, and S was significantly higher than that of CK by 3.5%, 17.6%, and 10.1% (p = 0.005), respectively. In 2021, compared with CK, the V_max of W, B, and S treatments significantly increased by 17.3%, 12.7%, and 12.7%, respectively, at the jointing stage (p = 0.066), and W and S significantly increased by 12.9% and 18.6%, respectively, at the maturity stage (p = 0.033). No significant difference in V_max among treatments was found in other measurement periods. There was no significant difference in V_pmax among the jointing, filling, and mature stages in 2020 and 2021 (p > 0.1).

3.5. Biological Factors of Leaves

The characteristics of biological factors of leaves affect the photosynthesis and yield formation of maize.

Table 4. Changes in biological factors of leaves among treatments.

Years	Parameters	Treatments	Jointing Stage	Filling Stage	Mature Stage
2020	LMA	CK	21.2 ± 3.29	33.8 ± 1.12	33.7 ± 0.8
	(g/m2)	W	23.3 ± 5.57	32.7 ± 0.94	35.5 ± 2.29
		B	24.8 ± 2.66	34.4 ± 1.72	35.1 ± 2.63
		S	21.0 ± 5.78	34.5 ± 2.31	35.0 ± 2.31
		p value	0.49	0.29	0.48
	Nmass	CK	37.3 ± 1.5ab	34.6 ± 3.0	30.9 ± 1.6
	(g/kg)	W	38.2 ± 0.9a	35.3 ± 1.9	31.3 ± 0.8
		B	37.5 ± 0.2ab	33.9 ±1.1	30.6 ± 1.2
		S	36.6 ± 0.6b	33.6 ± 2.8	30.6 ± 2.6
		p value	0.09	0.6	0.86
	SPAD	CK	—	54.4 ± 2.72	53.9 ± 3.19ab
		W	—	59.2 ± 3.28	57.4 ± 2.94c
		B	—	57.9 ± 3.58	51.1 ± 2.93a
		S	—	57.9 ± 2.73	55.0 ± 3.01bc
		p value		0.13	0.02
2021	LMA	CK	25.7 ± 1.86a	30.9 ± 1.32	29.4 ± 1.29a
	(g/m2)	W	28.2 ± 1.32b	30.4 ± 1.41	29.1 ± 1.62a
		B	26.1 ± 1.32ab	32.0 ± 1.01	30.1 ± 1.41ab
		S	27.0 ± 1.57b	30.7 ± 1.68	31.1 ± 1.77b
		p value	0.02	0.27	0.15
	Nmass	CK	36.7 ± 1.9	34.6 ± 0.8ab	29.6 ± 2.1
	(g/kg)	W	39.6 ± 4.1	35.7 ± 1.8b	30.5 ± 1.1
		B	37.3 ± 1.1	33.8 ± 1.4a	29.7 ± 1.7
		S	37.1 ± 0.9	34.0 ± 0.7a	29.5 ± 1.0
		p value	0.16	0.08	0.7
	SPAD	CK	52.6 ± 3.46	58.0 ± 1.30	49.2 ± 2.71a
		W	54.8 ± 1.21	58.4 ± 3.39	53.0 ± 2.01b
		B	53.2 ± 1.77	56.0 ± 2.35	49.7 ± 4.94ab
		S	54.7 ± 2.65	57.5 ± 2.48	53.0 ± 2.11b
		p value	0.21	0.63	0.09

Notes: Different lower-case letters (a, b, and c) in the table indicate statistically significant differences at 0.05 < p < 0.1.

showed the LMA, N_mass, and SPAD of all treatments at different growth stages. In 2021, at the jointing stage, W and S significantly increased LMA by 9.7% and 5.1% (p = 0.02), respectively, compared with CK. At the mature stage, the LMA of S was significantly higher than CK.

The N_mass was the highest at the jointing stage and gradually decreased with the advancement of the growth stage. In 2020, the N_mass of W was the highest in the whole growth period, and the N_mass of W at the jointing stage was significantly higher than that of S by 4.4% (p = 0.09). There was no significant difference in N_mass among other treatments. In 2021, the N_mass of W was significantly higher than that of B and S by 5.6% and 5.0% (p = 0.08), respectively, but no significant difference was found bet3ween the three treatments and CK.

SPAD did not show evident seasonal trends in the two seasons, and the SPAD of W was the highest. Compared with CK, the SPAD of W at the mature stage in 2020 increased significantly by 6.5% (p = 0.02), and the SPAD of the W and S at the mature stage in 2021 both increased significantly by 7.7% (p = 0.09). No significant difference in SPAD was found among treatments in other determination periods.

V_max was significantly correlated with V_pmax, N_mass, SPAD (Figure 5). In 2020, the correlation coefficients between V_max and V_pmax, LMA, N_mass, SPAD were 0.53 (p < 0.001), −0.47 (p < 0.001), 0.63 (p < 0.001), and 0.46 (p = 0.001), respectively; in 2021, the correlation coefficients between V_max and V_pmax, LMA, N_mass, and SPAD were 0.75 (p < 0.001), −0.23 (p = 0.07), 0.65 (p < 0.001), and 0.42 (p < 0.001), respectively.

The correlations between V_pmax and LMA, N_mass, and SPAD were significant. In 2020, the correlation coefficients between V_pmax and LMA, N_mass, and SPAD were −0.33 (p = 0.005), 0.37 (p = 0.002), and 0.29 (p = 0.056), respectively. In 2021, the correlation coefficients between V_pmax and LMA, N_mass, and SPAD were −0.40 (p = 0.001), 0.62 (p < 0.001), and 0.23 (p = 0.063), respectively.

The linear regressions between V_max, V_pmax, and N_mass were further analyzed. A significant linear correlation between V_max and N_mass was found, and the correlation among the treatments was affected by mulching. The intercept of the regression equation of W, B, and S was significantly higher than that of CK (Figure 6a, p < 0.001). There was a significant linear correlation between V_pmax and N_mass, but no significant treatment effects on the regressions (Figure 6b, p > 0.1).

3.6. A Comprehensive Evaluation of the Relationship between Factors Affecting Yield

Six index factors (

Table 5. Index system of yield influencing factors.

Indicator’s Category	Project	Factors	Unit
Field microenvironment	x1	θ	cm3cm−3
	x2	Soil temperature	℃
Growth index	x3	LAI
Physiological index	x4	Maximum rate of carboxylation	µmol m⁻2 s⁻¹
Yield components	x5	Ear diameter	mm
	x6	Kernels per ear

) of field microenvironment, growth index, physiological index, and yield components were selected for principal component analysis. The extraction conditions in the SPSS software environment were that the eigenvalue was greater than one, and the cumulative contribution rate was greater than 85%. Thus, according to the extraction conditions, the total variance decomposition table and the principal component load matrix (Tables S1 and S2) were obtained. According to Table S1, the eigenvalues of the two principal components were greater than one at 4.18 and 1.30. The cumulative variance contribution rate reached 91.27%, basically covering the information of all influencing factors. Therefore, this paper summarized two principal components. The first principal component was related to ear diameter (x5), LAI (x3), soil water content (x1), kernels per ear (x6), and soil temperature (x2). Among them, a strong positive correlation between ear diameter (x5) and LAI (x3) was found, mainly reflecting the impact of nonphysiological indicators such as yield components and growth indicators on yield. The second principal component was mainly positively correlated with the maximum rate of carboxylation (x4), which mainly reflected the positive effect of the physiological index on maize.

Through principal component analysis, we summarized six indexes as the effects of physiological and nonphysiological factors on maize yield, including both positive and negative correlation index factors. To better analyze the correlation degree between factors and crop yield, grey correlation analysis of maize yield was established (

Table 6. Grey correlation analysis of yield and influencing factors.

Factors	LAI	Soil Temperature	θ	Maximum Carboxylation Rate	Ear Diameter	Kernels Per Ear
Correlation coefficient	0.995	0.972	0.972	0.894	0.876	0.36
Correlation rank	1	2	3	4	5	6

), which was ranked as LAI (x3) > soil temperature (x2) > θ (x1) > maximum carboxylation rate (x4) > ear diameter (x5) > kernels per ear (x6).

4. Discussion

4.1. Growth and Yield

Compared with CK, W and B treatments significantly increased kernels per ear and hundred kernel mass, resulting in a significant increase of 7.2–12.4% in yield. At the same time, there was no significant difference in ear length, kernels per ear, hundred kernel mass, or yield between S and CK (Table 2). The changes in yield and components of yield could be contributed by the improvement in the θ, soil temperature, and growth indexes such as LAI. Transparent film mulching increased maize yield by providing appropriate soil temperature and water content, air temperature, and relative humidity [29]. Unlike Liu et al. in [29], we did not find a significant difference in air temperature and relative humidity at 20 cm above the soil surface among treatments (Figure S2). However, the θ and soil temperature were significantly improved in W and B treatments (Figure 1 and Figure 2), which may be the main reason for the better plant growth and higher yield and photosynthetic parameters of those treatments. In 2020 and 2021, W, B, and S increased the 0–40 cm θ by 7.4–13.6%, 9.1–11.1%, and 3.7–4.6%, respectively which was consistent with previous studies that found that transparent film, black film, and straw mulching can increase soil water content [9,30,31]. In addition, although there was sufficient rainfall during the growth period of summer maize in the experimental area, in 2020, because of the uneven distribution of rainfall, the θ of CK was less than 60%FC for seven days. Below 60%FC was generally considered as the beginning of crop drought stress [32]. W and B significantly reduced the number of such days to zero and four. Although S treatment did not shorten the drought stress period, the θ of S also increased significantly, indicating that all the mulching treatments in our study could reduce crop drought stress. The change in soil temperature and the increase in GDD₅ could also promote crop growth. Significant correlations between the growth rate of maize and the accumulated temperature at the depth of 5 cm below ground have been reported previously [33,34]. In this study, compared with CK, W and B increased the temperature and the GDD₅ (Figure 1 and Figure 2), so they probably would improve the growth rate of maize. Although S did not increase the GDD₅, the S treatment could cool down the soil temperature when the temperature was high and warm up it when low. Thus, the S treatment reduced the temperature variation experienced by crops (Figure 2), which may also have a positive effect on crop growth [21].

The increase in plant height, LAI, stem diameter, and dry matter also promoted the increase in yield. In particular, rising LAI is closely related to light interception. Studies have shown that 95% of radiation in crucial growth periods should be intercepted in order to maximize the growth rate and potential yield per plant [35,36]. Leaf area not only affects crop transpiration but also affects the size of sunlight exposure area and the ability of crop photosynthesis [37,38]. From the grey correlation analysis of yield and influencing factors in our study, LAI had a relatively significant impact on maize yield in this study (Table 6), and the increase in LAI under W and B treatments positively affected yield formation (Figure 3).

4.2. Photosynthetic Physiological Parameters

Leaf physiological and ecological indexes such as leaf nitrogen content, chlorophyll content, and LMA can significantly affect leaf photosynthetic capacity and then affect the formation of crop yield. V_pmax and V_max are two critical parameters for characterizing the photosynthetic capacity of maize as C4 plants. The variation range of V_max in this study was 27.9–50.2 µmol m⁻² s⁻¹, and V_pmax varied between 92.0 and 222.9 µmol m⁻² s⁻¹. This was consistent with previous research [39,40]. This study found that although there was no significant difference in V_pmax among different mulching treatments, the V_max of all the mulching treatments significantly increased compared with CK. The increase in V_max was significantly correlated with the N_mass, which was consistent with the results of previous studies [41,42]. Moreover, we found that different mulching treatments significantly affected the correlation between V_max and N_mass (Figure 6). The intercepts of the linear regression equation between V_max and N_mass in W, B, and S treatments were significantly higher than that in CK, suggesting that when N_mass was the same, the V_max of W, B, and S were higher than that of CK. Such results indicated that the mulched maize had higher leaf nitrogen utilization efficiency, which was consistent with the results of previous studies on the photosynthesis of maize under transparent film mulching [43] and straw mulching [21]. The improvement in leaf nitrogen utilization efficiency has been observed several times in maize with mulching treatments, though the internal reasons need to be further studied. In this study, both W and B increased SPAD, which was consistent with previous studies, but S also increased SPAD, which was contrary to a previous study in which S reduced SPAD [12]. The impact of straw mulching on crop growth is complex and affected by the interaction of farming measures, straw mulch amounts, soil texture, and other factors [44,45], which need to be further studied.

5. Conclusions

This study quantified the differences in field microenvironment, crop growth, yield, and photosynthetic physiological parameters among different mulching treatments with drip irrigation and tried to explain the contributions to yield improvement through correlation analysis. Our results partially supported our hypothesis. At least in North China Plain under the condition of full drip irrigation, transparent and black film mulching affected plant growth and yield to the same extent, because both mulching types increased soil water content in the root zone, increased the minimum soil temperature and accumulated temperature in 5 cm soil, and improved photosynthetic area (LAI) and photosynthetic capacity (V_max). However, straw mulching, even with the increase in soil water content, did not increase yield. Through principal component analysis and grey correlation analysis, it was found that LAI had the greatest impact on maize yield. This study also provided the biological factors of maize leaves in this area and clarified the environmental and physiological reasons for yield changes under different mulching treatments.