• Title/Summary/Keyword: estimation of damage

Search Result 930, Processing Time 0.034 seconds

Estimation of the Moisture Maximizing Rate based on the Moisture Inflow Direction : A Case Study of Typhoon Rusa in Gangneung Region (수분유입방향을 고려한 강릉지역 태풍 루사의 수분최대화비 산정)

  • Kim, Moon-Hyun;Jung, Il-Won;Im, Eun-Soon;Kwon, Won-Tae
    • Journal of Korea Water Resources Association
    • /
    • v.40 no.9
    • /
    • pp.697-707
    • /
    • 2007
  • In this study, we estimated the PMP(Probable Maximum Precipitation) and its transition in case of the typhoon Rusa which happened the biggest damage of all typhoons in the Korea. Specially, we analysed the moisture maximizing rate under the consideration of meteorological condition based on the orographic property when it hits in Gangneung region. The PMP is calculated by the rate of the maximum persisting 12 hours 1000 hPa dew points and representative persisting 12 hours 1000 hPa dew point. The former is influenced by the moisture inflow regions. These regions are determined by the surface wind direction, 850 hPa moisture flux and streamline, which are the critically different aspects compared to that of previous study. The latter is calculated using statistics program (FARD2002) provided by NIDP(National Institute for Disaster Prevention). In this program, the dew point is calculated by reappearance period 50-year frequency analysis from 5% of the level of significant when probability distribution type is applied extreme type I (Gumbel distribution) and parameter estimation method is used the Moment method. So this study indicated for small basin$(3.76km^2)$ the difference the PMP through new method and through existing result of established storm transposition and DAD(Depth-Area-Duration). Consequently, the moisture maximizing rate is calculated in the moisture inflow regions determined by meteorological fields is higher $0.20{\sim}0.40$ range than that of previous study. And the precipitation is increased $16{\sim}31%$ when this rate is applied for calculation.

Ecological Changes of Insect-damaged Pinus densiflora Stands in the Southern Temperate Forest Zone of Korea (I) (솔잎혹파리 피해적송림(被害赤松林)의 생태학적(生態学的) 연구(研究) (I))

  • Yim, Kyong Bin;Lee, Kyong Jae;Kim, Yong Shik
    • Journal of Korean Society of Forest Science
    • /
    • v.52 no.1
    • /
    • pp.58-71
    • /
    • 1981
  • Thecodiplosis japonesis is sweeping the Pinus densiflora forests from south-west to north-east direction, destroying almost all the aged large trees as well as even the young ones. The front line of infestation is moving slowly but ceaselessly norhwards as a long bottle front. Estimation is that more than 40 percent of the area of P. densiflora forest has been damaged already, however some individuals could escapes from the damage and contribute to restore the site to the previous vegetation composition. When the stands were attacked by this insect, the drastic openings of the upper story of tree canopy formed by exclusively P. densiflora are usually resulted and some environmental factors such as light, temperature, litter accumulation, soil moisture and offers were naturally modified. With these changes after insect invasion, as the time passes, phytosociologic changes of the vegetation are gradually proceeding. If we select the forest according to four categories concerning the history of the insect outbreak, namely, non-attacked (healthy forest), recently damaged (the outbreak occured about 1-2 years ago), severely damaged (occured 5-6 years ago), damage prolonged (occured 10 years ago) and restored (occured about 20 years ago), any directional changes of vegetation composition could be traced these in line with four progressive stages. To elucidate these changes, three survey districts; (1) "Gongju" where the damage was severe and it was outbroken in 1977, (2) "Buyeo" where damage prolonged and (3) "Gochang" as restored, were set, (See Tab. 1). All these were located in the south temperate forest zone which was delimited mainly due to the temporature factor and generally accepted without any opposition at present. In view of temperature, the amount and distribution of precipitation and various soil factor, the overall homogeneity of environmental conditions between survey districts might be accepted. However this did not mean that small changes of edaphic and topographic conditions and microclimates can induce any alteration of vegetation patterns. Again four survey plots were set in each district and inter plot distance was 3 to 4 km. And again four subplots were set within a survey plot. The size of a subplot was $10m{\times}10m$ for woody vegetation and $5m{\times}5m$ for ground cover vegetation which was less than 2 m high. The nested quadrat method was adopted. In sampling survey plots, the followings were taken into account: (1) Natural growth having more than 80 percent of crown density of upper canopy and more than 5 hectares of area. (2) Was not affected by both natural and artificial disturbances such as fire and thinning operation for the past three decades. (3) Lower than 500 m of altitude (4) Less than 20 degrees of slope, and (5) Northerly sited aspect. An intensive vegetation survey was undertaken during the summer of 1980. The vegetation was devided into 3 categories for sampling; the upper layer (dominated mainly by the pine trees), the middle layer composed by oak species and other broad-leaved trees as well as the pine, and the ground layer or the lower layer (shrubby form of woody plants). In this study our survey was concentrated on woody species only. For the vegetation analysis, calculated were values of intensity, frequency, covers, relative importance, species diversity, dominance and similarity and dissimilasity index when importance values were calculated, different relative weights as score were arbitrarily given to each layer, i.e., 3 points for the upper layer, 2 for the middle layer and 1 for the ground layer. Then the formula becomes as follows; $$R.I.V.=\frac{3(IV\;upper\;L.)+2(IV.\;middle\;L.)+1(IV.\;ground\;L.)}{6}$$ The values of Similarity Index were calculated on the basis of the Relative Importance Value of trees (sum of relative density, frequency and cover). The formula used is; $$S.I.=\frac{2C}{S_1+S_2}{\times}100=\frac{2C}{100+100}{\times}100=C(%)$$ Where: C = The sum of the lower of the two quantitative values for species shared by the two communities. $S_1$ = The sum of all values for the first community. $S_2$ = The sum of all values for the second community. In Tab. 3, the species composition of each plot by layer and by district is presented. Without exception, the species formed the upper layer of stands was Pinus densiflora. As seen from the table, the relative cover (%), density (number of tree per $500m^2$), the range of height and diameter at brest height and cone bearing tendency were given. For the middle layer, Quercus spp. (Q. aliena, serrata, mongolica, accutissina and variabilis) and Pinus densiflora were dominating ones. Genus Rhodedendron and Lespedeza were abundant in ground vegetation, but some oaks were involved also. (1) Gongju district The total of woody species appeared in this district was 26 and relative importance value of Pinus densiflora for the upper layer was 79.1%, but in the middle layer, the R.I.V. for Quercus acctissima, Pinus densiflora, and Quercus aliena, were 22.8%, 18.7% and 10.0%, respectively, and in ground vegetation Q. mongolica 17.0%, Q. serrata 16.8% Corylus heterophylla 11.8%, and Q. dentata 11.3% in order. (2) Buyeo district. The number of species enumerated in this district was 36 and the R.I.V. of Pinus densiflora for the uppper layer was 100%. In the middle layer, the R.I.V. of Q. variabilis and Q. serrata were 8.6% and 8.5% respectively. In the ground vegetative 24 species were counted which had no more than 5% of R.I.V. The mean R.I.V. of P.densiflora ( totaling three layers ) and averaging four plots was 57.7% in contrast to 46.9% for Gongju district. (3) Gochang-district The total number of woody species was 23 and the mean R.I.V. of Pinus densiflora was 66.0% showing greater value than those for two former districts. The next high value was 6.5% for Q. serrata. As the time passes since insect outbreak, the mean R.I.V. of P. densiflora increased as the following order, 46.9%, 57.7% and 66%. This implies that P. densiflora was getting back to its original dominat state again. The pooled importance of Genus Quercus was decreasing with the increase of that for Pinus densiflora. This trend was contradict to the facts which were surveyed at Kyonggi-do area (the central temperate forest zone) reported previously (Yim et al, 1980). Among Genus Quercus, Quercus acutissina, warm-loving species, was more abundant in the southern temperature zone to which the present research is concerned than the central temperate zone. But vice-versa was true with Q. mongolica, a cold-loving one. The species which are not common between the present survey and the previous report are Corpinus cordata, Beltala davurica, Wisturia floribunda, Weigela subsessilis, Gleditsia japonica var. koraiensis, Acer pseudosieboldianum, Euonymus japonica var. macrophylla, Ribes mandshuricum, Pyrus calleryana var. faruiei, Tilia amurensis and Pyrus pyrifolia. In Figure 4 and Table 5, Maximum species diversity (maximum H'), Species diversity (H') and Eveness (J') were presented. The Similarity indices between districts were shown in Tab. 5. Seeing Fig. 6, showing two-dimensional ordination of polts on the basis of X and Y coordinates, Ai plots aggregate at the left site, Bi plots at lower site, and Ci plots at upper-right site. The increasing and decreasing patterns as to Relative Density and Relative Importance Value by genus or species were given in Fig. 7. Some of the patterns presented here are not consistent with the previously reported ones (Yim, et al, 1980). The present authors would like to attribute this fact that two distinct types of the insect attack, one is the short war type occuring in the south temperate forest zone, which means that insect attack went for a few years only, the other one is a long-drawn was type observed at the temperate forest zone in which the insect damage went on continuously for several years. These different behaviours of infestation might have resulted the different ways of vegetational change. Analysing the similarity indices between districts, the very convincing results come out that the value of dissimilarity index between A and B was 30%, 27% between B and C and 35% between A and C (Table 6). The range of similarity index was obtained from the calculation of every possible combinations of plots between two districts. Longer time isolation between communities has brought the higher value of dissimilarity index. The main components of ground vegetation, 10 to 20 years after insect outbreak, become to be consisted of mainly Genus Lespedeza and Rhododendron. Genus Quercus which relate to the top dorminant state for a while after insect attack was giving its place to Pinus densiflora. It was implied that, provided that the soil fertility, soil moisture and soil depth were good enough, Genus Quercuss had never been so easily taken ever by the resistant speeies like Pinus densiflora which forms the edaphic climax at vast areas of forest land. Usually they refer Quercus to the representative component of the undisturbed natural forest in the central part of this country.

  • PDF

Detection of Wildfire Burned Areas in California Using Deep Learning and Landsat 8 Images (딥러닝과 Landsat 8 영상을 이용한 캘리포니아 산불 피해지 탐지)

  • Youngmin Seo;Youjeong Youn;Seoyeon Kim;Jonggu Kang;Yemin Jeong;Soyeon Choi;Yungyo Im;Yangwon Lee
    • Korean Journal of Remote Sensing
    • /
    • v.39 no.6_1
    • /
    • pp.1413-1425
    • /
    • 2023
  • The increasing frequency of wildfires due to climate change is causing extreme loss of life and property. They cause loss of vegetation and affect ecosystem changes depending on their intensity and occurrence. Ecosystem changes, in turn, affect wildfire occurrence, causing secondary damage. Thus, accurate estimation of the areas affected by wildfires is fundamental. Satellite remote sensing is used for forest fire detection because it can rapidly acquire topographic and meteorological information about the affected area after forest fires. In addition, deep learning algorithms such as convolutional neural networks (CNN) and transformer models show high performance for more accurate monitoring of fire-burnt regions. To date, the application of deep learning models has been limited, and there is a scarcity of reports providing quantitative performance evaluations for practical field utilization. Hence, this study emphasizes a comparative analysis, exploring performance enhancements achieved through both model selection and data design. This study examined deep learning models for detecting wildfire-damaged areas using Landsat 8 satellite images in California. Also, we conducted a comprehensive comparison and analysis of the detection performance of multiple models, such as U-Net and High-Resolution Network-Object Contextual Representation (HRNet-OCR). Wildfire-related spectral indices such as normalized difference vegetation index (NDVI) and normalized burn ratio (NBR) were used as input channels for the deep learning models to reflect the degree of vegetation cover and surface moisture content. As a result, the mean intersection over union (mIoU) was 0.831 for U-Net and 0.848 for HRNet-OCR, showing high segmentation performance. The inclusion of spectral indices alongside the base wavelength bands resulted in increased metric values for all combinations, affirming that the augmentation of input data with spectral indices contributes to the refinement of pixels. This study can be applied to other satellite images to build a recovery strategy for fire-burnt areas.

Freeze-thaw Resistance Estimation of Concrete using Surface Roughness and Image Analysis (콘크리트의 동결융해 저항성 추정을 위한 표면 거칠기 및 이미지 분석의 적용성)

  • Lee, Binna;Lee, Jong Suk
    • Journal of the Korea institute for structural maintenance and inspection
    • /
    • v.22 no.3
    • /
    • pp.1-7
    • /
    • 2018
  • As part of a research dedicated to the field evaluation of the durability of concrete subjected to freezing-thawing, this study analyzes the relationship between the surface roughness and the relative dynamic elastic modulus through image analysis. Four mix compositions with water-to-binder ratios (W/B) of 40%, 50%, 60% and 70% and without AE agent were considered to provoke early freezing. The basic physical properties of the mixes including the relative dynamic elastic modulus and the compressive strength were first evaluated experimentally according to W/B. Then, tests were performed to measure the surface roughness followed by photographs and SEM image analysis. The measured surface roughness tended to increase with larger number of freezing-thawing cycles regardless of W/B. The relative dynamic elastic modulus appeared to increase gradually with the number of cycles for the relatively denser mixes with W/B of 40% and 50%. Besides, the surface roughness increased only at rupture for the mixes with W/B of 60% and 70%. Moreover, the analysis of the photographs of the surface of the mixes with W/B of 40% and 50% revealed that the degradation progressed gradually from the surface with the freezing-thawing cycles. However, for the mixes with W/B of 60% and 70%, apparent change of the surface remained very insignificant until rupture at which damage like cracking could be observed. Consequently, the analysis of surface photograph or the measurement of the surface roughness presented some limitation in assessing the degree of freezing-thawing-induced degradation in case of relatively porous specimens. On the other hand, the photograph and surface roughness appeared to be sufficient for assessing such degradation for the mixes with W/B of 40% and 50%. Accordingly, the image of the surface and the surface roughness are potentially applicable on site for the assessment of freezing-thawing damages in relatively dense mixes.

A Study on the Effect of 2010 HNS Convention on Korean Industry (위험·유해물질 피해보상 국제협약의 우리 산업계에 대한 영향 고찰)

  • Kim, Ji-Hong
    • Journal of the Korean Society of Marine Environment & Safety
    • /
    • v.26 no.1
    • /
    • pp.57-64
    • /
    • 2020
  • The IOPC Fund general assembly reported that the International Convention on Liability and Compensation for Damage in Connection with the Carriage of Hazardous and Noxious Substances by Sea (known as the HNS Convention) will meet the requirement for the convention to take ef ect between 2021 and 2022. When the convention comes into effect, the liability-limit insurance of the HNS transport ship will be enforced and the shipper receiving the HNS will pay the share of the contribution from the International Fund for damages exceeding the limit of the ship's liability insurance. Korea is one of the major shipping and shipper countries in the world; thus, this study aimed to the need to analyze the effect of the convention on the related industries. The survey of ships and contribution targets analyzed the research data of the Ministry of Oceans and Fisheries. The P&I premium estimation was reviewed by the Korea Shipping Association and the K P&I as insured ships. In addition, the contribution of the HNS cargo volume was analyzed in an annual report by a representative international association for each cargo. About 1,500 ocean-going and domestic vessels have been identified as ships subject to the convention. The effect of changes in premiums under the convention was minimal for most ships. The effect of the shipping industry is expected, with about 150 domestic tankers expected to increase insurance premiums. In the case of shipper industries, 52 freight terminals were found to be eligible for the payment of the share of the international fund, as the proportion of freight volume in Korea was ranked second to fourth in the world by individual HNS accounts. This implies the obligation to pay contributions according to the convention. Considering the status of HNS transport ships entering and leaving ports and the quantity of HNS cargo, it can be concluded that the validity of Korea's convention is sufficient and that, it is necessary to coordinate with global major shipper countries.

Aggressive behavior of Male Rats following Hippocampal Ablation (뇌 해마를 떼어버린 흰쥐의 공격적 행동)

  • Park, Rho-Soon;Kim, Chul
    • The Korean Journal of Physiology
    • /
    • v.1 no.2
    • /
    • pp.169-175
    • /
    • 1967
  • An experiment was designed to see if the hippocampus exerts any influence upon the aggressive behavior of male rats. Fighting between rats was observed for the estimation of aggressiveness. Seventeen rats in which the hippocampus was almost totally removed through a small hole with a diameter around 3 mm made in the neocortex at the boundary between the parietal and occipital lobes (hippocampal group), 8 rats with similar neocortical damage alone (operated control group), and 17 normal control rats (normal group) were prepared and subjected to the experiment 3 months after the operation. Applying electric shock of short duration to the feet in a box with grid floor, a fight was provoked between an animal belonging to the hippocampal group and one belonging to the nor-mal group, between a rat of the hippocampal group and one of the operated control group, and also between a rat of the operated control group and one of the normal group. Three observers judged the performance of each animal independently and described it as winning, defeated, tied, or non-fighting. Fifteen shocked trials were administered to each pair of animals with around 2 minutes' interval between each trial. An animal received a 'judging score' of 3 when it won more frequently than was defeated, a judging score of 2 when it won as frequently as was defeated, when all fights were tied, or when no fighting occurred, while it received a judging score of 1 when it won less frequently than was defeated. Group differences in performances were analyzed in terms of judging score using Kolmogorov-Smirnov test for one sample. The results obtained were as follows: 1. In the fights between the hippocampal and the normal groups, the hippocampal animals made significantly better judging scores than the normal animals did (Table 1). 2. There was no significant difference between the hippocampal and the operated control group as to the judging scores they made in the course of fights between the two groups. However, the hippocampal animals tended to dominate over the operated control group as judged by comparing the total 'winning' of the former (30) to that of the latter (14) (Table 2). 3. The total judging score made by the operated control group in the course of the fights against the normal group was not significantly superior to that made by the normal group (Table 3). It was inferred from the above results that, though inconspicuous, the hippocampus tended to exert an inhibitory influence upon the aggressive behavior of male rats.

  • PDF

A Quantification Method for the Cold Pool Effect on Nocturnal Temperature in a Closed Catchment (폐쇄집수역의 냉기호 모의를 통한 일 최저기온 분포 추정)

  • Kim, Soo-Ock;Yun, Jin-I.
    • Korean Journal of Agricultural and Forest Meteorology
    • /
    • v.13 no.4
    • /
    • pp.176-184
    • /
    • 2011
  • Cold air on sloping surfaces flows down to the valley bottom in mountainous terrain at calm and clear nights. Based on the assumption that the cold air flow may be the same as the water flow, current models estimate temperature drop by regarding the cold air accumulation at a given location as the water-like free drainage. At a closed catchment whose outlet is blocked by man-made obstacles such as banks and roads, however, the water-like free drainage assumption is no longer valid because the cold air accumulates from the bottom first. We developed an empirical model to estimate quantitatively the effect of cold pool on nocturnal temperature in a closed catchment. In our model, a closed catchment is treated like a "vessel", and a digital elevation model (DEM) was used to calculate the maximum capacity of the cold pool formed in a closed catchment. We introduce a topographical variable named "shape factor", which is the ratio of the cold air accumulation potential across the whole catchment area to the maximum capacity of the cold pool to describe the relative size of temperature drop at a wider range of catchment shapes. The shape factor is then used to simulate the density profile of cold pool formed in a given catchment based on a hypsometric equation. The cold lake module was incorporated with the existing model (i.e., Chung et al., 2006), generating a new model and predicting distribution of minimum temperature over closed catchments. We applied this model to Akyang valley (i.e., a typical closed catchment of 53 $km^2$ area) in the southern skirt of Mt. Jiri National Park where 12 automated weather stations (AWS) are operational. The performance of the model was evaluated based on the feasibility of delineating the temperature pattern accurately at cold pool forming at night. Overall, the model's ability of simulating the spatial pattern of lower temperature were improved especially at the valley bottom, showing a similar pattern of the estimated temperature with that of thermal images obtained across the valley at dawn (0520 to 0600 local standard time) of 17 May 2011. Error in temperature estimation, calculated with the root mean square error using the 10 low-lying AWSs, was substantially decreased from $1.30^{\circ}C$ with the existing model to $0.71^{\circ}C$ with the new model. These results suggest the feasibility of the new method in predicting the site-specific freeze and frost warning at a closed catchment.

Estimation of fire Experiment Prediction by Utility Tunnels Fire Experiment and Simulation (지하공동구 화재 실험 및 시뮬레이션에 의한 화재 설칠 예측 평가)

  • 윤명오;고재선;박형주;박성은
    • Fire Science and Engineering
    • /
    • v.15 no.1
    • /
    • pp.23-33
    • /
    • 2001
  • The utility tunnels are the important facility as a mainstay of country because of the latest communication developments. However, the utilities tunnel is difficult to deal with in case of a fire accident. When a cable burns, the black smoke containing poisonous gas will be reduced. This black smoke goes into the tunnel, and makes it difficult to extinguish the fire. Therefore, when there was a fire in the utility tunnel, the central nerves of the country had been paralyzed, such as property damage, communication interruption, in addition to inconvenience for people. This paper is based on the fire occurred in the past, and reenacting the fire by making the real utilities tunnel model. The aim of this paper is the scientific analysis of the character image of the fire, and the verification of each fire protection system whether it works well after process of setting up a fire protection system in the utilities tunnel at a constant temperature. The fire experiment was equipped with the linear heat detector, the fire door, the connection water spray system and the ventilation system in the utilities tunnel. Fixed portion of an electric power supply cable was coated with a fire retardant coating, and a heating tube was covered with a fireproof. The result showed that the highest temperature was $932^{\circ}c$ and the linear heat detector was working at the constant temperature, and it pointed at the place of the fire on the receiving board, and Fixed portion of the electric power supply cable coated with the fire retardant coating did not work as the fireproof. The heating tube was covered with the fireproof about 30 minutes.

  • PDF

Extraction of Water Body Area using Micro Satellite SAR: A Case Study of the Daecheng Dam of South korea (초소형 SAR 위성을 활용한 수체면적 추출: 대청댐 유역 대상)

  • PARK, Jongsoo;KANG, Ki-Mook;HWANG, Eui-Ho
    • Journal of the Korean Association of Geographic Information Studies
    • /
    • v.24 no.4
    • /
    • pp.41-54
    • /
    • 2021
  • It is very essential to estimate the water body area using remote exploration for water resource management, analysis and prediction of water disaster damage. Hydrophysical detection using satellites has been mainly performed on large satellites equipped with optical and SAR sensors. However, due to the long repeat cycle, there is a limitation that timely utilization is impossible in the event of a disaster/disaster. With the recent active development of Micro satellites, it has served as an opportunity to overcome the limitations of time resolution centered on existing large satellites. The Micro satellites currently in active operation are ICEYE in Finland and Capella satellites in the United States, and are operated in the form of clusters for earth observation purposes. Due to clustering operation, it has a short revisit cycle and high resolution and has the advantage of being able to observe regardless of weather or day and night with the SAR sensor mounted. In this study, the operation status and characteristics of micro satellites were described, and the water area estimation technology optimized for micro SAR satellite images was applied to the Daecheong Dam basin on the Korean Peninsula. In addition, accuracy verification was performed based on the reference value of the water generated from the optical satellite Sentinel-2 satellite as a reference. In the case of the Capella satellite, the smallest difference in area was shown, and it was confirmed that all three images showed high correlation. Through the results of this study, it was confirmed that despite the low NESZ of Micro satellites, it is possible to estimate the water area, and it is believed that the limitations of water resource/water disaster monitoring using existing large SAR satellites can be overcome.

Development of Robotic Inspection System over Bridge Superstructure (교량 상판 하부 안전점검 로봇개발)

  • Nam Soon-Sung;Jang Jung-Whan;Yang Kyung-Taek
    • Proceedings of the Korean Institute Of Construction Engineering and Management
    • /
    • autumn
    • /
    • pp.180-185
    • /
    • 2003
  • The increase of traffic over a bridge has been emerged as one of the most severe problems in view of bridge maintenance, since the load effect caused by the vehicle passage over the bridge has brought out a long-term damage to bridge structure, and it is nearly impossible to maintain operational serviceability of bridge to user's satisfactory level without any concern on bridge maintenance at the phase of completion. Moreover, bridge maintenance operation should be performed by regular inspection over the bridge to prevent structural malfunction or unexpected accidents front breaking out by monitoring on cracks or deformations during service. Therefore, technical breakthrough related to this uninterested field of bridge maintenance leading the public to the turning point of recognition is desperately needed. This study has the aim of development on automated inspection system to lower surface of bridge superstructures to replace the conventional system of bridge inspection with the naked eye, where the monitoring staff is directly on board to refractive or other type of maintenance .vehicles, with which it is expected that we can solve the problems essentially where the results of inspection are varied to change with subjective manlier from monitoring staff, increase stabilities in safety during the inspection, and make contribution to construct data base by providing objective and quantitative data and materials through image processing method over data captured by cameras. By this system it is also expected that objective estimation over the right time of maintenance and reinforcement work will lead enormous decrease in maintenance cost.

  • PDF