• Journal of Korea Water Resources Association
• /
• v.55 no.12
• /
• pp.1041-1052
• /
• 2022

This study is to estimate the optimal ecological flow and analysis the spatial distribution of fish habitat for Andong dam downstream reach (4,565.7 km²) using PHABSIM (Physical Habiat Simulation System) and River2D. To establish habitat models, the cross-section informations and hydraulic input data were collected uisng the Nakdong river basic plan report. The establishment range of PHABSIM was set up about 410.0 m from Gudam streamflow gauging station (GD) and about 6.0 km including GD for River2D. To select representative fish species and construct HSI (Habitat Suitability Index), the fish survey was performed at Pungji bridge where showed well the physical characteristics of target stream located downstream of GD. As a result of the fish survey, Zacco platypus was showed highly relative abundance resulting in selecting as the representative fish species, and HSI was constructed using physical habitat characteristics of the Zacco platypus. The optimal range of HSI was 0.3~0.5 m/s at the velocity suitability index, 0.4~0.6 m at the depth suitability index, and the substrate was sand to fine gravel. As a result of estimating the optimal ecological flow by applying HSI to PHABSIM, the optimal ecological flow for target stream was 20.0 m³/sec. As a result of analysis two-dimensional spatial analysis of fish habitat using River2D, WUA (Weighted Usable Area) was estimated 107,392.0 m²/1000 m under the ecological flow condition and it showed the fish habitat was secured throughout the target stream compared with Q₃₅₅ condition.

• Korean Journal of Ichthyology
• /
• v.34 no.3
• /
• pp.160-171
• /
• 2022

This study investigated the characteristics of fish communities, habitat characteristics and the age of Lefua costata in the Jusucheon of Gangneung-si, Korea from January to December 2013. The survey collected 23 species belonging to ten families from 6 survey stations. The dominant and subdominant species were Zacco koreanus (relative abundance, 22.1%) and Rhynchocypris steindachneri (20.1%), respectively. The next most abundant species were Tridentiger brevispinis (10.7%), Orthrias nudus (9.7%), Zacco platypus (8.6%), Lefua costata (7.9%), Iksookimia koreensis (6.0%) and Pungitius sinensis (3.0%). Among the fish species collected, one species, P. sinensis, was class II endangered wildlife by the Ministry of Environment, and four species (Z. koreanus, I. koreensis, Silurus microdorsalis and Cottus koreanus), were endemic to Korea. Additionally, five migration fish species (Tribolodon hakonensis, Hypomesus nipponensis, Gasterosteus aculeatus, Oncorhynchus keta and Plecoglossus altivelis) were collected. The similarity index based on species composition and numbers separated fish communites in the Jusucheon according to main section, i.e., uppermost(St. 1), upper (St. 2~4), middle (St. 5), lower (St. 6). Ecosystem health of Jusucheon Stream (fish assessement index) evaluated using fish was assessed as very good (5 stations) and good (1 station). The main inhabit of L. costata was the middle-lower stream of aquatic plants habitats with mud bottoms, very slow water velocity and water depth 40~80 cm. The age groups for L. costata (female) estimated by the frequency distribution of total length in the spawning season (May) indicated that the 24~37 mm is 1-year old, the 38~51 mm group is 2-year old, the 52~63 mm is 3-year old, 64~77 mm is 4-year old and the 80~91 mm is more than 5-year old. Finally, characteristics of fish communities, habitat characteristics and the age of Lefua costata in the Jusucheon was discussed.

• Journal of the Korean Institute of Landscape Architecture
• /
• v.44 no.6
• /
• pp.84-97
• /
• 2016

This study was to investigate the user's thermal environments of the children's parks according to pavements and sunscreen types during periods of heat waves. The measurements were conducted at the sand pits, rubber chip pavement, shelters, and green shade ground of the two children's parks located in Jinju, Korea(Chilam: $N\;35^{\circ}11^{\prime}1.4{^{\prime}^{\prim}}$, $E\;128^{\circ}5^{\prime}31.7{^{\prime}^{\prime}}$, elevation 38m, Gaho: $N\;35^{\circ}09^{\prime}56.8{^{\prime}^{\prime}}$, $E\;128^{\circ}6^{\prime}41.1{^{\prime}^{\prime}}$, elevation 24m) over three days during 11-13, August, 2016. The highest ambient air temperatures at the Jinju Meteorological Office during the three measurement days were $35.9{\sim}36.8^{\circ}C$, which corresponded with the extremely hot weather. A series of experiments measured air temperature, relative humidity, wind velocity, black globe temperature, and long-wave and short-wave radiation of the six directions 0.6 m above ground level. The wet bulb globe temperature(WBGT) and the universal thermal climatic index(UTCI) were used to evaluate thermal stress. Surface temperature images of the play equipment were also taken using infrared thermography. Surface temperatures of the play equipment and grounds were used to evaluate burn risk through contact with playground materials. The results showed the following. The maximum air temperatures averaged over 1-hour period for three days were $36.6{\sim}39.4^{\circ}C$. The sun shades reduced those temperatures by up to $2.8^{\circ}C$(green shade) and $1.0^{\circ}C/2.3^{\circ}C$(shelters). The minimum relative humidity values averaged over 1-hour period for three days were 44~50%. The sun shades increased those humidity values by up to 6%(green shade) and 4%/6%(shelters). The risk of heat related illness at the measurement sites of the children's parks were extreme and high in the daytime hours. The maximum WBGT values averaged over a 30-minute period for three days were $31.2{\sim}33.6^{\circ}C$. The sun shades reduced those WBGT values by up to $2.4^{\circ}C$(green shade) and $0.5^{\circ}C/2.1^{\circ}C$(shelters) compared to sandpits, but would not block the risk of heat related illness in the daytime hours. The category of heat stress at the measurement sites of the children's parks were extreme and very strong in the daytime hours. The maximum UTCI values averaged over a 30-minute period for three days were $39.9{\sim}48.1^{\circ}C$. The sun shades reduced those UTCI values by up to $7.8^{\circ}C$(green shade) and $4.1^{\circ}C/8.2^{\circ}C$(shelters) compared to sandpits, but could not lower heat stress category from extreme and very strong to strong and moderate in the daytime hours. According to the burn threshold criteria when skin was in contact with playground materials, the maximum surface temperature of the stainless steels($70.8^{\circ}C$) surpassed three seconds $60^{\circ}C$ threshold for uncoated steel, that of the rubber chip($76.5^{\circ}C$) surpassed five seconds $74^{\circ}C$ threshold for the plastic, that of the plastic slide($68.5^{\circ}C$) and seats($71.0^{\circ}C$) surpassed the one min $60^{\circ}C$ threshold for plastic, respectively. The surface temperatures of shaded play equipment were lower approximately $20^{\circ}C$ than those of play equipment exposed to the sun. Therefore, sun shades can block the risk of burns in daytime hours. Because of the extreme and high risk of heat related illness and extreme and high heat stress at the children's parks during periods of heat waves, parents and administrators must protect children from the use of playgrounds. The risk of burn when contact with play equipments and grounds at the children's parks during periods of heat waves, was very high. The sun shades are essential to block the risk of burn from play equipments and grounds at the children's parks during heat waves.

• KOREAN JOURNAL OF CROP SCIENCE
• /
• v.34 no.s02
• /
• pp.45-65
• /
• 1989

The Korean peninsular having the complexity of the photography and variability of climate is located within passing area of a lots of typhoon occurring from the southern islands of Philippines. So, there are various patterns of wind and flooding damages in paddy field occuring by the strong wind and the heavy rain concentrated during the summer season of rice growing period in Korea. The wind damages to rice plants in Korea were mainly caused by saline wind, dry wind and strong wind when typhoon occurred. The saline wind damage having symptom of white head or dried leaves occurred by 1.1 to 17.2 mg of salt per dry weight stuck on the plant which was located at 2. 5km away from seashore of southern coastal area during the period(from 27th to 29th, August, 1986) of typhoon &Vera& accompanying 62-96％ of relative humidity, more than 6 m per second of wind velocity and 22.5 to 26.4$^{\circ}C$ of air temperature without rain. Most of the typhoons accompanying 4.0 to 8. 5m per second of wind and low humidity (lesp an 60％) with high temperature in the east coastal area and southen area of Korea. were changed to dry and hot wind by the foehn phenomenon. The dry wind damages with the symptom of the white head or the discolored brownish grain occurred at the rice heading stage. The strong wind caused the severe damages such as the broken leaves, cut-leaves and dried leaves before heading stage, lodging and shattering of grain at ripening stage mechanically during typhoon. To reduce the wind damages to rice plant, cultivation of resistant varieties to wind damages such as Sangpoongbyeo and Cheongcheongbyeo and the escape of heading stage during period of typhoon by accelerating of heading within 15th, August are effective. Though the flood disasters to rice plant such as earring away of field, burying of field, submerging and lodging damage are getting low by the construction of dam for multiple purpose and river bank, they are occasionally occurred by the regional heavy rain and water filled out in bank around the river. Paddy field were submerged for 2 to 4 days when typhoon and heavy rain occurred about the end of August. At this time, the rice plants that was in younger growing stage in the late transplanting field of southern area of Korea had the severe damages. Although panicles of rice plant which was in the meiotic growing stage and heading stage were died when flooded, they had 66％ of yield compensating ability by the upper tilling panicle produced from tiller with dead panicle in ordinary transplanting paddy field. It is effective for reduction of flooding damages to cultivate the resistant variety to flooding having the resistance to bacterial leaf blight, lodging and small brown planthopper simultaneously. Especially, Tongil type rice varieties are relatively resistant to flooding, compared to Japonica rice varieties. Tongil type rice varieties had high survivals, low elongation ability of leaf sheath and blade, high recovering ability by the high root activity and photosynthesis and high yield compensating ability by the upper tillering panicle when flooded. To minimize the flooding and wind damage to rice plants in future, following research have to be carried out; 1. Data analysis by telemetering and computerization of climate, actual conditions and growing diagnosis of crops damaged by disasters. 2. Development of tolerant varieties to poor natural conditions related to flooding and wind damages. 3. Improvement of the reasonable cropping system by introduction of other crops compensating the loss of the damaged rice. 4. Increament of utilization of rice plant which was damaged.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.16 no.2
• /
• pp.3361-3394
• /
• 1974

The purpose of this thesis is to disclose some characteristics of water consumption in relation to the quantities of dry matters through the growing period for two statured varieties of paddy rice which are a tall statured variety and a short one, including the water consumption during seedling period, and to find out the various coefficients of evapotranspiration that are applicable for the water use of an expected yield of the two varieties. PAL-TAL, a tall statured variety, and TONG-lL, a short statured variety were chosen for this investigation. Experiments were performed in two consecutive periods, a seedling period and a paddy field period, In the investigation of seedling period, rectangular galvanized iron evapotranspirometers (91cm${\times}$85cm${\times}$65cm) were set up in a way of two levels (PAL-TAL and TONG-lL varieties) with two replications. A standard fertilization method was applied to all plots. In the experiment of paddy field period, evapotanspiration and evaporation were measured separately. For PAL-TAL variety, the evapotranspiration measurements of 43 plots of rectangular galvanized iron evapotranspirometer (91cm${\times}$85cm${\times}$65cm) and the evaporation measurements of 25 plots of rectangular galvanized iron evaporimeter (91cm${\times}$85cm${\times}$15cm) have been taken for seven years (1966 through 1972), and for TONG-IL variety, the evapotranspiration measurements of 19 plots and the evaporation measurements of 12 plots have been collected for two years (1971 through 1972) with five different fertilization levels. The results obtained from this investigation are summarized as follows: 1. Seedling period 1) The pan evaporation and evapotranspiration during seedling period were proved to have a highly significant correlation to solar radiation, sun shine hours and relative humidity. But they had no significant correlation to average temperature, wind velocity and atmospheric pressure, and were appeared to be negatively correlative to average temperature and wind velocity, and positively correlative to the atmospheric pressure, in a certain period. There was the highest significant correlation between the evapotranspiration and the pan evaporation, beyond all other meteorological factors considered. 2) The evapotranpiration and its coefficient for PAL-TAL variety were 194.5mm and 0.94∼1.21(1.05 in average) respectively, while those for TONG-lL variety were 182.8mm and 0.90∼1.10(0.99 in average) respectively. This indicates that the evapotranspiration for TONG-IL variety was 6.2% less than that for PAL-TAL variety during a seedling period. 3) The evapotranspiration ratio (the ratio of the evapotranspiration to the weight of dry matters) during the seedling period was 599 in average for PAL-TAL variety and 643 for TONG-IL variety. Therefore the ratio for TONG-IL was larger by 44 than that for PAL-TAL variety. 4) The K-values of Blaney and Criddle formula for PAL-TAL variety were 0.78∼1.06 (0.92 in average) and for TONG-lL variety 0.75∼0.97 (0.86 in average). 5) The evapotranspiration coefficient and the K-value of B1aney and Criddle formular for both PAL-TAL and TONG-lL varieties showed a tendency to be increasing, but the evapotranspiration ratio decreasing, with the increase in the weight of dry matters. 2. Paddy field period 1) Correlation between the pan evaporation and the meteorological factors and that between the evapotranspiration and the meteorological factors during paddy field period were almost same as that in case of the seedling period (Ref. to table IV-4 and table IV-5). 2) The plant height, in the same level of the weight of dry matters, for PAL-TAL variety was much larger than that for TONG-IL variety, and also the number of tillers per hill for PAL-TAL variety showed a trend to be larger than that for TONG-IL variety from about 40 days after transplanting. 3) Although there was a tendency that peak of leaf-area-index for TONG-IL variety was a little retarded than that for PAL-TAL variety, it appeared about 60∼80 days after transplanting. The peaks of the evapotranspiration coefficient and the weight of dry matters at each growth stage were overlapped at about the same time and especially in the later stage of growth, the leaf-area-index, the evapotranspiration coefficient and the weight of dry matters for TONG-IL variety showed a tendency to be larger then those for PAL-TAL variety. 4) The evaporation coefficient at each growth stage for TONG-IL and PAL-TALvarieties was decreased and increased with the increase and decrease in the leaf-area-index, and the evaporation coefficient of TONG-IL variety had a little larger value than that of PAL-TAL variety. 5) Meteorological factors (especially pan evaporation) had a considerable influence to the evapotranspiration, the evaporation and the transpiration. Under the same meteorological conditions, the evapotranspiration (ET) showed a increasing logarithmic function of the weight of dry matters (x), while the evaporation (EV) a decreasing logarithmic function of the weight of dry matters; 800kg/10a x 2000kg/10a, ET=al+bl logl0x (bl>0) EV=a2+b2 log10x (a2>0 b2<0) At the base of the weight of total dry matters, the evapotranspiration and the evaporation for TONG-IL variety were larger as much as 0.3∼2.5% and 7.5∼8.3% respectively than those of PAL-TAL variety, while the transpiration for PAL-TAL variety was larger as much as 1.9∼2.4% than that for TONG-IL variety on the contrary. At the base of the weight of rough rices the evapotranspiration and the transpiration for TONG-IL variety were less as much as 3.5% and 8.l∼16.9% respectively than those for PAL-TAL variety and the evaporation for TONG-IL was much larger by 11.6∼14.8% than that for PAL-TAL variety. 6) The evapotranspiration coefficient, the evaporation coefficient and the transpiration coefficient and the transpiration coefficient were affected by the weight of dry matters much more than by the meteorological conditions. The evapotranspiratioa coefficient (ETC) and the evaporation coefficient (EVC) can be related to the weight of dry matters (x) by the following equations: 800kg/10a x 2000kg/10a, ETC=a3+b3 logl0x (b3>0) EVC=a4+b4 log10x (a4>0, b4>0) At the base of the weights of dry matters, 800kg/10a∼2000kg/10a, the evapotranspiration coefficients for TONG-IL variety were 0.968∼1.474 and those for PAL-TAL variety, 0.939∼1.470, the evaporation coefficients for TONG-IL variety were 0.504∼0.331 and those for PAL-TAL variety, 0.469∼0.308, and the transpiration coefficients for TONG-IL variety were 0.464∼1.143 and those for PAL-TAL variety, 0.470∼1.162. 7) The evapotranspiration ratio, the evaporation ratio (the ratio of the evaporation to the weight of dry matters) and the transpiration ratio were highly affected by the meteorological conditions. And under the same meteorological condition, both the evapotranspiration ratio (ETR) and the evaporation ratio (EVR) showed to be a decreasing logarithmic function of the weight of dry matters (x) as follows: 800kg/10a x 2000kg/10a, ETR=a5+b5 logl0x (a5>0, b5<0) EVR=a6+b6 log10x (a6>0 b6<0) In comparison between TONG-IL and PAL-TAL varieties, at the base of the pan evaporation of 343mm and the weight of dry matters of 800∼2000kg/10a, the evapotranspiration ratios for TONG-IL variety were 413∼247, while those for PAL-TAL variety, 404∼250, the evaporation ratios for TONG-IL variety were 197∼38 while those for PAL-TAL variety, 182∼34, and the transpiration ratios for TONG-IL variety were 216∼209 while those for PAL-TAL variety, 222∼216 (Ref. to table IV-23, table IV-25 and table IV-26) 8) The accumulative values of evapotranspiration intensity and transpiration intensity for both PAL-TAL and TONG-IL varieties were almost constant in every climatic year without the affection of the weight of dry matters. Furthermore the evapotranspiration intensity appeared to have more stable at each growth stage. The peaks of the evapotranspiration intensity and transpiration intensity, for both TONG-IL and PAL-TAL varieties, appeared about 60∼70 days after transplanting, and the peak value of the former was 128.8${\pm}$0.7, for TONG-IL variety while that for PAL-TAL variety, 122.8${\pm}$0.3, and the peak value of the latter was 152.2${\pm}$1.0 for TONG-IL variety while that for PAL-TAL variety, 152.7${\pm}$1.9 (Ref.to table IV-27 and table IV-28) 9) The K-value in Blaney & Criddle formula was changed considerably by the meteorological condition (pan evaporation) and related to be a increasing logarithmic function of the weight of dry matters (x) for both PAL-TAL and TONG-L varieties as follows; 800kg/10a x 2000kg/10a, K=a7+b7 logl0x (b7>0) The K-value for TONG-IL variety was a little larger than that for PAL-TAL variety. 10) The peak values of the evapotranspiration coefficient and k-value at each growth stage for both TONG-IL and PAL-TAL varieties showed up about 60∼70 days after transplanting. The peak values of the former at the base of the weights of total dry matters, 800∼2000kg/10a, were 1.14∼1.82 for TONG-IL variety and 1.12∼1.80, for PAL-TAL variety, and at the base of the weights of rough rices, 400∼1000 kg/10a, were 1.11∼1.79 for TONG-IL variety and 1.17∼1.85 for PAL-TAL variety. The peak values of the latter, at the base of the weights of total dry matters, 800∼2000kg/10a, were 0.83∼1.39 for TONG-IL variety and 0.86∼1.36 for PAL-TAL variety and at the base of the weights of rough rices, 400∼1000kg/10a, 0.85∼1.38 for TONG-IL variety and 0.87∼1.40 for PAL-TAL variety (Ref. to table IV-18 and table IV-32) 11) The reasonable and practicable methods that are applicable for calculating the evapotranspiration of paddy rice in our country are to be followed the following priority a) Using the evapotranspiration coefficients based on an expected yield (Ref. to table IV-13 and table IV-18 or Fig. IV-13). b) Making use of the combination method of seasonal evapotranspiration coefficient and evapotranspiration intensity (Ref. to table IV-13 and table IV-27) c) Adopting the combination method of evapotranspiration ratio and evapotranspiration intensity, under the conditions of paddy field having a higher level of expected yield (Ref. to table IV-23 and table IV-27). d) Applying the k-values calculated by Blaney-Criddle formula. only within the limits of the drought year having the pan evaporation of about 450mm during paddy field period as the design year (Ref. to table IV-32 or Fig. IV-22).

• Journal of Bio-Environment Control
• /
• v.5 no.2
• /
• pp.215-235
• /
• 1996

The thermal analysis by mathematical model simulation makes it possible to reasonably predict heating and/or cooling requirements of certain greenhouses located under various geographical and climatic environment. It is another advantages of model simulation technique to be able to make it possible to select appropriate heating system, to set up energy utilization strategy, to schedule seasonal crop pattern, as well as to determine new greenhouse ranges. In this study, the control pattern for greenhouse microclimate is categorized as cooling and heating. Dynamic model was adopted to simulate heating requirements and/or energy conservation effectiveness such as energy saving by night-time thermal curtain, estimation of Heating Degree-Hours(HDH), long time prediction of greenhouse thermal behavior, etc. On the other hand, the cooling effects of ventilation, shading, and pad ||||&|||| fan system were partly analyzed by static model. By the experimental work with small size model greenhouse of 1.2m$\times$2.4m, it was found that cooling the greenhouse by spraying cold water directly on greenhouse cover surface or by recirculating cold water through heat exchangers would be effective in greenhouse summer cooling. The mathematical model developed for greenhouse model simulation is highly applicable because it can reflects various climatic factors like temperature, humidity, beam and diffuse solar radiation, wind velocity, etc. This model was closely verified by various weather data obtained through long period greenhouse experiment. Most of the materials relating with greenhouse heating or cooling components were obtained from model greenhouse simulated mathematically by using typical year(1987) data of Jinju Gyeongnam. But some of the materials relating with greenhouse cooling was obtained by performing model experiments which include analyzing cooling effect of water sprayed directly on greenhouse roof surface. The results are summarized as follows : 1. The heating requirements of model greenhouse were highly related with the minimum temperature set for given greenhouse. The setting temperature at night-time is much more influential on heating energy requirement than that at day-time. Therefore It is highly recommended that night- time setting temperature should be carefully determined and controlled. 2. The HDH data obtained by conventional method were estimated on the basis of considerably long term average weather temperature together with the standard base temperature(usually 18.3$^{\circ}C$). This kind of data can merely be used as a relative comparison criteria about heating load, but is not applicable in the calculation of greenhouse heating requirements because of the limited consideration of climatic factors and inappropriate base temperature. By comparing the HDM data with the results of simulation, it is found that the heating system design by HDH data will probably overshoot the actual heating requirement. 3. The energy saving effect of night-time thermal curtain as well as estimated heating requirement is found to be sensitively related with weather condition: Thermal curtain adopted for simulation showed high effectiveness in energy saving which amounts to more than 50% of annual heating requirement. 4. The ventilation performances doting warm seasons are mainly influenced by air exchange rate even though there are some variations depending on greenhouse structural difference, weather and cropping conditions. For air exchanges above 1 volume per minute, the reduction rate of temperature rise on both types of considered greenhouse becomes modest with the additional increase of ventilation capacity. Therefore the desirable ventilation capacity is assumed to be 1 air change per minute, which is the recommended ventilation rate in common greenhouse. 5. In glass covered greenhouse with full production, under clear weather of 50% RH, and continuous 1 air change per minute, the temperature drop in 50% shaded greenhouse and pad ＆ fan systemed greenhouse is 2.6$^{\circ}C$ and.6.1$^{\circ}C$ respectively. The temperature in control greenhouse under continuous air change at this time was 36.6$^{\circ}C$ which was 5.3$^{\circ}C$ above ambient temperature. As a result the greenhouse temperature can be maintained 3$^{\circ}C$ below ambient temperature. But when RH is 80%, it was impossible to drop greenhouse temperature below ambient temperature because possible temperature reduction by pad ||||&|||| fan system at this time is not more than 2.4$^{\circ}C$. 6. During 3 months of hot summer season if the greenhouse is assumed to be cooled only when greenhouse temperature rise above 27$^{\circ}C$, the relationship between RH of ambient air and greenhouse temperature drop($\Delta$T) was formulated as follows : $\Delta$T= -0.077RH＋7.7 7. Time dependent cooling effects performed by operation of each or combination of ventilation, 50% shading, pad ＆ fan of 80% efficiency, were continuously predicted for one typical summer day long. When the greenhouse was cooled only by 1 air change per minute, greenhouse air temperature was 5$^{\circ}C$ above outdoor temperature. Either method alone can not drop greenhouse air temperature below outdoor temperature even under the fully cropped situations. But when both systems were operated together, greenhouse air temperature can be controlled to about 2.0-2.3$^{\circ}C$ below ambient temperature. 8. When the cool water of 6.5-8.5$^{\circ}C$ was sprayed on greenhouse roof surface with the water flow rate of 1.3 liter/min per unit greenhouse floor area, greenhouse air temperature could be dropped down to 16.5-18.$0^{\circ}C$, whlch is about 1$0^{\circ}C$ below the ambient temperature of 26.5-28.$0^{\circ}C$ at that time. The most important thing in cooling greenhouse air effectively with water spray may be obtaining plenty of cool water source like ground water itself or cold water produced by heat-pump. Future work is focused on not only analyzing the feasibility of heat pump operation but also finding the relationships between greenhouse air temperature(T$_{g}$ ), spraying water temperature(T$_{w}$ ), water flow rate(Q), and ambient temperature(T$_{o}$).