Effects of organic mulching on moisture and temperature of soil in greenhouse production of tomato under unheated greenhouse cultivation in the cold zone of China

Abstract In the cold zone of China, winter is cold and long and has a short duration of sunshine. Unheated earthen‐wall solar greenhouses are used for tomato production in winter in this region. This was an experimental investigation of different organic mulching materials (newspaper, bran, and grass) on the soil temperature, soil moisture, tomato yield, fruit quality, and water use efficiency. Organic mulching variously improved soil temperature, soil moisture, water use efficiency, and tomato yield, which is very important for greenhouse winter cultivation in this cold zone. Organic mulching regulated the soil temperature, with daily soil temperature ranges of bran, newspaper, and grass treatments being 1.6, 1.9, and 2.1°C lower than for bare land, respectively. Compared to bare land, newspaper mulching had little effect on soil temperature and fruit quality, but increased soil moisture (14.1%) and water use efficiency (WUE: WUE b , 31.3%; WUE y , 30.6%), and greatly increased yield (81.8%) and biomass (82.7%); bran mulching greatly increased soil temperature, moisture (16%), and WUE (WUE b , 60.1%; WUE y , 44.3%) and increased biomass (30.2%) and yield (17.3%); grass mulching greatly increased soil temperature and moisture (20.9%) and increased biomass (17.9%), yield (11.2%), and WUE (WUE b , 20.5%; WUE y , 13.6%). In addition, organic mulching had a good water retention effect on soil layer above 30 cm. The total soil water consumption during tomato growth was in the following order: newspaper (103 mm) > bare (74 mm) > grass (73 mm) > bran (60 mm). Soil water consumption mainly occurred in the 0‐ to 10‐cm soil layers.


| INTRODUC TI ON
Tomato (Solanum lycopersicum L.) is one of the most important vegetable plants in the world. Global production exceeds 180 million metric tons, with China and the USA as the leading producers in 2019, and production in China accounts for 20% of global output (FAOSTAT, 2021). Tomato is consumed fresh, cooked, or after processing by canning, or making juice, pulp, paste, or a variety of sauces (Zhang et al., 2016). Night temperature is the critical factor for tomato fruit setting, with the optimum range of 15-20°C. Fruit set is reduced markedly when average maximum temperature goes above 32°C and average minimum temperature goes above 21°C (Garg & Cheema, 2011). Wang et al. (2020) reported that daily water consumption ranged from 71.6 × 10 −3 to 91.9 × 10 −3 mm plant −1 under greenhouse soil cultivation.
In the cold zone (Beijing, Tianjin, Hebei, Shandong, Henan, Anhui, Shaanxi, Ningxia, southern Liaoning, east-central Gansu, southern Xinjiang, and southern Xizang) of China, winter is cold and long and with short durations of sunshine. The average temperature in January is between −1 and 5°C, the minimum temperature occasionally drops to −30°C, and the daily temperature range is 10-15°C (Wang, 2015). Greenhouses are the only facilities that allow mass production of vegetables in winter . Various engineering and horticultural measures are employed to provide an appropriate hydrothermal environment for vegetable growth in greenhouses in cold winter.
Engineering measures are mainly heating systems, such as hot water boilers, hot air furnaces, electric heaters, and heat pumps (Mellalou et al., 2021). However, added heating systems commensurately increase the capital cost and operating costs of greenhouses, and increase the release of greenhouse gases .
Unheated, single-slope solar greenhouses with a north wall (insulated and with high levels of sunlight) are mostly used to improve vegetable production in the northern part of China (Tong & Christopher, 2009;Wang et al., 2010). However, solar greenhouses are enclosed and are fragile production environments in winter due to it being colder outside with negligible ventilation. Soil evaporation and crop transpiration are the main pathways of water vapor and energy exchange between crops, soil, and atmosphere in greenhouses, which directly affect the hydrothermal and energy balance of soil and atmosphere inside greenhouses. In this fragile equilibrium environment, the change of one factor (especially soil) will inevitably lead to changes in other factors, which further affect crop traits (Li et al., 2013). The soil moisture influences the change in soil temperature, with higher soil moisture, there is less variation in soil temperature, and vice versa .
Mulching is the most efficient measure to regulate soil moisture and temperature (Kader et al., 2019), and plastic film mulches were widely used in greenhouses during winter in the cold zone of China.
However, most studies have reported that the residual plastic film accumulated in the soil damages the soil structure (Hu et al., 2020), soil physics (Wang et al., 2020), and soil hydraulic properties (Jiang et al., 2015). Meanwhile, more irrigation water or rainwater is retained in the topsoil layer due to the water-resistant layer formed by the contribution of residual plastic film, causing higher soil evaporation and shrinkage cracks in no-mulch soil surface (Wan et al., 2019).
Compared with plastic film mulches, organic material mulching has such advantages as accessibility , economy (Li et al., 2003), and environmental friendliness (Yin et al., 2016). In addition, some research has reported beneficial effects of organic mulching on soil temperature and moisture and crop yield. For example, Komariah et al. (2008) considered that organic mulches reduce heat conduction into the surface soil by retaining incoming solar radiation.
McMillen (2013) reported that soil moisture at a depth of 5-10 cm for wheat straw, grass clippings, and leaf debris mulching treatments was 10% higher than for bare soil. Kader, Senge, Mojid, and Nakamura (2017) mentioned that organic mulching had profound effects on microclimate, agronomic productivity, and yield of crops.
According to the above description, greenhouse cultivation is refined and sustainable agriculture, and proper cultivation management and techniques are very important, especially in the winter of cold zone of China. Meanwhile, selecting suitable mulching materials and methods to mulch the soil surface is crucial to cultivation in an unheated greenhouse when the daily average temperature is below zero. Refractory thick plastic film has been widely used in this area. To our knowledge, there are no previously published records and practice focusing on organic mulching materials for this cold environment.
It is envisaged that application of organic mulching and different organic mulches would modify the water uptake and physiological characteristics of tomato plants compared with traditional planting.
Field investigation is needed to explore the potential of this mulching system and different organic mulches, to enhance soil moisture availability, control soil temperature, and improve tomato physiological characteristics. The specific objectives of the study were to evaluate and compare soil temperature, soil moisture, soil water consumption, soil moisture extraction pattern (SMEP), water use efficiency (WUE), and tomato physiological characteristics of four mulching treatments (bare, newspaper, bran, and grass) under unheated greenhouse cultivation during winter of the cold zone of China.

| Experimental site
The experiment was conducted in an earthen-wall single-slope solar greenhouse (80 × 5.5 m, length and width) located in Yuci County, Shanxi Province, China (37°45′53″ N, 112°46′44″ E). Figure 1 shows the greenhouse cross-section. Heat preservation quilts were unrolled after sunset and rolled up after sunrise to preserve the inside temperature in the cold season. The experimental period was from 11 October 2020 to 23 February 2021. The temperature, relative humidity, and solar radiation inside and outside the greenhouse were measured by standard meteorological screen (RS-BYH-M, LvBo Instrument, Co. Ltd; Figure 2) for 10-s intervals during the experimental period.

| Experimental design and treatments
In the experiment, the middle position of greenhouse was selected to avoid errors due to environmental difference among plots. In order to reduce the effect of freezing plants far from earth wall in winter, the experiment used the plot near the earth wall ( Figure 1). Four plots were used with different organic mulching treatments. The large plot size for each mulching treatment was 2.1 m 2 (3.0 × 0.7 m) with a buffer zone of 0.7 m surrounding the treatments. To prevent the irrigation water from spilling out, 15 cm ridges were built around the plot. Each large plot was divided into three small plots (1.0 × 0.7 m). A no-mulching treatment was used as control (bare).
The permeability, water retention, and water loss of each treatment were measured before experiment (Table 1). The permeability rate of each treatment was measured by double ring infiltrometer (Touma & Albergel, 1992). The double rings were driven vertically into the ground to a depth of 10 cm, and the organic mulching material was then spread into the double rings. Ten centimeters of water F I G U R E 1 Cross-section of an earthen-wall unheated solar greenhouse.
F I G U R E 2 Field layout of the experiment and setup of various devices in the greenhouse.
was injected into the outer and the inner ring simultaneously, and subsequent injections were required to meet the benchmark of 10 cm. Then, 500 mL of water was added to the inner ring, and the time required for water infiltration to 10 cm was recorded. After repeating the process multiple times, water infiltration in the inner ring gradually slowed and began to stabilize. The water infiltration rate was calculated with the following formula: where V delineates the water infiltration rate for a certain duration (mm s −1 ). Q n refers to the cumulative amount of water injected within the nth measurement time (mL). S indicates the infiltration area of the inner ring (mm 2 ). T n is the interval among the nth measurements (s). The result showed that the initial (permeability rate before steady permeability) and steady (permeability rate after steady permeability) permeability rates among treatments did not significantly differ ( Table 1).
The water retention and saturated water loss of each mulching material were measured by the soaking-extracting method (Dong et al., 2021). The naturally air-dried organic mulches were placed in a special diamond mesh (Mesh diameter 0.5 mm) cage (0.2 × 0.1 × 0.1 m; Figure 3) Then, the diamond mesh cage was completely soaked in water for 24 h, suspended (air drying) for 24 h, and weighed. The calculation formulas are as follows: where NWC is the natural water content of mulching material (%), MWR is the maximum water retention of mulching material (%), AWR is the water retention of mulching material after irrigation (kg m −2 ), SWL is the water loss rate of mulching material after adequate water absorption (g kg −1 h −1 ), G n is the weight of the natural mulching material (kg), G d is the weight of mulching material after drying (kg), G a24 is the weight of mulching material after soaking for 24 h (kg), G l24 is the weight of mulching material after water loss for 24 h (kg), and S m is an area of natural mulching material (m 2 ).
The water content of natural organic mulching material was in the following order: grass (38.8%) > bran (24.5%) > newspaper (5.7%). The grass had the highest maximum water retention (410%), followed by newspaper (364%) and bran (236%) after soaking for 24 h. The highest actual water retention after adequate irrigation was for bran treatment (6.7 kg m −2 ), followed by grass (6.3 kg m −2 ) and newspaper (0.83 kg m −2 ). During the 24-h water loss experiment, the organic mulching water loss rate after adequate water absorption (soaking for 24 h) was in the following order: newspaper (137 g kg −1 h −1 ) > grass (91.3 g kg −1 h −1 ) > bran ( The hydrological characteristics of organic mulching material. Initial permeability rate (IPR) is the permeability rate before steady permeability. Note: Steady permeability rate (SPR) is the permeability rate after steady permeability. Natural water content (NWC) is the natural water content of mulching material. MWR is the maximum water retention of mulching material. Actual water retention (AWR) is the water retention of mulching material after irrigation. Saturated water loss (SWL) is the water loss rate of mulching material after adequate water absorption. Values followed by different lowercase letters (a-c) in the same column significantly differ (p < .05) based on Duncan's test.

F I G U R E 3
The schematic diagram of water retention and saturated water loss (a is a special diamond mesh cage; b is the schematic of the 5-cm bran; c is the schematic of the cage fitted with organic mulch immersed in water).

| Soil environment
The depth of the main root zone of the tomato crop was <30 cm from field observation. Undisturbed soil samples were collected from the tomato root zone: 0-10, 10-20, and 20-30 cm soil profile under bare land treatment. The soil texture of the study site was silty loam; soil bulk density was 1.4 g cm −3 in the 0-to 30-cm soil layer.
The field capacity and permanent wilting point moisture contents were 0.36 and 0.14 cm 3 cm −3 , respectively.
The soil temperature and absolute soil moisture content were

| Agronomic, physiological, and food chemistry management
Spray irrigation (simulated rainfall) was carried out randomly when plant dehydration or soil drying was observed.
Relative leaf chlorophyll (soil-plant analysis development: SPAD values) and relative leaf nitrogen (N) content (mg g −1 ) were measured once a week using a plant nutrition meter (LYS-4N, LvBo Instrument, Co. Ltd). There is an exact linear relationship between SPAD values and extracted leaf chlorophyll (Markwell et al., 1995;Yadava, 1986); therefore, this is used to represent leaf chlorophyll in this study.
Red fruits were harvested once a week during harvest season. The yield (g plant −1 ) on each plant was measured during the harvest season.
The WUE was calculated according to biomass (WUE b , g mm −1 : grams of biomass produced per mm of irrigation water) and yield (WUE y , g mm −1 : grams of yield produced per mm of irrigation water).
Fully mature tomato fruits were used for the fruit quality mea-

| Data analysis
Statistical analysis was performed using one-way analysis of variance, and the results were compared using Duncan's test with a confidence level of 5% using the R programming language.

| Meteorological parameters
Inside the greenhouse, the mean monthly air temperature during    Figure 5). For grass mulching treatment, the DMT at 10-cm depth was 0.32°C higher than for bare treatment, followed by 5 cm (0.31°C), 20 cm (0.17°C), and 30 cm (0.06°C). However, for newspaper mulching treatment, the DMT only at 5-cm depth was 0.15°C higher than for bare treatment; in other soil layers temperatures were 0.3-0.54°C lower than for bare treatment ( Figure 5). The difference in DMT at a 5-cm depth between the bare and mulching treatments (newspaper, bran, and grass) decreased through the tomato growth period (from October 2020 to February 2021); notably, newspaper and grass treatments had lower temperatures than the bare treatment from January ( Figure 6).

| Soil moisture and SWC
The soil moisture content in the tomato plots under the mulching treatments (newspaper, bran, and grass) fluctuated with TA B L E 3 The daily mean soil temperature and daily mean soil moisture at four depths (5, 10, 20, and 30 cm) in the experimental tomato plots under various mulching treatments.
The SWC of mulch treatments (newspaper, bran, and grass) during the tomato growth period showed that the SWC was low in the early part of the growth season (~60 DAT, days after transplants), but increased over time until tomato maturity ( Figure 8).

| SMEP
The SMEP among the treatments is shown in Table 4. For each treatment, the SMEP at 5-cm depth was significantly higher than for other soil layers, and decreased with increasing soil depth except for the bran treatment. The SMEP under bare, newspaper, and grass treatments did not significantly differ except at 10-cm depth where SMEP was slightly lower for bare than the other treatments.

| Biomass, yield, WUE, and fruit quality
Plant biomass, relative leaf chlorophyll, relative leaf N content, and yield under bare, bran, and grass treatments did not significantly differ, but these were significantly lower than newspaper treatments F I G U R E 7 Daily mean soil moisture (%) in the experimental tomato plots under various mulching treatments along with irrigation during the crop growing period.  (Table 5). The WUEs for plant biomass and yield significantly differed among all treatments. The WUE b and WUE y were in the following order: bran > newspaper > grass > bare (Table 5).
The fruit sugar content under bare treatment was significantly higher than for other treatments, and bran treatment had the lowest value. The sugar content among all treatments was in the following order: bare (9.3%) > grass (8.4%) > newspaper (7.6%) > bran (7.3%; Table 6). The values of fruit acid content among all treatments showed the same order as sugar content, but differences were not significant.

| DISCUSS ION
The mulching treatments (newspaper, bran, and grass) reduced the daily maximum soil temperature but raised the daily minimum soil temperature compared to the bare treatment. The reason might be because organic mulching (e.g., straw, grass, and wood chips) increases albedo and decreases thermal conductivity of the soil surface, reducing solar energy that penetrates the soil (Li et al., 2013), meanwhile reducing the soil energy dissipating, resulting in a significant reduction in daily soil temperature range.
Mulch acts as a buffer to insulate soil from heat and cold temperatures (Kader et al., 2019). Olasantan (1999) also reported that compared to nonmulched soil, mulched soil showed higher soil temperature during colder weather and lower soil temperature during warmer weather. In addition, Wang et al. (2021) summarized that organic mulch helps moderate temperature and protect soil and roots from extreme temperatures.
For DMT, bran and newspaper treatments significantly increased and decreased soil temperature compared to bare treatment, respectively, but grass treatment had no effect on soil temperature.
This may be because covering the soil surface with dense and poorly breathable bran seriously hinders the exchange of heat capacity between soil and atmosphere, thus maintaining soil temperature TA B L E 4 Soil moisture extraction pattern at 5-, 10-, 20-, and 30-cm depths over the 0-to 30-cm soil profile during the growing period under four different treatments. Note: Data are the means ± standard deviations. Values followed by different lowercase letters (a-c) in the same column significantly differ (P < .05) based on the Duncan's test.

TA B L E 6
Effect of different mulching on fruit quality. (Kodešova et al., 2014). However, the smoothness and lighter color of newspaper reflect sunlight (Haapala et al., 2014;. The heat was exchanged quickly between soil covered by fluffy grass and the atmosphere, and soil heat capacity was lost more easily for grass mulching in a cooler external environment (Al-Shammary et al., 2020).
The difference in DMT between bare and mulching treatments at 5-cm depth (decreased through the tomato growth period, Figure 6) may be because mulching warms the topsoil quickly during the day and dissipates heat slowly at night compared to bare land under direct sunlight (Zhou et al., 2009). However, as the canopy grows, the ground receives progressively less sunlight and the heat harbored by the mulch layer does not transfer to the soil (Dai et al., 2021). In addition, organic mulch degrades on the soil surface and this reduces the action of the barrier during later stages of plant growth.
The DMSM under organic mulching treatments was significantly higher than bare treatment, and grass treatment had the highest value (20.8%), followed by bran (16%) and newspaper (14.1%). The SWC for bran and grass treatments did not significantly differ from that for bare treatment but was significantly higher for newspaper than bare treatment. This may be because mulching alters the moisture and temperature environment of the soil by affecting the hydrothermal regime . Each type of mulching material has a particular set of characteristics. Mulching materials affect the soil moisture environment mainly through their water retention, water permeability, and decomposability. In this study, the influence of water permeability of mulching materials on soil moisture can be ignored due to the flood irrigation method. The DMSM contents of mulching treatments were higher than for bare treatment. The reason may be that mulching physically prevents vertical evaporation from the soil surface by blocking direct sunlight reaching the soil surface, and mulching measures also reduce water runoff on the soil surface (Saglam et al., 2017), thus preserving soil water. In addition, some organic mulch acts as a sponge and retains rainfall and irrigation water, thus preventing runoff and providing water at the time of crop requirement . Zhao et al. (2014) also reported that mulching treatments stored more soil moisture compared to bare soil.
Among the mulching treatments (newspaper, bran, and grass), the newspaper had the lowest DMSM and highest SWC.
Newspaper mulching with its smooth surface, low porosity, negligible actual water retention (0.83 kg m −2 ), and low water loss characteristics (137 g kg −1 h −1 ) greatly affected the soil and resulted in high evaporation efficiency after irrigation (Haapala et al., 2014).
The dynamic of hourly mean soil moisture among the treatments in the earlier ( Figure 9A) and later ( Figure 9B) growth stages after

F I G U R E 9
The dynamic of hourly mean soil moisture among the treatments in the earlier (3-8 November 2020, a) and later (2-7 February 2021, b) growth stages after irrigation. irrigation further substantiates these findings. Soil water under newspaper and bare treatments was rapidly lost in the earlier growth stages compared to grass and bran treatments; however, there was no difference in the rate of water loss among treatments in the later growth stages. This is because water loss at the earlier stage of plant growth is mainly by surface evaporation, yet the smooth and lighter color of the newspaper surface speeds up evaporation of irrigation water. However, with plant canopy growth, water loss is mainly dominated by plant consumption because sunlight does not reach the ground directly. In general, with plant growth and newspaper degradation, the function of newspaper on temperature preservation and water insulation gradually decreases. However, the organic matter produced by the degradation of newspapers is good fertilizer for later growth of plants.
In this study, the coverage of newspaper was two layers, and the actual water retention after irrigation was only 0.83 kg m −2 , far below that for bran (6.7 kg m −2 ) and grass (6.3 kg m −2 ). Increasing the coverage layers of newspaper would theoretically increase the soil temperature, water retention of newspaper, and organic matter production, but would reduce water permeability.
Therefore, appropriate coverage of newspaper requires further experimentation.
Bran had great water retention (6.7 kg m −2 ) and slow water loss characteristics (35.2 g kg −1 h −1 ) due to its inherently fine porous structure, which results in irrigation water not only reaching the soil surface but also some being retained in the bran layer. As the soil water infiltrates and is absorbed by the crops, the water in bran will be slowly released into the soil and to maintain the soil moisture. In addition, when irrigation water reaches the soil surface, some infiltrates and some evaporates. Yuan et al. (2009) also reported that porous materials such as camellia shell, straw, and straw and ecological pad, and slow water vapor transmission to the atmosphere.
However, evaporated water vapor reaches the mulch and then falls back to the soil as droplets . In particular, the bran layer with fine pores seriously hinders soil water evaporation.
Meanwhile, compared with newspaper and grass treatments, after irrigation, the bran treatment more easily forms a crust due to temperature increase from direct sunlight, which helps to preserve the water in the bran layer as well as the soil moisture.
Grass mulching greatly increased DMSM, consistent with the conclusion that mulching increased soil moisture. However, the SWC was close to bare treatment and much higher than that of bran treatment. This may be due to fluffy grass mulch macropores that block direct sunlight reaching the soil surface, but not preventing the gas exchange between soil and atmosphere. Water contained in grass was readily transformed into vapor to satisfy evaporative demand from the soil and grass following the energy balance at the soil surface. In this process, most of the water vapor is lost through evaporation although some evaporated water vapor reaches the mulch and then falls back. However, grass mulching treatment resulted in the highest SWC during the whole growth period. This may be for two reasons, grass mulching blocks direct sunlight reaching the soil surface, resulting in reduced evaporation and a large proportion of evaporated water vapor condenses into grass, then falls back to the soil surface and evaporates again after irrigation, this repeated water cycle delays soil moisture loss until the next irrigation.
The SMEP mainly focuses on the soil layer of 0-10 cm, especially 0-5 cm. This result is consistent with the findings of Kader, Senge, Mojid, Onishi, and Ito (2017). Soil moisture in the upper surface layer (0-10 cm) was highly dynamic due to water vapor fluxes across the soil-atmosphere interface, thus increasing the evaporation loss at the surface layer (Bittellia et al., 2008), and most of the water consumed by plants is in this zone due to high root mass (Lenka et al., 2009). However, the bran treatment at 5-cm depth showed significantly lowest SMEP among treatments. This is because the thickness of bran cover was about 5 cm, and the soft, thick layer of small bran particles had poor porosity, thus hindering the water vapor exchange between atmosphere and soil surface. Especially, irrigation water filled the pores between bran particles (high water retention) after irrigation, which almost cut off the circulation channel between soil and atmosphere, and so reducing the surface soil moisture evaporation loss. In addition, crusting of the bran layer surface further seriously affected water vapor circulation between soil surface and atmosphere.
Plant growth is a very complex process, and is affected by many factors such as external environment (e.g., soil, water, and climate; Chakma et al., 2021), plant differences (e.g., variety, disease resistance, and growth phenology; , and cultivation environment (e.g., field cultivation, greenhouse cultivation, and hydroponics; Wang et al., 2020;Zhang et al., 2016). In this study, the tomato biomass, relative leaf  et al. (2015).
Tomato fruit quality is susceptible to various biotic and abiotic stresses, with drought and salinity being major abiotic stresses that can positively affect fruit quality (Chakma et al., 2021). In this study, the bare treatment had the highest fruit quality, possibly due to the lowest root zone moisture during the fruit ripening period. Water deficit decreases the photo-assimilate dilution by decreasing the water level within fruits, thus promoting the accumulation of assimilates within fruits and improving quality parameters (Chakma et al., 2021;Kuscu et al., 2014). In addition, root zone temperature influences transport from source to sink during fruit development and the biochemical transformation during fruit ripening. High temperature increases the fruit's sugar content while decreasing the acid content (Mesa et al., 2022). Meanwhile, large diurnal temperature difference increases tomato fruit's sugar contents while organic acid contents are diminished (Ruiz-Nieves et al., 2020). In this study, the bare treatment had highest maximum temperature in the root zone and largest diurnal temperature difference. The superimposed effects of the above multiple factors promoted the best fruit quality in the bare treatment.

| CON CLUS ION
Overall, organic mulching tended to improve soil temperature, soil moisture, WUE, and tomato yield, which is very important for greenhouse cultivation in winter in the cold zone of China. In particular, the water retention of organic mulching was clear, and this influence extended to a depth of 30 cm.
Based on the study results, the three organic mulching materials had advantages and disadvantages. Newspaper mulching showed poor soil temperature and water protection, good biomass and yield, and good fruit quality. Bran mulching showed good soil temperature and water protection, middle biomass and yield, and poor fruit quality. Grass mulching had middle soil temperature and water protection, poor biomass and yield, and middle fruit quality. Newspaper mulching had no effect on tomato fruit quality, and bran and grass mulching reduced fruit quality. In addition, mulching regulated the soil temperature; the daily soil temperature ranges of bran, newspaper, and grass were 1.6, 1.9, and 2.1°C lower than for bare land, respectively.
However, the amount of organic mulch used in this study may not have reached the necessary threshold, and further research is needed to determine the threshold of organic mulch suitable for greenhouse cultivation in the cold zone in winter. writing -review and editing (equal).

ACK N O WLE D G E M ENTS
All the institutions and individuals who provided support for the paper are listed as funding and authors.

CO N FLI C T O F I NTER E S T S TATEM ENT
The authors declare that there are no conflicts of interest regarding the publication of this paper.

DATA AVA I L A B I L I T Y S TAT E M E N T
All data generated or analyzed during this study are included in this published article.

E TH I C S S TATEM ENT
Not applicable.

CO N S ENT TO PA RTI CI PATE
Not applicable.

CO N S E NT FO R PU B LI C ATI O N
Not applicable.

CO D E AVA I L A B I LIT Y
Not applicable.