Characteristic Air Temperature Distributions Observed in Summer and Winter in Urban Area in Japan

The air temperature distributions in August (summer) and December (winter) were measured in an approximately 15 × 15-km urban area in Hyogo Prefecture, Japan, in order to study the spatial distribution of the air temperature and to propose effective measures against the heat island phenomenon. The air temperature was measured mainly by using thermometer shelters installed in an elementary school and a junior high school. The characteristic air temperature distribution depended on the season. The air temperature was higher inland than in the coastal region in August but was higher in the coastal region in December. The air temperature index indicated that the area where higher air temperatures would most likely appear was 5 to 10 km inland from the coast in August and around the coast in December. The seasonal air temperature distribution was presumably due to the strength of solar radiation and anthropogenic exhaust heat.     

Land use change and landslide characteristics analysis for community-based disaster mitigation

On August 8, 2009, Typhoon Morakot brought heavy rain to Taiwan, causing numerous landslides and debris flows in the Taihe village area of Meishan Township, Chiayi County, in south-central Taiwan. In the Taihe land is primary used for agriculture and land use management may be a factor in the area’s landslides. This study explores Typhoon Morakot-induced landslides and land use changes between 1999 and 2009 using GIS with the aid of field investigation. Spot 5 satellite images with a resolution of 2.5 m are used for landslide interpretation and manually digitalized in GIS. A statistical analysis for landslide frequency–area distribution was used to identify the landslide characteristics associated with different types of land use. There were 243 landslides with a total area of 2.75 km² in the study area. The area is located in intrinsically fragile combinations of sandstone and shale. Typhoon Morakot-induced landslides show a power–law distribution in the study area. Landslides were mainly located in steep slope areas containing natural forest and in areas planted with bamboo, tea, and betel nut. Land covered with natural forest shows the highest landslide ratio, followed by bamboo, betel nut, and tea. Landslides thus show a higher ratio in areas planted with shallow root vegetation such as bamboo, betel nut, and tea. Furthermore, the degree of basin development is proportional to the landslide ratio. The results show that a change in vegetation cover results in a modified landslide area and frequency and changed land use areas have higher landslide ratios than non-changed. Land use management and community-based disaster prevention are needed in mountainous areas of Taiwan for hazard mitigation.     

锦州湾海域生态系统健康状况评价

基于结构功能指标体系评价法,利用层次分析法构建锦州湾海域生态系统健康评价指标体系,并选择2005—2009年8月生态调查数据对其进行初步的定量评价。结果表明,可以通过海洋化学指标、生态学指标和社会经济学指标三大类33个指标来建立锦州湾生态系统健康评价指标体系;2005—2009年8月锦州湾生态系统的整合健康指数依次分别为1.01、1.90、2.00、1.45和1.10,均低于系统健康指数值2,锦州湾海域生态系统连续5年一直处于亚健康状态,且从目前海域的污染现状来看,这一海域的水质、沉积物和生物体污染均比较严重,要修复其生态功能需要投入巨大的物质和能量。     

Columnar aerosol optical and radiative properties according to season and air mass transport pattern over East Asia

The column-integrated optical and radiative properties of aerosols in the downwind area of East Asia were investigated based on sun/sky radiometer measurements performed from February 2004 to June 2005 at Gwangju (35.23° N, 126.84° E) and Anmyeon (36.54° N, 126.33° E), Korea. The observed aerosol data were analyzed for differences among three seasons: spring (March-May), summer (June-August), and autumn/winter (September-February). The data were also categorized into five types depending on the air mass origin in arriving in the measurement sites: (a) from a northerly direction in spring (S(N)), (b) from a westerly direction in spring (S(W)), (c) cases with a low ?ngstr?m exponent (<0.8) in spring (dust), (d) from a northerly direction in autumn/winter (AW(N)), and (e) from a westerly direction during other seasons (AW(W)). The highest ?ngstr?m exponents (α) at Gwangju and Anmyeon were 1.43?±?0.30 and 1.49?±?0.20, respectively, observed in summer. The lowest column-mean single-scattering albedo (ω) at 440 nm observed at Gwangju and Anmyeon were 0.89?±?0.02 and 0.88?±?0.02, respectively, during a period marked by the advection of dust from the Asian continent. The highest ω values at Gwangju and Anmyeon were 0.95?±?0.02 and 0.96?±?0.02, respectively, observed in summer. Variations in the aerosol radiative-forcing efficiency (β) were related to the conditions of the air mass origin. The forcing efficiency in summer was -131.7 and -125.6 W?m(-2) at the surface in Gwangju and Anmyeon, respectively. These values are lower than those under the atmospheric conditions of spring and autumn/winter. The highest forcing efficiencies in autumn/winter were -214.3 and -255.9 W?m(-2) at the surface in Gwangju and Anmyeon, respectively, when the air mass was transported from westerly directions.     

Pool of mobile and immobile phosphorus in sediments of the large, shallow Lake Peipsi over the last 100?years

Temporal variations in sediment phosphorus (P) composition and mobility were estimated in surface sediments of accumulation (core PS509 44?cm) and erosion (core PS2009 30?cm) areas of the shallow, large Lake Peipsi sensu stricto (s.s.; 2611?km(2), unregulated water level). The P pool in sediments including buried and mobile P is evaluated for the first time, which will provide baseline data for the future modelling of internal loading in L. Peipsi. Five sedimentary P fractions (including inorganic and nonreactive P) were separately quantified: loosely adsorbed and pore-water P (NaCl-P); redox-sensitive fraction P (NaBD-P); P bound to oxides of non-reducible Fe and Al (NaOH-P); calcium-bound P, mainly from apatite minerals (HCl-P) and refractory, mainly organic P (Res-P). Concentrations of P fractions varied during the 100?years with the highest values around 2007-2008 and 1923-1935. The P in "active" layers that are available for bacteria and algae or those undergoing changes and diagenetic transformations in the sediment could follow sediments with the water content of ～88%. Potentially mobile P is not decreased in the sediments deposited 50-100?years ago and makes up ca 63?mg?m(-2)?y(-1) (with range 8.3-23.7% of the total P (TP)) in the accumulation area, and ca 0.047?mg?cm(-3) (with range 1.3-22.4% of TP) in the erosion area. The result shows that 13-60% of TP contained in the surface sediment (from 34 to 398?mg P m(-2) y(-1)) has been remobilised during accumulation and could be exported to the overlying water.     

Spatial and temporal distribution of gaseous elemental mercury in Chongqing, China

In order to investigate the spatial and temporal variability of atmospheric mercury (Hg) in Chongqing, China, gaseous elemental mercury (GEM) was measured from August 2006 to September 2007, using Lumex multifunctional mercury analyzer RA-915(+) (Lumex Ltd., Russia). The mean GEM concentration was 6.74 +/- 0.37 ng m(-3) in Chongqing, much higher than the accepted global background values (1.5-2 ng m(-3)). The GEM concentrations were different in different function areas. GEM in transport, industrial and commercial areas were 7.07 +/- 1.04, 7.05 +/- 0.96 and 6.71 +/- 1.10 ng m(-3), respectively, while GEM was 6.14 +/- 1.30 and 4.32 +/- 1.04 in the educational/recreational and nature conservation areas, suggesting that Hg emissions from mobile vehicles and industrial sources (specially coal combustion) were the most important contributors to atmospheric Hg in Chongqing. Mean Hg concentrations also had monthly variations with highest in November (8.24 +/- 0.50 ng m(-3)) and lowest values in August (5.36 +/- 0.70 ng m(-3)). Additionally, the diurnal variation of GEM concentrations was dependent on the local/regional atmospheric conditions. At Jinyun Mountain site (natural conservation area), hourly GEM concentrations had much higher values in daytime than at night. At Power Plant site, however, the hourly GEM concentrations were lower in daytime than at night. GEM concentrations in the air were correlated significantly with meteorological parameters except for barometric pressure.     

Changes in the area of inland lakes in arid regions of central Asia during the past 30 years

Inland lakes are major surface water resource in arid regions of Central Asia. The area changes in these lakes have been proved to be the results of regional climate changes and recent human activities. This study aimed at investigating the area variations of the nine major lakes in Central Asia over the last 30 years. Firstly, multi-temporal Landsat imagery in 1975, 1990, 1999, and 2007 were used to delineate lake extents automatically based on Normalized Difference Water Index (NDWI) threshold segmentation, then lake area variations were detailed in three decades and the mechanism of these changes was analyzed with meteorological data and hydrological data. The results indicated that the total surface areas of these nine lakes had decreased from 91,402.06 km² to 46,049.23 km² during 1975?C2007, accounting for 49.62% of their original area of 1975. Tail-end lakes in flat areas had shrunk dramatically as they were induced by both climate changes and human impacts, while alpine lakes remained relatively stable due to the small precipitation variations. With different water usage of river outlets, the variations of open lakes were more flexible than those of other two types. According to comprehensive analyses, different types of inland lakes presented different trends of area changes under the background of global warming effects in Central Asia, which showed that the increased human activities had broken the balance of water cycles in this region.     

Analysis of spatial and seasonal distributions of MODIS aerosol optical properties and ground-based measurements of mass concentrations in the Yellow Sea region in 2009

Satellite-retrieved data on aerosol optical depth (AOD) and Ångström exponent (AE) using a moderate resolution imaging spectrometer (MODIS) were used to analyze large-scale distributions of atmospheric aerosols in East Asia. AOD was relatively high in March (0.44?±?0.25) and low in September (0.24?±?0.21) in the East Asian region in 2009. Sandstorms originating from the deserts and dry areas in northern China and Mongolia were transported on a massive scale during the springtime, thus contributing to the high AOD in East Asia. However, whereas PM10 with diameters ≤10 μm was the highest in February at Anmyon, Cheongwon, and Ulleung, located leeward about halfway through the Korean Peninsula, AOD rose to its highest in May. The growth of hygroscopic aerosols attendant on increases in relative humidity prior to the Asian monsoon season contributed to a high AOD level in May. AE typically appears at high levels (1.30?±?0.37) in August due to anthropogenic aerosols originating from the industrial areas in eastern China, while AOD stays low in summer due to the removal process caused by rainfall. The linear correlation coefficients of the MODIS AOD and ground-based mass concentrations of PM10 at Anmyon, Cheongwon, and Ulleung were measured at 0.4~0.6. Four cases (6 days) of mineral dustfall from sandstorms and six cases (12 days) of anthropogenically polluted particles were observed in the central area of the Korean Peninsula in 2009. PM10 mass concentrations increased at both Anmyon and Cheongwon in the cases of mineral dustfall and anthropogenically polluted particles. Cases of dustfall from sandstorms and anthropogenic polluted particles, with increasing PM10 mass concentrations, showed higher AOD values in the Yellow Sea region.     

Determination of carbonyls and their sources in three sites of the metropolitan area of Costa Rica,Central America

Ambient levels of carbonyl compounds and their possible sources were studied at three places in the metropolitan area of Costa Rica, including a residential, an industrial, and a commercial downtown area with high vehicular flow, during the periods of April–May and September–December 2009. Fifteen carbonyl compounds were identified in the ambient air, of which acetone was the most abundant carbonyl, followed by formaldehyde and acetaldehyde. Concentrations were highest in rainy season at all sites and lower in dry season. These decreases in concentration are explained by the influences of both photochemical reactions and local meteorological conditions. The strong correlation between C1–C2 and C3 indicated a common origin for these carbonyls. The C1/C2 ratios varied between 0.49 to 1.05, values which can be considered typical of an urban area.     

盐城市空气微生物调查与评价

为了解盐城市空气微生物污染状况,于2009—2013年在盐城市布设了4个城市主要功能区点位和1个市郊清洁参照点,采用自然沉降法采集空气中细菌和马丁霉菌样品进行分析。结果表明,盐城市空气中细菌数量最多的为交通区,其次依次为文教区、工业区、居民区,清洁参照点最少;马丁霉菌数量最多的是交通区,其次依次为工业区、居民区、文教区,清洁参照点最少。盐城市空气微生物污染级别除交通区外,其他功能区均处于较清洁或轻微污染,污染程度呈逐年波动下降趋势。     

Interaction effects of elevated CO2 and temperature on microbial biomass and enzyme activities in tropical rice soils

The impacts of elevated CO(2) and temperature on microbial biomass and soil enzyme activities in four physicochemically different types of tropical rice soils (Aeric Endoaquept, Aeric Tropoaquept, Ultic Haplustalf and Udic Rhodostalf) were investigated in a laboratory incubation study. Soil samples were incubated under 400, 500 and 600 μmol mol(-1) CO(2) concentration at 25°C, 35°C and 45°C for 2 months. Elevated CO(2) significantly increased the mean microbial biomass carbon (MBC) content, across the soils, over control by 6.2%, 38.0% and 49.2% at 400, 500 and 600 μmol mol(-1) CO(2) concentration, respectively. Soil enzyme activities (fluorescein diacetate hydrolase, dehydrogenase, β-glucosidase, urease, alkaline and acid phosphatases) also increased significantly ranging from 1.3% (urease) to 53.2% (alkaline phosphatase) under high CO(2) in the soils studied. Both MBC and soil enzyme activities were further stimulated at high temperatures suggesting elevated CO(2) and high temperature interaction accelerated the general turnover of the organic C fractions of the soil and through increase in microbially mediated processes.     

胶州湾大型底栖动物的次级生产力

基于国外大型底栖动物次级生产力的评估方法,结合国内研究现状,选取了3种Brey模型对胶州湾大型底栖动物的次级生产力进行估算,通过2015—2017年连续3个季度对胶州湾海域大型底栖动物的调查,最终结果表明:胶州湾大型底栖动物的次级生产力与其他海域相比具备较高的水平,次级生产力的季节性变化明显,其中夏季为最高;湾内高值区位于湾顶中部红岛以南到向西偏移至大沽河口和红石崖附近。分析发现历年来大型底栖动物群落的物种组成基本保持稳定,菲律宾蛤仔(Ruditapes philippinarum)资源的变化是造成次级生产力变化的最主要因素。研究推荐使用Brey(2012)模型来估算胶州湾大型底栖动物次级生产力,春、夏、秋季的估算结果分别为880.17、1 537.64、458.33 kJ/(m~2·a)。     

Assessing the changes in land use and ecosystem services in an oasis agricultural region of Yanqi Basin,Northwest China

The Yanqi Basin, one of the most productive agricultural areas, has a high population density in Xinjiang, Northwest China. Land use changes, mainly driven by oasis expansion, significantly impact ecosystem services and functions, but these effects are difficult to quantify. The valuation of ecosystem services is important to clarify the ecological and environmental changes caused by agriculturalization of oasis. This study aimed to investigate variations in ecosystem services in response to land use changes during oasis agricultural expansion activities in the Yanqi Basin from 1964 to 2009. The methods used were based on formula of ecosystem service value (ESV) and ESV coefficients. Satellite data were combined with the ESV coefficients to quantify land use changes and ecosystem service changes in the study area. Sensitivity analysis determined the effect of manipulating the coefficients on the estimated values. The results show that the total ESVs in the Yanqi Basin were $1,674, $1,692, $1,471, $1,732, and $1,603 million in 1964, 1973, 1989, 1999, and 2009, respectively. The net deline in ESV was $71 million in the past 46 years, but the ESVs of each types of landscape changed significantly. The aggregated ESVs of water areas and wetlands were approximately 80 % of the total ESV. Water supply and waste treatment were the two largest service functions and contributed approximately 65 % of the total ESV. The estimated ESVs in this study were elastic with respect to the value coefficients. Therefore, the estimations were robust in spite of uncertainties on the value coefficients. These significant changes in land use occur within the entire basin over the study period. These changes cause environmental problems, such as land degradation, vegetation degeneracy, and changes in aquatic environment.     

Evaluation of river pollution of neonicotinoids in Osaka City (Japan) by LC/MS with dopant-assisted photoionisation

An atmospheric pressure photoionisation (APPI) source for liquid chromatography/mass spectrometry (LC/MS) was applied to determine neonicotinoid pesticides in the aquatic environment. Dopant-assisted APPI was very effective in the ionisation of neonicotinoids. Neonicotinoids generated protonated molecules in APPI with high sensitivity, while adduct ions, such as sodiated molecules, were predominantly generated in conventional electrospray ionisation. The ionisation of neonicotinoids was confirmed by ultra-high-resolution MS. An analytical method coupled with solid phase extraction was developed for acetamiprid, clothianidin, dinotefuran, imidacloprid, nitenpyram, and thiamethoxam. Method detection limits were 0.47 to 2.1 ng L(-1) for six neonicotinoids. Dinotefuran was the most frequent and highest among the neonicotinoids examined in the aquatic environment in Osaka, Japan. The maximum concentration of dinotefuran was 220 ng L(-1). Given the toxicity of neonicotinoids for aquatic creatures, the concentrations observed here were substantially low. The change in concentrations was temporally coincident with the period of the neonicotinoid application. Although rapid photodegradation and some degradation products have been elucidated, the degradation products in the aquatic environment were not identified in the present study.     

Effects of temperature on the acute toxicity of cadmium to Corophium Insidiosum

This study aimed to optimise the methodology for the use of Corophium insidiosum in a bioassay. Taking into account that it would be suitable to execute the bioassay with organisms having a good sensitivity during the year and low mortality in control sediment, the influence of different temperatures (10–15–20 and 25°C) has been examined. C. insidiosum was collected during August, November 2005 and January 2006, in Mar Piccolo basin (Ionian sea). The results obtained show that this species mortality in the negative control sediment, ranged from 2.6 ± 0.6% at 10°C in August to 17 ± 2.2% at 20°C in November, at different temperatures tested. At 20°C there were significant differences in mortality among different months examined. Indeed no relationship among months was found at 15°C. Significant differences between August and November at 25°C, between November and January were not found at 10°C. The 96-h LC50 values found for cadmium at all temperature experimental conditions ranged from 2.11 mg/l (1.57–2.82) to 0.70 mg/l (0.54–0.93). The highest values were found at 10°C in November and January. The results showed that the optimal temperature for the bioassay seems to be between 15°C and 20°C. Even if, at 20°C the mortality differs significantly among organisms sampled.     

Regulation of the dissolved phosphate concentration of a mountainous stream, Kitakyushu, southwestern Japan

The phosphate concentration in mountainous stream water can be a measure of the forest condition, because its concentration will be low when the biomass in the forest is increasing and vice versa when the forest is declining. To investigate the seasonal change in the dissolved phosphate concentration of the mountainous stream water of the Yamakami River, Kitakyushu, from June 2009 to August 2010, and the regulation mechanism of the phosphate concentration, solid-phase spectrophotometry, which can be applicable to natural water without any pretreatment procedures, was employed for the determination of phosphate at μg P L(-1) levels in natural water. The phosphate concentrations in the mountainous stream waters at 6 sites ranged from 2.2 to 13 μg P L(-1), and those from the catchment area of the steady state forest were 5.3 ± 1.6 (±1 SD) μg P L(-1). Changes in the concentration were fairly small even during a storm runoff. The average phosphate concentration of rain was 2.8 ± 0.7 μg P L(-1), about half of the concentration in the stream water. The rate of runoff in forest areas is generally considered to be about 50% of the total precipitation. For a forest under a climax condition, the phosphate concentration is estimated to be regulated by the fallout and evapotranspiration (α = 0.05). At one of the sites, an upstream tributary, where a fairly big landslide occurred before July in 2009, the phosphate concentration was the highest, suggesting that the biomass may still be decreasing. For all of the six sites examined, a characteristic seasonal change in phosphate concentration was observed, reflecting the local budget between the biological decomposition of plant matter and the consumption by the biomass. The increase in the phosphate concentration during late spring and early summer may result from the extensive decomposition of plant litter mainly supplied in autumn and of plant matter relating to spring blooming such as fallen flowers, pollen and immature fruits. The proposed method using the phosphate concentration in surface stream waters without the period of the seasonal change mentioned above is expected to be very helpful in diagnosing the condition of forests.     

Accumulation of endocrine disrupting chemicals (EDCs) in the polychaete Paraprionospio sp. from the Yodo River mouth,Osaka Bay,Japan

This study presents the levels of endocrine disrupting chemicals (EDCs) accumulated by Paraprionospio sp. from the Yodo River mouth, Osaka Bay. Since high concentrations of nonylphenol (NP), bisphenol A (BP), octylphenol (OP), 17β-estradiol (E2), and estrone (E1) have been measured in sediment from Osaka Bay, some bioaccumulation could be expected particularly in benthic animals. EDCs were analysed in Paraprionospio sp., a dominant benthic species in Osaka Bay. The results showed that Paraprionospio sp. had accumulated varying concentrations (wet weight; w.w.) of NP at 1,460–4,410 ng/g; BP at 22.5–39.6 ng/g; OP at 18.9–45.4 ng/g; E2 at 0.89–4.35 ng/g; and E1 at 0.06–2.50 ng/g. Accumulation of NP and OP were highest among the samples gathered in summer (July 2008), while concentrations of BP, E2, and E1 did not much differs within 3 years. EDC levels in Paraprionospio sp. were apparently greater than those in sediments showing bioaccumulation.     

Natural regeneration of the herbaceous community in a semiarid region in Northeastern Brazil

This study aimed to describe and compare the interannual changes in the diversity and population structure of herbaceous plants in an anthropogenic area that has been regenerating for 15 years and to identify the similarities and differences in the biological attributes of the community compared with the characteristics of a regenerating conserved area. In total, 105 plots measuring 1 m² were established. In each plot, the herbaceous plants were identified, and their height and stem diameter were measured for two consecutive years. The herbaceous flora of the anthropogenic area was represented by 86 species in 70 genera and 27 families, and there were no significant differences in the average richness between years. The conserved area was represented by 71 species in 63 genera and 35 families, and there was a significant difference in the total richness between areas and between years, except when comparing the richness between the conserved area and the anthropogenic area during the second year. Considering both the anthropogenic and conserved areas, 123 herbaceous species were listed, and the similarity between areas was 60 %. For the anthropogenic area, the floristic similarity between years was 95 %, and in the fragment of the conserved area, the similarity was 74 %. The diversity and density were significantly different between years and between areas. Given these results, this study suggests that 15 years of natural regeneration for the caatinga is not sufficient to reestablish its native flora with respect to its herbaceous component.     

Recent Changes in Soil Total Phosphorus in the Everglades: Water Conservation Area 3

We assessed recent changes in the distribution of soil total phosphorus (TP) in Water Conservation Area 3 (WCA-3) of the Everglades. Soil cores were collected in 1992 and 2003 at 176 sites. To reflect hydrologic boundaries within the system, WCA-3 was divided into three zones (3AN, 3AS, and 3B). Total P was mapped on both a mass (TPm) and a volumetric basis (TPv) to determine if spatial distributions varied depending on the choice of units. Interpolated maps for both years showed that the highest levels of TPm were located in 3AN and in boundary areas of all zones that received surface water inputs of P from canals. Increases in TPm were greatest in central 3AN in an area adjacent to the Miami Canal that received inputs from a water control structure. Interpolated maps for TPv illustrated that a hotspot present in 1992 had disappeared by 2003. The highest levels of TPv in 2003 were located in northwestern 3AN, a region of WCA-3 that has been chronically overdrained and burned in 1999. From 1992 to 2003, increases in TPm were observed for 53% of the area of WCA-3, while only 16% of WCA-3 exhibited increases in TPv. In 1992, approximately 21% of WCA-3 had TPm concentrations in the 0–10 cm layer >500 mg kg⁻¹, indicating P enrichment beyond historic levels. Eleven years later, 30% of the area of WCA-3 had TPm >500 mg kg⁻¹. This indicated that during this period, the area of WCA-3 with enriched TPm concentrations increased about one % year⁻¹.     

Spatio-temporal changes in agrochemical inputs and the risk assessment before and after the grain-for-green policy in China

China’s Grain-For-Green Policy (GFGP) of returning marginal cropland to forest or grassland is one of the most important large-scale initiatives to combat land degradation in its ecologically vulnerable regions. In order to maintain and increase crop production from decreasing areas of cropland, substantial spatio-temporal changes in agrochemical inputs have occurred, which have strongly influenced the ecological and environmental status of land in China. Based on the agrochemical inputs (chemical fertilizer, pesticide, plastic sheeting, and agricultural diesel oil) at the provincial level between 1993 and 2009, cluster analysis and gravity center modeling were used to trace these spatio-temporal changes. A regional comparative study was also undertaken to investigate the changes in the relative size of agrochemical inputs in the eastern, central, and western regions of China. It was found that the agrochemical inputs increased considerably at the nation level after the GFGP, which in order of increasing rate were: plastic sheeting > agricultural diesel oil > pesticide > chemical fertilizer. The gravity centers of agrochemical inputs moved substantially towards the northwest or west during the latter period of GFGP and regional comparative analysis showed that the agrochemical inputs increased substantially in the western region between 2004 and 2009. The ecological degradation caused by the expansion of the area devoted to crop production in the western region and the potential risk of agricultural non-point pollution caused by the increasing agrochemical inputs are the main factors restricting this area’s sustainable development.