• Title/Summary/Keyword: Heating-cooling

Search Result 2,022, Processing Time 0.028 seconds

Effect of supercooling on the storage stability of rapidly frozen-thawed pork loins (과냉각 온도가 급속냉동-해동 처리된 돈육 등심의 저장성에 미치는 영향)

  • Choi, Eun Ji;Park, Hae Woong;Chung, Young Bae;Kim, Jin Se;Park, Seok Ho;Chun, Ho Hyun
    • Food Science and Preservation
    • /
    • v.24 no.2
    • /
    • pp.168-180
    • /
    • 2017
  • This study was performed to determine the rapid thawing method for reducing the thawing time of frozen pork loins and to examine the effects of supercooling on the microbiological, physicochemical, and sensory qualities of fresh and frozen-thawed pork during storage at -1.5, 4, and $15^{\circ}C$. Forced-air thawing at $4^{\circ}C$ was the most time-consuming process, whereas radio frequency thawing time was the shortest by dielectric heating. The supercooling storage temperature was chosen to be $-1.5^{\circ}C$ because microstructural damages were not observed in the pork sample after cooling at $-1.5^{\circ}C$ for 24 h. Fresh or frozen-thawed pork loins stored at $-1.5^{\circ}C$ had lower drip loss and total volatile base nitrogen, thiobarbituric acid-reactive substance, and Hunter b* levels than loins stored at 4 and $15^{\circ}C$. In addition, the least degree of increase in preexisting microorganisms counts of the fresh or frozen-thawed pork loin samples was obtained during supercooled storage at $-1.5^{\circ}C$. Sensory quality results of fresh and frozen-thawed pork loin samples stored at $-1.5^{\circ}C$ showed higher scores than the samples stored at 4 and $15^{\circ}C$. These data indicate that supercooling at $-1.5^{\circ}C$ in the meat processing industry would be effective for maintaining the quality of pork meats without ice crystal nucleation and formation.

Actual State of Structures and Environmental Control Facilities for Tomato Greenhouses in Chungnam Region (충남지역 토마토 재배온실의 구조와 환경조절설비 실태분석)

  • Nam, Sang-Woon;Kim, Young-Shik
    • Korean Journal of Agricultural Science
    • /
    • v.36 no.1
    • /
    • pp.73-85
    • /
    • 2009
  • An investigation was conducted to get the basic data for establishing structural safety and environmental management of tomato greenhouses in Chungnam region. The contents of the investigation consisted of actual state of greenhouse structures and environmental control facilities. Most of greenhouses were arch type single-span plastic houses and they had too low height for growing tomatoes. Frameworks of multi-span greenhouses were suitable, but those of single-span were mostly insufficient. Every greenhouse had thermal curtain movable or covering fixed inside the greenhouse for energy saving, and heating facilities were mostly warm air heater. Irrigation facilities were mostly drip tube and controlled by manual operation or timer. Almost all of the greenhouses didn't install high level of environmental control facilities such as ventilator, air circulation fan, $CO_2$ fertilizer, insect screen, supplemental light, and cooling device.

  • PDF

Strategies about Optimal Measurement Matrix of Environment Factors Inside Plastic Greenhouse (플라스틱온실 내부 환경 인자 다중센서 설치 위치 최적화 전략)

  • Lee, JungKyu;Kang, DongHyun;Oh, SangHoon;Lee, DongHoon
    • Journal of Bio-Environment Control
    • /
    • v.29 no.2
    • /
    • pp.161-170
    • /
    • 2020
  • There is systematic spatial variations in environmental properties due to sensitive reaction to external conditions at plastic greenhouse occupied 99.2% of domestic agricultural facilities. In order to construct 3 dimensional distribution of temperature, relative humidity, CO2 and illuminance, measurement matrix as 3 by 3 by 5 in direction of width, height and length, respectively, dividing indoor space of greenhouse was designed and tested at experimental site. Linear regression analysis was conducted to evaluate optimal estimation method in terms with horizontal and vertical variations. Even though sole measurement point for temperature and relative humidity could be feasible to assess indoor condition, multiple measurement matrix is inevitably required to improve spatial precision at certain time domain such as period of sunrise and sunset. In case with CO2, multiple measurement matrix could not successfully improve the spatial predictability during a whole experimental period. In case with illuminance, prediction performance was getting smaller after a time period of sunrise due to systematic interference such as indoor structure. Thus, multiple sensing methodology was proposed in direction of length at higher height than growing bed, which could compensate estimation error in spatial domain. Appropriate measurement matrix could be constructed considering the transition of stability in indoor environmental properties due to external variations. As a result, optimal measurement matrix should be carefully designed considering flexibility of construction relevant with the type of property, indoor structure, the purpose of crop and the period of growth. For an instance, partial cooling and heating system to save a consumption of energy supplement could be successfully accomplished by the deployment of multiple measurement matrix.

Time Resolution Improvement of MRI Temperature Monitoring Using Keyhole Method (Keyhole 방법을 이용한 MR 온도감시영상의 시간해상도 향상기법)

  • Han, Yong-Hee;Kim, Tae-Hyung;Chun, Song-I;Kim, Dong-Hyeuk;Lee, Kwang-Sig;Eun, Choong-Ki;Jun, Jae-Ryang;Mun, Chi-Woong
    • Investigative Magnetic Resonance Imaging
    • /
    • v.13 no.1
    • /
    • pp.31-39
    • /
    • 2009
  • Purpose : This study proposes the keyhole method in order to improve the time resolution of the proton resonance frequency(PRF) MR temperature monitoring technique. The values of Root Mean Square (RMS) error of measured temperature value and Signal-to-Noise Ratio(SNR) obtained from the keyhole and full phase encoded temperature images were compared. Materials and Methods : The PRF method combined with GRE sequence was used to get MR temperature images using a clinical 1.5T MR scanner. It was conducted on the tissue-mimic 2% agarose gel phantom and swine's hock tissue. A MR compatible coaxial slot antenna driven by microwave power generator at 2.45GHz was used to heat the object in the magnetic bore for 5 minutes followed by a sequential acquisition of MR raw data during 10 minutes of cooling period. The acquired raw data were transferred to PC after then the keyhole images were reconstructed by taking the central part of K-space data with 128, 64, 32 and 16 phase encoding lines while the remaining peripheral parts were taken from the 1st reference raw data. The RMS errors were compared with the 256 full encoded self-reference temperature image while the SNR values were compared with the zero filling images. Results : As phase encoding number at the center part on the keyhole temperature images decreased to 128, 64, 32 and 16, the RMS errors of the measured temperature increased to 0.538, 0.712, 0.768 and 0.845$^{\circ}C$, meanwhile SNR values were maintained as the phase encoding number of keyhole part is reduced. Conclusion : This study shows that the keyhole technique is successfully applied to temperature monitoring procedure to increases the temporal resolution by standardizing the matrix size, thus maintained the SNR values. In future, it is expected to implement the MR real time thermal imaging using keyhole method which is able to reduce the scan time with minimal thermal variations.

  • PDF

Development of Summer Leaf Vegetable Crop Energy Model for Rooftop Greenhouse (옥상온실에서의 여름철 엽채류 작물에너지 교환 모델 개발)

  • Cho, Jeong-Hwa;Lee, In-Bok;Lee, Sang-Yeon;Kim, Jun-Gyu;Decano, Cristina;Choi, Young-Bae;Lee, Min-Hyung;Jeong, Hyo-Hyeog;Jeong, Deuk-Young
    • Journal of Bio-Environment Control
    • /
    • v.31 no.3
    • /
    • pp.246-254
    • /
    • 2022
  • Domestic facility agriculture grows rapidly, such as modernization and large-scale. And the production scale increases significantly compared to the area, accounting for about 60% of the total agricultural production. Greenhouses require energy input to create an appropriate environment for stable mass production throughout the year, but the energy load per unit area is large because of low insulation properties. Through the rooftop greenhouse, one of the types of urban agriculture, energy that is not discarded or utilized in the building can be used in the rooftop greenhouse. And the cooling and heating load of the building can be reduced through optimal greenhouse operation. Dynamic energy analysis for various environmental conditions should be preceded for efficient operation of rooftop greenhouses, and about 40% of the solar energy introduced in the greenhouse is energy exchange for crops, so it should be considered essential. A major analysis is needed for each sensible heat and latent heat load by leaf surface temperature and evapotranspiration, dominant in energy flow. Therefore, an experiment was conducted in a rooftop greenhouse located at the Korea Institute of Machinery and Materials to analyze the energy exchange according to the growth stage of crops. A micro-meteorological and nutrient solution environment and growth survey were conducted around the crops. Finally, a regression model of leaf temperature and evapotranspiration according to the growth stage of leafy vegetables was developed, and using this, the dynamic energy model of the rooftop greenhouse considering heat transfer between crops and the surrounding air can be analyzed.

Comparing Net CO2 Uptake of Schlumbergera truncata 'Pink Dew' Phylloclades in a Growth Chamber and a Greenhouse (생육상과 온실에서 게발선인장 '핑크듀'의 엽상경별 CO2 흡수율 비교)

  • Seo Hee Jung;Ah Ram Cho;Yoon Jin Kim
    • Journal of Bio-Environment Control
    • /
    • v.32 no.1
    • /
    • pp.64-71
    • /
    • 2023
  • Crassulacean acid metabolism (CAM) plants use surplus CO2 generated by cooling and heating at night when ventilation is not needed in a greenhouse. Schlumbergera truncata 'Pink Dew' is a multi-flowering cactus that needs more phylloclades for high-quality production. This study examined photosynthetic characteristics by the phylloclade levels of S. truncata in a growth chamber and a greenhouse for use of night CO2 enrichment. The CO2 uptake rate of the S. truncata's top phylloclade in a growth chamber exhibited a C3 pattern, and the second phylloclade exhibited a C3-CAM pattern. The CO2 uptake rate of the top phylloclade in a greenhouse showed a negative value both day and night, but those of the second phylloclade exhibited a CAM pattern. The stomatal conductance and water-use efficiency (WUE) of S. truncata at both the top and second phylloclades were higher in a growth chamber than in a greenhouse. The WUE of S. truncata in a growth chamber and a greenhouse was higher at the second phylloclade, which is a CAM pattern compared with those of the top phylloclade. The daily total net CO2 uptake of S. truncata was higher in a growth chamber than in a greenhouse. The daily total net CO2 uptake of S. truncata at the second phylloclade had the highest value of 155 mmol·m-2·d-1 in a growth chamber. The night total CO2 uptake of S. truncate at the second phylloclade was 3-fold higher in a growth chamber than in a greenhouse. S. truncata's second phylloclade exhibited a CAM pattern that uptake CO2 at night, and the second phylloclade, was more mature than the top phylloclade. A multi-flowering cactus S. truncata 'Pink Dew' efficiently uptake night surplus CO2 in the proper environmental condition with matured phylloclade.

Study on the Controlling Mechaniques of the Environmental Factors in the Mushroom Growing House in Chonnam Province (전남 지방에 있어서의 양송이 재배에 최적한 환경조건 조절법 분석에 관한 연구)

  • Chung, Byung-Jae;Lee, Eun-Chol
    • Journal of the Korean Wood Science and Technology
    • /
    • v.2 no.2
    • /
    • pp.32-34
    • /
    • 1974
  • The important results which have been obtained in the investigation can be recapitulated as follows. 1. As demonstrated by the experimental results and analyses concerning their effects in the on-ground type mushroom house, the constructions in relation to the side wall and ceiling of the experimental house showed a sufficient heat insulation on effect to protect insides of the house from outside climatic conditions. 2. As the effect on the solar type experimental mushroom house which was constructed in a half basement has been shown by the experimental results and analyses, it has been proved to be effective for making use of solar heat. However there were found two problems to be improved for putting solar house to practical use in the farm mushroom growing: (1) the construction of the roof and ceiling should be the same as for the on ground type house, and (2) the solar heat generating system should be reconstructed properly. 3. Among several ventilation systems which have been studied in the experiments, the underground earthen pipe and ceiling ventilation, and vertical side wall and ceiling ventilation systems have been proved to be most effective for natural ventilation. 4. The experimental results have shown that ventilation systems such as the vertical side wall and underground ventilation systems are suitable to put to practical use as natural ventilation systems for farm mushroom house. These ventilation systems can remarkably improve the temperature of fresh air which is introduced into the house by heat transfers within the ventilation passages, so as to approach to the desired temperature of the house without any cooling or heating operation. For example, if it is assuming that X is the outside temperature and Y is the amount of temperature adjustment made by the influence of the ventilation system, the relationships that exist between X and Y can be expressed by the following regression lines. Underground iron pipe ventilation system. Y=0.9X-12.8 Underground earthen pipe ventilation system. Y=0.96X-15.11 Vertical side wall ventilation system. Y=0.94X-17.57 5. The experimental results have 8hown that the relationships existing between the admitted and expelled air and the $CO_2$ concentration can be described with experimental regression lines or an exponent equation as follows: 5.1 If it is assumed that X is an air speed cm/sec. and Y is an expelled air speed in cm/sec. in a natural ventilation system, since the Y is a function of the X, the relationships that exist between X and Y can be expressed by the regression lines shown below: 5.2 If it IS assumed that X is an admitted volume of air in $m^3$/hr. and Y is an expelled volume of air in $m^3$/hr. in a natural ventilation system, since the Y is a function of the X, the relationships that exist between X and Y can be expressed by the regression lines shown below. 5.3 If it is assumed that expelled air speed in emisec. and replacement air speed in cm/sec. at the bed surface in a natural ventilation system are shown as X and Y. respectively, since the Y is a function of the X. the relationships that exist between X and Y can be expressed by the following regression line: GE(100%)-CV (50%) ventilation system. Y=-0.54X+0.84 5.4 If it is assumed that the replacement air speed in cm/sec. at the bed surface is shown as X, and $CO_2$ concentration which is expressed by multiplying 1000 times the actual value of $CO_2$ % is shown as Y, in a natural ventilation system, since the Y is a function of the X, the relationships that exist between X and Y can be expressed by the following regression line: GE(100%)-CV(50%) ventilation system. Y=114.53-6.42X 5.5 If it is assumed that the expelled volume of air is shown as X and the $CO_2$ concencration which is expressed by multiplying 1000 times the actual of $CO_2$% is shown as Y in a natural ventilation system, since the Y is a function of the X, the relationships that exist between X and Y can be expressed by the following exponent equation: GE(100%)-CV(50%) ventilation system. Y=$127.18{\times}1.0093^{-x}$ 5.6 The experimental results have shown that the ratios of the cross sectional area of the GE and CV vent to the total cubic capacity of the house, required for providing an adequate amount of air in a natural ventilation system, can be estimated as follows: GE(admitting vent of the underground ventilation) 0.3-0.5% (controllable) CV(expelling vent of the ceiling ventilation) 0.8-1.0% (controllable) 6. Among several heating devices which were studied in the experiments, the hot-water boilor which wasmodified to be fitted both as hot-water boiler and as a pressureless steam-water was found most suitable for farm mushroom growing.

  • PDF

A STUDY ON THE IONOSPHERE AND THERMOSPHERE INTERACTION BASED ON NCAR-TIEGCM: DEPENDENCE OF THE INTERPLANETARY MAGNETIC FIELD (IMF) ON THE MOMENTUM FORCING IN THE HIGH-LATITUDE LOWER THERMOSPHERE (NCAR-TIEGCM을 이용한 이온권과 열권의 상호작용 연구: 행성간 자기장(IMF)에 따른 고위도 하부 열권의 운동량 강제에 대한 연구)

  • Kwak, Young-Sil;Richmond, Arthur D.;Ahn, Byung-Ho;Won, Young-In
    • Journal of Astronomy and Space Sciences
    • /
    • v.22 no.2
    • /
    • pp.147-174
    • /
    • 2005
  • To understand the physical processes that control the high-latitude lower thermospheric dynamics, we quantify the forces that are mainly responsible for maintaining the high-latitude lower thermospheric wind system with the aid of the National Center for Atmospheric Research Thermosphere-Ionosphere Electrodynamics General Circulation Model (NCAR-TIEGCM). Momentum forcing is statistically analyzed in magnetic coordinates, and its behavior with respect to the magnitude and orientation of the interplanetary magnetic field (IMF) is further examined. By subtracting the values with zero IMF from those with non-zero IMF, we obtained the difference winds and forces in the high-latitude 1ower thermosphere(<180 km). They show a simple structure over the polar cap and auroral regions for positive($B_y$ > 0.8|$\overline{B}_z$ |) or negative($B_y$ < -0.8|$\overline{B}_z$|) IMF-$\overline{B}_y$ conditions, with maximum values appearing around -80$^{\circ}$ magnetic latitude. Difference winds and difference forces for negative and positive $\overline{B}_y$ have an opposite sign and similar strength each other. For positive($B_z$ > 0.3125|$\overline{B}_y$|) or negative($B_z$ < -0.3125|$\overline{B}_y$|) IMF-$\overline{B}_z$ conditions the difference winds and difference forces are noted to subauroral latitudes. Difference winds and difference forces for negative $\overline{B}_z$ have an opposite sign to positive $\overline{B}_z$ condition. Those for negative $\overline{B}_z$ are stronger than those for positive indicating that negative $\overline{B}_z$ has a stronger effect on the winds and momentum forces than does positive $\overline{B}_z$ At higher altitudes(>125 km) the primary forces that determine the variations of tile neutral winds are the pressure gradient, Coriolis and rotational Pedersen ion drag forces; however, at various locations and times significant contributions can be made by the horizontal advection force. On the other hand, at lower altitudes(108-125 km) the pressure gradient, Coriolis and non-rotational Hall ion drag forces determine the variations of the neutral winds. At lower altitudes(<108 km) it tends to generate a geostrophic motion with the balance between the pressure gradient and Coriolis forces. The northward component of IMF By-dependent average momentum forces act more significantly on the neutral motion except for the ion drag. At lower altitudes(108-425 km) for negative IMF-$\overline{B}_y$ condition the ion drag force tends to generate a warm clockwise circulation with downward vertical motion associated with the adiabatic compress heating in the polar cap region. For positive IMF-$\overline{B}_y$ condition it tends to generate a cold anticlockwise circulation with upward vertical motion associated with the adiabatic expansion cooling in the polar cap region. For negative IMF-$\overline{B}_z$ the ion drag force tends to generate a cold anticlockwise circulation with upward vertical motion in the dawn sector. For positive IMF-$\overline{B}_z$ it tends to generate a warm clockwise circulation with downward vertical motion in the dawn sector.

The Influence of Ventilation and Shade on the Mean Radiant Temperature of Summer Outdoor (통풍과 차양이 하절기 옥외공간의 평균복사온도에 미치는 영향)

  • Lee, Chun-Seok;Ryu, Nam-Hyung
    • Journal of the Korean Institute of Landscape Architecture
    • /
    • v.40 no.5
    • /
    • pp.100-108
    • /
    • 2012
  • The purpose of the study was to evaluate the influence of shading and ventilation on Mean Radiant Temperature(MRT) of the outdoor space at a summer outdoor. The Wind Speed(WS), Air Temperature(AT) and Globe Temperature(GT) were recorded every minute from $1^{st}$ of May to the $30^{th}$ of September 2011 at a height of 1.2m above in four experimental plots with different shading and ventilating conditions, with a measuring system consisting of a vane type anemometer(Barini Design's BDTH), Resistance Temperature Detector(RTD, Pt-100), standard black globe(${\O}$ 150mm) and data acquisition systems(National Instrument's Labview and Compfile Techs' Moacon). To implement four different ventilating and shading conditions, three hexahedral steel frames, and one natural plot were established in the open grass field. Two of the steel frames had a dimension of $3m(W){\times}3m(L){\times}1.5m(H)$ and every vertical side covered with transparent polyethylene film to prevent lateral ventilation(Ventilation Blocking Plot: VP), and an additional shading curtain was applied on the top side of a frame(Shading and Ventilation Blocking Plot: SVP). The third was $1.5m(W){\times}1.5m(L){\times}1.5m(H)$, only the top side of which was covered by the shading curtain without the lateral film(Shading Plot: SP). The last plot was natural condition without any kind of shading and wind blocking material(Natural Open Plot: NP). Based on the 13,262 records of 44 sunny days, the time serial difference of AT and GT for 24 hour were analyzed and compared, and statistical analysis was done based on the 7,172 records of daytime period from 7 A.M. to 8 P.M., while the relation between the MRT and solar radiation and wind speed was analyzed based on the records of the hottest period from 11 A.M. to 4 P.M.. The major findings were as follows: 1. The peak AT was $40.8^{\circ}C$ at VP and $35.6^{\circ}C$ at SP showing the difference about $5^{\circ}C$, but the difference of average AT was very small within${\pm}1^{\circ}C$. 2. The difference of the peak GT was $12^{\circ}C$ showing $52.5^{\circ}C$ at VP and $40.6^{\circ}C$ at SP, while the gap of average GT between the two plots was $6^{\circ}C$. Comparing all four plots including NP and SVP, it can be said that the shading decrease $6^{\circ}C$ GT while the wind blocking increase $3^{\circ}C$ GT. 3. According to the calculated MRT, the shading has a cooling effect in reducing a maximum of $13^{\circ}C$ and average $9^{\circ}C$ MRT, while the wind blocking has heating effect of increasing average $3^{\circ}C$ MRT. In other words, the MRT of the shaded area with natural ventilation could be cooler than the wind blocking the sunny site to about $16^{\circ}C$ MRT maximum. 4. The regression and correlation tests showed that the shading is more important than the ventilation in reducing the MRT, while both of them do an important role in improving the outdoor thermal comfort. In summary, the results of this study showed that the shade is the first and the ventilation is the second important factor in terms of improving outdoor thermal comfort in summer daylight hours. Therefore, it can be apparently said that the more shade by the forest, shading trees etc., the more effective in conditioning the microclimate of an outdoor space reducing the useless or even harmful heat energy for human activities. Furthermore, the delicately designed wind corridor or outdoor ventilation system can improve even the thermal environment of urban area.

Application of LCA on Lettuce Cropping System by Bottom-up Methodology in Protected Cultivation (시설상추 농가를 대상으로 하는 bottom-up 방식 LCA 방법론의 농업적 적용)

  • Ryu, Jong-Hee;Kim, Kye-Hoon;Kim, Gun-Yeob;So, Kyu-Ho;Kang, Kee-Kyung
    • Korean Journal of Soil Science and Fertilizer
    • /
    • v.44 no.6
    • /
    • pp.1195-1206
    • /
    • 2011
  • This study was conducted to apply LCA (Life cycle assessment) methodology to lettuce (Lactuca sativa L.) production systems in Namyang-ju as a case study. Five lettuce growing farms with three different farming systems (two farms with organic farming system, one farm with a system without agricultural chemicals and two farms with conventional farming system) were selected at Namyangju city of Gyeonggi-province in Korea. The input data for LCA were collected by interviewing with the farmers. The system boundary was set at a cropping season without heating and cooling system for reducing uncertainties in data collection and calculation. Sensitivity analysis was carried out to find out the effect of type and amount of fertilizer and energy use on GHG (Greenhouse Gas) emission. The results of establishing GTG (Gate-to-Gate) inventory revealed that the quantity of fertilizer and energy input had the largest value in producing 1 kg lettuce, the amount of pesticide input the smallest. The amount of electricity input was the largest in all farms except farm 1 which purchased seedlings from outside. The quantity of direct field emission of $CO_2$, $CH_4$ and $N_2O$ from farm 1 to farm 5 were 6.79E-03 (farm 1), 8.10E-03 (farm 2), 1.82E-02 (farm 3), 7.51E-02 (farm 4) and 1.61E-02 (farm 5) kg $kg^{-1}$ lettuce, respectively. According to the result of LCI analysis focused on GHG, it was observed that $CO_2$ emission was 2.92E-01 (farm 1), 3.76E-01 (farm 2), 4.11E-01 (farm 3), 9.40E-01 (farm 4) and $5.37E-01kg\;CO_2\;kg^{-1}\;lettuce$ (farm 5), respectively. Carbon dioxide contribute to the most GHG emission. Carbon dioxide was mainly emitted in the process of energy production, which occupied 67~91% of $CO_2$ emission from every production process from 5 farms. Due to higher proportion of $CO_2$ emission from production of compound fertilizer in conventional crop system, conventional crop system had lower proportion of $CO_2$ emission from energy production than organic crop system did. With increasing inorganic fertilizer input, the process of lettuce cultivation covered higher proportion in $N_2O$ emission. Therefore, farms 1 and 2 covered 87% of total $N_2O$ emission; and farm 3 covered 64%. The carbon footprints from farm 1 to farm 5 were 3.40E-01 (farm 1), 4.31E-01 (farm 2), 5.32E-01 (farm 3), 1.08E+00 (farm 4) and 6.14E-01 (farm 5) kg $CO_2$-eq. $kg^{-1}$ lettuce, respectively. Results of sensitivity analysis revealed the soybean meal was the most sensitive among 4 types of fertilizer. The value of compound fertilizer was the least sensitive among every fertilizer imput. Electricity showed the largest sensitivity on $CO_2$ emission. However, the value of $N_2O$ variation was almost zero.