Next Article in Journal
Ca2+/Na+ Ratio as a Critical Marker for Field Evaluation of Saline-Alkaline Tolerance in Alfalfa (Medicago sativa L.)
Previous Article in Journal
Changes in Biochemistry and Yield in Response to Biostimulants Applied in Bean (Phaseolus vulgaris L.)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling the Interaction of Plastic Film Mulch and Potato Canopy Growth with Soil Heat Transport in a Semiarid Area

1
College of Hydraulic Science and Engineering, Yangzhou University, Yangzhou 225009, China
2
Center for Agricultural Water Research in China, China Agricultural University, No. 17 Qinghua East Road, Haidian, Beijing 100083, China
3
Oregon State University, Malheur Experiment Station, 595 Onion Ave., Ontario, OR 97914, USA
*
Author to whom correspondence should be addressed.
Agronomy 2020, 10(2), 190; https://doi.org/10.3390/agronomy10020190
Submission received: 9 December 2019 / Revised: 18 January 2020 / Accepted: 29 January 2020 / Published: 31 January 2020

Abstract

:
Plastic film mulch is an important agricultural technology to reduce water evaporation and modify the soil thermal conditions for crop production. The optical properties of plastic film mulch and the crop canopy growth are both key factors impacting soil heat transport in the soil-film-canopy-atmosphere ecosystem. In this study, a process-oriented model was developed to better understand the interaction among the plastic film mulch, potato (Solanum tuberosum L.) canopy growth, and soil thermal conditions. Canopy growth, photosynthetically active radiation transmittance, net radiation, soil heat flux, and temperature were monitored in a two-year plastic mulch field experiment in Wuwei (Gansu Province, China). Results showed that the simulation of daily soil surface temperature had a good performance with 2.8 and 1.5 °C of root mean square error (RMSE) for the transparent film mulch (TM) and black film mulch (BM), respectively. Moreover, the simulation of the daily net radiation and soil heat flux model indicated reasonable fluctuations with potato phenological development with the daily R2 ranging from 0.89 to 0.98 in 2014 and 2015 for the TM and BM treatments. It was shown that the canopy temperature under BM was greater than that in TM treatment, and the maximum value difference could be up to 7 °C during the early potato growing period, which implied that the BM may perform better in modifying the canopy thermal condition. The model could provide heat distribution information for plastic film choosing in potato field to avoid heat stress.

1. Introduction

Plastic film mulch has been widely used in agriculture, especially in arid and semiarid regions, since it can save soil water by reducing evaporation [1] and optimize the light and thermal environments for crops [2,3,4,5,6,7,8,9,10,11]. Accordingly, plastic film mulching can improve seed emergence [12,13], root growth [14,15], nutrient uptake [16], stem elongation [5,11], final yield and water use efficiency [17,18,19,20,21].
The two widely used plastic films are black plastic film and transparent plastic film. Although both black plastic film (with high shortwave absorptance) and transparent plastic film (with high shortwave transmittance) can raise soil temperature [22], the effects of these two plastic films on plant growth can be contrary. For example, some researchers report that the use of black plastic film mulch leads to higher crop yield than transparent plastic film mulch [12,23,24]. Others report that transparent plastic film mulch results in higher crop yield than black plastic film mulch [25,26]. Such contradictory results are possible due to the complexity of heat transport in the soil-film-canopy-atmosphere system. Heat transport can be affected by the extent of contact between mulch and soil, the optical properties of the plastic film, and the geometry of the mulched surface [27].
Models have been developed to study heat transport considering different influencing factors since field experiments are time consuming and laborious. For example, some models have been developed to describe heat transport between soil, transparent plastic film mulch, and the atmosphere with energy balance equations considering the impacts of soil characteristics and the optical properties of mulch [28,29,30]. One validated heat transport model considers the effects of the optical properties of different plastic film mulches and the air gap between plastic film mulch and soil surface [6]. Since these models have limitations in simulating heat transport in field crops with plastic film mulch, other models have been developed. For example, Wu et al. [31] proposed a model to simulate soil temperature and water content in winter wheat (Triticum aestivum L.) with transparent plastic film mulch. This model mainly focuses on the distribution of soil temperature and water content. Yang et al. [32] added a plastic film mulch layer model into the SHAW model, constructed by Flerchinger et al. [33] and Flerchinger and Pierson [34], to consider the effects of plastic film mulch on heat transport in maize (Zea mays L.).
Although these models have been developed to study the effects of plastic film mulch on heat transport in fields with or without crops, some problems still need to be solved with new models. Heat stress, affecting the above-ground and below-ground parts growth of potato [35,36], can cause the low potato yield [37]. It is necessary to know the heat distribution for plastic film mulch choosing to provide optimal heat conditions for potato growth. One important problem is how to simulate the heat transport as the potato plant canopy changes. Zhang et al. [11] observed that the effects of transparent and black plastic film mulch on soil temperature in potato fields can be different in two growing seasons because of the different canopy growth. How can the interaction of canopy growth and the optical properties of plastic film mulch on heat transport be considered in the model?
The objective of this study was to develop a model to simulate the radiation, soil heat flux, and mulch surface and soil surface temperature at different potato growth stages considering the interaction of canopy growth and the optical properties of plastic film mulch. Leaf area index, canopy coverage, and photosynthetically active radiation transmittance were introduced in the model to account for the effects of plant canopy growth on heat transport. The model can be used to predict radiation, temperature, and soil heat flux in potato field and explain the effects of different plastic film mulches on potato growth during different growth stages.

2. Materials and Methods

2.1. Study Area and Experimental Design

The study area was located at the Shiyanghe Experimental Station (N 37°52′, E 102°50′, altitude 1581 m), China Agricultural University, Wuwei, Gansu Province, on the border of the Tenger Desert with a typical continental temperate climate. The area has sufficient heat resources and limited water resources with mean annual temperature 8 °C, annual accumulated temperature (>0 °C) 3550 °C, annual sunshine duration 3000 h, annual precipitation 164 mm, pan evaporation 2000 mm, and groundwater table 5–30 m below ground surface [38]. The soil texture is sandy loam with mean soil bulk density 1.53 g/cm3 at 0–1.0 m depth [11].
To calibrate and verify the model, experiments were conducted in potato fields from April to August in 2014 and 2015. Two soil surface treatments were compared: transparent plastic film mulch (TM) and black plastic film mulch (BM), each with a thickness of 0.008 mm. The soil surface of each treatment was fully covered with polyethylene (PE) plastic film. Each treatment was replicated three times in a randomized complete block design.

2.2. Agronomic and Irrigation Practices

As the agronomic and irrigation practices have been described in detail by Zhang et al. [11], only a general introduction is presented here to avoid repetition. In the first year the potatoes were planted on 22 April and harvested on 21 August 2014. Since the weather was warm earlier in 2015, the potatoes were planted on 15 April and harvested on 20 August. The plot size was 6 m long and 5.6 m wide with 7 beds (0.8 m wide and 0.2 m high). The potato crop was irrigated using a drip irrigation system with drip tape placed under the plastic mulch on the center of the beds. Irrigation was initiated when the soil matric potentials measured with tensiometers reached −25 kPa. The specific irrigation information has been shown by Zhang et al. [11].

2.3. Measurements

Ten plants in each plot were labeled for height measurement every 7 days (fewer than 7 days during vigorous plant growth). The leaves of five plants in two plots were collected to measure leaf area index in 2014. Total leaf weight and the weight of 30 typical leaves were measured. To reduce labor intensity only the area of the 30 typical leaves was measured with the AM300 (ADC BioScientific Ltd., Hertfordshire, UK) a portable leaf area meter. The total leaf area per plant was calculated using the ratio between the weight of the 30 leaves and the total leaf weight.
The fractional canopy coverage, proportion of the ground covered by potato canopy (Pcc), was determined with digital images which were taken at a height of 2 m above ground using a Canon IXUS 132 digital camera (Canon Inc., Tokyo, Japan) [39,40]. The images were taken at 12:00–13:00 to minimize shadows. A typical region (1.6 m × 0.9 m) containing two rows (row spacing 0.8 m) with 6 plants (plant spacing 0.3 m) was selected in each plot to take photographs. The plants were captured in the images using quick selection tool, magic wand tool, and lasso tool of Adobe Photoshop CS5 (Adobe Systems Inc., San Jose, CA, USA). Then, the images were processed into black and white binary image (black pixels denote plants) with the Adobe Photoshop CS5. The proportion of black pixels in the image, read from the Adobe Photoshop CS5 histogram represented the fractional canopy coverage. The example of images before and after processing is shown in Figure 1.
The photosynthetically active radiation transmittance of the potato canopy (Tpar) was measured with SunScan Canopy Analysis System (Delta-T Devices Ltd., Cambridge, UK) at 12:00–13:00 on a sunny day. Given the inhomogeneity of the potato canopy, the SunScan probe was placed on the bed (perpendicular to the row direction) at 7 random positions for each plot. The average of the transmittance was recorded.
Temperatures on the plastic film mulch surface, on the soil surface, and at 0.05 and 0.1 m soil depths were measured with sensors in the middle of the beds in one replication of each treatment. The sensors on the plastic film mulch surface, on the soil surface, and at 0.05 m soil depth were 0.25-mm copper-constantan thermocouples (ST10, Beijing Unism Technologies, Inc., Beijing, China). The sensors at 0.1 m soil depth were soil temperature/water sensors (FDS120, Beijing Unism Technologies, Inc.).
Weather variables were measured with a standard automatic weather station (HOBO H21-001, Onset Computer Corp., Cape Cod, MA, USA), 2 m aboveground. The measurable weather variables consisted of global radiation, wind speed, relative humidity, air temperature, atmospheric pressure, and precipitation. Global radiation and air temperature has been shown in previous report of Zhang et al. [11]. Net radiation above the canopy was measured with net radiometers (NR Lite2, Kipp&Zonen, Delft, Netherlands), installed in the middle of one replication of each treatment.
Soil heat fluxes at 0.05 m soil depth in the middle of the bed were monitored in one replication of each treatment using soil heat flux plates (HFP01, Hukseflux, Delft, The Netherlands). According to Giambelluca et al. [41] the soil heat flux on the soil surface can be estimated as follows:
G = F + M
M = Δ T s t d ( ρ b C s + ρ w θ C w )
where G is the soil heat flux on the soil surface (W/m2) (W is the amount of energy received per second), F is the soil heat flux at 0.05 m soil depth (W/m2), M is the soil heat flux change between soil surface and soil heat flux sensor (W/m2), t is the time interval (s) (1800 s), ΔTs is the soil temperature change in the 0 to 0.05 m soil layer during the time interval t (°C), d is the distance between the soil surface and soil heat flux sensor (m) (0.05 m), ρb is the soil dry bulk density (kg/m3) (1530 kg/m3), ρw is the water density (kg/m3), Cs is the soil specific heat [J/(kg °C)] [840 J/(kg °C)], Cw is the water specific heat [J/(kg °C)], and θ is the volumetric soil water content in the 0 to 0.05 m soil layer (m3/m3) (where the volumetric soil water content at 0.1 m soil depth is used as an approximation).
The temperature sensors were connected to a data logger (SMC6108, Beijing Unism Technologies, Inc.). The net radiometers and soil heat flux plates were connected to another data logger (CR1000, Campbell Scientific Inc., Logan, UT, USA). The measurements were sampled at 10 s intervals and the averages were stored every 10 min.

3. Description of the Model

3.1. Energy Balance Equations

The energy transport in potato fields with plastic film mulch is illustrated in Figure 2. According to Ham and Kluitenberg [6], these assumptions are made for the model: (1) Soil evaporation is neglected because of the full plastic film mulch. Although soil evaporation still occurs at the hole of plant emergence, the evaporation is quite small compared with plant transpiration; (2) There is no heat stored in the plastic film mulch and plant; and (3) Heat transfer in the process of evaporation/condensation between plastic film mulch and soil is neglected. Thus, the energy balance of the plant canopy can be expressed as follows:
R n c H c a λ T + H m c = 0
where Rnc is the net radiation absorbed by the plant canopy (W/m2), Hca is the sensible heat flux from the plant canopy to the atmosphere (W/m2), λT is the latent heat flux of the plant canopy (W/m2), and Hmc is the sensible heat flux from the plastic film mulch to the plant canopy (W/m2).
According to Ham and Kluitenberg [6], the energy balance of the plastic film mulch can be defined as follows:
R n m H m c + C s m = 0
where Rnm is the net radiation absorbed by the plastic film mulch (W/m2), and Csm is the conduction from the soil surface to the plastic film mulch (W/m2). The energy balance of the soil surface can be defined as follows:
R n s C s m G = 0
where Rns is the net radiation absorbed by the soil surface (W/m2), and G is the soil heat flux at the soil surface (W/m2).

3.1.1. Governing Equations for Net Radiation of Each Layer

The radiation transfer is more complicated in the soil-film-canopy-atmosphere system than the soil-film-atmosphere system. The physical representation of the shortwave radiation transfer can be described in Figure 3.
The shortwave radiation Rsc, Rsm, and Rss absorbed by the plant canopy, the plastic film mulch, and the soil surface can be described as follows:
R s c = α s c R s + α s c ρ s m τ s c R s
R s m = α s m τ s c R s + α s m ρ s s τ s m τ s c R s
R s s = ( 1 ρ s s ) τ s m τ s c R s
where Rs is the global radiation (W/m2); αsc and αsm are the shortwave absorptances of the plant canopy and plastic film mulch, respectively; ρsm and ρss are the shortwave reflectances of the plastic film mulch and soil, respectively; and τsc and τsm are the shortwave transmittances of the plant canopy and plastic film mulch, respectively.
As shown in Figure 4a–d, the longwave radiation Rlc, Rlm, and Rls absorbed by the plant canopy, the plastic film mulch, and the soil surface can be described as follows:
R l c = ε c ε s k y σ T s k y 4 ( 1 + ρ l m τ l c ) + ε c σ T c 4 ( ρ l m ε c 2 ) + ε c ε m σ T m 4 + ε c τ l m ε s σ T s 4
R l m = ε m τ l c ε s k y σ T s k y 4 [ 1 + ( 1 ε s ) τ l m ] + ε m ε c σ T c 4 [ 1 + ( 1 ε s ) τ l m ] + ε m σ T m 4 [ ρ l c ε m + ( 1 ε s ) ε m 2 ] + ε m ε s σ T s 4 ( 1 + ρ l c τ l m )
R l s = ε s τ l m τ l c ε s k y σ T s k y 4 + ε s τ l m ε c σ T c 4 + ε s ε m σ T m 4 + ε s σ T s 4 ( ε s ρ l m 1 )
where σ is the Stefan-Boltzmann constant [W/(m2 °C4)]; Tsky (set equal to the air temperature) [6], Tc, Tm, and Ts are the temperatures of the sky, plant canopy, plastic film mulch, and soil surface (°C), respectively; εsky, εc, εm, and εs are the emissivities (or infrared absorptances) of the sky, plant canopy, plastic film mulch, and soil surface, respectively; ρlm and ρlc are the longwave reflectances of the plastic film mulch and plant canopy, respectively; and τlc and τlm are the longwave transmittances of the plant canopy and plastic film mulch, respectively.
Then, the net radiation Rnc, Rnm, and Rns absorbed by the plant canopy, the plastic film mulch, and the soil surface can be expressed as:
R n c = α s c R s ( 1 + ρ s m τ s c ) + ε c ε s k y σ T s k y 4 ( 1 + ρ l m τ l c ) + ε c σ T c 4 ( ρ l m ε c 2 ) + ε c ε m σ T m 4 + ε c τ l m ε s σ T s 4
R n m = α s m τ s c R s ( 1 + ρ s s τ s m ) + ε m τ l c ε s k y σ T s k y 4 [ 1 + ( 1 ε s ) τ l m ] + ε m ε c σ T c 4 [ 1 + ( 1 ε s ) τ l m ] + ε m σ T m 4 [ ρ l c ε m + ( 1 ε s ) ε m 2 ] + ε m ε s σ T s 4 ( 1 + ρ l c τ l m )
R n s = ( 1 ρ s s ) τ s m τ s c R s + ε s τ l m τ l c ε s k y σ T s k y 4 + ε s τ l m ε c σ T c 4 + ε s ε m σ T m 4 + ε s σ T s 4 ( ε s ρ l m 1 )

3.1.2. Governing Equations for Latent Heat Flux

According to Kang et al. (page 39) [42] the latent heat flux from the plant canopy can be defined as follows:
λ T = ρ a c p γ [ ( e s e a ) + Δ ( T c T a ) ] / ( r s c + r a c )
where ρa is the density of air (kg/m3), cp is the specific heat at constant pressure [J/(kg °C)], γ is the psychrometric constant (kPa/°C), Tc is the temperature of the plant canopy (°C), Ta is the air temperature (°C), es is the saturation vapour pressure (kPa), ea is the actual vapour pressure (kPa), Δ is the slope of the vapour pressure curve (kPa/°C), r s c is the bulk stomatal resistance of the plant canopy (s/m), and r a c is the bulk boundary layer resistance of the vegetative elements in the plant canopy (s/m).
The bulk stomatal resistance r s c and r a c are suggested by Shuttleworth and Wallace [43] as follows:
r s c = r S T 2 L A I
r a c = r b 2 L A I
where rST is the mean stomatal resistance (s/m), rb is the mean boundary layer resistance (s/m), and LAI is the leaf area index.

3.1.3. Governing Equations for Sensible Heat Flux

The sensible heat flux can be defined as follows:
H c a = ρ a c p ( T c T a ) r a a
H m c = ρ a c p ( T m T c ) r a m
where Tm is the temperature of the plastic film mulch (°C), r a a is the aerodynamic resistance between the plant canopy and the reference height (s/m), and r a m is the aerodynamic resistance between the plastic film mulch and the plant canopy (s/m).
The aerodynamic resistances r a a and r a m are given by Shuttleworth and Wallace [43] as follows:
r a a = 1 4 L A I   r a a ( α ) + 1 4 ( 4 L A I ) r a a ( 0 ) r a m = 1 4 L A I   r a m ( α ) + 1 4 ( 4 L A I ) r a m ( 0 ) }     0 L A I 4
r a a = r a a ( α ) r a m = r a m ( α ) }     L A I > 4
where r a a (α) and r a m (α) are the values of r a a and r a m , respectively, when the ground is covered fully by a crop canopy (s/m); and r a a (0) and r a m (0) are the values of r a a and r a m , respectively, when there is no crop (s/m).
The r a a (α) and r a m (α) can be expressed as follows:
r a a ( α ) = ln [ ( h r d ) / z 0 ] k 2 u ln [ ( h r d ) / ( h p d ) ] + h p n ( h p d ) × { exp [ n ( 1 d + z 0 h p ) ] 1 } r a m ( α ) = ln [ ( h r d ) / z 0 ] k 2 u h p n ( h p d ) exp n exp { n [ 1 ( d + z 0 ) / h p ] } }
where hr is the reference height (m), u is the wind speed at the reference height (m/s), n is the eddy diffusivity decay constant when the ground is fully covered by crop, k is von Kármán’s constant, hp is the plant height (m), d is the zero plane displacement when the ground is covered by crop fully (m), and z0 is the roughness height when the ground is covered fully by a crop canopy (m). The d and z0 can be given as follows:
d = 0.63 h p
z 0 = 0.13 h p
The r a a (0) and r a m (0) can be expressed as follows:
r a a ( 0 ) = ln 2 ( h r / z 0 ) k 2 u r a m ( 0 ) r a m ( 0 ) = ln ( h r / z 0 ) ln [ ( d + z 0 ) / z 0 ] k 2 u }
where z0′ is the roughness height when there is no crop (m).

3.1.4. Governing Equations for Soil Heat Flux

The soil heat flux G suggested by Horton and Chung [44] is as follows:
G = λ T l n 1 T s n Δ z + C ( T s n T s n 1 ) Δ z 2 Δ t
where Ts is the temperature of the soil surface (°C), Tl is the temperature of the first vertical node below the soil surface (°C), Δz is the distance between the soil surface and the first vertical node (m), Δt is the time step increment (s), n is the time step number, C is the volumetric soil heat capacity [J/(m3 °C)], λ is the soil thermal conductivity [W/(m °C)].
The volumetric soil heat capacity C is given by de Vries [45] as follows:
C ( 1.92 θ n + 2.51 θ o + 4.18 θ ) × 10 6
where θn is the solid phase content, θo is the organic matter content, and θ is the volumetric soil water content (m3/m3).
According to Chung and Horton [46] the soil thermal conductivity λ can be estimated as follows:
λ = b 1 + b 2 θ + b 3 θ 0.5
where b1, b2, and b3 are empirical parameters [W/(m °C)].

3.1.5. Governing Equations for Conduction From the Soil Surface to the Plastic Film Mulch

The heat conduction from the soil surface to the plastic film mulch Csm is given by Ham Kluitenberg [6] as follows:
C s m = T s T m r c
where rc is the thermal contact resistance between the plastic mulch and soil (m2 °C/W). According to Ham Kluitenberg [6] the conduction is dominated by heat transfer when the distance between the plastic film mulch and soil is less than 0.01 m. Then the rc can be defined as follows:
r c = z g k a N u d
where zg is the distance between the plastic film mulch and the soil (m), ka is the air thermal conductivity [W/(m °C)], and Nud is the Nusselt number reflecting the regime of conduction.
In summary, the energy balance Equations (3)–(5) can be expressed as follows:
α s c R s ( 1 + ρ s m τ s c ) + ε c ε s k y σ T s k y 4 ( 1 + ρ l m τ l c ) + ε c σ T c 4 ( ρ l m ε c 2 ) + ε c ε m σ T m 4 + ε c τ l m ε s σ T s 4 ρ a c p ( T c T a ) / r a a ρ a c p γ [ ( e s e a ) + Δ ( T c T a ) ] / ( r s c + r a c ) + ρ a c p ( T m T c ) / r a m = 0
α s m τ s c R s ( 1 + ρ s s τ s m ) + ε m τ l c ε s k y σ T s k y 4 [ 1 + ( 1 ε s ) τ l m ] + ε m ε c σ T c 4 [ 1 + ( 1 ε s ) τ l m ] + ε m σ T m 4 [ ρ l c ε m + ( 1 ε s ) ε m 2 ] + ε m ε s σ T s 4 ( 1 + ρ l c τ l m ) ρ a c p ( T m T c ) / r a m + ( T s T m ) / r c = 0
( 1 ρ s s ) τ s m τ s c R s + ε s τ l m τ l c ε s k y σ T s k y 4 + ε s τ l m ε c σ T c 4 + ε s ε m σ T m 4 + ε s σ T s 4 ( ε s ρ l m 1 ) ( T s T m ) / r c + λ ( T l n 1 T s n ) / Δ z C Δ z ( T s n T s n 1 ) / ( 2 Δ t ) = 0

3.2. Parameterization for the Energy Balance Equations

As many parameters were involved in the energy balance equations, the parameterization for the net radiation, latent heat flux, and heat conduction were introduced separately.

3.2.1. Parameterization for Net Radiation

Net Radiation of the Plant Canopy

The global radiation (Rs) was measured with the weather station. The shortwave transmittance of the potato canopy (τsc) is supposed to be equal to the photosynthetically active radiation transmittance (Tpar). According to Haverkort et al. [47] the shortwave reflectance of potato canopy can be approximated as follows:
ρ s c = 0.314 P c c + 0.203
where ρsc is the shortwave reflectance of the potato canopy and Pcc is the canopy coverage. Then the shortwave absorptance of the potato canopy (αsc) is:
α s c = 1 τ s c ρ s c
The emissivity (or infrared absorptance) of the potato canopy (εc) is 0.97 when the ground is covered fully by the canopy [27]. The relationship between εc and Pcc can be assumed as follows:
ε c = 0.97 P c c
The longwave reflectance of the potato canopy (ρlc) is about 0.01 [27]. According to Ham Kluitenberg [6] the longwave transmittance of potato canopy (τlc) is approximated as follows:
τ l c = 1 0.01 0.97 P c c

Net Radiation of the Plastic Film Mulch

The shortwave absorptance (αsm) and the shortwave transmittance (τsm) of the plastic film mulch were determined through parameter calibration with the field experiment data. The shortwave reflectance of the plastic film mulch (ρsm) can be obtained by the equation:
ρ s m = 1 α s m τ s m
The emissivity (or infrared absorptance) (εm) and the longwave transmittance (τlm) of the plastic film mulch were determined by parameter calibration. The longwave reflectance (ρlm) of the black film mulch and transparent film mulch are about 0.01 and 0.13, respectively [27].

Net Radiation of the Soil Surface

The shortwave reflectance (ρss) and emissivity (or infrared absorptance) (εs) of soil were determined by parameter calibration. According to Campbell and Norman (page 164) [48] the emissivity of sky (εsky) is defined as follows:
ε s k y = 9.2 × 10 6 T a 2

3.2.2. Parameterization for Latent Heat Flux

The air temperature (Ta), actual atmospheric pressure (P), and relative humidity (RH) were monitored with weather station. According to Maiti et al. [49] the density of air (ρa) can be calculated as follows:
ρ a = ρ 0 273 273 + T a
where ρ0 (the air density at 0 °C and 0.1013 Mpa) is about 1.29 kg/m3. According to Allen et al. [50], the psychrometric constant (γ), the slope of the vapour pressure curve (Δ), the saturation vapour pressure (es), and the actual vapour pressure (ea) can be expressed as follows:
γ = 0.665 × 10 3 P
Δ = 4098 [ 0.6108 exp ( 17.27 T a T a + 237.3 ) ] ( T a + 237.3 ) 2
e s ( T a ) = 0.6108 exp ( 17.27 T a T a + 237.3 )
e a = e s ( T min ) R H max 100 + e s ( T max ) R H min 100 2
where Tmin and Tmax are the minimum and maximum air temperature at a certain period of time (°C), respectively; es(Tmin) and es(Tmax) are the saturation vapour pressure at Tmin and Tmax (kPa), respectively; RHmin and RHmax are the minimum and maximum relative humidity (%), respectively.
The mean stomatal resistance (rST) and the mean boundary layer resistance (rb) of the potato canopy were determined through parameter calibration with field experiment data.
As the direct measurement of leaf area index (LAI) is destructive and time consuming, the daily LAI was determined by the fraction of intercepted photosynthetically active radiation (Ipar) indirectly. According to Haverkort et al. [47] and Boyd et al. [51] the relationship between Ipar and LAI can be expressed as follows:
I p a r = 1 e κ L A I
where κ is the extinction coefficient. The κ can be determined through the measured Tpar and LAI. In this experiment κ is 0.92, close to the result of Spitters [52]. The relationship between the photosynthetically active radiation transmittance (Tpar) and Ipar is:
T p a r = 1 I p a r

3.2.3. Parameterization for Sensible Heat Flux

The wind speed at the reference height (u) was monitored with the weather station. The plant height (hp) was determined using interpolation method with plant height measured in the field.
The reference height (hr) was 2 m. The roughness height (z0′) is 0.01 m [43]. The eddy diffusivity decay constant (n) is 2.5 and the von Kármán’s constant (k) is 0.41 [43].

3.2.4. Parameterization for Soil Heat Flux

The temperature (Tl) of the first vertical node above the soil surface was monitored with sensors. The distance (Δz) between the soil surface and the first vertical node was assumed to be 0.1 m and the time step increment (Δt) was 3600 s.

3.2.5. Parameterization for Heat Conduction

The air thermal conductivity (ka) is about 0.025 W/(m °C) ([48]; page 118, Table 8.2). The Nusselt number (Nud) is 1 [6]. The distance (zg) between the plastic film mulch and soil was determined through parameter calibration or measurement with field experiment data.
The parameterization for the net radiation is shown in Table 1. The parameterizations for the latent heat flux, sensible heat flux, soil heat flux, and heat conduction are shown in Table 2.

3.3. Numerical Solution of the Energy Balance Equations

After the parameterization, the energy balance equations were functions of Tc, Tm, and Ts. Equations (31)–(33) can be expressed as a nonlinear equation group F(T) as follows:
F ( T ) = { f 1 ( T c , T m , T s ) = 0 f 2 ( T c , T m , T s ) = 0 f 3 ( T c , T m , T s ) = 0
The equation group (47) can be solving by Newton-Raphson iteration technique [6]. The partial derivatives of the iteration technique can be described as a Jacobian matrix F’(T) as follows [6]:
F ( T ) = [ f 1 ( T ) T c f 1 ( T ) T m f 1 ( T ) T s f 2 ( T ) T c f 2 ( T ) T m f 2 ( T ) T s f 3 ( T ) T c f 3 ( T ) T m f 3 ( T ) T s ]
where:
f 1 ( T ) T c = 4 ε c σ T c 3 ( ρ l m ε c 2 ) ρ a c p / r a a ρ a c p γ Δ / ( r s c + r a c ) ρ a c p / r a m
f 1 ( T ) T m = 4 ε c ε m σ T m 3 + ρ a c p / r a m
f 1 ( T ) T s = 4 ε c τ l m ε s σ T s 3
f 2 ( T ) T c = 4 ε m ε c σ T c 3 [ 1 + ( 1 ε s ) τ l m ] + ρ a c p / r a m
f 2 ( T ) T m = 4 ε m σ T m 3 [ ρ l c ε m + ( 1 ε s ) ε m 2 ] ρ a c p / r a m 1 / r c
f 2 ( T ) T s = 4 ε m ε s σ T s 3 ( 1 + ρ l c τ l m ) + 1 / r c
f 3 ( T ) T c = 4 ε s τ l m ε c σ T c 3
f 3 ( T ) T m = 4 ε s ε m σ T m 3 + 1 / r c
f 3 ( T ) T s = 4 ε s σ T s 3 ( ε s ρ l m 1 ) 1 / r c λ ( θ ) / Δ z C ( θ ) Δ z / ( 2 Δ t )

3.4. Statistical Analysis

To evaluate the model performance three statistical parameters, including root mean square errors (RMSE), coefficient of determination (R2), and mean errors (ME) were applied during model calibration and validation. These statistical parameters were calculated from simulated and measured temperatures, soil heat fluxes, and net radiations using the equations as follows:
R M S E = 1 N i = 1 N ( O i P i ) 2
R 2 = [ i = 1 N ( O i O ¯ ) ( P i P ¯ ) i = 1 N ( O i O ¯ ) 2 i = 1 N ( P i P ¯ ) 2 ] 2
M E = 1 N i = 1 N ( O i P i )
where Oi and Pi are the observed and simulated values, respectively; O ¯ and P ¯ are the averages of the observed and simulated values, respectively; and N is the number of observations.

4. Results and Discussion

The Pcc (fractional canopy coverage) and Tpar (photosynthetically active radiation transmittance of the potato canopy) were measured to reflect the potato canopy growth (Figure 5; Figure 6). The black plastic film mulch had greater plant canopy in the BM treatment than in the TM treatment in 2014 and 2015. It might be caused by the higher air temperature near the black mulch surface [11].

4.1. Model Calibration and Validation

Model parameters were calibrated with the temperature of the plant canopy (Tc), the temperature of the plastic film mulch (Tm), the temperature of the soil surface (Ts), the net radiation (Rn), and the soil heat flux (G) which were measured in the potato field in 2014. The calibrated parameters for the transparent and black plastic film mulch treatments are shown in Table 3. The calibrated model was validated with Tc, Tm, Ts, Rn, and G measured in the potato field in 2015.

4.2. Model Calibration and Validation for Temperature

Generally, the model can reflect the daily fluctuations of mulch surface temperature and soil surface temperature for the TM and BM treatments by comparing the simulated and observed temperatures (Figure 7 and Figure 8; Table 4; Table 5). In 2015, the daily RMSE of soil surface temperature for the TM and BM treatments, which were smaller than the hourly RMSE, were 2.8 and 1.5 °C, respectively. This result was similar with Ham and Kluitenberg [6] who reported that the RMSE of soil temperature was 2.4 °C for the TM treatment and 1.8 °C for the BM treatment. In 2014 and 2015, the hourly and daily R2 of temperatures ranged from 0.65 to 0.84 and 0.76 to 0.91, respectively.
During the early plant growth stage (23 May to 27 May), the simulations were not satisfactory (Figure 7 and Figure 8). The larger simulation deviation might be caused by the neglect of soil evaporation from the hole (8 cm in diameter) of the plant emergence. Moreover, the temperature sensors on the mulch surface were affected by solar radiation during the early plant growth stage, which might also cause the deviation. As the plant canopy got larger, the simulation got better during the later plant growth stages. The variations of mulch surface temperature and soil surface temperature were smaller during later potato growth because of the larger area shaded by the plant canopy. Generally, the model could simulate the daily temperature variations at different potato growth stages for different plastic film mulch treatments.

4.3. Model Calibration and Validation for Net Radiation and Soil Heat Flux

The hourly simulated and observed net radiation for the TM and BM treatments were compared and the statistical hourly and daily RMSE, R2, and ME were calculated (Figure 9 and Table 6; Table 7). The daily simulation of net radiation was reasonably good with daily RMSE for the TM and BM treatments 20.30 and 32.79 W/m2 during validation, respectively. The ME of daily net radiation for the two treatments ranged from 12.96 to 36.14 W/m2 in 2014 and 2015. The R2 ranged from 0.89 to 0.98, which was greater than 0.8, a typical value for model performance according to Baldocchi and Wilson [53]. An under-estimation of net radiation for the TM treatment and over-estimation for the BM treatment were observed at noon in 2015 (Figure 9). Despite the simulation deviations, the model reflected the daily fluctuations of net radiation during different growth stages of potato growth.
Generally, the model reflected the daily variation of soil heat flux with daily RMSE of the TM and BM treatments ranging from 17.03 to 25.27 W/m2 and ME ranging from −11.30 to 23.35 W/m2 during the calibration and validation presses (Figure 10 and Table 6; Table 7). However, the hourly and daily R2 of soil heat flux were rather small. This might be caused by some abnormal values during soil heat flux measurement.

4.4. Canopy Temperature Analysis With the Model

The canopy temperature can affect the aboveground plant growth, especially during the early plant growth stage [11,54]. The canopy temperature was estimated with the calibrated model for the TM and BM treatments in 2015 (Figure 11). The canopy temperature of the BM treatment was greater than for the TM treatment. The maximum temperature difference between the TM and BM treatments could be as much as 7 °C during early potato growth. The higher canopy temperature in the BM treatment could be the reason why the aboveground growth in the BM treatment was greater than in the TM treatment [5,11]. The model could be a useful tool to study the spatiotemporal distribution of temperature, net radiation, and soil heat flux in potato field with different plastic film mulches and their effects on potato growth during different growth stages.

5. Summary and Conclusions

A numerical model was developed to simulate the heat transport in potato field with full plastic film mulch in an arid area. This model predicted the interaction of potato plant growth and the optical properties of plastic film mulch on heat transport by utilizing the leaf area index, canopy coverage, photosynthetically active radiation transmittance, and considering the radiation transfer in atmosphere-canopy-mulch-soil system. The heat transfer parameters (net radiation, mulch surface temperature, soil surface temperature, and soil heat flux) were measured in the field experiments in 2014 and 2015 growing seasons in potato field with two conventional plastic film mulch treatments (TM and BM treatments) to calibrate (2014) and validate (2015) the model.
Generally, the temperature simulation was reasonably good with daily RMSE of soil surface temperature 2.8 and 1.5 °C for the TM and BM treatments, respectively. The model reflected the daily temperature variations during different potato growth stages except during early plant growth. Moreover, the model simulated the daily fluctuations of net radiation with daily RMSE for the TM and BM treatments 20.30 and 32.79 W/m2 during validation, respectively. The daily R2 of net radiation ranged from 0.89 to 0.98, while the daily R2 of soil heat flux were small because of some abnormal values during soil heat flux measurement. The simulation showed that the BM treatment had higher canopy temperature than the TM treatment. The maximum canopy temperature difference between TM and BM treatments could be as high as 7 °C during early potato growth.
In conclusion, the model could give an explanation for the interaction of plastic film mulch and potato canopy growth on heat transport. The model can be used to simulate heat transport in potato fields with different plastic film mulches in semiarid areas. The heat conditions simulation can provide guidance in plastic film choosing for potato cultivation to avoid heat stress. One of the limitations to study in the future is the cause of the simulation deviations, especially during early crop growth.

Author Contributions

Data curation, Y.-L.Z.; Formal analysis, Y.-L.Z.; Funding acquisition, F.-X.W. and Y.-L.Z.; Methodology, Y.-L.Z. and F.-X.W.; Software, Y.-L.Z.; Supervision, F.-X.W.; Writing—original draft, Y.-L.Z.; Writing—review & editing, F.-X.W., C.C.S. and S.-Y.F. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported financially by the National Natural Science Foundation of China (51579240, 51809225), the Ministry of Water Resources of China (201501017), and the Natural Science Foundation of Jiangsu Province, China (BK20180929).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, F.M.; Wang, J.; Xu, J.Z.; Xu, H.L. Productivity and soil response to plastic film mulching durations for spring wheat on entisols in the semiarid loess plateau of China. Soil Tillage Res. 2004, 78, 9–20. [Google Scholar] [CrossRef]
  2. Bonachela, S.; Granados, M.R.; López, J.C.; Hernández, J.; Magán, J.J.; Baeza, E.J.; Baille, A. How plastic mulches affect the thermal and radiative microclimate in an unheated low-cost greenhouse. Agric. For. Meteorol. 2012, 152, 65–72. [Google Scholar] [CrossRef]
  3. Cavero, J.; Ortega, R.G.; Zaragoza, C. Clear plastic mulch improved seedling emergence of direct-seeded pepper. Hortscience 1996, 31, 70–73. [Google Scholar] [CrossRef]
  4. Chellemi, D.O.; Olson, S.M.; Mitchell, D.J.; Secker, I.; McSorley, R. Adaptation of soil solarization to the integrated management of soilborne pests of tomato under humid conditions. Phytopathology 1997, 87, 250–258. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Decoteau, D.R.; Kasperbauer, M.J.; Daniels, D.D.; Hunt, P.G. Plastic mulch color effects on reflected light and tomato plant growth. Sci. Hortic. 1988, 34, 169–175. [Google Scholar] [CrossRef]
  6. Ham, J.M.; Kluitenberg, G.J. Modeling the effect of mulch optical properties and mulch-soil contact resistance on soil heating under plastic mulch culture. Agric. For. Meteorol. 1994, 71, 403–424. [Google Scholar] [CrossRef]
  7. Hou, X.Y.; Wang, F.X.; Han, J.J.; Kang, S.Z.; Feng, S.Y. Duration of plastic mulch for potato growth under drip irrigation in an arid region of northwest China. Agric. For. Meteorol. 2010, 150, 115–121. [Google Scholar] [CrossRef]
  8. Liakatas, A.; Clark, J.A.; Monteith, J.L. Measurements of the heat balance under plastic mulches. Part I. Radiation balance and soil heat flux. Agric. For. Meteorol. 1986, 36, 227–239. [Google Scholar] [CrossRef]
  9. Liu, C.A.; Jin, S.L.; Zhou, L.M.; Jia, Y.; Li, F.M.; Xiong, Y.C.; Li, X.G. Effects of plastic film mulch and tillage on maize productivity and soil parameters. Eur. J. Agron. 2009, 31, 241–249. [Google Scholar] [CrossRef]
  10. Qin, S.J.; Li, S.E.; Yang, K.; Hu, K.L. Can plastic mulch save water at night in irrigated croplands? J. Hydrol. 2018, 564, 667–681. [Google Scholar] [CrossRef]
  11. Zhang, Y.L.; Wang, F.X.; Shock, C.C.; Yang, K.J.; Kang, S.Z.; Qin, J.T.; Li, S.E. Effects of plastic mulch on the radiative and thermal conditions and potato growth under drip irrigation in arid northwest China. Soil Tillage Res. 2017, 172, 1–11. [Google Scholar] [CrossRef]
  12. Anikwe, M.A.N.; Mbah, C.N.; Ezeaku, P.I.; Onyia, V.N. Tillage and plastic mulch effects on soil properties and growth and yield of cocoyam (colocasia esculenta) on an ultisol in southeastern Nigeria. Soil Till. Res. 2007, 93, 264–272. [Google Scholar] [CrossRef]
  13. Liang, S.M.; Ren, C.; Wang, P.J.; Wang, X.T.; Li, Y.S.; Xu, F.H.; Wang, Y.; Dai, Y.Q.; Zhang, L.; Li, X.P.; et al. Improvements of emergence and tuber yield of potato in a seasonal spring arid region using plastic film mulching only on the ridge. Field Crops Res. 2018, 223, 57–65. [Google Scholar] [CrossRef]
  14. Gao, Y.H.; Xie, Y.P.; Jiang, H.Y.; Wu, B.; Niu, J.Y. Soil water status and root distribution across the rooting zone in maize with plastic film mulching. Field Crops Res. 2014, 156, 40–47. [Google Scholar] [CrossRef]
  15. Zhang, Y.L.; Wang, F.X.; Shock, C.C.; Yang, K.J.; Kang, S.Z.; Qin, J.T.; Li, S.E. Influence of different plastic film mulches and wetted soil percentages on potato grown under drip irrigation. Agric. Water Manag. 2017, 180, 160–171. [Google Scholar] [CrossRef]
  16. Kasirajan, S.; Ngouajio, M. Polyethylene and biodegradable mulches for agricultural applications: A review. Agron. Sustain. Dev. 2012, 32, 501–529. [Google Scholar] [CrossRef]
  17. Dong, B.D.; Liu, M.Y.; Jiang, J.W.; Shi, C.H.; Wang, X.M.; Qiao, Y.Z.; Liu, Y.Y.; Zhao, Z.H.; Si, F.Y. Growth, grain yield, and water use efficiency of rain-fed spring hybrid millet (Setaria italica) in plastic-mulched and unmulched fields. Agric. Water Manag. 2014, 143, 93–101. [Google Scholar] [CrossRef]
  18. Liu, X.E.; Li, X.G.; Hai, L.; Wang, Y.P.; Li, F.M. How efficient is film fully-mulched ridge–furrow cropping to conserve rainfall in soil at a rainfed site? Field Crops Res. 2014, 169, 107–115. [Google Scholar] [CrossRef]
  19. Qin, S.H.; Zhang, J.L.; Dai, H.L.; Wang, D.; Li, D.M. Effect of ridge–furrow and plastic-mulching planting patterns on yield formation and water movement of potato in a semi-arid area. Agric. Water Manag. 2014, 131, 87–94. [Google Scholar] [CrossRef]
  20. Zhang, F.; Zhang, W.J.; Qi, J.G.; Li, F.M. A regional evaluation of plastic film mulching for improving crop yields on the Loess Plateau of China. Agric. For. Meteorol. 2018, 248, 458–468. [Google Scholar] [CrossRef]
  21. Zhao, H.; Xiong, Y.C.; Li, F.M.; Wang, R.Y.; Qiang, S.C.; Yao, T.F.; Mo, F. Plastic film mulch for half growing-season maximized WUE and yield of potato via moisture-temperature improvement in a semi-arid agroecosystem. Agric. Water Manag. 2012, 104, 68–78. [Google Scholar] [CrossRef]
  22. Kader, M.A.; Senge, M.; Mojid, M.A.; Ito, K. Recent advances in mulching materials and methods for modifying soil environment. Soil Tillage Res. 2017, 168, 155–166. [Google Scholar] [CrossRef]
  23. Qin, X.L.; Li, Y.Z.; Han, Y.L.; Hu, Y.C.; Li, Y.J.; Wen, X.X.; Liao, Y.C.; Siddique, K.H.M. Ridge-furrow mulching with black plastic film improves maize yield more than white plastic film in dry areas with adequate accumulated temperature. Agric. For. Meteorol. 2018, 262, 206–214. [Google Scholar] [CrossRef]
  24. Ruiz, J.M.; Hernandez, J.; Castilla, N.; Romero, L. Potato performance in response to different mulches. 1. Nitrogen metabolism and yield. J. Agric. Food Chem. 1999, 47, 2660–2665. [Google Scholar] [CrossRef]
  25. Farias-Larios, J.; Orozco-Santos, M. Effect of polyethylene mulch colour on aphid populations, soil temperature, fruit quality, and yield of watermelon under tropical conditions. N. Z. J. Crop Hortic. Sci. 1997, 25, 369–374. [Google Scholar] [CrossRef]
  26. Yaghi, T.; Arslan, A.; Naoum, F. Cucumber (cucumis sativus, L.) water use efficiency (WUE) under plastic mulch and drip irrigation. Agric. Water Manag. 2013, 128, 149–157. [Google Scholar] [CrossRef] [Green Version]
  27. Tarara, J.M. Microclimate modification with plastic mulch. Hortscience 2000, 35, 169–180. [Google Scholar] [CrossRef] [Green Version]
  28. Mahrer, Y. Prediction of soil temperatures of a soil mulched with transparent polyethylene. J. Appl. Meteorol. 1973, 18, 1263–1267. [Google Scholar] [CrossRef] [Green Version]
  29. Mahrer, Y.; Katan, J. Spatial soil temperature regime under transparent polyethylene mulch: Numerical and experimental studies. Soil Sci. 1981, 131, 82–87. [Google Scholar] [CrossRef]
  30. Mahrer, Y.; Naot, O.; Rawitz, E.; Katan, J. Temperature and moisture regimes in soils mulched with transparent polyethylene. Soil Sci. Soc. Am. J. 1984, 48, 362–367. [Google Scholar] [CrossRef]
  31. Wu, C.L.; Chau, K.W.; Huang, J.S. Modelling coupled water and heat transport in a soil–mulch–plant–atmosphere continuum (SMPAC) system. Appl. Math. Model. 2007, 31, 152–169. [Google Scholar] [CrossRef] [Green Version]
  32. Yang, Q.D.; Zuo, H.C.; Xiao, X.; Wang, S.J.; Chen, B.L.; Chen, J.W. Modelling the effects of plastic mulch on water, heat and CO2 fluxes over cropland in an arid region. J. Hydrol. 2012, 452–453, 102–118. [Google Scholar] [CrossRef]
  33. Flerchinger, G.N.; Sauer, T.J.; Aiken, R.A. Effects of crop residue cover and architecture on heat and water transfer at the soil surface. Geoderma 2003, 116, 217–233. [Google Scholar] [CrossRef]
  34. Flerchinger, G.N.; Pierson, F.B. Modeling plant canopy effects on variability of soil temperature and water. Agric. For. Meteorol. 1991, 56, 227–246. [Google Scholar] [CrossRef]
  35. Hijmans, R.J. The effect of climate change on global potato production. Am. J. Potato Res. 2003, 80, 271–280. [Google Scholar] [CrossRef]
  36. Van Dam, J.; Kooman, P.L.; Struik, P.C. Effects of temperature and photoperiod on early growth and final number of tubers in potato (Solanum tuberosum L.). Potato Res. 1996, 39, 51–62. [Google Scholar] [CrossRef]
  37. Rykaczewska, K. The effect of high temperature occurring in subsequent stages of plant development on potato yield and tuber physiological defects. Am. J. Potato Res. 2015, 92, 339–349. [Google Scholar] [CrossRef] [Green Version]
  38. Zhang, Y.L.; Feng, S.Y.; Wang, F.X.; Binley, A. Simulation of soil water flow and heat transport in drip irrigated potato field with raised beds and full plastic-film mulch in a semiarid area. Agric. Water Manag. 2018, 209, 178–187. [Google Scholar] [CrossRef] [Green Version]
  39. Bojacá, C.R.; García, S.J.; Schrevens, E. Analysis of potato canopy coverage as assessed through digital imagery by nonlinear mixed effects models. Potato Res. 2011, 54, 237–252. [Google Scholar] [CrossRef]
  40. Purcell, L.C. Soybean canopy coverage and light interception measurements using digital imagery. Crop Sci. 2000, 40, 834–837. [Google Scholar] [CrossRef]
  41. Giambelluca, T.W.; Scholz, F.G.; Bucci, S.J.; Meinzer, F.C.; Goldstein, G.; Hoffmann, W.A.; Franco, A.C.; Buchert, M.P. Evapotranspiration and energy balance of Brazilian savannas with contrasting tree density. Agric. For. Meteorol. 2009, 149, 1365–1376. [Google Scholar] [CrossRef]
  42. Kang, S.Z.; Liu, X.M.; Xiong, Y.Z. Theory and Application of Water Transport in Soil-Plant-Atmosphere Continuum; Water Conservancy and Electric Power Press: Beijing, China, 1994. [Google Scholar]
  43. Shuttleworth, W.J.; Wallace, J.S. Evaporation from sparse crops‒an energy combination theory. Q. J. R. Meteor. Soc. 1985, 111, 839–855. [Google Scholar] [CrossRef]
  44. Horton, R.; Chung, S.O. Soil Heat Flow. In Modeling Plant and Soil Systems; Hanks, J., Ritchie, J., Eds.; American Society of Agronomists: Madison, WI, USA, 1991. [Google Scholar]
  45. de Vries, D.A. The thermal properties of soils. In Physics of Plant Environment; van Wijk, R.W., Ed.; North Holland: Amsterdam, The Netherlands, 1963; pp. 210–235. [Google Scholar]
  46. Chung, S.O.; Horton, R. Soil heat and water flow with a partial surface mulch. Water Resour. Res. 1987, 23, 2175–2186. [Google Scholar] [CrossRef] [Green Version]
  47. Haverkort, A.J.; Uenk, D.; Veroude, H.; van de Waart, M. Relationships between ground cover, intercepted solar radiation, leaf area index and infrared reflectance of potato crops. Potato Res. 1991, 34, 113–121. [Google Scholar] [CrossRef]
  48. Campbell, G.S.; Norman, J.M. An Introduction to Environmental Biophysics, 2nd ed.; Springer-Verlag: New York, NY, USA, 1998. [Google Scholar]
  49. Maiti, S.; Patel, P.; Vyas, K.; Eswaran, K.; Ghosh, P.K. Performance evaluation of a small scale indirect solar dryer with static reflectors during non-summer months in the Saurashtra region of western India. Sol. Energy 2011, 85, 2686–2696. [Google Scholar] [CrossRef]
  50. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration: Guidelines for Computing Crop Water Requirements; FAO Irrigation and Drainage Paper No. 56; FAO: Rome, Italy, 1998. [Google Scholar]
  51. Boyd, N.S.; Gordon, R.; Martin, R.C. Relationship between leaf area index and ground cover in potato under different management conditions. Potato Res. 2002, 45, 117–129. [Google Scholar] [CrossRef]
  52. Spitters, C.J.T. An analysis of variation in yield among potato cultivars in terms of light absorption, light utilization and dry matter partitioning. Acta Hortic. 1987, 214, 71–84. [Google Scholar] [CrossRef]
  53. Baldocchi, D.D.; Wilson, K.B. Modeling CO2 and water vapor exchange of a temperate broadleaved forest across hourly to decadal time scales. Ecol. Model. 2001, 142, 155–184. [Google Scholar] [CrossRef]
  54. Ham, J.M.; Kluitenberg, G.J.; Lamont, W.J. Optical properties of plastic mulches affect the field temperature regime. J. Am. Soc. Hortic. Sci. 1993, 118, 188–193. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Example of potato canopy cover images taken at Wuwei, China on 2 June 2015 before and after processing for obtaining fractional canopy coverage.
Figure 1. Example of potato canopy cover images taken at Wuwei, China on 2 June 2015 before and after processing for obtaining fractional canopy coverage.
Agronomy 10 00190 g001
Figure 2. Schematic diagram of energy transport in a potato field with plastic film mulch (Rnc, Rnm, and Rns are the net radiations absorbed by the plant canopy, the plastic film mulch, and the soil surface, respectively; Hca is the sensible heat flux from the plant canopy to the atmosphere; Hmc is the sensible heat flux from the plastic film mulch to the plant canopy; Csm is the conduction from the soil surface to the plastic film mulch; G is the soil heat flux at the soil surface; λT is the latent heat flux of the plant canopy; λET is the latent heat flux of complete canopy; r s c is the bulk stomatal resistance of the plant canopy; r a c is the bulk boundary layer resistance of the vegetative elements in the plant canopy; r a a is the aerodynamic resistance between the plant canopy and the reference height; r a m is the aerodynamic resistance between the plastic film mulch and the plant canopy; rc is the thermal contact resistance between the plastic mulch and soil).
Figure 2. Schematic diagram of energy transport in a potato field with plastic film mulch (Rnc, Rnm, and Rns are the net radiations absorbed by the plant canopy, the plastic film mulch, and the soil surface, respectively; Hca is the sensible heat flux from the plant canopy to the atmosphere; Hmc is the sensible heat flux from the plastic film mulch to the plant canopy; Csm is the conduction from the soil surface to the plastic film mulch; G is the soil heat flux at the soil surface; λT is the latent heat flux of the plant canopy; λET is the latent heat flux of complete canopy; r s c is the bulk stomatal resistance of the plant canopy; r a c is the bulk boundary layer resistance of the vegetative elements in the plant canopy; r a a is the aerodynamic resistance between the plant canopy and the reference height; r a m is the aerodynamic resistance between the plastic film mulch and the plant canopy; rc is the thermal contact resistance between the plastic mulch and soil).
Agronomy 10 00190 g002
Figure 3. Schematic diagram of shortwave radiation transfer in the potato canopy, plastic film mulch, and soil surface (Rs is the global radiation; αsc and αsm are the shortwave absorptances of the plant canopy and plastic film mulch, respectively; ρsm and ρss are the shortwave reflectances of the plastic film mulch and soil, respectively; and τsc and τsm are the shortwave transmittances of the plant canopy and plastic film mulch, respectively).
Figure 3. Schematic diagram of shortwave radiation transfer in the potato canopy, plastic film mulch, and soil surface (Rs is the global radiation; αsc and αsm are the shortwave absorptances of the plant canopy and plastic film mulch, respectively; ρsm and ρss are the shortwave reflectances of the plastic film mulch and soil, respectively; and τsc and τsm are the shortwave transmittances of the plant canopy and plastic film mulch, respectively).
Agronomy 10 00190 g003
Figure 4. Schematic diagram of emitted longwave radiation transfer between the potato canopy, plastic film mulch, and soil surface from: (a) the atmosphere, (b) the potato canopy, (c) the plastic film mulch, and (d) the soil surface (σ is the Stefan-Boltzmann constant; Tsky, Tc, Tm, and Ts are the temperatures of the sky, plant canopy, plastic film mulch, and soil surface, respectively; εsky, εc, εm, and εs are the emissivities (or infrared absorptances) of the sky, plant canopy, plastic film mulch, and soil surface, respectively; ρlm and ρlc are the longwave reflectances of the plastic film mulch and plant canopy, respectively; and τlc and τlm are the longwave transmittances of the plant canopy and plastic film mulch, respectively).
Figure 4. Schematic diagram of emitted longwave radiation transfer between the potato canopy, plastic film mulch, and soil surface from: (a) the atmosphere, (b) the potato canopy, (c) the plastic film mulch, and (d) the soil surface (σ is the Stefan-Boltzmann constant; Tsky, Tc, Tm, and Ts are the temperatures of the sky, plant canopy, plastic film mulch, and soil surface, respectively; εsky, εc, εm, and εs are the emissivities (or infrared absorptances) of the sky, plant canopy, plastic film mulch, and soil surface, respectively; ρlm and ρlc are the longwave reflectances of the plastic film mulch and plant canopy, respectively; and τlc and τlm are the longwave transmittances of the plant canopy and plastic film mulch, respectively).
Agronomy 10 00190 g004aAgronomy 10 00190 g004b
Figure 5. Fractional canopy coverage by potato grown with transparent plastic film (TM) and black plastic film (BM) treatments at Wuwei, China, in 2014 and 2015.
Figure 5. Fractional canopy coverage by potato grown with transparent plastic film (TM) and black plastic film (BM) treatments at Wuwei, China, in 2014 and 2015.
Agronomy 10 00190 g005
Figure 6. Photosynthetically active radiation transmittance of the transparent plastic film mulch (TM) treatment and black plastic film mulch (BM) treatment for potato production at Wuwei, China, in 2014 and 2015.
Figure 6. Photosynthetically active radiation transmittance of the transparent plastic film mulch (TM) treatment and black plastic film mulch (BM) treatment for potato production at Wuwei, China, in 2014 and 2015.
Agronomy 10 00190 g006
Figure 7. Comparison of the observed and simulated hourly temperatures on the transparent plastic film mulch (TM) surface and soil surface at Wuwei, China, in 2015.
Figure 7. Comparison of the observed and simulated hourly temperatures on the transparent plastic film mulch (TM) surface and soil surface at Wuwei, China, in 2015.
Agronomy 10 00190 g007
Figure 8. Comparison of the observed and simulated hourly temperature on the black plastic film mulch (BM) surface and soil surface at Wuwei, China, in 2015.
Figure 8. Comparison of the observed and simulated hourly temperature on the black plastic film mulch (BM) surface and soil surface at Wuwei, China, in 2015.
Agronomy 10 00190 g008
Figure 9. Comparison of the observed and simulated hourly net radiation of the transparent plastic film mulch (TM) treatment and black plastic film mulch (BM) treatment at Wuwei, China, in 2015.
Figure 9. Comparison of the observed and simulated hourly net radiation of the transparent plastic film mulch (TM) treatment and black plastic film mulch (BM) treatment at Wuwei, China, in 2015.
Agronomy 10 00190 g009
Figure 10. Comparison of the observed and simulated hourly soil heat flux of the transparent plastic film mulch (TM) treatment and black plastic film mulch (BM) treatment at Wuwei, China, in 2015.
Figure 10. Comparison of the observed and simulated hourly soil heat flux of the transparent plastic film mulch (TM) treatment and black plastic film mulch (BM) treatment at Wuwei, China, in 2015.
Agronomy 10 00190 g010
Figure 11. Simulated hourly canopy temperature of potato plants grown with transparent plastic film mulch (TM) and black plastic film mulch (BM) at Wuwei, China, in 2015.
Figure 11. Simulated hourly canopy temperature of potato plants grown with transparent plastic film mulch (TM) and black plastic film mulch (BM) at Wuwei, China, in 2015.
Agronomy 10 00190 g011
Table 1. Input parameters for the net radiation calculation using in modeling the interaction of plastic film mulch and potato canopy growth on soil heat transport.
Table 1. Input parameters for the net radiation calculation using in modeling the interaction of plastic film mulch and potato canopy growth on soil heat transport.
SymbolDescriptionValueSource
Canopy
PccCanopy coverageMeasured in the fieldDigital camera
RsGlobal radiationMeasured in the fieldMeteorological station
αscShortwave absorptance of plant canopyCalculated
εcEmissivity (or infrared absorptance) of plant canopyCalculated
ρlcLongwave reflectance of plant canopy0.01Tarara [27]
ρscShortwave reflectance of plant canopyCalculated
τlcLongwave transmittance of plant canopyCalculated
τscShortwave transmittance of plant canopyMeasured in the fieldSunscan
Plastic film mulch
αsmShortwave absorptance of transparent plastic film mulch0.01~0.1 *Tarara [27]
Shortwave absorptance of black plastic film mulch0.8~0.98 *Tarara [27]
εmEmissivity (or infrared absorptance) of transparent plastic film mulch0.05~0.97 *Tarara [27]
Emissivity (or infrared absorptance) of black plastic film mulch0.8~0.98 *Tarara [27]
ρlmLongwave reflectance of transparent plastic film mulch0.13Tarara [27]
Longwave reflectance of black plastic film mulch0.01Tarara [27]
ρsmShortwave reflectance of transparent plastic film mulchCalculatedTarara [27]
Shortwave reflectance of black plastic film mulchCalculatedTarara [27]
τlmLongwave transmittance of transparent plastic film mulch0.7~0.9 *Tarara [27]
Longwave transmittance of black plastic film mulch0.1~0.67 *Tarara [27]
τsmShortwave transmittance of transparent plastic film mulch0.8~0.98 *Tarara [27]
Shortwave transmittance of black plastic film mulch0.01~0.1 *Tarara [27]
Soil
εsEmissivity (or infrared absorptance) of soil surface0.78~0.98 *Tarara [27]
ρssShortwave reflectance of soil0.15~0.3 *Tarara [27]
* The value was calibrated with field experiments.
Table 2. Input parameters for the latent heat flux, sensible heat flux, soil heat flux, and heat conduction calculations using in modeling the interaction of plastic film mulch and potato canopy growth on soil heat transport.
Table 2. Input parameters for the latent heat flux, sensible heat flux, soil heat flux, and heat conduction calculations using in modeling the interaction of plastic film mulch and potato canopy growth on soil heat transport.
SymbolDescriptionValueUnitsSource
Latent heat flux
PActual atmospheric pressureMeasured in the fieldkPaMeteorological station
RHmaxMaximum relative humidityMeasured in the field%Meteorological station
RHminMinimum relative humidityMeasured in the field%Meteorological station
TaAir temperatureMeasured in the field °CMeteorological station
TmaxMaximum air temperatureMeasured in the field °CMeteorological station
TminMinimum air temperatureMeasured in the field °CMeteorological station
TparPhotosynthetically active radiation transmittanceMeasured in the field-Sunscan
rbMean boundary layer resistance10~50 *s/mShuttleworth and Wallace [43]
rSTMean stomatal resistance60~240 *s/mShuttleworth and Wallace [43]
ρ0Air density at 0 oC and 0.1013 Mpa1.29kg/m3Allen et al. [50]
Sensible heat flux
hpPlant heightMeasured in the fieldmSteel tape
hrReference height2mShuttleworth and Wallace [43]
kVon Kármán’s constant0.41-Shuttleworth and Wallace [43]
nEddy diffusivity decay constant2.5-Shuttleworth and Wallace [43]
uWind speed at the reference heightMeasured in the fieldm/sMeteorological station
z0Roughness height0.01mShuttleworth and Wallace [43]
Soil heat flux
b1Empirical parameter0.243W/(m °C)Chung and Horton [46]
b2Empirical parameter0.393W/(m °C)Chung and Horton [46]
b3Empirical parameter1.534W/(m °C)Chung and Horton [46]
TlSoil temperature at Δz depthMeasured in the field°CSensors
ΔtTime step increment3600sSet in field experiments
ΔzDistance between the soil surface and the first vertical node0.1mSet in field experiments
θVolumetric soil water content at Δz depthMeasured in the fieldm3/m3Sensors
θnSolid phase content at Δz depth0.66-Set in field experiments
θoOrganic matter content at Δz depth0-Set in field experiments
Heat conduction
kaAir thermal conductivity0.025W/(m °C)Campbell and Norman ([48]
NudNusselt number reflecting regime of conduction1-Ham and Kluitenberg [6]
zgDistance between the plastic film mulch and soil0.001~0.01 *mHam and Kluitenberg [6]
* The value was calibrated with field experiments.
Table 3. The model parameters (shortwave transmittance τsm, shortwave absorptance αsm, emissivity εm, and longwave transmittance τlm of the plastic film mulch; emissivity εs and shortwave reflectance ρss of soil; mean stomatal resistance rST, mean boundary layer resistance rb; and distance between the plastic film mulch and soil zg) calibrated with hourly temperature, net radiation, and soil heat flux measured for the transparent plastic film mulch (TM) and black plastic film mulch (BM) in the potato field at Wuwei, China in 2014.
Table 3. The model parameters (shortwave transmittance τsm, shortwave absorptance αsm, emissivity εm, and longwave transmittance τlm of the plastic film mulch; emissivity εs and shortwave reflectance ρss of soil; mean stomatal resistance rST, mean boundary layer resistance rb; and distance between the plastic film mulch and soil zg) calibrated with hourly temperature, net radiation, and soil heat flux measured for the transparent plastic film mulch (TM) and black plastic film mulch (BM) in the potato field at Wuwei, China in 2014.
TreatmentτsmαsmεmτlmεsρssrST (s/m)rb (s/m)zg (m)
TM0.930.050.380.720.860.1790350.004
BM0.030.930.820.110.860.1790350.004
Table 4. The root mean square errors (RMSE), coefficient of determination (R2), and mean errors (ME) in simulated and observed hourly temperatures on the mulch surface and soil surface for transparent plastic film mulch and black plastic film mulch treatments at Wuwei, China.
Table 4. The root mean square errors (RMSE), coefficient of determination (R2), and mean errors (ME) in simulated and observed hourly temperatures on the mulch surface and soil surface for transparent plastic film mulch and black plastic film mulch treatments at Wuwei, China.
SeasonParametersTransparent Plastic Film MulchBlack Plastic Film Mulch
Mulch SurfaceSoil SurfaceMulch SurfaceSoil Surface
2014RMSE (°C)3.52.54.43.8
R20.800.790.840.82
ME (°C)2.2−0.94.03.6
2015RMSE (°C)3.83.54.12.7
R20.820.790.730.65
ME (°C)2.2−2.2−1.1−1.1
Table 5. The root mean square errors (RMSE) and coefficient of determination (R2) in simulated and observed daily temperatures (on mulch surface and soil surface) for the transparent plastic film mulch and black plastic film mulch treatments at Wuwei, China.
Table 5. The root mean square errors (RMSE) and coefficient of determination (R2) in simulated and observed daily temperatures (on mulch surface and soil surface) for the transparent plastic film mulch and black plastic film mulch treatments at Wuwei, China.
SeasonParametersTransparent Plastic Film MulchBlack Plastic Film Mulch
Mulch SurfaceSoil SurfaceMulch SurfaceSoil Surface
2014RMSE (°C)2.41.44.23.7
R20.830.840.760.86
2015RMSE (°C)2.62.81.71.5
R20.780.830.870.91
Table 6. The root mean square errors (RMSE), coefficient of determination (R2), and mean errors (ME) in simulated and observed hourly net radiation and soil heat flux for the transparent plastic film mulch and black plastic film mulch treatments at Wuwei, China.
Table 6. The root mean square errors (RMSE), coefficient of determination (R2), and mean errors (ME) in simulated and observed hourly net radiation and soil heat flux for the transparent plastic film mulch and black plastic film mulch treatments at Wuwei, China.
SeasonParametersTransparent Plastic Film MulchBlack Plastic Film Mulch
Net RadiationSoil Heat FluxNet RadiationSoil Heat Flux
2014RMSE (W/m2)66.4728.7364.9323.47
R20.980.890.970.67
ME (W/m2)20.84−11.3036.1420.22
2015RMSE (W/m2)71.8442.8383.8339.67
R20.920.630.920.56
ME (W/m2)12.96−18.4320.1223.35
Table 7. The root mean square errors (RMSE) and coefficient of determination (R2) in simulated and observed daily net radiation and soil heat flux for the transparent plastic film mulch and black plastic film mulch treatments at Wuwei, China.
Table 7. The root mean square errors (RMSE) and coefficient of determination (R2) in simulated and observed daily net radiation and soil heat flux for the transparent plastic film mulch and black plastic film mulch treatments at Wuwei, China.
SeasonParametersTransparent Plastic Film MulchBlack Plastic Film Mulch
Net RadiationSoil Heat FluxNet RadiationSoil Heat Flux
2014RMSE (W/m2)26.5717.0341.1320.63
R20.920.520.960.62
2015RMSE (W/m2)20.3023.4232.7925.27
R20.890.230.930.36

Share and Cite

MDPI and ACS Style

Zhang, Y.-L.; Wang, F.-X.; Shock, C.C.; Feng, S.-Y. Modeling the Interaction of Plastic Film Mulch and Potato Canopy Growth with Soil Heat Transport in a Semiarid Area. Agronomy 2020, 10, 190. https://doi.org/10.3390/agronomy10020190

AMA Style

Zhang Y-L, Wang F-X, Shock CC, Feng S-Y. Modeling the Interaction of Plastic Film Mulch and Potato Canopy Growth with Soil Heat Transport in a Semiarid Area. Agronomy. 2020; 10(2):190. https://doi.org/10.3390/agronomy10020190

Chicago/Turabian Style

Zhang, You-Liang, Feng-Xin Wang, Clinton C. Shock, and Shao-Yuan Feng. 2020. "Modeling the Interaction of Plastic Film Mulch and Potato Canopy Growth with Soil Heat Transport in a Semiarid Area" Agronomy 10, no. 2: 190. https://doi.org/10.3390/agronomy10020190

APA Style

Zhang, Y. -L., Wang, F. -X., Shock, C. C., & Feng, S. -Y. (2020). Modeling the Interaction of Plastic Film Mulch and Potato Canopy Growth with Soil Heat Transport in a Semiarid Area. Agronomy, 10(2), 190. https://doi.org/10.3390/agronomy10020190

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop