Accepted in Journal of Atmospheric and Oceanic Technology Corresponding Author Address
- 格式:pdf
- 大小:1.37 MB
- 文档页数:41
2A review of research progress on CO capture,storage,and utilization in 34Q156789111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576with 58%of the sources being located within the east and south 77central regions.The contributions of large point sources in each 78sector to total CO 2emissions in China are listed in Fig.2[13].With 79rapid development of energy technologies in the 21st century,fos-80sil fuels,especially coal,will still remain the dominant energy0016-2361/$-see front matter Ó2011Elsevier Ltd.All rights reserved.doi:10.1016/j.fuel.2011.08.022⇑Corresponding authors.Q2Address:State Key Laboratory of Coal Conversion,Institute of Coal Chemistry,Chinese Academy of Sciences,Taiyuan 030001,China (Y.Sun).Tel.:+863514049612;fax:+863514041153.E-mail addresses:zhaoning@ (N.Zhao),weiwei@ (W.Wei),yhsun@ (Y.Sun).Please cite this article in press as:L Q1i L et al.A review of research progress on CO 2capture,storage,and utilization in Chinese Academy of Sciences.Fuel(2011),doi:10.1016/j.fuel.2011.08.022source in China for decades to come.Chinese government recognized the huge challenge of CO 2abatement while satisfying ever-increasing energy demand.In the light of this situation,November 26,2009,China officially announced action to control emissions per unit of GDP by 40–45%by 2020,based on 86levels [14].To address this,China is undertaking a range of techni-87cal research and development projects on CCSU,including the na-88tional fundamental research and high-tech programs,as well as a 89large number of international programs.The CCS projects,fun-90dings,and research institutes in China is shown in Table 1.91Since 1990,China had carried out a series of climate change 92projects under framework of national programs,such as China’s 93National Climate Change Program (CNCCP),National Hi-tech R&D949596979899100101102103104105106107108109110111112113114115116117118119120121cooperation with international on CCSU.On the basis of the finical 122support of both Chinese government and CAS,a lot of progresses 123were obtained in several academic institutes in CAS including 124CO 2capture;enhanced oil recovery (EOR)and enhanced coal bed 125methane (ECBM)projects as well as CO 2chemical utilizations.This 126brief review has covered the research progress in CO 2capture,stor-127age,and utilization in CAS.2.The contributions of large point sources in each sector to overall total emissions in China [12].Please cite this article in press as:L Q1i L et al.A review of research progress on CO 2capture,storage,and utilization in Chinese Academy of Sciences.Fuel(2011),doi:10.1016/j.fuel.2011.08.022128129130131132133134135136137aration,and cryogenic fractionation.138 2.1.Amine-based scrubbing solvent139Amine scrubbing is a well known technology for capturing CO 2140from flue gas,which has been widely deployed on a large scale 141across several industries [25–28].The industrially most important142143144145146147148149150151152the environment [30].1532.2.Ionic liquids154Therefore,a nonvolatile solvent that could facilitate CO 2capture 155without the loss of solvent into the gas stream would be advanta-156geous.Ionic liquids (ILs)are commonly defined as liquids whichTable 1CCS projects,fundings,and research institutes in China.China and international cooperation on CCS projects and fundingsResearch institutes aBNLMS–CAS IET–CAS RCEES–CAS ICC–CAS IPE–CAS LICP–CAS CIAC–CAS National High Technology Research and Development Program of China (863)SIC–CAS National Key Basic Research and Development Program of China (973)IGG–CAS China’s National Climate Change Program (CNCCP)IAP–CAS L.L Q1i et al./Fuel xxx (2011)xxx–xxx3Please cite this article in press as:L Q1i L et al.A review of research progress on CO 2capture,storage,and utilization in Chinese Academy of Sciences.Fuel(2011),doi:10.1016/j.fuel.2011.08.022157are composed entirely of ions with a melting point of less than 158100°C.ILs have many unique properties in comparison to other 159solvents as extremely low volatility,broad range of liquid temper-160ature,high thermal and chemical stability,and tunable physico-161chemical characteristics and as a result,ILs have been considered 162as a potential substitute of aqueous amine solutions for CO 2cap-163ture [31–34].164In Changchun Institute of Applied Chemistry,Chinese Academy 165of Sciences (CIAC–CAS),a novel dissolving process for chitin and 166chitosan was developed by using the ionic liquid 1-butyl-3-167methyl-imidazolium chloride ([Bmim]Cl)as a solvent for capturing 168and releasing CO 2.The results showed that the chitin/IL and chito-169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201tant to consider the maximum mass loading when considering the 202support of ILs on inert substrates –yes,these can enhance the ILs’203ability to take up CO 2,but at the expense of cycling an inert sorbent 204round a thermal cycle.The ILs are also monstrously expensive,the 205complex structure,and high cost for preparation,compared to sim-206pler solvents such as MEA or ammonia.Thus,the potential for 207improving CO 2solubilities and reducing cost of the ILs still needs 208to be studied for future applications.In 2010,in Beijing National 209Laboratory for Molecular,Chinese Academy of Sciences (BNLMS–210CAS),Zhang and co-workers first reported on CO 2capture by 211hydrocarbon surfactant liquids.It also found that CO 2had high sol-212ubility in low-cost hydrocarbon surfactant liquids,and the ab-213and the 214215216such as corrosion,at 217regenerable solid sor-218concept for CO 2recov-219into amine-based and 220221222with various so-223as silica gels,activated 224have been shown to phys-225enhance the sorp-226of many amine-based 227In Dalian Institute of 228(DICP–CAS),Zhang 229silica foam (MCF)materi-230with polyethyl-231The results showed that 232having large window 2333.45mmol CO 2/g sorbent 234In Institute of Coal Chem-235Zhao et al.studied 236materials derived 237sorption capabili-238°C and 1bar.The as-pre-239selectivity for CO 2over 240prepared a series of CO 2241pentamine (TEPA)was 242(PMHS)based mesopor-243The highest absorption 24475°C with the 10vol.%245was higher than most 246Desorption could 247in 1h [46].248of amine-based sorbent 249capacity of solid sorbent.250have poor mechanical 251amine-based sorbents 252and require signifi-253processes.254255sorbents for CO 2capture 256Alkali earth metal,such 257form alkali earth metal-258vapor at high tempera-259and post-combus-260simplified process flow 261the calcium looping 262vessel (the carbonator)4L.L Q1i et al./Fuel xxx (2011)xxx–xxxPlease cite this article in press as:L Q1i L et al.A review of research progress on CO 2capture,storage,and utilization in Chinese Academy of Sciences.Fuel(2011),doi:10.1016/j.fuel.2011.08.022the carbonation reaction between CO 2and solid CaO separates CO from coal-combustion flue gas at a temperature between 600°and 650°C.The CaCO 3formed is then passed to another vessel (the calciner),where it is heated to reverse the reaction (900–950°C),releasing the CO 2suitable for sequestration,and regener-ating the CaO-sorbent which is then return to the carbonator.The carbonation process is exothermic,which is matched with the temperature of a steam cycle,allowing recuperation of the heat.In IPE–CAS,the decomposition conditions of CaCO 3particles for CO 2capture in a steam dilution atmosphere (20–100%steam 307308309310311312313314315316317318319320321322323324325326327328329330332333334335336337338339340341342343344345346347348349350351352353354355Fig.3.The process flow diagram of post-combustion capture using the calcium looping cycle [47,48].Please cite this article in press as:L Q1i L et al.A review of research progress on CO 2capture,storage,and utilization in Chinese Academy of Sciences.Fuel(2011),doi:10.1016/j.fuel.2011.08.022strength (998N/cm 2)and exhibited good stability in multiple cy-cles [74,75].Furthermore,the application of a conceptual CO 2cap-ture process using this sorbent was proposed for an existing coal fired power plant [75].However,to optimize CO 2sorption capacity,understand of the interaction between CO 2and the sorbent need to be studied in the further work.Moreover,much work remains before the technology fluidized bed CO 2capture can be commercialized.Simulta-neously,the numerical simulation based on the computational fluid 365dynamics (CFD)method will become a research focus in the future.366 3.CO 2storage367Following the capture and transport process,CO 2can be dis-368posed of in natural sites such as deep geological sequestration,369mineral carbonation,or ocean storage [76].There are three geolog-370ical formations that have also been recognized as major potential 371CO 2sinks:deep saline-filled sedimentary (DSFs),depleted oil 372natural gas reservoirs,and unmineable coal-seams.The geology 373also suggests possibilities for CO 2enhanced oil recovery (CO 374EOR),CO 2enhanced gas recovery (CO 2–EGR)and CO 2enhanced 375coal-bed methane recovery (CO 2–ECBM)projects [12].3763.1.Geological sequestration377Geological storage involves injecting CO 2at depths greater than 3781000m into porous sedimentary formations using technologies 379derived from the oil and gas industry [77].CO 2can be stored in 380supercritical state at depth below 800–1000m,which provides 381the potential for efficient utilization of the space,due to the li-382quid-like density of supercritical CO 2.The point at which CO 2Table 2Performance summary of K-based sorbents capturing CO 2.Material Temperature (°C)CO 2partial pressure (bar)e Total capacity (mmol CO 2/g sorbent)Method f Regenerature temperature (°C)Ref.K 2CO 3/AC a 600.01 1.95TCD g 150[62]K 2CO 3/SiO 2600.010.23TCD g –[62]K 2CO 3/USY 600.010.43TCD g –[62]K 2CO 3/CsNaX 600.01 1.35TCD g –[62]K 2CO 3/Al 2O 3600.01 1.93TCD g 350[62]K 2CO 3/CaO 600.01 1.11TCD g –[62]K 2CO 3/MgO 600.01 2.70TCD g 400[62]K 2CO 3/TiO 2600.01 1.89TCD g 150[62]K 2CO 3/Al 2O 3600.01 1.96TCD g >300[63]Re-KAl(I)30b 600.01 1.86TCD g <200[63]g Fig.5.The schematic diagram of experimental apparatus for the fluidized bed [74,75].6L.L Q1i et al./Fuel xxx (2011)xxx–xxxPlease cite this article in press as:L Q1i L et al.A review of research progress on CO 2capture,storage,and utilization in Chinese Academy of Sciences.Fuel(2011),doi:10.1016/j.fuel.2011.08.022383transforms from critical to supercritical point is 31.1°C and 3847.38MPa [78].CO 2is injected usually in the supercritical form into 385the saline aquifer or depleted oil or gas reservoir.Four major clas-386ses of deep geologic reservoirs present within China have been 387identified and evaluated as candidates for the long-term storage 388of anthropogenic CO 2:deep saline-filled sedimentary (DSFs)for-389mations,depleted gas basins,depleted oil basins with potential 390for CO 2–EOR,and deep unmineable coal seams with potential for 391CO 2–ECBM.Fig.6shows the map of the combined location and ex-392tent of candidate geologic CO 2storage formations in China [13].393Because the CO 2industry is not mature,there are few active CO 2394storage projects which can provide site specific information;hence417China are also potential storage candidates.Recently,other re-418search has also focused on estimating the distance between CO 2419sources and potential sinks.Zheng et al.superimposed the loca-420tions of these 27facilities onto maps of sedimentary basins in each 421of the five regions of China (Huabei,Ordos,Dongbei,Yuwan,and 422Xinjiang).The majority of the candidate CO 2sources are found in 423the Ordos,Huabei and Dongbei regions [85].424The China–UK Near Zero Emissions Coal (NZEC)Initiative exam-425ined options for carbon (CO 2)capture,transport and geological 426storage in China,which was developed under the 2005EU–China 427NZEC Agreement that aims to demonstrate CCS in China and the 428EU [16,86–88].The NZEC Initiative has evaluated the potential to 429430431432433434435436437438439440441442443444445446447448449450451L.L Q1i et al./Fuel xxx (2011)xxx–xxx7Please cite this article in press as:L Q1i L et al.A review of research progress on CO 2capture,storage,and utilization in Chinese Academy of Sciences.Fuel(2011),doi:10.1016/j.fuel.2011.08.022452Kailuan mining area (Hebei Province)and deep saline aquifers in 453the Jiyang Depression (Shandong province)[89,90].The results 454show that the Dagang oilfield is not suitable for large-scale storage,455though could be considered for EOR pilots.The Shengli oilfield was 456considered more promising for storage (472Mt in eight selected 457fields).Storage potential in the Kailuan mining area is 504,000458Mt adsorbed onto the coal and 38,100Mt void storage capacity.459However,the coals have low porosity and permeability that will af-460fect future energy resources [90].The Institute of Geology and Geo-461physics,Chinese Academy of Sciences (IGG–CAS)studied the 462potential for storage in the Jiyang Depression.The results revealed 463that Guantao Formation in the Jiyang Depression has good porosity 464and permeability 465areas was 4662010,in South China 467of Sciences 468storage capacity in 469the Pearl River 470CCS-related 471China [78].There 472saline formations 473tive storage 474including 60Mt in 475large for storaging 476in Guangdong in 477In a word,these 478age of CO 2in deep 479Although this is 480countries,it will 481and there was 482characteristics.483 3.1.2.CO 2–EOR484Although CO 2485oil recovery (EOR)486this process can be 487oil,the cost of CO 2,488the CO 2source [92].489the production of 490be an ideal option 49184commercial or 492tion worldwide [1].493been implemented 494Oil Corporation 495in the Daqing,496ernments of Japan 497out a project to 498plant in China into a 499duced from the 500that between 270501ered by using CO 2502including IGG–CAS 503three large oil fields 504oil reservoirs in the 505were suitable both 506found suitable for 507showed that the 508CO 2storage 509the oil recovery by steam injection has been already applied at Lia-510ohe oil field.Each single well,in average,had conducted 7.6times 511of steam injection-oil recovery processing for EOR propose.The to-512tal recovered oil amounts were 12.06Mt [92].Active oil producing 513fields where CO 2–EOR is technically possible provide credible 514opportunities to initiate CO 2storage demonstration projects.How-515ever,significant further investigations,including detailed site516appraisals would be necessary before such fields can be considered 517as technically and economically suitable for CO 2storage.5183.1.3.CO 2–ECBM519In a similar manner,ECBM recovery can be used to store CO 2520while improving methane recovery.A bright prospect of gas injec-521tion technology for ECBM production has been suggested by Chi-522nese engineers since the late 1990s [94].More recently,a joint 523venture was formed between the China United Coal Bed Methane 524Corporation and the Alberta Research Council of Canada to develop 525a project entitled ‘‘Development of China’s coalbed methane tech-526nology/CO 2sequestration’’[12].This project was initiated in March 527project was performed 528in the anthracitic coals of 529China [95],which is 530in China up to now 531at ICC–CAS in 2005532were investigated based 533An equipment simulated 534middle pressures was 535in coal seam,536behaviors were studied.537coal mine and salt-water 538four coals of various rank 539China were tested for 540one of the most impor-541process.The result 542capacities for methane 543>Bulianta coal >Zhangji 544adsorption isotherms 545lattice model [100].546given to estimate the 547[101],which was 548underground stress is so 5492and CH 4respectively 550the mechanical sta-551of the impact factors of 552stress could be obviously 553of the casing or by using 554and Young’s modulus 555China was also estimated 556prospecting data of coal 557and the replacement ratio 558different ranks,it is esti-559methane resources will 560technology is uti-561in coalbeds is about 562as the total CO 2emission 563also developed 564simulation of the CO 2–565is a lack of knowledge 566due to the complexity 567fluid transport processes.568will be the next 569570571Large amounts of CO 2can also be fixed by a process called min-572eral carbonation,which is natural or artificial fixation of CO 2into 573carbonates.It has been proposed as a promising CO 2sequestration 574technology e.g.the silicate rocks (calcium or magnesium)could be 575turned into carbonates by reacting with CO 2following this mech-576anism [8,105]:577ðMg ;Ca Þx Si y O x þ2y þx CO 2!x ðMg ;Ca ÞCO 3þy SiO 25798L.L Q1i et al./Fuel xxx (2011)xxx–xxxPlease cite this article in press as:L Q1i L et al.A review of research progress on CO 2capture,storage,and utilization in Chinese Academy of Sciences.Fuel(2011),doi:10.1016/j.fuel.2011.08.022580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639hanced by the fact that this method of storage is highly verifiable 640and unquestionably permanent,the grinding energy required to 641produce particles of the size required to react rapidly with the 642acids is large,and the residence times on the order of hours re-643quired to allow carbonation of the solids,via either route,is so long 644that immense reactors would be required,associating environmen-645tal concerns.Furthermore,mineral carbonation will always be646expensive than most applications of geological storage 647important gap in mineral carbonation is the lack of 648onstration plant.649Ocean storage650Captured CO 2also could help reduce the atmospheric 6516526536546556566576586596606616626636646656666676686696706716726736746756766776786796806816826836842685carbonic acid,which would be likely harmful to ocean organisms 686and ecosystems [17].Additionally,it is not known whether the 687public will accept the deliberate storage of CO 2in the ocean as part 688of a climate change mitigation strategy.The development of ocean 689storage technology is generally at a conceptual stage;thus,further 690research and development would be needed to make technologies 691available.3.Reaction mechanism for enhanced carbonation crystallization Q1Please cite this article in press as:L Q1i L et al.A review of research progress on CO 2capture,storage,and utilization in Chinese Academy of Sciences.Fuel(2011),doi:10.1016/j.fuel.2011.08.022692693694695696697698699700701702703704705706707708709710711712713714715716717718719721722723the reaction [131–135],such as in ICC–CAS,Zhao et al.had reported 724a catalyst system composed of KI supported on metal oxides for 725cycloaddition of propylene oxide with CO 2.It was found that the 726activity of KI for cycloaddition was greatly enhanced by ZnO as both 727support and promoter,resulting in a high yield of propylene 728carbonate within a short reaction time.The mechanism is also pro-729posed (Scheme 4)[133].Recently,a large number of catalytic730systems,such as metal oxides,transition metal,ammonium 731well as main group complexes,were reported to be active 732reactions [136–139].In ICC–CAS,the efficient ultrasonic tech-733nique was used for the preparation of amine-functionalized porous 734catalysts for CO 2coupling with epoxide.According to 735study by Zhang and co-workers [140],the reaction conditions 736great influence on the performance and the silanols on the surface 737played an important role in the chemical fixation of CO 2.In addi-738they also proposed the possible reaction mechanism for 739coupling with epoxide over such type of catalysts (Scheme 5).740In recent years,ionic liquids as environmentally benign media 741organic synthesis and catalytic reaction significant progress 7427437447457467477487492750catalyst system without using additional organic solvents was 751achieved in excellent selectivity and TOF (5410h À1)[144].In 752IPE–CAS,an efficient Lewis acid/base catalyst composed of ZnCl 2/753PPh 3C 6H 13Br was developed and showed high activity and selectiv-754ity for the coupling reaction of CO 2and epoxide under the mild 755conditions [145].Sun et al.prepared a series of hydroxyl-function-756alized ionic liquids (HFILs)which showed efficient reactivity andScheme 4.The proposed diagram of reaction mechanism [133].Q1Please cite this article in press as:L Q1i L et al.A review of research progress on CO 2capture,storage,and utilization in Chinese Academy of Sciences.Fuel(2011),doi:10.1016/j.fuel.2011.08.022780781782783784785786787788789790791792793794795796797798799800802803804805806807Please cite this article in press as:L Q1i L et al.A review of research progress on CO 2capture,storage,and utilization in Chinese Academy of Sciences.Fuel(2011),doi:10.1016/j.fuel.2011.08.022CO 2þCH 4¼2CO þ2H 2841842In the past decade,a lot of researches have been devoted to the 843catalytic performance of noble metals,including Pt,Ru,Rh,Pd,844and Ir for this reaction [165–169].It showed that Rh and Ru845exhibited both high activity and stability in CH 4dry reforming,846while Pd,Pt and Ir were less active and prone to deactivation.847Nevertheless,considering the aspects of high cost and limitedPlease cite this article in press as:L Q1i L et al.A review of research progress on CO 2capture,storage,and utilization in Chinese Academy of Sciences.Fuel(2011),doi:10.1016/j.fuel.2011.08.022848availability of noble metals,it is more practical to develop non-849noble metal catalysts which exhibited both high activity and pared with noble metals,Ni-based catalysts have been 851widely investigated because of their high activity and relatively 852low price [170–172].Nevertheless,application of Ni-based cata-853lysts in a large scale process is not so straightforward due to rapid 854carbon deposition,resulting in the deactivation of the catalyst 855[173].It was found that when Ni is supported on a alkaline earth 856metal oxide such as MgO,CaO,and BaO with strong Lewis basi-857city,carbon deposition can be attenuated or even suppressed 858[174]which is because that the support could promote chemi-859sorption of CO 2and thus,accelerated the reaction of CO 2and C 8608618628638648658668678688698708718728738748758768778788798808818828838848858868878888898908918928938948958968974.4.Reaction of CO 2with ethane and propane898Ethylene and propylene are basic raw material in the petrol-899chemical industry.Thermal cracking of hydrocarbons (such as eth-900ane)in the presence of steam is currently the main source of eth-901ylene [181,182].Nevertheless,steam cracking of ethane to 902ethylene is a highly endothermic process that must be performed 903at high temperatures,which means the consumption of a large 904amount of energy.The introduction of CO 2could reduce the extent 905of deep oxidation which results in many byproducts whereas eth-906ylene selectivity drops when oxygen is used as oxidant [183].907Thermodynamics analysis and experimental results have indi-908909910911912913915916917918919920921922923924925926927928929930931932933934935936938939940941L.L Q1i et al./Fuel xxx (2011)xxx–xxx13Please cite this article in press as:L Q1i L et al.A review of research progress on CO 2capture,storage,and utilization in Chinese Academy of Sciences.Fuel(2011),doi:10.1016/j.fuel.2011.08.022969970971972973974975976977978979980981982983984985986987988990991992993994995996997998999100010011002100310041005100610071008100910101011101210131014Please cite this article in press as:L Q1i L et al.A review of research progress on CO 2capture,storage,and utilization in Chinese Academy of Sciences.Fuel(2011),doi:10.1016/j.fuel.2011.08.0221015DMC in supercritical phase was proposed in Scheme 11.Recently,developed a supported Cu-Ni/V 2O 5-SiO 2heterogeneous because the reaction can be carried out in a fixed-bed the side production of water molecules 10341035103610371038103910401041104210431044and theoretical approaches,which enable the development of 1045CO 2selective sorbents.Besides,the sorbent performance,lifetime,10461047104810491050105110521053105410551056105710581059106010611062106310641065Scheme 11.The proposed catalytic reaction mechanism Please cite this article in press as:L Q1i L et al.A review of research progress on CO 2capture,storage,and utilization in Chinese Academy of Sciences.Fuel(2011),doi:10.1016/j.fuel.2011.08.0221066of component costs,specific Chinese market conditions,and other 1067factors impacting costs of deployment in China will be important 1068to consider in greater detail.To propose the storage mechanism,1069monitoring,simulation,risk assessment,control methods as well 1070as engineering design will be studied in future.1071The utilization of CO 2to chemicals has attracted considerable 1072attention as a possible way to manufacture useful commercial 1073chemicals from CO 2in some specific locations (Scheme 12)[222].1074The utilization of CO 2as a raw material in the synthesis of chemi-1075cals was also conducted by CAS,including synthesis of cyclic car-1076bonate from CO 2and epoxide,reaction of CO 2and propylene 1077glycol (PG),CO 2reforming of CH 4,reaction of CO 2and ethane 1078and propane,CO 21079methyl carbonate 1080amount of CO 2can 1081to the order of 1082ent the typical 1083cations is only 1084laminates are used 1085of the materials can 1086and overall net 1087ation.1088and fundamental 1089lue-added chemicals 1090tive net carbon 1091moderate reaction 1092the energy or 1093clear,wind,1094also important to 1095 6.Uncited 1096[201].Q31097Acknowledgments1098This work was 1099vation Programme 1100323);the Ministry of 1101lic of China 1102Climate Change:1103Academy of 1104the Natural Science 1105References1106[1]Mikkelsen M,1107fixation and 11082010;3(1):43–81.1109[2]Xu XC,Song CS,1110separation from 111111121113[3]Shukla R,Ranjith P,1114and caprock 1115[4]Bredesen R,Jordal 1116generation with CO 1117[5]Barelli L,Bidini 1118sorption-enhanced 1119review.Energy 1120[6]Drage TC,1121adsorbents for CO 2capture in gasification.Energy Fuels 2009;23:2790–6.1122[7]Stewart C,Hessami M.A study of methods of carbon dioxide capture and1123sequestration –the sustainability of a photosynthetic bioreactor approach.1124Energy Convers Manage 2005;46:403–20.1125[8]Yang HQ,Xu ZH,Fan MH,Gupta R,Slimane RB,Bland AE,et al.Progress in1126carbon dioxide separation and capture:a review.J Environ Sci 11272008;20(1):14–27.1128[9]Khatri RA,Chuang SSC,Soong Y,Gray M.Thermal and chemical stability of1129regenerable solid amine sorbent for CO 2capture.Energy Fuels 11302006;20(4):1514–20.1131[10]Jin HG,Gao L,Han W,Hong H.Prospect options of CO 2capture technology1132suitable for China.Energy 2009:1–8.1133[11]Hashim H,Douglas P,Elkamel A,Croiset E.Optimization model for energy1134planning with CO 2emission considerations.Ind Eng Chem Res 11352005;44(4):879–90.1136[12]Meng KC,Williams RH,Celia MA.Opportunities for low-cost CO 2storage1137demonstration projects in China.Energy Policy 2007;35(4):2368–78.1138[13]Dahowski RT,Li X,Davidson CL,Wei N,Dooley JJ,Gentiled RH.A preliminary1139cost curve assessment of carbon dioxide capture and storage potential in 1140China.Energy Procedia 2008;1(1):2849–56.1141[14]Wang R,Liu WJ,Xiao LS,Liu J,Kao W.Path towards achieving of China’s 20201142carbon emission reduction target –a discussion of low-carbon energy policies 1143at province level.Energy Policy 2011;39(5):2740–7.1144[15]China’s scientific &technological actions on climate change;2007.<http:///WebSite/CCChina/UpFile/File199.pdf >.1146Pearce J,et al.Carbon capture1147near zero emissions coal 11481149capture and storage for CO 211502010;38(9):5281–9.1151R,Mezghani K,Imashuku S,1152capture utilizing oxy-fuel 1153membrane systems.Int J 11541155and size effects of activated1156Ind Eng Chem Res 11571158to produce hydrogen with1159combustion.Int J Hydrogen 11601161combustion CO 2capture by1162Prog Energy Combust Sci 11631164DA,McMichael WJ.Carbon1165sorbents.Energy Fuels 11661167Park AHA,et al.High efficiency1168based on amine-functionalized 11691170Gimenez A,Sanchez-Biezma A,1171for low cost CO 2capture in 11722008;49(10):2809–14.1173organogels via ‘‘latent’’gelators.1174their ammonium carbamates.11751176and prediction of the solubility1177and mixed alkanolamine 11781179of the structural features1180of CO 2and regeneration in 11811182Svendsen HF.Modeling and1183in aqueous alkanolamine 1184Ind Eng Chem Res 11851186of products of the oxidative1187aqueous monoethanolamine 1188gases.Ind Eng Chem Res 11891190Galindo A,Jackson G,et al.An1191Energy Environ Sci 11921193and thermodynamic11941-n-butyl-3-methylimidazolium 11951196High-pressure phase behavior1197ionic liquids.J Phys Chem B 11981199JF.Anion effects on gas12001201liquids for post-combustion CO 212021203dissolved in ionic liquids as12041205W,Zhang XP.Dual amino-1206for CO 2capture.Chem Eur J 12071208[37]Zhang JM,Zhang SJ,Dong K,Zhang YQ,Shen YQ,Lv XM.Supported absorption1209of CO 2by tetrabutylphosphonium amino acid ionic liquids.Chem Eur J 12102006;12(15):4021–6.1211[38]Sairi NA,Yusoff R,Alias Y,Aroua MK.Solubilities of CO 2in aqueous N-1212methyldiethanolamine and guanidinium trifluoromethanesulfonate ionic 1213liquid systems at elevated pressures.Fluid Phase Equilibria 2011;300(1–12142):89–94.1215[39]Zhang JL,Han BX,Zhao YJ,Li JS,Hou MQ,Yang GY.CO 2capture by1216hydrocarbon surfactant liquids.Chem Commun 2011;47(3):1033–5.16L.L Q1i et al./Fuel xxx (2011)xxx–xxxPlease cite this article in press as:L Q1i L et al.A review of research progress on CO 2capture,storage,and utilization in Chinese Academy of Sciences.Fuel(2011),doi:10.1016/j.fuel.2011.08.022。
V OL .56,N O .215J ANUARY 1999J O U R N A L O F T H E A T M O S P H E R I C S C I E N C E S ᭧1999American Meteorological Society127SCIAMACHY:Mission Objectives and Measurement ModesH.B OVENSMANN ,J.P .B URROWS ,M.B UCHWITZ ,J.F RERICK ,S.N OE¨L ,ANDV .V .R OZANOVInstitute of Environmental Physics,University of Bremen,Bremen,GermanyK.V .C HANCEHarvard–Smithsonian Center for Astrophysics,Cambridge,MassachusettsA.P .H.G OEDESRON Ruimetonderzoek,Utrecht,the Netherlands(Manuscript received 5September 1997,in final form 16June 1998)ABSTRACTSCIAMACHY (Scanning Imaging Absorption Spectrometer for Atmospheric Chartography)is a spectrometerdesigned to measure sunlight transmitted,reflected,and scattered by the earth’s atmosphere or surface in the ultraviolet,visible,and near-infrared wavelength region (240–2380nm)at moderate spectral resolution (0.2–1.5nm,/⌬ഠ1000–10000).SCIAMACHY will measure the earthshine radiance in limb and nadir viewing geometries and solar or lunar light transmitted through the atmosphere observed in occultation.The extraterrestrial solar irradiance and lunar radiance will be determined from observations of the sun and the moon above the atmosphere.The absorption,reflection,and scattering behavior of the atmosphere and the earth’s surface is determined from comparison of earthshine radiance and solar irradiance.Inversion of the ratio of earthshine radiance and solar irradiance yields information about the amounts and distribution of important atmospheric constituents and the spectral reflectance (or albedo)of the earth’s surface.SCIAMACHY was conceived to improve our knowledge and understanding of a variety of issues of importance for the chemistry and physics of the earth’s atmosphere (troposphere,stratosphere,and mesosphere)and potential changes resulting from either increasing anthropogenic activity or the variability of natural phenomena.Topics of relevance for SCIAMACHY areR tropospheric pollution arising from industrial activity and biomass burning,R troposphere–stratosphere exchange processes,R stratospheric ozone chemistry focusing on the understanding of the ozone depletion in polar regions as well as in midlatitudes,andR solar variability and special events such as volcanic eruptions,and related regional and global phenomena.Inversion of the SCIAMACHY measurements enables the amounts and distribution of the atmospheric con-stituents O 3,O 2,O 2(1⌬),O 4,BrO,OClO,ClO,SO 2,H 2CO,NO,NO 2,NO 3,CO,CO 2,CH 4,H 2O,N 2O,and aerosol,as well as knowledge about the parameters pressure p,temperature T,radiation field,cloud cover,cloud-top height,and surface spectral reflectance to be determined.A special feature of SCIAMACHY is the combined limb–nadir measurement mode.The inversion of the combination of limb and nadir measurements will enable tropospheric column amounts of O 3,NO 2,BrO,CO,CH 4,H 2O,N 2O,SO 2,and H 2CO to be determined.1.IntroductionLarge and significant changes in the composition and behavior of the global atmosphere have emphasized the need for global measurements of atmospheric constit-uents.Examples are (i)the precipitous loss of Antarctic (WMO 1995)and Arctic stratospheric ozone (O 3)(New-Corresponding author address:Dr.Heinrich Bovensmann,Institute of Environmental Physics,University of Bremen (FB1),P .O.Box 330440,D-28334Bremen,Germany.E-mail:bov@gome5.physik.uni-bremen.deman et al.1997;Mu ¨ller et al.1997)resulting from the tropospheric emission of chlorofluorocarbon com-pounds (CFCs,halones,and HFCs)(WMO 95);(ii)the global increase of tropospheric O 3(WMO 1995);(iii)the observed increase of tropospheric ‘‘greenhouse gas-es’’such as CO 2,CH 4,N 2O,and O 3(IPCC 1996);and (iv)the potential coupling between polar stratospheric ozone loss and increased greenhouse gas concentrations (Shindell et al.1998).To assess the significance of such changes a detailed understanding of the physical and chemical processes controlling the global atmosphere is required.Similarly knowledge about the variability and temporal behavior128V OLUME56J O U R N A L O F T H E A T M O S P H E R I C S C I E N C E Sof atmospheric trace gases is necessary to test the pre-dictive ability of the theories currently used to model the atmosphere.Consequently,the accurate assessment of the impact of current and future anthropogenic ac-tivity or natural phenomena on the behavior of the at-mosphere needs detailed knowledge about the temporal and spatial behavior of several atmospheric trace con-stituents(gases,aerosol,clouds)on a global scale,in-cluding the troposphere.Over the past two decades pioneering efforts have been made by the scientific community to establish both ground-based networks and satellite projects that will eventually result in an adequate global observing sys-tem.Examples of satellite borne elements of such pro-grams are the Solar Backscatter Ultraviolet(SBUV)and Total Ozone Mapping Spectrometer(TOMS)on NASA’s Nimbus-7satellite(Heath et al.1975);the Stratospheric Aerosol and Gas Experiment(SAGE)(McCormick et al.1979);the Upper Atmosphere Research Satellite (UARS)(Reber et al.1993)with the Microwave Limb Sounder(MLS),the Halogen Occultation Experiment (HALOE),the Cryogenic Limb Array Etalon Spectrom-eter(CLAES),and the Improved Stratospheric and Me-sospheric Sounder(ISAMS)instruments on board;and the Second European Remote Sensing satellite(ERS-2), which carries the Global Ozone Monitoring Experiment (GOME)(Burrows et al.1999).In the near future,sev-eral new missions will be launched and will contribute significantly to research in thefields of atmospheric chemistry and physics:NASA’s Earth Observing System (EOS)satellites EOS-AM and EOS-CHEM,the Japa-nese Advanced Earth Observing System(ADEOS),and the European Space Agency’s(ESA)Environmental Satellite(ENVISAT).The Scanning Imaging Absorption Spectrometer for Atmospheric Chartography(SCIAMACHY)is part of the atmospheric chemistry payload onboard ENVISAT being prepared by ESA.Following the call for earth observation instrumentation in the Announcement of Opportunity for the Polar Platform issued by ESA,the SCIAMACHY proposal(Burrows et al.1988)was sub-mitted to ESA by an international team of scientists led by Principal Investigator J.P.Burrows.After peer re-view SCIAMACHY was selected as part of the payload for the satellite now known as ENVISAT,which is planned to be launched in2000.The heritage of SCIAMACHY(Burrows et al.1988) lies in both the ground-based measurements using Dif-ferential Optical Absorption Spectroscopy(DOAS) (Brewer et al.1973;Platt and Perner1980;Solomon et al.1987)and previous satellite atmospheric remote sensing missions.SCIAMACHY combines and extends the measurement principles and observational modes of the nadir scattered sunlight measuring instruments SBUV and TOMS(Heath et al.1975),the solar occul-tation instrument SAGE(McCormick et al.1979;Maul-din et al.1985),and the limb scattered sunlight mea-suring instrument Solar Mesospheric Explorer(SME)(Barth et al.1983)within one instrument.SCIAMA-CHY measures in the wavelength range from240nm to2380nm the following:R The scattered and reflected spectral radiance in nadir and limb geometry,R the spectral radiance transmitted through the atmo-sphere in solar and lunar occultation geometry,and R the extraterrestrial solar irradiance and the lunar ra-diance.Limb,nadir,and occultation measurements are planned to be made during every orbit.Trace gases,aerosols, clouds,and the surface of the earth modify the light observed by SCIAMACHY via absorption,emission, and scattering processes.Inversion of the radiance and irradiance measurements enables the amounts and dis-tributions of a significant number of constituents to be retrieved from their spectral signatures and is discussed in section4.Figure1shows the wavelength range to be observed by SCIAMACHY and the position of spec-tral windows where atmospheric constituents are to be retrieved.SCIAMACHY and GOME,which is a small-scale version of SCIAMACHY(see Burrows et al.1999and references therein),represent a new generation of space-based remote sounding sensors,which rely on and uti-lize the simultaneous spectrally resolved measurement of light upwelling from the atmosphere to determine amounts of atmospheric constituents.Using data from GOME,which was launched on board the European Remote Sensing satellite ERS-2in April1995,the feasibility of the instrument and retrieval concepts have been successfully demonstrated for nadir observations.The trace gases O3,NO2,BrO,OClO, SO2,and H2CO have been observed as predicted(Bur-rows et al.1999),and studies of ClO,NO,and aerosol retrieval are proceeding.The determination of O3profile information,including tropospheric O3,from GOME measurements(Burrows et al.1999;Munro et al.1998; Rozanov et al.1998)has a large number of potential applications.In addition,the retrieval of tropospheric column information of SO2,H2CO,NO2,and BrO from GOME measurements was demonstrated(Burrows et al. 1999).The goal of this paper is to provide a comprehensive overview of the SCIAMACHY mission and instrument, to summarize the retrieval strategies,to report on planned data products and expected data quality,and to demonstrate the range of applications and the potential that lies in the concept of this new generation of hy-perspectral UV–VIS–NIR sensors.Section2provides details about the targeted constituents.In section3the instrument design and observational modes are pre-sented.The proposed retrieval strategies are summa-rized in section4.Section5focuses on the expected data precision and section6summarizes the current sta-tus of operational data products.15J ANUARY 1999129B O V E N S M A N N E T A L.F IG .1.Wavelength range covered by SCIAMACHY and absorption windows of the targeted constituents.2.Scientific objectives and targeted constituents The main objective of the SCIAMACHY mission is to improve our knowledge of global atmospheric change and related issues of importance to the chemistry and physics of our atmosphere (cf.WMO 1995and IPCC 1996)such asR the impact of tropospheric pollution arising from in-dustrial activity and biomass burning,R exchange processes between the stratosphere and tro-posphere,R stratospheric chemistry in the polar regions (e.g.,un-der ‘‘ozone hole’’conditions)and at midlatitudes,and R modulations of atmospheric composition resulting from natural phenomena such as volcanic eruptions,solar output variations (e.g.,solar cycle),or solar pro-ton events.Figure 2lists the constituents targeted by SCIA-MACHY and shows the altitude where measurements are to be made.In Fig.2,the combined use of nadir and limb measurements is assumed to yield tropospheric amounts of the constituents down to the ground or the cloud top,depending on cloud cover.a.Tropospheric chemistrySCIAMACHY will measure the backscattered sun-light that reaches the earth’s surface (Ն280nm).The retrieval of tropospheric constituents is influenced and limited by clouds.SCIAMACHY is the only atmo-spheric chemistry sensor on ENVISAT capable of de-termining trace gases and aerosol abundances in the lower troposphere including the planetary boundary lay-er under cloud-free conditions.From the SCIAMACHY nadir and limb measurements tropospheric columns of O 3,NO 2,BrO,CO,CH 4,H 2O,N 2O,SO 2,and H 2CO (cf.Fig.2)will be retrieved.In addition,surface spectral reflectance,aerosol and cloud parameters (cover and cloud-top height),and the tropospheric flux from 280to 2380nm will be retrieved.These data are required for studies of the oxidizing capacity of the troposphere,photochemical O 3production and destruction,and tro-pospheric pollution (biomass burning,industrial activ-ities,aircraft).b.Stratosphere–troposphere exchangeFor the investigation of stratosphere–troposphere ex-change (Holton et al.1995)SCIAMACHY measure-130V OLUME 56J O U R N A L O F T H E A T M O S P H E R I C S C I E N C ES F IG .2.Altitude ranges of atmospheric constituents targeted by SCIAMACHY.Retrieval from the occultation measurements yields infor-mation over a wider altitude range than the limb measurements,due to its higher S/N ratio.ments of the height-resolved profiles of the tracers O 3,H 2O,N 2O,CH 4,and aerosol will be of primary sig-nificance.These measurements enable investigations of the downward transport of stratospheric O 3and upward transport of important species (e.g.,aerosol,CH 4,H 2O,and N 2O).The CH 4and N 2O molecules are emitted into the planetary boundary layer.Their long tropospheric lifetime results in being transported to the stratosphere,where they are the dominant source of the ozone-de-stroying HO x and NO x radicals.Studies of relatively small-scale features such as tropopause folding at mid-latitudes require a high spatial resolution and are un-likely to be unambiguously observed by SCIAMACHY .However,larger-scale stratosphere–troposphere ex-change as envisaged by Holton et al.(1995)will be readily observed.In the neighborhood of the tropopause the different measurements modes of SCIAMACHY will have dif-ferent vertical and horizontal resolutions.Solar and lu-nar occultation modes yield measurements with a ver-tical resolution of 2.5km and a horizontal resolution of 30km across track,determined by the solar diameter,and extending roughly 400km along track.For the limb measurements the geometrical spatial resolution is ap-proximately 3km vertically and typically 240km hor-izontally across track,determined by scan speed and integration time,and extending roughly 400km along track (see Table 3).More details about the geometricalresolution of the different measurement modes will be given in section 3b.c.Stratospheric chemistry and dynamicsThe study of the stratospheric chemistry and dynam-ics will utilize the simultaneous retrieval of total col-umns from nadir measurements and vertical stratospher-ic profiles from limb and occultation measurements of O 3,NO 2,BrO,H 2O,CO,CH 4,and N 2O (and OClO and possibly ClO under ozone hole conditions),as well as aerosol and stratospheric cloud information.Tem-perature and pressure profiles can be determined from limb and occultation observations of the well-mixed gases CO 2and O 2assuming local thermal equilibrium.SCIAMACHY will be making measurements when halogen loading of the stratosphere maximizes around the turn of the century (WMO 1995).It has recently been pointed out by Hofmann (1996)that the springtime polar lower-stratospheric O 3,specifically the layer from 12to 20km,will be the first region to show a response to the international control measures on chlorofluoro-carbon compounds (CFCs)defined in the Montreal Pro-tocol of 1987and its Copenhagen and London amend-ments.SCIAMACHY will enable this preposition to be studied in detail.In general,SCIAMACHY measurements will yield detailed information about the development of strato-15J ANUARY 1999131B O V E N S M A N N E T A L .spheric O 3above the Arctic and Antarctica,the global stratospheric active halogen species (BrO,ClO,OClO),and the global O 3budget as a function of the height in the atmosphere.As SCIAMACHY measures simulta-neously the backscattered radiation field and constituent profiles,an important objective is to test the accuracy of current stratospheric photochemical models and their predictive capability.d.Mesospheric chemistry and dynamicsIn the upper stratosphere and lower mesosphere SCIAMACHY measurements yield profiles of O 3,H 2O,N 2O,NO,O 2,and O 2(1⌬).These measurements will be used to study the distribution of H 2O and O 3and their global circulation.There has recently been much dis-cussion of upper-stratospheric and mesospheric chem-istry in the context of the ‘‘ozone deficit problem’’(Crutzen at al.1995;Summers et al.1997).It has also been suggested that monitoring of H 2O in the lower mesosphere may offer an opportunity for the early de-tection of climate change (Chandra et al.1997).The O 3destruction by mesospheric and upper-stratospheric NO will be investigated.Finally,the mesospheric source of stratospheric NO x will be quantified.In contrast to the retrieval of the majority of trace gases from SCIAMACHY data,NO and O 2(1⌬)profiles are to be determined from their emission features rather than their absorptions.Satellite measurements of NO via the ␥-band emission had been demonstrated by SME to determine profile information from the limb scan (Barth et al.1983,1988)and by SBUV to determine column amounts above 45km from nadir measurements (McPeters 1989).NO can be detected above 40km via the emission from the excited A 2⌺ϩstate into the ground state X 2⌸1/2,3/2(NO ␥-band transitions,200–300nm)as determined in a model sensitivity study by Frederick and Abrams (1982).SCIAMACHY will be able to detect several bands in the 240–300-nm spectral region of the ␥-band emissions of NO in limb as well as in nadir observation mode.O 2(1⌬)can be detected using its emission around 1.27m as shown by results from the SME (Thomas et al.1984).The combination of height-resolved O 3,O 2(1⌬),and UV radiance products from SCIAMACHY provides detailed information about the photolysis of O 3in the upper stratosphere and mesosphere.This will provide an excellent opportunity to test our current photochem-ical knowledge of the mesosphere.e.Climate researchFor use in climate research,SCIAMACHY measure-ments will provide the distributions of several important greenhouse gases (O 3,H 2O,CH 4,N 2O,and CO 2),aero-sol and cloud data,surface spectral reflectance (280–2380nm),the incoming solar spectral irradiance and the outgoing spectral radiance (240–2380nm),and pro-files of p and T (via O 2and CO 2).As it is intended that SCIAMACHY observations are to be made for many years,this long-term dataset will provide much unique information useful for the study of the earth–atmosphere system and variations of the solar output and its impact on climate change.To reach continuity with other spec-trometers measuring solar spectral irradiance such as SBUV or GOME,it is foreseen that SCIAMACHY will be calibrated with standard methods also applied to the GOME or SBUV calibration (Weber et al.1998).3.The instrumentDetails of the instrument concept and design have been given by Burrows and Chance (1991),Goede et al.(1994),Burrows et al.(1995),and Mager et al.(1997).The design is summarized in the following sub-sections.Since the development of the design of SCIA-MACHY two significant changes have occurred.1)The original concept (Burrows and Chance 1991;Burrows et al.1995)used an active Stirling cooler to maintain the infrared detectors of SCIAMACHY at their operational temperature of 150K.During the development phase it was found that a passive cooler could be used for this purpose.This has the advantage of reducing the electrical power con-sumption and potentially extending the lifetime of the mission.2)As an outcome of phase B studies an additional sev-enth polarization measurement device (PMD),mea-suring the 45Њcomponent of the incoming radiance,was added to the spectrometer,to improve the ra-diometric accuracy for the limb mode.a.Design and performanceThe SCIAMACHY instrument is a passive remote sensing moderate-resolution imaging spectrometer.It comprises a mirror system,a telescope,a spectrometer,and thermal and electronic subsystems.A schematic view of the light path within the instrument is depicted in Fig.3.The incoming radiation enters the instrument via one of three ports.1)For nadir measurements the radiation from the earth’s scene is directed by the nadir mirror into a telescope (off-axis parabolic mirror),which focuses the beam onto the entrance slit of the spectrometer.2)For limb and solar/lunar occultation measurements the radiation is reflected by the limb (elevation)mir-ror to the nadir (azimuth)mirror and then into the telescope,which focuses the beam onto the entrance slit of the spectrometer.3)For internal and subsolar calibration measurements the radiation of internal calibration light sources or the solar radiation is directed by the nadir mirror into the telescope.Except for the scan mirrors,all spectrometer parts are132V OLUME 56J O U R N A L O F T H E A T M O S P H E R I C S C I E N C ES F IG .3.Schematic view of the SCIAMACHY optical layout.All imaging optical components (mirrors,redirecting prisms,lenses,etc.)areomitted.All used gratings are in a fixed position.Each detector contains a 1024-pixel photo diode array.fixed and the spectra are recorded simultaneously from 240to 1750nm and in two smaller windows,1940–2040nm and 2265–2380nm,in the near-infrared.The solar radiance varies by a factor of about 100between 240and 400nm.In comparison,the earthshine radiance varies approximately four orders of magnitude over the same spectral range.Spectrometers that measure these quantities therefore need to suppress well any stray light within the instrument.The SCIAMACHY spectrometer achieves this by the combination of a predispersing prism and gratings.This is equivalent in principle to a double spectrometer design.Initially light from the spectrometer slit is collimated and directed onto the pre-dispersing prism.The main beam of light leaving the predispersing prism forms a spectrum in the middle of the instrument.Reflective optics are used to separate the spectrum into four parts.The shorter wavelengths of the spectrum are directed to channel 1(240–314nm)and channel 2(314–405nm)respectively.The majority of the light in the spectrum (405–1750nm)passes without reflection to channels 3–6.The infrared part of the spec-trum (1940–2380nm)is reflected toward channels 7and 8.Dichroic mirrors are used to select the wavelength ranges for channels 3,4,5,and 6,and to separate light for channel 7from that for channel 8.Each individual channel comprises a grating,transmission optics,and a diode array detector.The grating further disperses the light,which is then focused onto eight linear 1024pixel detector arrays.To minimize detector noise and dark current,the diode arrays are cooled:the detector for channels 1and 2to 200K,those for channels 3–5to 235,that for channel 6to 200K,and those for channels 7and 8to 150K.The entire instrument is cooled to 253K in order to minimize the infrared emission from the instrument that might influence the detectors of channels 6–8.In channels 1–5the detectors are silicon monolithic diode arrays (EG&G Reticon RL 1024SR).For the NIR channels 6to 8InGaAs detectors were developed by Epitaxx,Inc.(Joshi et al.1992),and space qualified specifically for SCIAMACHY (see, e.g.,Goede et al.1993;van der A et al.1997).The spectral and radiometric characteristics of the SCIAMACHY spectrometer are summarized in Table 1.The spectral resolution of the spectrometer varies be-tween 0.24and 1.48nm depending on channel number (see Table 1).For DOAS retrieval (see section 4)a high spectral stability is required.The instrument is designed to have a spectral stability of 1/50of a detector pixel,which requires a temperature stability of the spectrom-eter of better than 250mK over one orbit in combination with dedicated calibration measurements.The second relevant retrieval strategy (see section 4),the Full Re-trieval Method (FURM)based on optimal estimation (Rodgers 1976),requires in addition to high spectral stability a high radiometric accuracy of the SCIAMA-CHY measurements.Knowledge of the state of polar-ization of the incoming light and the polarization re-sponse of the instrument determines the radiometric ac-curacy of the radiance,irradiance,and higher-level data products.To achieve the required radiometric accuracy15J ANUARY1999133B O V E N S M A N N E T A L.T ABLE1.Optical parameters of the spectrometer from the designanalysis.ChannelSpectralrange(nm)Resol-ution(nm)Stability(nm)High-resolution channels 1234240–314309–405394–620604–8050.240.260.440.480.0030.0030.0040.005 5678785–10501000–17501940–20402265–23800.541.480.220.260.0050.0150.0030.003Polarization measurement devices PMD1PMD2PMD3PMD4310–377450–525617–705805–900broadbandbroadbandbroadbandbroadband PMD5PMD6PMD71508–16452265–2380802–905broadbandbroadbandbroadbandRadiometric accuracy2–4%Ͻ1%absoluterelativeof2%–4%(depending on the spectral region),dedicated on-ground and in-flight radiometric calibration mea-surements have to be performed in combination with measurements of the polarization properties of the at-mosphere.For the latter purpose SCIAMACHY is equipped with seven polarization measurement devices. Six of these devices(PMD1–6)measure light polarized perpendicular to the SCIAMACHY optical plane,gen-erated by a Brewster angle reflection at the second face of the predispersing prism.This polarized beam is split into six different spectral bands,as described in Table 1.The spectral bands are quite broad and overlap with spectral regions of channels2,3,4,5,6,and8.The PMDs and the light path to the array detectors(including the detectors)have different polarization responses. Consequently,the appropriate combination of PMD data,array detector data,and on-ground polarization calibration data enables the polarization of the incoming light for the nadir measurements(Kruizinga et al.1994; Frerick et al.1997)to be determined.For atmospheric limb measurements,where both limb and nadir mirrors are used,the light is off the optical plane of the spec-trometer.This requires the measurement of additional polarization information of the incoming light.A sev-enth PMD(PMD7)will therefore measure the45Њcom-ponent of the light extracted from the channels3–6light path,as depicted in Fig.3.All PMDs are read out every 1/40s and they observe the same atmospheric volume as channels1–8.In addition to these PMD data being used for the determination of the polarization charac-teristics of the incoming light,they are also planned to be used to determine the fractional cloud cover of the observed ground scene.Additional information about the polarization of the incoming light can be obtained from the diode array overlap regions1/2(309–314nm),2/3(394–405nm), 3/4(604–620nm),4/5(785–805nm),and5/6(1000–1050nm).The polarization efficiency is different for the measurements of the same wavelength in the dif-ferent channels.Inversion of these measurements yields the ratio of plane to parallel polarization components of the incoming light in a manner similar to that used for the array and PMD detectors.The advantage of the over-lap regions is that they are in small wavelength bands, having the same spectral resolution as the corresponding channel.SCIAMACHY aims to retrieve trace gas amounts of relatively weak absorbers.For example,the dif-ferential optical density due to the BrO absorption around350nm detected with GOME(Burrows et al. 1999)is in the order of10Ϫ3and below.Therefore, to achieve a high retrieval precision,a high signal-to-noise ratio(S/N)is required for the scattered ra-diance as well as for the solar irradiance and lunar radiance from the UV to the NIR.The predicted in-strumental S/N values as a function of wavelength are depicted in Fig.4.These S/N values are calculated for an individual detector pixel,for example,of nadir, limb,and occultation measurements.In most cases the predicted S/N is well above103.Exceptions are found in channels1,7,and8.In channel1S/N de-creases toward the UV primarily because the sun is weaker and ozone absorption increases strongly from 320to250nm.In the IR channels7and8the lower S/N values arise from the higher noise of the InGaAs detectors.For these channels the S/N is limited by the detector noise.The apparent missing S/N in Fig. 4c for channel1is the result of the almost complete absorption of the solar photons by the ozone layer when observing the tangent height of15km.In gen-eral,higher S/N values can be obtained by averaging measurements either temporally or spectrally at the cost of losing temporal(and consequently spatial)or spectral resolution.This strategy enables the optimal set of radiance and irradiance data to be generated for a given inversion.Summation of succeeding mea-surements on board(so-called onboard co-adding)is to be used to match optimally the amount of down-linked data to the ENVISAT data rate allowed for SCIAMACHY.In order to cope with the large dy-namic range of the input signals(limb scattered ra-diance vs solar irradiance),which is of six to eight orders of magnitude,the exposure time of each chan-nel can be selected independently over a wide range of values from0.03125to80s.In addition,an ar-rangement involving an aperture stop and a neutral densityfilter is used to limit the intensity of the in-coming light during solar occultation measurements. To optimize S/N over the orbit,exposure times are varied as a function of the solar zenith angle.To calibrate the instrument inflight and to monitor the instrument performance,SCIAMACHY is equipped with a Pt/Cr/Ne hollow cathode(spectral calibration),a。
Atmospheric Chemistry and PhysicsAtmospheric chemistry and physics are two branches of science that are closely related to each other. The atmosphere, which consists of different gases, plays a vital role in regulating the Earth's climate. The study of atmospheric chemistry and physics helps us understand the chemical and physical processes that occur in the atmosphere. This, in turn, helps us understand how human activities affect the Earth's climate.Atmospheric ChemistryAtmospheric chemistry deals with the chemical reactions that take place in the atmosphere. These chemical reactions involve different gases, such as carbon dioxide, methane, and ozone. Many of these gases are greenhouse gases, which means they contribute to the warming of the Earth's atmosphere.One of the most important chemical reactions that occur in the atmosphere is the formation of ozone. Ozone is a molecule made up of three oxygen atoms. It is present in two layers of the Earth's atmosphere: the stratosphere and the troposphere. In the stratosphere, ozone plays a vital role in protecting the Earth's surface from harmful ultraviolet radiation. In the troposphere, however, ozone is a pollutant that can harm human health and the environment.Other important chemical reactions that occur in the atmosphere include the production of acid rain and the formation of smog. Acid rain is formed when sulfur dioxide and nitrogen oxides react with water vapor in the atmosphere. This can damage ecosystems and harm human health. Smog is formed when pollutants, such as nitrogen oxides and volatile organic compounds, react with sunlight. This can cause respiratory problems for humans and damage plant life.Atmospheric PhysicsAtmospheric physics deals with the physical properties of the atmosphere. This includes the study of atmospheric dynamics, which is the study of how air moves in theatmosphere. Atmospheric dynamics is important because it affects weather patterns and climate.Another important area of atmospheric physics is the study of radiation and energy transfer in the atmosphere. The Earth's atmosphere absorbs and reflects both incoming and outgoing radiation. This plays a crucial role in regulating the Earth's temperature. Atmospheric physics helps us understand how this radiation and energy transfer occurs.The study of atmospheric physics also includes the study of clouds. Clouds are an important part of the Earth's climate system because they reflect sunlight and trap heat. The study of clouds helps us understand how they form, how they affect the Earth's climate, and how they are affected by climate change.ConclusionIn conclusion, atmospheric chemistry and physics are two branches of science that are crucial for understanding the Earth's climate. The study of atmospheric chemistry helps us understand the chemical reactions that occur in the atmosphere, which in turn helps us understand how human activities affect the Earth's climate. The study of atmospheric physics helps us understand the physical properties of the atmosphere, including atmospheric dynamics, radiation and energy transfer, and the formation and behavior of clouds. Together, atmospheric chemistry and physics provide a comprehensive understanding of the Earth's climate system and the role that humans play in affecting it.。
ClickHereforFullArticleLong GPS coordinate time series:Multipath and geometry effectsMatt A.King1and Christopher S.Watson2Received16April2009;revised1October2009;accepted21October2009;published3April2010.[1]Within analyses of Global Positioning System(GPS)observations,unmodeledsubdaily signals propagate into long‐period signals via a number of different mechanisms.In this paper,we investigate the effects of time‐variable satellite geometry andthe propagation of a time‐constant unmodeled multipath signal.Multipath reflectors atH=0.1m,0.2m,and1.5m below the antenna are modeled,and their effects on GPScoordinate time series are examined.Simulated time series at20global IGS sites for2000.0–2008.0were derived using the satellite geometry as defined by daily broadcastorbits.We observe the introduction of time‐variable biases in the time series of up toseveral millimeters.The frequency and magnitude of the signal is dependent on sitelocation and multipath source.When adopting realistic GPS observation geometriesobtained from real data(e.g.,including the influence of local obstructions and hardwarespecific tracking),we observe generally larger levels of coordinate variation.In thesecases,we observe spurious signals across the frequency domain,including very highfrequency abrupt changes(offsets)in addition to secular trends.Velocity biases of morethan0.5mm/yr are evident at some sites.The propagated signal has noise characteristicsthat fall between flicker and random walk and shows spectral peaks at harmonics ofthe draconitic year for a GPS satellite(∼351days).When a perfectly repeating synthetic constellation is used,the simulations show near‐negligible time correlated noisehighlighting that subtle variations in the GPS constellation can propagate multipathsignals differently over time,producing significant temporal variations in time series.Citation:King,M.A.,and C.S.Watson(2010),Long GPS coordinate time series:Multipath and geometry effects,J.Geophys. Res.,115,B04403,doi:10.1029/2009JB006543.1.Introduction[2]Geophysical interpretation of GPS coordinate time series most commonly involves the determination of rates of crustal motion and amplitudes and phases of periodic motions caused by seasonal mass loads(for example, atmospheric,hydrological,or oceanic)[Dong et al.,2002; van Dam and Wahr,1998].Each geophysical parameter of interest derived from long GPS time series is biased at some level by residual systematic error,particular those that manifest as long‐period spurious trends or(quasi)periodic signals.In the case of GPS,these systematic errors and their propagation into GPS time series are not yet all well understood.Recent literature has highlighted that long‐period systematic errors may occur in coordinate time series due to one of two primary mechanisms.[3]First,spurious signals may occur directly due to unmodeled long‐period signals,such as satellite antenna modeling errors that propagate differently as the satellite constellation changes[Ge et al.,2005].The second origi-nates in the presumption that all subdaily signals are mod-eled at the observation level,either completely within the functional model or through partial mitigation from the sto-chastic model(e.g.,elevation‐dependent weighting).How-ever,this is not yet the case and residual subdaily(systematic) errors remain which have been shown to propagate into time series[e.g.,Penna et al.,2007]due to a combination of dif-ferent mechanisms[e.g.,Stewart et al.,2005].[4]One well‐studied example of how subdaily signals propagate into longer‐period signals is the case of an un-modeled subdaily tidal signal.Unmodeled semidiurnal and diurnal signals at certain tidal periods have been shown to propagate into fortnightly,semiannual and annual periods [Penna and Stewart,2003;Stewart et al.,2005]with admit-tances exceeding100%in some cases[Penna et al.,2007]. The propagated signal was shown to be highly sensitive to input signal frequency,coordinate component of the un-modeled signal and site location.King et al.[2008]showed that,even after modeling for solid earth tides and ocean tide loading displacements,substantial signal at subdaily periods remained in GPS coordinate time series and these propagate into annual and semiannual periods in24h solutions with median amplitudes of∼0.5mm,but reaching several milli-meters at several sites.As indicated,studies of seasonal geophysical loading phenomena[e.g.,Blewitt et al.,2001;1School of Civil Engineering and Geosciences,Newcastle University,Newcastle upon Tyne,UK.2Surveying and Spatial Science Group,School of Geography andEnvironmental Studies,University of Tasmania,Hobart,Tasmania,Australia.Copyright2010by the American Geophysical Union.0148‐0227/10/2009JB006543JOURNAL OF GEOPHYSICAL RESEARCH,VOL.115,B04403,doi:10.1029/2009JB006543,2010B044031of23Wu et al.,2003]are therefore adversely affected,as are esti-mates of linear velocity,by the presence of such systematic errors.[5]Systematic errors in least squares solutions,such as used in GPS data analysis,propagate according to the least squares design matrix.Therefore,the propagation mecha-nism is controlled by the observation geometry(as defined by the receiver location(s),satellite constellation and local obstructions)plus the chosen parameterization of the solu-tion.Parameters typically include(but are not limited to) site coordinates,adjustments to tropospheric zenith delay, atmospheric gradients,phase ambiguity terms(when not fixed to integers)and receiver clocks(depending on the data differencing approach adopted).Temporal variations in the observation geometry or the number or type of parameters estimated will therefore likely produce temporal variations in the propagated signal.[6]It is notable,therefore,that even in the trivial yet unrealistic case of the GPS antenna location and number of estimated parameters remaining constant with time,the GPS satellite constellation is constantly evolving as satellites are commissioned and decommissioned,or removed from the solution due to satellite eclipse or maneuver.Furthermore, site specific obstructions such as vegetation or man‐made structures may change with time,producing a further change in the observation geometry.The consequence is that,even if an unmodeled signal remains completely constant in time, the way in which it will propagate is likely to change with time.If these temporal variations are suitably large,then the ensuing systematic error will likely bias the GPS time series significantly,resulting in erroneous interpretation of geo-physical signals such as tectonic velocity,glacial isostatic adjustment,vertical motion of tide gauges or seasonal geophysical loading signals.Furthermore,such errors would degrade the GPS contribution to the International Terrestrial Reference Frame[Altamimi et al.,2007].[7]In this paper we test this hypothesis,using one source of presently unmodeled subdaily signal:carrier phase mul-tipath,taking“multipath”to mean signal reflections from planar surfaces not part of the antenna itself,as adopted by Georgiadou and Kleusberg[1988].Errors of this type are similar in nature to antenna phase center variation mis-modeling and antenna imaging(changes in the antenna phase pattern induced by conducting material in the vicinity of the antenna)[Georgiadou and Kleusberg,1988],in that they exhibit an elevation dependency.We do not consider these additional sources of systematic error here,yet simply note that the propagation mechanism is likely to be some-what similar.Early studies investigating multipath found very small effects on geodetic time series,leading to the thought that multipath effects are mitigated by averaging when sufficient observational time spans are used[e.g., Davis et al.,1989;Lau and Cross,2007a].However,years of experience with continuous GPS time series,combined with analysis advances leading to improved signal‐to‐noise ratio of long time series,have shown that GPS time series are highly sensitive to GPS hardware changes(including receiver,receiver firmware and antenna),suggesting that multipath may be playing an important role;yet the mech-anism for this is not well understood.As a contribution to understanding the effect of carrier phase multipath on multiyear GPS coordinate time series,we perform several trials using simulated and real data,using various simulated carrier phase multipath signals.[8]We commence in section2by introducing the adopted model of carrier phase multipath that we use throughout this paper to perturb a multiyear,simulated GPS coordinate time series.The simulation approach is introduced(section2.2) before detailing the different satellite constellation config-urations adopted to assess the different characteristics of the multipath propagation(section2.3).First,in order to show the influence of time variability on the propagation,we start with a theoretical constellation that has a fixed orbit repeat time,a constant number of satellites and no obstructions above the elevation mask at each site(section 2.3.1). Second,we adopt the same clear horizon but with a realistic time variable satellite constellation taken from the broadcast orbits,(section2.3.2).Finally,we detail the most realistic constellation that includes site specific time variable changes to the observation tracking as observed in real data (for example hardware changes and physical obstructions on the horizon,section2.3.3).[9]In section3we compare and discuss the simulated time series generated using the three constellation config-urations,in addition to investigating the effect of changes to the adopted functional and stochastic models within the evolving constellation scenario.In section4,we pro-vide a comparison against time series computed using real data in a PPP approach with GIPSY5software[Webb and Zumberge,1995],and in section5we present two possible mitigation strategies that involve novel weighting strategies of the input observations in order to minimize the time‐and geometry‐dependent propagation of the multipath signal.2.Simulations2.1.Signal Multipath[10]The space surrounding an antenna may be subdivided into three regions including the reactive near field in the region nearest the antenna,the radiating near field and the far field out to infinity[e.g.,Balanis,2005].The boundaries of these regions are not sharply defined although criteria have been developed in order to delineate them in practice. Given an antenna with maximum dimension D,and signal with wavelength l,and for D>l[Balanis,2005]defines the first and second boundaries occurring at distancesR1¼0:62ffiffiffiffiffiffiD2rand R2¼2D2from the antenna surface,respectively.For a GPS choke ring antenna with D=0.38m,l L1=0.19m and l L2=0.24m, then R1L1=0.03m,R1L2=0.02m,R2L1=1.52m and R2L2=1.20m.In this paper we consider multipath sources solely within the radiating near field,or distances in the range ∼0.03m to∼1.2–1.5m.For an examination of the effect of reactive near field sources on GPS time series,see Dilssner et al.[2008],and for the effect of phase center modeling errors on site velocities,see Steigenberger et al.[2009]. [11]To date there is no widely accepted model for multi-path,partly due to the complexity of real world GPS antenna environments.One relatively simple model shown to approx-imate observed multipath has been described by EloseguiKING AND WATSON:MULTIPATH AND GEOMETRY EFFECTS ON GPS B04403 B044032of23et al.[1995],based on the earlier work of Georgiadou and Kleusberg [1988]and Young et al.[1985].This model is based on the assumption that multipath is caused solely by a horizontal reflector at some height,H ,below the GPS antenna phase center causing an attenuation,a ,of the signal voltage amplitude.So,for a satellite with elevation angle,",phase bias d L (in units of meters)of the L1or L2phase due to multipath may be modeled as [Elosegui et al.,1995]L ¼ 2 tan À1 sin 4 Hsin "!1þ cos 4 Hsin "!0B B @1C C A ð1Þwhere l is the carrier phase wavelength for the L1or L2carrier phase signal.Since most geodetic GPS positioning is performed using the ionosphere ‐free linear combination (LC)of the raw carrier phase observables (L1and L2),the LC phase delay can be computed asLC L "; ;H ðÞ%2:5457Â L "; ;H ; L 1ðÞÀ1:5457Â L "; ;H ; L 2ðÞð2Þ[12]We show in Figure 1(left)d LLC for a range of heights above the reflector (H =0.1,0.2,and 1.5m)and for two different attenuation values (a =0.05and 0.1).As noted by Elosegui et al.[1995],the effect of changing the height above the reflector is to change the rate at which d L LC varies with elevation (increased rate of variation with increasing height).Within a reasonable range of attenuation values (0.01to 0.1),a change in attenuation has the effect ofapproximately uniformly scaling the bias as seen by com-paring Figure 1(left).[13]Despite this model providing a useful approximation,it is based on geometric ray optics which is not appropriate in the radiating near field.In addition,values of d L LC will be modified by the antenna gain pattern.To improve on this,we adopt a model as developed by T.A.Herring (personal communication,2009,hereafter denoted HMM)that extends the Elosegui et al.[1995]model through the use of Fresnel equations relating to electric field amplitudes,and attempting to take into account antenna gain properties:L ¼2 tan À1a sin 4 H sin "ÂÃg d þa cos 4 Hsin "ÂÃ!ð3Þwith the antenna gain pattern,consisting of the direct (g d )andreflected (g r )gain,modeled using a simple modified dipole model expressed as a function of rate of change (G )of the antenna gain with signal zenith angle,such thatg d ¼cos z =G ðÞg r ¼cos 90=G ðÞ1Àsin "ðÞðÞalso,a =Sg r R a is the amplitude of the reflected signal for agiven surface roughness (S ),withR a ¼n 1cos z Àffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffin 22Àn 1sin z ðÞ2q n 1cos z þffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffin 22Àn 1sin z ðÞ2q 264375being the Fresnel equation for an electric field perpendicularto the plane ofincidence.Figure 1.Carrier phase bias due to reflectors using two different multipath models.(left)The Elosegui et al.[1995]model with two different attenuation values,a ,at H =0.1m (gray),H =0.2m (blue),and H =1.5m (magenta),sampled every 1°.(right)The HMM for the same values of H and with two different surface reflectivity terms,S ,sampled every 1°.Np is the refractive index of the reflective medium.The dotted black line (bottom right)is for S =1.0for H =0.16m.KING AND WATSON:MULTIPATH AND GEOMETRY EFFECTS ON GPS B04403B044033of 23[14]This depends on refractive indices n 1and n 2,with n 1=1for air and n 2appropriate for the reflective medium.In practice this can be derived from the square root of tabulated dielectric constants,and for concrete this is typi-cally taken as 4(n 2=2).[15]Here d LLC may be formed analogously to equation (2).[16]For the purposes of this paper we adopt values which approximate the amplitude of the signal in GPS phase residuals (T.A.Herring,personal communication,2009),namely S =0.3,G =1.1and n 2=2,although these are not definitive,and indeed some sites may require a different value of S ,as we shall demonstrate.We show in Figure 1(right)d L LC for S =0.3and 0.5based on HMM.The ampli-tude of the modeled signals scales approximately linearly with variations in S .Increasing either G to 1.2or n 2to 100result in increased model signal amplitude of about 50%.Only small changes in frequency or phase occur.Two dif-ferences to the model of Elosegui et al.[1995]can be identified.First,the signal in HMM is substantially reduced at high elevations,which is in general agreement with the pattern seen in GPS carrier phase residuals.Second,both amplitude and frequency of the HMM is proportional to H ,whereas the Elosegui et al.[1995]model does not show the same sensitivity of amplitude to variations in H (compare Figure 1,left versus right).2.2.Simulation Approach[17]Our GPS simulator is based on the undifferenced GPS observable,and is conceptually similar to that of Santerre [1991].Observation ‐level biases,such as d L LC ,are entered via the “observed minus computed ”term (b )of the batch least squares adjustment:^x ¼A T WA ÀÁÀ1A T Wbð4Þ[18]The design matrix,A ,is defined by the partialderivatives of the functional model with respect to the parameters.The effect of the unmodeled signal on theseparameters is reflected in the estimated values ^x .W is the inverse of the observation variance ‐covariance matrix.[19]Initial station coordinates are taken from the ITRF2005coordinate set and satellite positions are taken from one of three orbit configurations (section 2.3).Satellite orbits and clocks are assumed known and fixed.Unless otherwise specified below,estimated parameters are the corrections to initial station coordinates,tropospheric zenith delays,receiver clocks (normally every epoch)and realvalued ambiguity terms (normally one per satellite pass).In these simulations,nonzero values of the estimated parameters represent parameter bias.Correct integer ambiguity fixing may be simulated by simply removing these terms from the design matrix [Santerre ,1991].[20]This simulator has previously been used by King et al.[2003]to study propagation of tidal signals in sub-daily coordinate estimates;the observed biases in output time series were accurately reproduced in the simulator.We have also verified the simulator by reproducing the propa-gated periodic biases observed in the GIPSY solutions of Penna et al.[2007]to within very small errors.We therefore assume the simulator is capable of reproducing the effects of systematic errors on real GPS solutions (that use an undifferenced observation strategy),but with the advantage of controlling all systematic error sources.[21]For the tests described here,we used a typical 24h “observation ”session and estimated adjustments to topo-centric station coordinates (north,east and up)once per day (i.e.,once per session),tropospheric zenith delay parameters once per hour,and used an elevation cutoff angle of 7°.By default we applied uniform weighting to all observations and include a single real value ambiguity parameter per satellite pass (we assess the impact of elevation ‐dependent weighting and ambiguity fixing throughout section 3).We used one measurement epoch every 300s (to reduce the computational burden)and considered data over the years 2000–2007inclusive for a sample of 20sites (Figure 2)in the International GNSS Service (IGS)[Dow et al.,2005].[22]By fixing the satellite orbits,clocks and earth orien-tation parameters to predetermined values we are potentially oversimplifying the problem.However,since multipath isFigure 2.Site locations from the IGS network used in this study.KING AND WATSON:MULTIPATH AND GEOMETRY EFFECTS ON GPS B04403B044034of 23heights above reflectors and reflector conditions are all site specific)they are unlikely to systematically bias parameters with large spatial scales,and hence we argue that a site ‐by ‐site analysis is reasonable.This remains to be verified in further work.[23]Our simulations represent a substantial advance on the work of Elosegui et al.[1995],which was limited in part by the available computational power at the time.First,they performed a simplified adjustment,considering mainly a one ‐dimensional (vertical)coordinate with,at most,one tropospheric zenith delay term per day and they did not examine the effects of ambiguity fixing.Second,they report on simulations for only a small number of (noncontinuous)days.Third,they reported on the propagation at only a single midlatitude site.Fourth,they used a multipath model not entirely appropriate to the near field.Each of these dif-ferences affects the least squares design matrix (or b in the last case),and hence could alter the propagation into long ‐term GPS time series.2.3.Satellite Orbit Configurations2.3.1.Clear Horizon:Constant Constellation[24]In order to provide a non time variable reference orbit,we adopted an artificial configuration that used a fixed 24satellite constellation,each with a fixed orbit repeat period of 24h minus 246s,close to the average GPS satellite orbital repeat time [Agnew and Larson ,2007].To achieve this constellation,we adapted the “perfect GPS orbit ”scenario as developed in a simulator implemented by Penna and Stewart [2003].The horizon at each site was assumed clear for this configuration;that is,satellites are observed without ob-struction down to the elevation cutoff angle of 7°.2.3.2.Clear Horizon:Evolving Constellation[25]To investigate the propagation impacts of a more realistic time variable orbit in our simulations,this “evolving constellation ”configuration adopted orbits from the daily broadcast ephemerides,with satellites set as unhealthy eliminated from the simulations.Again,this configuration assumed no obstructions above the elevation mask at each site.[26]The time variability of this configuration can be clearly seen in Figure 3(blue)where we show constellation statistics as a function of elevation over time.The minimum,maximum and standard deviation of epoch by epoch satellite elevations (computed in daily bins)are shown for a repre-sentative set of sites.Changes in the blue lines reflect the changing constellation geometry over time.The mini-mum elevation angle is governed by the 7°elevation cut-off used in the simulations.The maximum elevation angle and standard deviation of all elevation angles are governed by the actual satellite constellation as observed from each site location.[27]The temporal variability of the satellite geometry in Figure 3is striking.MAW1exhibits a strong secular term in maximum satellite elevation over time,although this is a feature of all high ‐latitude sites to some extent.More equatorial sites exhibit high ‐frequency variation in maxi-mum elevation and low ‐frequency fluctuations in elevation standard deviation.A clear trend in elevation standard de-viation between 2000and 2002is also a common feature to many sites.2.3.3.Obstructed Horizon:Evolving Constellation [28]The final orbit configuration assessed within the simulations takes into consideration time variable changes that occur due to the influence of local tracking issues,hardware changes and site obstructions above the elevation cutoff angle that are also potentially time variable and highly site specific.We derived this constellation by using only those satellite observations that are used in analyses of real world GPS data in a conventional GIPSY PPP analysis [Webb and Zumberge ,1995;Zumberge et al.,1997].Con-stellation statistics from this constellation (shown in brown,Figure 3)show considerable variation from the clear hori-zon,constant constellation configuration discussed previ-ously.Notable differences are observed at sites such as MCM4and DUBO which exhibit both secular and clear quasi ‐periodic characteristics.3.Simulation Output3.1.Clear Horizon:Constant Versus Evolving Constellation[29]Figure 4shows the simulation output highlighting the propagation of the HMM d L LC into the height compo-nent of three representative sites using H =0.1m,H =0.2m,and H =1.5m,for both the ambiguity free and fixed scenarios.The effect of a time ‐varying constellation may be seen by comparing the constant and evolving constellation scenarios (Figures 4(left)and 4(right),respectively).The three sites are NTUS,an equatorial site;POTS,a midlatitude site;and MCM4,a high ‐latitude site.Summary statistics including the offset,slope and RMS (white noise)are pro-vided for all sites in Table 1.[30]On first inspection,the most noticeable difference between these simulations is a clear shift from a stationary time series with high ‐frequency periodic variability for the constant constellation,toward clear low ‐frequency vari-ability present within the evolving constellation output.Also present in the later output are occasional transient events.For example,the H =1.5m NTUS ambiguity fixed example (Figure 4)shows what could be interpreted as a significant offset in the coordinate time series during late 2000.Closer examination of the time series shows that the offset is ∼1mm,and occurs over no more than 2–3days,and pos-sibly as quickly as from one day to the next.This offset is purely an artifact of the propagation of the systematic error introduced into the simulations as a function of the time variable satellite constellation.[31]In relation to the north and east coordinate compo-nents,we observe a similar pattern of variability yet ap-proximately an order of magnitude smaller in magnitude.This finding agrees with intuition given the simulated multipath signal is azimuthally symmetrical and hence one would expect the effects on north and east components to be relatively small.In the main,we show below only Figures 4–11pertaining to the height component,including the results for north and east components in the supple-mentary material where warranted.[32]Returning to the vertical component,and now con-sidering the height bias (offset in Table 1),we observe biases reaching just over 7mm for all tested values of H ,and note that the offset is approximately equal for the constant and evolving constellation cases.Biases for H =0.2m tendMULTIPATH AND GEOMETRY EFFECTS ON GPS B044035of 23to be smaller,but the mean bias will change depending on elevation cutoff angle (see below).Table 1confirms the output for the constant constellation has a zero trend com-puted over the entire period (2000.0–2008.0).The evolving constellation case shows trends reaching ±0.16mm yr −1insome instances;the effect over shorter periods could be sub-stantially larger,as is suggested by Figure 4,and it will also scale with input multipath amplitude.[33]Considering next the variability of the simulation output,(RMS in Table 1),we observe a progressiveincreaseFigure 3.Daily elevation angle statistics generated using the broadcast constellation (blue)and those observations actually used in a GIPSY PPP analysis (brown).Units are degrees.Note the different scales for each plot.KING AND WATSON:MULTIPATH AND GEOMETRY EFFECTS ON GPS B04403B044036of 23in scatter from H =0.1m to H =1.5m,and from constant to evolving constellations (mean values of 0.01,0.07and 0.15mm for the constant case,and 0.04,0.21and 0.49mm for the evolving case,H =0.1,0.2,and 1.5,respectively).The RMS shows a small yet interesting correlation with latitude for the evolving constellation case,with the maxi-mum RMS present near the equator,reducing to a minimum at approximately ±50°,and subsequently increasing toward the poles.This pattern is expressed across all values of H ,increasing significantly from H =0.1to H =1.5.[34]The difference between the ambiguity free and fixed cases is interesting for both the constant and evolving constellation scenarios.In the time domain (Figure 4),there appears to be a clear reduction in the energy of the periodic components within the time series when moving from the ambiguity free to fixed cases.To assist in the analysis of these periodic components and potentially time correlated noise structures,we compute global stacks of power spectra using the Lomb ‐Scargle periodogram [Scargle ,1982]as implemented according to Press et al.[1992]with an over-sampling factor of 4.We repeated the analysis with an oversampling factor of 1and found noin the expression of various harmonic peaks present in the spectra.In order to aid comparison between different simulations,prior to stacking in the frequency domain,the power spectra for each time series are normalized by an arbitrary variance of 1mm 2for each coordinate component.We generate our stacked spectra by computing the median normalized power estimate across each frequency.[35]The stacked spectra for the constant constellation (Figure 5,left)confirm that these time series are composed of a stationary series with a complex harmonic comb of energy superimposed.The energy within this harmonic comb is maximum for the H =0.2m case,and a minimum for the H =0.1m scenario.It is interesting that dominant peaks occur at harmonics of a period close to one year.On closer inspection,these harmonics are in fact multiples of 351.4days,equivalent to the so ‐called GPS satellite draconitic year (abbreviated here as dy)which defines the average repeat period of the GPS satellites in inertial space relative to the sun.In a study of stacked spectra taken from the weekly solutions of the IGS,Ray et al.[2008]noticedFigure 4.Height bias time series for three sites when propagating d L LC with H =0.1m (gray/black),H =0.2m (blue)and H =1.5m (magenta)for (left)constant constellation and (right)evolving constellation.Colors match those in Figure 1.Ambiguity ‐free (lighter shades)and ambiguity ‐fixed (darker shades)solutions are shown.KING AND WATSON:MULTIPATH AND GEOMETRY EFFECTS ON GPS B04403B044037of 23。
郑州市2023 年高中毕业年级第二次质量预测英语试题卷本试卷分四部分,考试时间120 分钟,满分150 分(听力成绩算作参考分)。
考生应首先阅读答题卡上的文字信息,然后在答题卡上作答,在试题卷上作答无效。
第二部分阅读理解(共两节,满分40 分)第一节(共15 小题;每小题2 分,满分30 分)阅读下列短文,从每题所给的A、B、C 和D 四个选项中,选出最佳选项,并在答题卡上将该项涂黑。
AHigh school programs at the National Gallery of Art value depth over breadth, exploring original works of art through a single specific question or theme.Studio WorkshopsSingle museum visit, 2.5 hoursGrades 9-12These half-day art workshops include an in-depth examination and discussion of works of art in the galleries, followed by behind-the-scenes access to the Education Studio, where students create a related art project.National Gallery of Art educators will encourage students to look carefully at works of art and then share their responses and come up with theories based on their observations.Students will have the chance to create a work of art in the studio inspired by what they have seen in the galleries.Museum MakersExploring Art and MuseumsGrades 11-12The Museum Makers program explains how museums operate and what they have to offer. It gives upper-level high school students the tools to experience, understandand interpret art. Participants will gain an insider's view of how an art museum works.Students meet for seven Saturday sessions from 10:00 a.m. to 3:00 p.m. Completion of the program requires attendance of all seven sessions.Creative WritingGrades 7-12, 90 minutesStudents will provide a voice for their personal responses to art through creative writing while they are looking at a selection of artworks in the galleries. Using close observation, group discussion and personal reflection, they will be guided through exercises that use different writing forms, including free-form poetry.A maximum of 30 students ( minimum of 15) will be accepted at each session.21.What can students do at Studio Workshops?A.Get basic training as an artist.B.Put forward their own art theories.C.Discuss with artists about their works.D.Learn about artists, inspiration for their works.22.What can students get from Museum Makers?A.Methods of how artworks are created.B.Experience in running an art museum.C.Knowledge about how an art museum works.D.Academic credits for completing the program.23.Which group can attend Creative Writing at a session?A.15 college students.B. 10 Grade 9 students.C. 25 Grade 10 students.D. 40 high school students.BI never imagined that someone telling me I looked skinny would anger me. And yet, I was made very angry when a colleague pinched (捏)my waist and screamed, “Rosa, you've lost weight. You look great In The truth is that I was tired and not taking care of myself. I decided to start a proper weight-loss program.The first to go would be road rage (路怒).I am in far less control of this weight than any other. Every time something gets in my path, I fly off the handle. I need to lose the road rage, and fast! No, no more speed. Instead, I now repeat the words:"I am not in a hurry. " This year, I will drive safely, allowing "stupid" to happen all aroundme. From that, I hope to gain patience.Next is guilt. When guilt drives my conscience to do better, it's functional. But when it presents itself as an internal dialogue that goes nowhere, it's useless. This year, I want to stop feeling guilty for not keeping a cleaner house, for spending time away from my children to be with friends, for not attending every party because I’d rather be at home, or for watching TV when I should be reading. My image and performance are not at the front of anyone else’s mind but my own. From this, I hope to gain freedom to be myself.The last is fear. Fear has held me back. Fear of failure has prevented me from being a writer. Fear of embarrassment has prevented me from giving an opinion. Fear of rejection has stopped me from aiming higher in my life. Fear of regret has led me into situations that made me uncomfortable. If I can lose any one of these fears, I stand to gain experience.So, if I can lose the rage, shake off some guilt, and take fear off my plate, I stand to gain patience, freedom, and experience. Pound for pound I have not lost a thing but I will be much lighter. Next time, I hope my colleague looks me in the eye to see my glow instead of pinching part of me that has nothing to do with how great I really look.24.What does the underlined part “fly off the handle” in paragraph 2 probably mean?A.Pick up speed.B. Drive off.C. Feel quite nervous.D. Become very angry.25.What has made the author feel guilty before?A.Attending too many parties.B.Reading much with her children.C.Wasting her time in watching TV.D.Spending little time with her friends.26.How has fear affected the author?A.It has prevented her achieving her goals.B.It has stopped her furthering her education.C.It has made it difficult for her to make friends.D.It has caused her to quit her dream of being a teacher.27.What does the author expect to gain from her weight-loss program?A.Respect.B. Independence.C. Optimism.D. Tolerance.CSix months before she died, my grandmother moved into an old people's home and I visited her there. The room was clean and warm, and the care assistants were kind and cheerful. A general knowledge quiz show was on the television, and the only other sound was snoring (打呼噜).People moved only when they needed to be helped to the bathroom. It was disappointing. Grandmother talked a lot about how much she missed seeing her grandchildren, but I knew from my sister that they hated going to visit her there.So I was interested to read a newspaper article about a new concept in old people's homes in France. The idea is simple, but revolutionary - combining a residential home for the elderly with a nursery school in the same building. The children and the residents eat lunch together and share activities. In the afternoons, the residents enjoy reading or telling stories to the children, and if a child is feeling sad or tired, there is always a kind lap to sit on.The advantages are huge for everyone concerned. The children are happy, because they get a lot more individual attention. The residents are happy because they feel useful and needed. And the staff are happy because they see an improvement in the physical and psychological health of the residents and have an army of assistants to help with the children.Nowadays there is less and less contact between the old and the young in an increasing number of countries. There are many reasons for this, including the breakdown of the extended family, working parents with no time to care for ageing relations, families that have moved away, and smaller flats with no room for grandparents. But the result is the same - increasing numbers of children without grandparents and old people who have no contact with children, and more old people who are lonely and feel useless, along with more and more families with young children who desperately need more support. Ifs a major problem in many societies.That's why intergenerational programs, designed to bring the old and the young together, are growing in popularity all over the world.28.What does the underlined word "residents" in paragraph 2 probably refer to?A.Old people.B. School teachers.C. Assistants.D. Staff.29.How were the old people at the home the author's grandmother was in?A.They felt lonely and useless.B.They weren't allowed to be visited.C.They weren't looked after properly.D.They lived in a dirty and uncomfortable room.30.What does the author think is a major problem in many societies today?A.The extended family is broken down.B.There isn't much room for grandparents.C.Working parents have no time to care for their children.D.There isn't much contact between the old and the young.31.What will be probably talked about later in the passage?A.Advice on how to communicate with children.B.Plans for setting up more homes for old people,C.Examples of successful intergenerational programs.D.Ways of teaching entertainment skills to old people.DMicroplastic pollution is increasing greatly around the globe, according to a study of plastic particles (微粒)carried in the air.People are already known to breathe, drink and eat microplastics, and research suggests that pollution levels will continue to rise rapidly. The researchers said that breathing in these particles can be harmful to lung tissue and lead to serious diseases.Professor Natalie Mahowald, at Cornell University in the US and part of the research team, said, " But maybe we could solve this before it becomes a huge problem, if we manage our plastics better, before they accumulate in the environment and move around everywhere."The research, published in the journal Proceedings of the National Acadenry of Sciences, examined airborne (空气传播的)microplastics, which have been far less studied than plastics in oceans and rivers.The team gathered more than 300 samples of airborne microplastics from 11 s ites across the western US. These were the basis for atmospheric modeling that estimated the contribution from different sources (来源),and it was the first such study to do so.They found that roads were the main factor (因素)in the western US, linked to about 85% of the microplastics in the air. These are likely to include particles from tiresand brake pads on vehicles, and plastics from litter that had been broken down.The researchers extended their modeling work to a global level and this suggested that while roads are also likely to be the major driver of airborne plastics in Europe, South America and Australia, plastic particles blown up from fields may be a much bigger factor in Africa and Asia.Professor Andreas Stohl of the University of Vienna's Faculty of Earth Sciences said, “The study confirms the global-scale (全球规模的)nature of microplastic transport in the atmosphere and does a good job in highlighting highly relevant and concerning possibilities, but more measurement data is needed to get a better idea of the sources.”32.What can be known about microplastic pollution from this text?A.The particles can do great harm to our lungs.B.Airborne microplastics have been widely studied.C.It has become the most pressing environmental problem.D.There is less plastic in the air than in oceans and rivers.33.What did the researchers find out about microplastic pollution?A.Its results differ across many continents.B.Africa and Asia are suffering most from it.C.Roads and fields are largely to blame for it.D.It spreads fast from one continent to another.34.What should the researchers do next according to Professor Andreas Stohl?A.To predict the potential damage of microplastics.B.To understand the nature of microplastic pollution.C.To improve the method of collecting samples of microplastics.D.To collect more data to understand the sources of microplastics.35.What can be the best title for this text?A.Effects of microplastics on human healthB.Microplastic pollution on the global scaleC.Possible solutions to microplastic pollutionD.Microplastic pollution rising quickly in the air第二节(共5 小题;每小题2 分,满分10 分)根据短文内容,从短文后的选项中选出能填入空白处的最佳选项,并在答题卡上将该项涂黑。
在核心期刊发表的英文Publishing in core journals is a significant achievement for researchers and scholars. It not only enhances their academic reputation but also contributes to the advancement of knowledge in their respective fields. In this document, we will discuss the importance of publishing in core journals and provide some practical tips for achieving this goal.First and foremost, publishing in core journals allows researchers to reach a wider audience and have a greater impact on their field of study. Core journals are highly regarded within the academic community and are often the first choice for scholars seeking high-quality research articles. By publishing in these journals, researchers can ensure that their work is read and cited by other experts in their field, thereby increasing their visibility and influence.Furthermore, publishing in core journals is a mark of academic excellence. It demonstrates that the research has undergone rigorous peer review and has been deemed to be of high quality and significance. This can enhance the researcher's reputation and credibility within the academic community, leading to increased opportunities for collaboration, funding, and career advancement.In order to publish in core journals, researchers must adhere to high standards of scholarship and research integrity. This includes conducting original and rigorous research, presenting their findings clearly and concisely, and adhering to ethical guidelines and best practices in their field. Additionally, researchers must be familiar with the specific requirements and preferences of the core journals they are targeting, such as formatting guidelines, citation styles, and submission processes.To increase the likelihood of publishing in core journals, researchers should carefully select the journals to which they submit their work. It is important to target journals that are well-regarded in the researcher's field and have a high impact factor. Researchers should also consider the scope and focus of the journal, ensuring that their work aligns with the journal's editorial priorities and audience.In addition, researchers should pay close attention to the feedback and suggestions provided by peer reviewers and editors. Constructive criticism and guidance from experts in the field can help researchers improve their work and increase its chances of acceptance in core journals. Researchers should be open to revising and refining their work based on the feedback received, demonstrating their commitment to producinghigh-quality research.In conclusion, publishing in core journals is a significant achievement for researchers and scholars. It provides an opportunity to reach a wider audience, enhance academic reputation, and contribute to the advancement of knowledge in their field. By adhering to high standards of scholarship and research integrity, carefully selecting target journals, and responding to feedback from peer reviewers and editors, researchers can increase their chances of publishing in core journals and making a meaningful impact in their field.。
advances in earth sciences影响因子Advances in Earth Sciences (AES) is a reputable, peer-reviewed journal that publishes high-quality research articles spanning various topics in the field of Earth sciences. As the advancement of Earth sciences are crucial for a better understanding of our planet and addressing the challenges it faces, the impact factor of AES plays a significant role in determining the influence and reach of the journal. In this article, we will delve into the factors that contribute to the impact factor of AES, starting from its inception to the present day.The impact factor of a scientific journal is a measure of the average number of citations received per article published in that journal within a specific time period. It offers a quantitative assessment of the quality and influence of a journal within its field. For AES, the impact factor is a reflection of the importance and significance of the research articles it publishes, as well as the wide readership and citation rate it attracts.In the early years of AES, the impact factor was relatively low as the journal was still in its infancy. However, over time, with the increasing recognition and credibility of the journal, the impactfactor experienced a steady rise. This can be attributed to various factors, which we will now explore in more detail.Firstly, the quality of the research articles published in AES is a crucial determinant of its impact factor. The journal follows a rigorous peer-review process, ensuring that only high-quality and groundbreaking research is accepted for publication. By maintaining high standards, AES attracts renowned scientists and researchers to submit their work, thereby increasing the overall caliber and impact of the journal.Secondly, the relevance and timeliness of the research topics covered in AES also contribute to its impact factor. As the field of Earth sciences constantly evolves and addresses emerging global challenges, AES strives to publish cutting-edge research in areas such as climate change, natural disasters, geological processes, and environmental management. By focusing on these pressing issues, AES attracts the attention of researchers, policymakers, and other scientists who actively cite the published articles, thereby increasing the impact factor.Furthermore, the accessibility and international reach of the journalplay a vital role in determining its impact. AES aims to ensure that its publications are widely available and easily accessible to the global scientific community. By utilizing online platforms and open-access options, AES maximizes the visibility and dissemination of its articles. This broad accessibility increases the chances of articles being cited by researchers from various disciplines and geographical locations, positively influencing the impact factor.Another factor that enhances the impact factor of AES is its commitment to fostering collaborations and interdisciplinary research. Earth sciences inherently encompass a wide range of disciplines like geology, geography, atmospheric sciences, and hydrology. By encouraging interdisciplinary approaches and publishing articles that bridge different fields, AES attracts a diverse readership and consequently increases the likelihood of citations, thereby raising its impact factor.Lastly, the reputation and standing of the editorial board and reviewers of AES contribute significantly to its impact factor. AES is comprised of renowned scientists and experts in the field, who bring their expertise and knowledge to the review and selectionprocess. The reputation of the editorial board as well as the stringent review process assures the scientific community of the reliability and credibility of the articles published in AES. This confidence leads to higher citations and ultimately increases the impact factor of the journal.In conclusion, the impact factor of Advances in Earth Sciences (AES) is influenced by various factors that contribute to its recognition and influence within the scientific community. These factors include the quality of research articles, the relevance of the topics covered, the accessibility and reach of the journal, its interdisciplinary approach, and the reputation of the editorial board and reviewers. By continuously nurturing these factors, AES maintains a high impact factor, thereby promoting the advancement and dissemination of knowledge in the field of Earth sciences.。
Distribution of sulfur and pyrite in coal seams from Kutai Basin(East Kalimantan,Indonesia):Implications for paleoenvironmental conditionsSri Widodo a ,⁎,Wolfgang Oschmann b ,Achim Bechtel c ,Reinhard F.Sachsenhofer c ,Komang Anggayana d ,Wilhelm Puettmann eaDepartment of Mining Engineering,Moslem University of Indonesia,Jln.Urip Sumoharjo,Makassar,Indonesia bInstitute of Geosciece,J.W.Goethe-University,Altenhöferallee 1,D-60438Frankfurt a.M.,Germany cDepartment of Applied Geoscience and Geophysics,University of Leoben,Peter-Tunner-Str.5,A-8700Leoben,Austria dDepartment of Mining Engineering,Bandung Institute of Technology,Jln.Ganesa 10,I-40132Bandung,Indonesia eInstitute of Atmospheric and Environmental Sciences,Dapartment of Analytical Enviromental Chemistry,J.W.Goethe-University,Altenhöferallee 1,D-60438Frankfurt a.M.,Germanya b s t r a c ta r t i c l e i n f o Article history:Received 12August 2009Received in revised form 29November 2009Accepted 3December 2009Available online 13December 2009Keywords:Kutai Basin Pyrite SulfurFramboidal Ombrogenous TopogenousThirteen Miocene coal samples from three active open pit and underground coal mines in the Kutai Basin (East Kalimantan,Indonesia)were collected.According to our microscopical and geochemical investigations,coal samples from Sebulu and Centra Busang coal mines yield high sulfur and pyrite contents as compared to the Embalut coal mine.The latter being characterized by very low sulfur (b 1%)and pyrite contents.The ash,mineral,total sulfur,iron (Fe)and pyrite contents of most of the coal samples from the Sebulu and Centra Busang coal mines are high and positively related in these samples.Low contents of ash,mineral,total sulfur,iron (Fe)and pyrite have been found only in sample TNT-32from Centra Busang coal mine.Pyrite was the only sulfur form that we could recognize under re flected light microscope (oil immersion).Pyrite occurred in the coal as framboidal,euhedral,massive,anhedral and epigenetic pyrite in cleats/fractures.High concentration of pyrite argues for the availability of iron (Fe)in the coal samples.Most coal samples from the Embalut coal mine show lower sulfur (b 1wt.%)and pyrite contents as found within Centra Busang and Sebulu coals.One exception is the coal sample KTD-38from Embalut mine with total sulfur content of 1.41wt.%.The rich ash,mineral,sulfur and pyrite contents of coals in the Kutai Basin (especially Centra Busang and Sebulu coals)can be related to the volcanic activity (Nyaan volcanic)during Tertiary whereby aeolian material was transported to the mire during or after the peati fication process.Moreover,the adjacent early Tertiary deep marine sediment,ma fic igneous rocks and melange in the center of Kalimantan Island might have provided mineral to the coal by uplift and erosion.The inorganic matter in the mire might also originate from the ground and surface water from the highland of central Kalimantan.©2009Elsevier B.V.All rights reserved.1.IntroductionMost of the inorganic matter in coal is present as minerals which are dispersed throughout the coal macerals.Individual grains of minerals vary largely in size from less than one micrometre to tens or hundreds of micrometres.Sometimes mineral-rich layers are even thick enough to be visible on the coal surface (Taylor et al.,1998).Mineral components in the coals were classi fied in three groups according to their origin (Stach et al.,1975):(1)Mineral from the original plants;(2)mineral that formed during the first stage of the coali fication process or which was introduced by water and wind into the later coal deposits;and (3)mineral deposited during the second phase of the coali fication process,after consolidation of the coal,byascending or descending solutions in cracks,fissures,or cavities or by alteration of primarily deposited minerals.The dominant mineral of coals is usually composed of sul fides,clay,carbonates,and quartz and sometimes additional phosphates,heavy minerals,and salts as minor contributions to inorganic matter of coal.In most coals,sul fides are preferentially composed of pyrite and marcasite but pyrite is in general dominating by far (Balme,1956;Mackowsky,1943).Sul fides can be categorized as either syngenetic (primary),early-diagenetic or epigenetic (secondary)in origin.During peati fication,syngenetic or early-diagenetic fine-crystalline or fine-concretionary pyrite appears,commonly in the form of framboids.Syngenetic pyrite formed during accumulation of the peat and/or during early (humi fication)processes,and is usually small in size,and intimately dispersed throughout the coal (Renton and Cecil,1979;Reyes-Navarro and Davis,1976).Occasionally,the cell walls of plant material have been replaced by pyrite (Taylor et al.,1998).Falcon and Snyman (1986)suggest that the accumulation of pyrite in coal might also ariseInternational Journal of Coal Geology 81(2010)151–162⁎Corresponding author.Tel.:+62411454775;fax:+62411453009.E-mail address:srwd007@ (S.Widodo).0166-5162/$–see front matter ©2009Elsevier B.V.All rights reserved.doi:10.1016/j.coal.2009.12.003Contents lists available at ScienceDirectInternational Journal of Coal Geologyj o u r n a l ho m e p a g e :w w w.e l s ev i e r.c o m /l o c a t e /i j c o a l g e ofrom the aeolian andfluviatile import of iron-rich mineral at the time of peat accumulation followed by in-situ precipitation.Epigenetic pyrite is incorporated in the coal after compaction or partial consolidation(Reyes-Navarro and Davis,1976)and is generally much larger(coarse grained)andfills cracks,cleats,and cavities (Renton and Cecil,1979).The formation of epigenetic pyrite is dependent primarily on the availability of reduced sulfur,dissolved cation(ferrous iron)and a suitable site for formation i.e.,cleat (Casagrande et al.,1977;Spears and Caswell,1986;Demchuk,1992). Moreover,epigenetic pyrite might be precipitated from water percolating into fractures,cavities and pores present in coal seams long after accumulation of the peat(Falcon and Snyman,1986).In general,coals deposited in paralic basins contain more pyrite than those in limnic basins.Among the paralic deposits,coal seams which have been influenced by marine transgressions are consistently characterized by a particularly high content of pyrite and sometimes also of organic sulfur,especially in the upper part of the seams(Balme, 1956;Dai et al.,2002;Mackowsky,1943).In sulfur-rich humic coals, pyrite in the form offine grains orfine concretions is particularly common in microlithotypes containing a high proportion of vitrinite; these forms also tend to be common in sapropelic coals.In the absence of other criteria(such as marine fossils or coal balls),a relatively high proportion of synsedimentary or early-diagenetic pyrite can be useful for seam correlation.Many previous investigations(Anggayana et al.,2003;Baruah,1995; Dai et al.,2002,2003,2006,2007,2008;Dai and Chou,2007;Elswick et al.,2007;Frankie and Hower,1987;Kortenski and Kostova,1996; Lόpez-Buendía et al.,2007;Querol et al.,1989;Renton and Bird,1991; Strauss and Schieber,1989;Turner and Richardson,2004;Wiese and Fyfe,1986)have described the characteristics,type,morphology, genesis,and distribution of pyrite in coal seams from different deposits. Our investigation deals with sulfur and pyrite occurrences in the Miocene coal seams from Centra Busang,Sebulu and Embalut coal mines,Kutai Basins,East Kalimantan,Indonesia.The primary purpose of the study is to explain why most Miocene coal seams from Kutai Basin have very low sulfur contents,whereas in some coal seams higher sulfur contents are observed.The second purpose is to identify the types of pyrite and factors affecting their appearance and it's relation with paleoenvironmental conditions during deposition of the coals.2.Geological settingKalimantan is surrounded by marginal basins of South China,Celebes and Sulu seas,microcontinental fragments of south China in the north,and mainland SE Asia(Indochina and peninsular Malaysia)in the west(Moss and Chambers,1999).Kalimantan was interpreted as the product of Mesozoic accretion of oceanic crustal material(ophiolite),marginal basin fill,island arc material and microcontinental fragments onto the Paleozoic continental core of the Schwaner Mountain in the SW of the island(Fig.1; Hall and Nichols,2002;Hutchison,1989;Moss and Wilson,1998;Widodo et al.,2009).During early Tertiary times,Kalimantan formed a promontory of the Sundaland Craton:the stable eastern margin of the Eurasian plate (Hall,1996;Metcalfe,1998).In the east,Kalimantan is separated from Sulawesi by the deep Makassar Basins(Fig.1),formed during Paleogene times(Situmorang,1982).The main areas of eastern, central and northern part of Kalimantan are coated by Tertiary sediments(Fig.1)which were deposited influvial,marginal-marine or marine environments.Tertiary sedimentation in these areas occurred at the same time with,and subsequent to,a period of widespread Paleogene extension and subsidence,which may have commenced in the middle Eocene or earlier(Moss and Wilson,1998).A number of Tertiary sedimentary basins,of which the Kutai Basin is the largest,are identified across Kalimantan(Moss et al.,1997).The island of Kalimantan and in particular the Kutai Basin has experienced a complex tectonic history from the Paleogene to the present day.The Kutai Basin was formed during early Tertiary times and wasfilled-up with clastic sediments progressing from the western to the eastern part of the basin.This basin was subdivided into the Upper Kutai Basin,consisting of Paleogene outcrops with Cenozoic volcanics possessing a strong northwest–southeast struc-tural grain,and the Lower Kutai Basin with Miocene strata cropping out in a north–northeast-trending structure.The Meratus Mountains to the southwest and the Central Kalimantan Mountains to the north of the Kutai Basin have an ophiolitic basement together with Paleogene strata striking dominantly in a north–northeast direction (Clay et al.,2000).The coal mining companies are located in the vicinity of the Mahakam River,Kutai Basin,East Kalimantan Province(Fig.2).The precise geographic position of Sebulu coal mine is S00°26′40.4″/E116°52′54.1″and Centra Busang coal mine is S00°44′22.2″/E116°89′16.6″,whereas the Embalut coal mine is situated S00°33′34.9″/E117°12′15.5″.Centra Busang is located in the Busang village,East Kutai regency and Sebulu coal mine in the Sebulu village,Kutai Kertanegara regency,East Kalimantan province. The Embalut coal mine is located in Embalut village,Kutai Kertanegara regency,East Kalimantan province.Coal seams in the Centra Busang and Sebulu mines were found in the Pulau Balang Formation with Middle Miocene age,and coal seams in the Embalut mine was found in the Pulau Balang(Middle Miocene age)and Balikpapan Formation with Upper Miocene age(Fig.3).Previous studies of the sedimentary evolution of the Kutai Basin, based onfield survey and oil wells,have shown that the Tertiary sequences are broadly regressive in general with a(dominantly) offshore marine argillaceous sequence of Palaeocene age followed by a coal bearing deltaic and coastal plain succession of Miocene age. Shoreline progradation was generally towards the east(Samuel and Muchsin,1976;Rose and Hartono,1978in Land and Jones,1987).According to Supriatna and Rustandi(1986)the Neogene succession in the Kutai Basin includes from bottom to top the following formations:Pamaluan,Bebuluh,Pulau Balang,Balikpapan, Kampung Baru Formation and alluvial.The Pamaluan Formation(up to1500m thick)consists of sandstones with insertion of claystone, shale,limestones,and siltstone.It formed in a deep marine environment during Late Oligocene and Early Miocene times.The Lower Miocene Bebuluh Formation(up to900m thick)consists of limestones with insertion of sandy limestones and argillaceous shales.The deposition of the formation occurred in a shallow sea.The Bebuluh Formation interfingers with the Pamaluan Formation.The Pulau Balang Formation of Early to Middle Miocene age overlies the Bebuluh Beds concordantly. It is composed of graywackes,quartz sandstones,limestones,claystones, dacitic tuff and coal insertion.The thickness of coal seams ranges from 3to4m.The depositional environment can be characterized as deltaic to shallow marine according to Supriatna and Rustandi(1986).The formation is approximately900m thick.The Middle to Upper Miocene Balikpapan Formation(1000–1500m thick)uniformly overlies the Pulau Balang Formation and consists of quartz sandstone,clay with insertion of shale,and coal seams5to10m thick.The deposition of the Balikpapan Formation occurred in a delta environment.The Upper Miocene to Pliocene Kampung Baru Formation disconcordantly overlies the Balikpapan Formation.It is composed of quartz sandstone with insertion of clay,shale,silt and approximately3m thick coal(lignite). The deposition of the Kampung Baru Formation up to900m thick, sedimentary succession occurred in a delta.Alluvium deposit(alluvial) consists of sandy clay and clayey sands.3.Samples and methodCoal samples were collected in-situ with channel sampling method from three active surface mines in the Centra Busang(3samples),152S.Widodo et al./International Journal of Coal Geology81(2010)151–162Sebulu (2samples),and Embalut (8samples)coal mines,Kutai Basin,East Kalimantan,Indonesia.The sample preparation and microscopic examination generally followed the procedures described by Taylor et al.(1998).Coal particles of about 1mm in diameter were used for preparation of polished sections,which were embedded in a silicone mould (diameter 40mm)using epoxy resin as an embedding medium.After hardening,the samples were ground and polished.MicroscopicFig.1.Simpli fied geology of Kalimantan (modi fied from Moss and Wilson 1998;Hall and Nichols,2002).153S.Widodo et al./International Journal of Coal Geology 81(2010)151–162analyses were performed by a single-scan method with a Leica MPV microscope using re flected white and fluorescent light.At least 300points were counted for coal macerals and mineral.Pyrites were counted separately from the other minerals (for example,clay,carbonate,quartz)which were counted together in one group.Total sulfur contents were determined using an automated Leco SC-344carbon sulfur analyzer.A weighed coal sample is mixed with iron chips and a tungsten accelerator and is then combusted in an oxygen atmosphere at 1370°C.The moisture and dust are removed from the combustion product and the SO 2gas is measured by a solid-state infrared detector.Ash yields were determined following standard procedure DIN 51719using dry coal samples.One gram of each coal sample is heated 2h to 815°C (±15°C)in a muffel furnace,the residue was then cooled to room temperature and weighed.For trace element analyses,carried out by ACTLABS,Ancaster,Canada,the coal ash (0.5g)was dissolved in aqua regia (0.5ml H 2O,0.6ml concentrated HNO 3and 1.8ml concentrated HCl).After cooling,samples were diluted to 10ml with deionized water and homogenized.The digestion is near total for base metals,but will only be partial for silicates and oxides.The solutions were then analyzed using a Perkin Elmer OPTIMA 3000Radial ICP-MS for the 30elementsFig.2.Geographic map showing the locations of Sebulu,Centra Busang and Embalut coal mines in the Kutai Basin,East Kalimantan,Indonesia.Fig.3.Sedimentary sequences and the distribution of some fundamental parameters of the Embalut,Centra Busang and Sebulu coal field,Kutai Basin,East Kalimantan,Indonesia.154S.Widodo et al./International Journal of Coal Geology 81(2010)151–162suite.A series of USGS-geochemical standards were used as controls. The iron(Fe)is one of the elements to be discussed in this paper. 4.Results and discussion4.1.Ash,mineral,total sulfur and iron(Fe)contentsAsh yields,minerals,total sulfur,and Fe contents in the coal samples from Centra Busang,Sebulu and Embalut coal mines are summarized in Table1and cross-correlation are shown in Fig.4.The ash yield of the Centra Busang and Sebulu coals varies from1.40to 5.80(wt.%,db).The Embalut coal samples show higher variability in ash yields(1.33to8.98wt.%,db).Four coal samples from Embalut yield more than2wt.%(db)ash(Table1)while ash contents lower than2wt.%(db)are found in the samples KTD-36,KTD-35,KTD-43, and KTD-37(Table1).The ash contents of Cetra Busang and Sebulu coals show a positive correlation to the total sulfur+Fe(r2=0.83). Embalut coals show also a positive correlation between ash contents to the total sulfur+Fe(r2=0.77).The mineral contents of Centra Busang and Sebulu coals vary from 0.70to5.60(vol.%).Mineral is predominantly composed of pyrite.Clay and carbonates are also found in a significant quantity.Quartz is observed only in a trace proportion.In the coal samples from Centra Busang coal mine,mineral contents vary from3.70vol.%(TNT-30,seam BL-4)over1.60vol.%(TNT-31BL-7)to0.70vol.%(TNT-32,BL-7.1).The coal samples from Sebulu coal mine show mineral contents between 2.3vol.%(TNT-33)and5.6vol.%(TNT-34).The mineral contents of Centra Busang and Sebulu coals show a strong positive correlation to the ash contents(r2=0.96;Fig.4).Unlike to the Centra Busang and Sebulu coals,mineral in the Embalut coals is dominated by clay minerals. Carbonate,pyrite and quartz are found in small amounts.The mineral contents of Embalut coals vary from0.30to2.0(vol.%).Most of mineral contents are lower than2vol.%.The mineral contents(from micro-scopical analysis)of Embalut coal mine show a negative correlation to the ash contents in the coals(r2=0.05,Fig.5).The studied coal samples from the Centra Busang and Sebulu coal mines have relatively high total sulfur contents(up to3.15wt.%; Table1).Total sulfur contents show a positive correlation to the ash contents in the samples(r2=0.67;Fig.4).On the other hand,most samples from the Embalut coal mine show very low total sulfur contents(≤0.2wt.%).Sole exception is the coal sample KTD-38from Pulau Balang Formation with a total sulfur content of1.41wt.%.The correlation of total sulfur to the ash contents in the Embalut coal mine also shows a positive correlation(r2=0.84,Fig.5).Iron(Fe)contents in the samples from the Centra Busang and Sebulu coal mines fall within the range of0.17to1.51wt.%.As shown in Table1 and Fig.6,the iron(Fe)contents in the coal samples from Centra Busang and Sebulu coal mines show a strong positive correlation to the pyrite contents and total sulfur.The correlation coefficients are r2=0.96(Fe vs pyrite contents),and r2=0.86(Fe vs total sulfur).Whereas,iron(Fe) contents of Embalut coal samples are low and range from0.15to1.32%. Highest content are found in the sample KTD-38with a value of1.32%. Fig.7also shows a positive correlation of iron(Fe)contents to the pyrite contents and total sulfur in the coal samples from Embalut coal mine. The correlation coefficients are r2=0.97(Fe vs pyrite contents)and r2=0.99(Fe vs total sulfur).The sum of total sulfur and iron(Fe)is close to the amount of ash for Centra Busang and Sebulu coals(Table1).This indicates that a large proportion of the mineral in the coal samples must be composed of pyrite.In the two samples(TNT-30and TNT-33)the sum of the total sulfur content and the Fe content exceeds the ash content.In these samples some of the sulfur must be organically bond sulfur, which evaporates together with the organic matter during heating. The lowest total sulfur content(0.14wt.%)of the Centra Busang and Sebulu coal mines is observed in TNT-32,which provided a very low ash yield of1.40wt.%,db.This allows to estimate the non-pyritic mineral content to approximately1wt.%.4.2.Pyrite content and pyrite types in the coal seamsBased on microscopical analyses,the amount of pyrite in the Centra Busang and Sebulu coal samples varies from0.4to4.3vol.% (Table1).In comparison,most of pyrite contents of Embalut coal samples are lower than in Centra Busang and Sebulu coals and pyrite is only observed in three samples(KTD-36,KTD-43and KTD-38). Differences in the pyrite contents of coals from Centra Busang,Sebulu and Embalut coal mines are reflected by the type of pyrite found(i.e. framboidal,euhedral,massive,anhedral,and epigenetic/syngenetic pyrite in cleats and fractures).4.2.1.Framboidal pyriteFramboidal forms of pyrite are categorized as syngenetic(Dai et al., 2007).Some authors proposed that the type of this pyrite originates from pyritization of sulfur bacteria(Casagrande et al.,1977; Casagrande et al.,1980;Kortenski and Kostova,1996;Lόpez-Buendía et al.,2007;Querol et al.,1989;Renton and Cecil,1979).Kortenski and Kostova(1996)also proposed the possibility of pyritization of other kinds of bacteria which might have coexisted along with the sulfur metabolizing bacteria and supported the decomposition and assim-ilation of the plant tissue.Moreover,it has been suggested that framboidal pyrite might be generated from mineral solutions in inorganic material(Dai et al., 2002,2003;Kortenski and Kostova,1996).Other theories(Wilkin and Barnes,1997)suggested for the formation of framboidal pyrite is in first step the activity of biogenic processes i.e.,pyritic fossilization of bacterial colonies,and in a second step further growing of framboidalTable1Ash yield,minerals(from microscopical analysis),total sulfur,pyrite(from microscopical analysis)and iron(Fe)contents in the studied coal samples.Samples Seams Formations Coal mines Ash(wt.%,db)Mineral(vol.%)Total sulfur(wt.%,db)Pyrite(vol.%)Fe(wt.%)KTD-4221Balikpapan Embalut 2.66 1.30.050.00.18 KTD-3618 1.33 1.60.050.20.15 KTD-3517 1.650.30.060.00.21 KTD-4312 1.61 1.00.070.30.17 KTD-4010 5.050.70.200.00.20 KTD-399 2.81 1.00.100.00.22 KTD-388Pulau Balang8.98 1.0 1.410.7 1.32 KTD-377 1.77 2.00.180.00.20 TNT-30BL-4Pulau Balang Centra Busang 3.50 3.7 3.15 3.0 1.15 TNT-31BL-7.9 1.90 1.6 1.10 1.00.34 TNT-32BL-7.10 1.400.70.140.40.17 TNT-33ST Pulau Balang Sebulu 2.10 2.3 1.69 2.00.51 TNT-34STU 5.80 5.6 2.92 4.3 1.51155S.Widodo et al./International Journal of Coal Geology81(2010)151–162pyrite by organic processes,based on laboratory syntheses over a wide range of thermal conditions.Bacterial framboidal pyrite preserves the independence of the separate globules even when they form aggregates (Kortenski and Kostova,1996;Querol et al.,1989;Renton and Cecil,1979).During peati fication,syngenetic or early-diagenetic fine-crystal-line or fine-concretionary pyrite appears,commonly in the form of framboids.The Centra Busang and Sebulu coal samples contain bacterial framboidal pyrite in high abundance.An example is shown in sample TNT-34(Fig.8a).In some samples such as TNT-30from the Centra Busang coal mine the primary framboidal pyrite shows an overgrowth by secondary pyrite which is generally associated with clay minerals (Fig.8b).In most of the framboidal pyrite globules the crystals are densely intergrown and consist of some aggregates as previously described by Skripchenko and Berberian (1975;in Kortenski and Kostova,1996).Most of bacterial framboidal pyrite in the sample TNT-30appears as single bodies or solitary.In Embalut coal sample framboidal pyrite is not observed.4.2.2.Euhedral pyriteEuhedral pyrite is recognizable as well shaped pyrite crystals (Kortenski and Kostova,1996).Querol et al.(1989)described euhedral pyrite in the coal samples from Maestrazgo Basin,northeastern Spain.Kortenski and Kostova (1996)observed this type of pyrite in the coal samples from Bulgaria and divided euhedral pyrite into isolated and clustered varieties and isolated anhedral crystals and aggregates of euhedral crystals.Most of the euhedral pyrite is syngenetic and is generated during deposition of peat and/or during early humi fication.In general,the crystals of euhedral pyrite are small in size and intimately dispersed throughout the coal (Dai et al.,2007,2008;Renton and Cecil,1979;Reyes-Navarro and Davis,1976;Turner and Richardson,2004).Isolated euhedral pyrite was found only in small amounts insampleFig.4.The variation and cross correlation of ash to the mineral matter,total sulfur,Fe and pyrite contents in the coal samples from Centra Busang and Sebulu coalmines.Fig.5.The variation and cross correlation of ash to the mineral matter,total sulfur,Fe and pyrite contents in the coal samples from Embalut coal mine.156S.Widodo et al./International Journal of Coal Geology 81(2010)151–162TNT-34from Sebulu coal mine (Fig.8c).Clustered euhedral pyrite could not be detected in all analyzed samples.4.2.3.Massive pyriteMassive pyrite is usually found as cleat-/cell-fillings,cementing or coating framboids,euhedral or detrital minerals (Querol et al.,1989).Massive pyrite has also been found as a replacement of organic matter in different macerals (Querol et al.,1989).Many authors denoted pyrite grains with irregular shapes and different sizes by the term massive pyrite (Dai and Chou,2007;Grady,1977;Kortenski and Kostova,1996;Wiese and Fyfe,1986).Renton and Bird (1991)described this type of pyrite as irregular.Massive pyrite was found in most coal samples from Centra Busang and Sebulu.The homogeneous massive pyrite was generally porous and not compact,which is due to the inclusion of relict organic matter and clay minerals during the crystallization processes.Homogeneous massive pyrite was present in lenticular or irregular form.Fig.8d and e show the homogeneous massive pyrite in the coal samples of the Sebulu mine (sample TNT-34).This type of pyrite was not observed in the Embalut coal samples.4.2.4.Anhedral pyriteAnhedral pyrite corresponds to pyrite forms whose shape depends on the shape of the plant debris in which they were deposited.Theanhedral pyrite was divided into two types,the replacement anhedral pyrite and the in filling anhedral pyrite (Kortenski and Kostova,1996;Wiese and Fyfe,1986),which are of late syngenetic and epigenetic origin,respectively.The anhedral pyrite in the sample from Centra Busang and Sebulu coal mines was found in small amounts.Replacement anhedral pyrite was deposited in the lumens of densinite maceral (Fig.8f).The replacement anhedral pyrite was a result of mineralization of cell wall and described to originate from replacement of plant material or massive pyrite replacement of organic matter (Kortenski and Kostova,1996;Querol et al.,1989;Wiese and Fyfe,1986).Anhedral pyrite was not found in the Embalut coal samples.4.2.5.Epigenetic pyrite in cleats and fracturesThe term epigenetic pyrite in cleats and fractures is used for pyrite deposited in fractures or cleats which determine the path of solutions penetrating a coal seam.There are two types of epigenetic pyrite in cleat and fracture:in filling and replacing epigenetic pyrite.The in filling epigenetic pyrite in cleat and fracture has again been divided in two types:fracture and cleat filling (Kortenski and Kostova,1996;Querol et al.,1989;Renton and Cecil,1979).Epigenetic pyrite in cleats and fractures is observed only in a very small quantity in the studied coal samples.In filling epigenetic pyrite in cleats and fractures was observed in the Centra Busang coal samples (e.g.,TNT-30,Fig.8f).TheFig.6.Cross correlation of the concentration ratios of Fe (wt.%)to the pyrite content (vol.%)and total sulfur (wt.%)in the coal samples from Centra Busang and Sebulu coalmines.Fig.7.Cross correlation of the concentration ratios of Fe (wt.%)to the pyrite content (vol.%)and total sulfur (wt.%)in the coal samples from Embalut coal mine.157S.Widodo et al./International Journal of Coal Geology 81(2010)151–162pyrite grains were deposited in the fractures and cleats of the coals and are of epigenetic origin.Replacing epigenetic pyrite from Sebulu coal samples is shown in the Fig.8g.The grains of pyrite replaced (filled)the cell lumens.Replacing epigenetic in fracture pyrite was also found in the Embalut coal samples (Fig.8h).The pyrite from the Embalut coal sample can be characterized as massive fracture filling pyrite with density-fill in the lumens or fractures (Fig.8i).4.3.Interpretation of coal depositional environment and the difference of pyrite in coalsEnrichment of minerals in coal suggests that the peat was deposited under topogenous conditions resulting in increased (richer)nutrient supply (minerals/detrital in flux).In contrast,coals of low mineral content are generated in general from ombrogenous peat bogs with low nutrient supply (minerals).Supply of mineral to the topogenous peat occurs preferentially by surface water and ground water,while in ombrogenous peat bogs the supply takes place through atmospheric deposition.Pyrite is sometimes the dominating inorganic matter present in the coal and the morphology of pyrite can help to reconstruct the environmental condition during and after peatformation (Casagrande et al.,1977,1980;Dai et al.,2008;Demchuk,1992;Taylor et al.,1998;Wiese and Fyfe,1986).The amount and type of pyrite identi fied in the coals from Centra Busang,Sebulu and Embalut differs signi ficantly.In the Centra Busang and Sebulu coals both syngenetic and epigenetic pyrites appear.Syngenetic framboidal pyrite is partly overgrown by epigenetic pyrite.Typically,primary framboidal pyrite occurs in coals and carbonaceous shales overlain by marine strata (Taylor et al.,1998).In case of the Centra Busang and Sebulu coals,the high abundance of mineral might originate from the erosion of Early Tertiary marine sediments of the Central Kalimantan ridge (Fig.9),delivering suf ficient iron and sulphate for pyrite formation under subaquatic conditions (Fig.9).It is assumed that the Centra Busang and Sebulu coals have been evolved from a topogenous peat deposited under wet forest swamp conditions (Fig.9).This interpretation is veri fied by the high proportion of ash,mineral,sulfur,pyrite and iron contents in the Centra Busang and Sebulu coals.In the Embalut coals,pyrite is found only in epigenetic form resulting from post depositional processes.The absence of framboidal pyrite is consistent with the formation of the coal from an ombrogenous peat.This type of peat receives the water through heavy rainfall and the groundwater table lying below the peatformingFig.8.All photos taken under oil immersion.(a)Bacterial framboidal pyrite (Py)from the Sebulu coal mine (TNT-34),(b)overgrowth framboidal pyrite (Py)from the Centra Busang coal mine (TNT-30),(c)euhedral pyrite (Py)from the Sebulu coal mine (TNT-34),(d)massive pyrite (Py)from the Sebulu coal mine (TNT-34),(e)massive pyrite (Py)from Sebulu coal mine (TNT-34),(f)epigenetic pyrite (Py)in cleats/fractures and anhedral pyrite (Py)from the Sebulu coal mine (TNT-30),(g)replacing epigenetic pyrite (Py)in cleats/fractures from the Sebulu coal mine (TNT-30),(h)replacing epigenetic pyrite (Py)in fractures from Embalut coal mine (KTD-38)(i)and replacement massive pyrite (Py)from Embalut coal mine (sample KTD-38).158S.Widodo et al./International Journal of Coal Geology 81(2010)151–162。
文章主要亮点有:①“科学界应重视发表那些负向或失败的研究结果”观点不同寻常,有助于科学界深思;②论证手法多样(对比,因果,引用等);③结构明确,按照“提出观点(第一段)——论证观点(第二至五段)——分析观点所涉错误现象产生的原因(第六段)——总结收篇,重申观点(第七段)”的脉络展开论述。
【原文】ⅠHypothesis-driven research is at the heart of scientific endeavor, and it is often the positive,confirmatory data that get the most attention and guide further research. But many studies produce non-confirmatory data—observations that refute current ideas and carefully constructed hypotheses. And it can be argue d that these “negative data,” far from having little value in science, are actually an integral part of scientific progress that deserve more attention.ⅡAt first glance, this may seem a little nonsensical; after all, how can non-confirmatory results help science to progress when they fail to substantiate anything? But in fact, in a philosophical sense, only negative data resulting in rejection of a hypothesis represent real progress. As philosopher of science Karl Popper stated: “Every refutation should b e regarded as a great success; not merely a success of the scientist who refuted the theory, but also of the scientist who created the refuted theory and who thus in the first instance suggested, if only indirectly, the refuting experiment.”ⅢOn a more practical level, Journal of Negative Results in Biomedicine (JNRBM) was launched on the premise that scientific progress depends not only on the accomplishments of individuals but requires teamwork and open communication of all results—positive and negative. After all, the scientific community can only learn from negative results if the data are published.ⅣThough not every negative result will turn out to be of groundbreaking significance, it is imperative to be aware of the more balanced perspective that can result from the publication of non-confirmatory findings. The first and most obvious benefits of publishing negative results are a reduction in the duplication of effort between researchers, leading to the acceleration of scientific progress, and greater transparency and openness.ⅤMore broadly, publication of negative data might also contribute to a more realistic appreciation of the “messy” nature of science. Scientific endeavors rarely result in perfect discoveries of elements of “truth” about the world. This is largely because they are frequently based on methods with real limitations and hypotheses based on uncertain premises.ⅥIt is perhaps this “messy” aspect of science that contributes to a hesitation within the scientific community to publish negative data. In an ever more competitive environment, it may be that scientific journals prefer to publish studies with clear and specific conclusions. Indeed, Daniele Fanelli of the University of Edinburgh suggests that results may be distorted by a “p ublish or perish” culture in which the progress of scientific careers depends on the frequency and quality of citations. This leads to a situation in which data that support a hypothesis may be perceived in a more positive light and receive more citations than data that only generate more questions and uncertainty.ⅦDespite the effects of this competitive environment, however, a willingness to publish negative data is emerging among researchers. Publications that emphasize positive findings are of course useful, but a more balanced presentation of all the data, including negative or failed experiments, would also make a significant contribution to scientific progress.【词汇短语】1.hypothesis [haɪˈpɒθəsɪs] n.假设2.confirmatory [kən'fɜ:məˌtərɪ] a.证实的,确实的3.refute [rɪˈfju:t] v.驳斥4.negative [ˈnegətɪv] a.否定的5.integral [ˈɪntɪgrəl] a.构成整体所必须的6.nonsensical [nɒnˈsensɪkl] a.荒谬的7.substantiate [səbˈstænʃieɪt] v.证明8.rejection [rɪ'dʒekʃn] n.拒绝,驳回unch [lɔ:ntʃ] v.发射,发起,推出10.premise ['premɪs] n.前提11.accomplishment [əˈkʌmplɪʃmənt] n.成就munity [kəˈmju:nəti] n.团体,界13.groundbreaking [ˈgraʊndbreɪkɪŋ] a.创新的14.significance [sɪgˈnɪfɪkəns] n.意义15.imperative [ɪmˈperətɪv] a.必要的16.duplication [ˌdju:plɪ'keɪʃn] n.双重,重复17.acceleration [əkˌseləˈreɪʃn] n.加速18.transparency [trænsˈpærənsi] n.透明度19.contribute to 增益,有助于20.appreciation [əˌpri:ʃiˈeɪʃn] n.欣赏,评论21.messy [ˈmesi] a.散乱的22.endeavor [ɪn'devə] n.努力23.hesitation [ˌhezɪ'teɪʃn] n.犹豫24.distort [dɪˈstɔ:t] v.扭曲25.perish [ˈperɪʃ] v.毁灭26.frequency [ˈf ri:kwənsi] n.频率27.emerge [iˈmɜ:dʒ] v.出现,显露【翻译点评】Ⅰ①Hypothesis-driven research is at the heart of scientific endeavor, and it is often the positive, confirmatory data that get the most attention and guide further research. ②But many studies produce non-confirmatory data—observations that refute current ideas and carefully constructed hypotheses. ③And it can be argued that these “negative data,” far from having little value in science, are actually an integral part of scientific progress that deserve more attention.翻译:假设驱动型研究对于科学事业至关重要, 且经常是那些正向的、验证性数据最受关注并引导深入研究。
running head是什么是你对论文标题的总结,用最简短的几个单词表述出来..给自己题目用一个很简短的句子或单词表达出来。
就是出版时,在页眉上出现的东西,不能太长,可出用简称,能体现你文章的价值的句子female reproductive toxicity and mechanism of HDovarian toxicity and relevant mechanism of HDcover letterHere within enclosed is our paper for consideration to be published on " Journal of Applied Toxicology ****". The further information about the paper is in the following:The Title:The Authors:This paper is about the ******The authors claim that none of the material in the paper has been published or is under consideration for publication elsewhere. Correspondence and phone calls about the paper should be directed to****at the following address,phone and fax number,and e-mail address:Corresponding author:Tel:Fax:E-mail:Thanks very much for your attention to our paper.Sincerely yours,Dear Editors:On behalf of my co-authors,I am submitting the enclosed material “Apoptosis of Rat Ovarian Granulosa Cells by 2,5-hexanedione in vitro and its Relevant Gene Regulation” for possible publication on “Journal of Applied Toxicology”. The authors claim that none of the material in the paper has been published or is under consideration for publication elsewhere.Thank you very much for your considering our manuscript for potential publication. I'm looking forward to hearing from you soon..Here is our Contact information:Corresponding author: Wenchang ZhangTel.: 86-0591-********;Fax:0591-********;E-mail address: wenchang2008@Thanks very much for your attention to our paper.Best Regards,Yours Sincerely.COVER LETTER 实用指南1、什么是cover letter?指的是投稿信2、cover letter的内容主要包括那些?应该简述所投稿件的核心内容、主要发现和意义,拟投期刊,对稿件处理有无特殊要求等(如“not to review” list)。
scientific reports收稿范围Scientific Reports is an open-access journal that publishes original research articles in all areas of natural sciences, including physics, chemistry, biology, environmental sciences, and earth sciences. With an aim to provide a platform for scientists to share their findings and contribute to the scientific community, Scientific Reports accepts submissions from researchers worldwide.The journal encourages the submission of research articles that present significant and novel findings in their respective fields. These articles should provide thorough experimental evidence and analysis, as well as clearly articulate the implications and impact of the research. The scope of the journal is broad, encompassing a wide range of sub-disciplines within the natural sciences.In physics, for example, Scientific Reports welcomes studies related to classical mechanics, quantum mechanics, astrophysics, cosmology, condensed matter physics, optics, and electromagnetism. Research articles in this field may focus on theoretical advancements, experimental investigations, or the application of physics principles to solve real-world problems.Chemistry research articles accepted by Scientific Reports cover diverse areas such as organic chemistry, inorganic chemistry, analytical chemistry, physical chemistry, and materials chemistry. Contributions may range from the synthesis and characterization of new compounds to the development of novel chemical processes. The journal values studies that have broad implications in various scientific and industrial fields.In the field of biology, Scientific Reports is interested in receiving submissions that explore topics like cellular and molecular biology, microbiology, genetics, ecology, and evolutionary biology. Highly impactful research articles may investigate the mechanism of disease or provide insights into the functioning of biological systems. Studies that advance our understanding of biodiversity, conservation, and sustainable practices are also welcome.Among the environmental sciences, Scientific Reports publishes research articles that investigate topics such as climate change, environmental pollution, ecology, atmospheric sciences, and oceanography. The journal seeks to promote studies that contribute to the development of sustainable practices and environmental management.Lastly, the field of earth sciences encompasses research articles on geology, geophysics, paleontology, seismology, and geochemistry. Scientific Reports is interested in studies that explore the Earth's processes, history, and resources, as well as those that contribute to the understanding of natural hazards and geological phenomena.In summary, Scientific Reports provides a platform for researchers to share their original findings in various sub-disciplines within the natural sciences. The journal welcomes high-quality, impactful research articles that present significant advancements, thorough experimental evidence, and clear implications.。
2022考研英语阅读捕获希格斯粒子Looking for the Higgs捕获希格斯粒子Enemy in sight?敌军现身?The search for the Higgs boson is closing in on its quarry希格斯玻色子的讨论接近其目标ON JULY 22nd two teams of researchers based at CERN, Europe s main particle-physicslaboratory, near Geneva, told a meeting of the European Physical Society in Grenoble thatthey had found the strongest hints yet that the Higgs boson does, in fact, exist.7月22日,驻欧洲粒子物理讨论所的两组讨论人员在格勒诺布尔欧洲物理协会的一次会议上声称,他们已经得到迄今为止最有力的线索,将力证希格斯玻色子的确真实存在。
The Higgs is thelast unobserved part of the Standard Model, a 40-year-old theory which successfullydescribes the behaviour of all the fundamental particles and forces of nature bar gravity.希格斯粒子是基础模型中最终一个尚未观测到的组件,基础模型已有40年的历史,它胜利地描述了全部基础粒子的行为及除重力以外的全部自然力。
Mathematically, the Higgs is needed to complete the modelbecause, otherwise, none of theother particles would have any mass.在数学层面上,希格斯粒子对于完成模型是必不行少的,这是由于,一旦缺少它,全部的其它粒子都将会失去质量。
生态毒理学报Asian Journal of Ecotoxicology第15卷第6期2020年12月V ol.15,No.6Dec.2020㊀㊀基金项目:国家自然科学基金资助项目(91543202,91843301);山西省自然科学基金资助项目(201801D121260)㊀㊀第一作者:白丽荣(1995 ),女,硕士研究生,研究方向为环境污染与健康,E -mail:*****************㊀㊀*通讯作者(Corresponding author ),E -mail:************.cnDOI:10.7524/AJE.1673-5897.20191121001白丽荣,谭子康,龚航远,等.大气PM 2.5全身暴露对大鼠多脏器病理㊁氧化应激指标和炎症因子的影响[J].生态毒理学报,2020,15(6):132-140Bai L R,Tan Z K,Gong H Y ,et al.Effects of atmospheric fine particulate matter whole body exposure on pathology of multiple organs,oxidative stress indicators and inflammatory factors in rats [J].Asian Journal of Ecotoxicology,2020,15(6):132-140(in Chinese)大气PM 2.5全身暴露对大鼠多脏器病理㊁氧化应激指标和炎症因子的影响白丽荣,谭子康,龚航远,耿红,董川,李瑞金*山西大学环境科学研究所,太原030006收稿日期:2019-11-21㊀㊀录用日期:2020-03-16摘要:探讨太原市冬季大气细颗粒物(PM 2.5)全身系统暴露对雄性SD 大鼠5种脏器病理学㊁氧化应激及炎症反应的影响㊂于2018年11月至2019年1月于太原市采用IVC 通气笼对大鼠进行大气PM 2.5的全身暴露染毒,分别染毒1个月和2个月;同时设立相应空气过滤对照组㊂采用HE 染色方法观察大鼠肝㊁肾㊁脾㊁胃和小肠的组织病理学变化,并测定肝㊁肾㊁脾和胃的脏器系数;采用生化方法测定大鼠组织匀浆液中的超氧化物歧化酶(SOD)活性㊁丙二醛(MDA)含量㊁白细胞介素-6(IL -6)和肿瘤坏死因子-α(TNF -α)的水平㊂PM 2.5暴露后大鼠5种组织均体现出不同程度的病理学损伤,MDA 含量与对照组相比显著增加㊂PM 2.5染毒使大鼠肝和肾的脏器系数比对照组显著增高;肾㊁脾㊁胃和小肠的SOD 活性比相应对照组明显上升;肝和小肠的IL -6和TNF -α水平显著高于对照组㊂2个月暴露的损伤效应比1个月暴露更为严重㊂大气PM 2.5全身系统暴露(结合气态污染物)可导致大鼠多器官损害,且随着暴露时间延长而损害加重;肝㊁肾和小肠对PM 2.5刺激表现敏感㊂其机制可能与PM 2.5引起的组织MDA 含量和炎症因子水平增加有关㊂关键词:PM 2.5;大鼠;病理学;SOD ;丙二醛;炎症因子;通气笼暴露系统文章编号:1673-5897(2020)6-132-09㊀㊀中图分类号:X171.5㊀㊀文献标识码:AEffects of Atmospheric Fine Particulate Matter Whole Body Exposure on Pathology of Multiple Organs ,Oxidative Stress Indicators and Inflamma-tory Factors in RatsBai Lirong,Tan Zikang,Gong Hangyuan,Geng Hong,Dong Chuan,Li Ruijin *Institute of Environmental Science,Shanxi University,Taiyuan 030006,ChinaReceived 21November 2019㊀㊀accepted 16March 2020Abstract :To investigate pathology,oxidative stress and inflammatory response of whole body exposure of ambient fine particulate matter (PM 2.5)in winter of Taiyuan in China in five organs of male SD rats in this study,rats were exposed to airborne PM 2.5using a whole -body inhalation system (individual ventilated cage (IVC)exposure system)for 1month and 2months respectively from November 2018to January 2019in Taiyuan.Meanwhile,the corre -sponding control groups (filter air groups)were respectively established.The histopathological changes of liver,第6期白丽荣等:大气PM2.5全身暴露对大鼠多脏器病理㊁氧化应激指标和炎症因子的影响133㊀kidney,spleen,stomach and small intestine in rats were observed by hematoxylin-eosin(HE)staining,and the or-gan coefficients of liver,kidney,spleen and stomach were measured.The levels of superoxide dismutase(SOD), malondialdehyde(MDA),interleukin-6(IL-6)and tumor necrosis factor-α(TNF-α)in rat tissue homogenates were determined by biochemical methods.PM2.5exposure caused pathological damage in varying degrees in five organs in rats,and significantly increased MDA contents compared with the corresponding an coefficients ofliver and kidney in rats exposed to PM2.5were raised relative to the controls.The SOD activities of the kidney, spleen,stomach and small intestine were significantly higher than that of the controls.The levels of IL-6and TNF-αin the liver and small intestine were obviously elevated compared to the controls.The injury effects of2-monthexposure were more severe than that of1-month exposure.Whole body exposure to atmospheric PM2.5(combined with gaseous pollutants)caused multiple organ damage in rats,and such injury aggravated with the prolonged ex-posure time.Liver,kidney and small intestine were sensitive to PM2.5stimulation.Increases of MDA contents andinflammatory factor levels might implicate in the impairment of multiple murine organs induced by PM2.5.Keywords:PM2.5;rats;pathology;SOD;malondialdehyde;inflammatory factor;Individual Ventilated Cage(IVC) exposure system㊀㊀细颗粒物(PM2.5)是我国主要的大气污染物及雾霾的主要成因㊂太原市大气全年PM2.5质量浓度在冬季最高,超标严重,太原市10月超标污染物主要是PM2.5和PM10,超标时间均为14d,超标率为64%,而11月除了颗粒物外,SO2也超标,PM2.5㊁PM10和SO2超标时间分别为10㊁9和4d,超标率分别为56%㊁50%和22%㊂采暖开始后的约2周内(11月1 16日),大气PM2.5和PM10平均质量浓度分别为(57.7+31.3)μg㊃m-3和(108.2ʃ48.8)μg㊃m-3,较采暖前有所下降,均未超标;11月17 23日其平均质量浓度分别为(133.0ʃ44.5)μg㊃m-3和(213.5+ 55.5)μg㊃m-3,大气PM2.5和PM10浓度大幅上升,超标倍数增加㊂太原市大气中PM2.5的主要成分是有机碳(OC)㊁元素碳(EC)㊁硫酸盐㊁硝酸盐和铵盐,金属元素成分包括Na㊁K㊁Zn㊁Pb㊁As㊁Cd和Cu[1]㊂化石燃料燃烧㊁工业生产排放㊁汽车尾气㊁生物质燃烧和道路扬尘是我国PM2.5的主要来源㊂国际癌症研究机构确定空气颗粒物为人类一类致癌物[2]㊂许多流行病学研究显示,PM2.5暴露导致呼吸系统疾病㊁心血管系统疾病㊁神经系统疾病及肺癌的发病率和死亡率上升[3-4],其毒理机制主要涉及氧化应激损伤机制和炎症介质机制[5]㊂PM2.5对肺㊁心脏和脑的健康效应和毒性研究是人们关注的主要课题㊂胃肠道和肝脏属于机体的消化系统,脾脏属于免疫系统,肾脏属于泌尿系统㊂研究发现,PM2.5暴露与胃癌㊁慢性肾病和肝癌的发生有关[6-8]㊂毒理学研究也表明,PM2.5可引起肝脏㊁肾脏㊁肠道和脾脏损伤效应[9-14]㊂这说明,人们不仅非常关注PM2.5对心肺和神经系统损伤的研究,同时也注意到PM2.5对全身其他器官的影响,PM2.5的全身毒性效应研究正被逐渐重视㊂在PM2.5毒理学相关研究中,实验动物暴露方式有滴注(气道滴注和鼻滴注)[10,12,15-16]和全身暴露[9,11,17]㊂动物全身PM2.5吸入系统主要有浓缩式和实时暴露2种,目前已成为科学的PM2.5暴露研究手段㊂采用PM2.5混悬液人工气道滴注的方法与动物吸入暴露的方式有差别,全身吸入暴露系统更能模拟动物在现实大气实际暴露的情境,已成为人们更认可的暴露模式㊂太原市工业以煤炭㊁焦化和冶金为主,大气污染以煤烟型为主,特别是其盆地结构和严重的冬季辐射逆温,加剧了大气污染㊂本研究于2018年11月中旬开始采用全身PM2.5真实暴露模型(通气笼暴露系统)进行SD大鼠暴露实验,分别暴露1个月和2个月,研究了太原市冬季PM2.5亚慢性暴露对大鼠肝㊁肾㊁胃㊁脾和小肠的病理损伤㊁脂质过氧化效应和炎症因子变化,揭示PM2.5对消化系统㊁免疫系统和泌尿系统的损伤效应,为研究PM2.5全身性毒性机制提供实验数据㊂这种PM2.5真实暴露模型模拟现实生活中PM2.5亚慢性暴露情境,笼内动物吸入PM2.5,符合实际情况,更有助于揭示PM2.5真实致病机制㊂1㊀材料与方法(Materials and methods)1.1㊀暴露系统本研究采用了一个独特新颖的PM2.5暴露系统(IVC通气笼),购自苏州市君圣实验动物设备有限公司㊂该系统的特点是实验动物可在通气笼内全身134㊀生态毒理学报第15卷暴露㊁真实吸入大气环境PM2.5㊂对照笼使用3层HEAP过滤膜,有效阻挡了颗粒物,使得笼内PM2.5浓度低于TSI颗粒物测定仪的检出限(1μg㊃m-3);暴露笼使用特殊的过滤膜,只允许PM2.5颗粒进入笼内㊂经测定,笼内PM2.5浓度约为室外大气浓度的70%,即61μg㊃m-3㊂实验过程中严密监测笼内气象条件,保持相对恒定的温度(22~27ħ)㊁湿度(30%~40%),温/湿度通过仪器显示屏控制,在显示屏上设定适宜的温湿度后仪器自动调整,并保证PM2.5暴露笼内的物理条件与对照笼一致㊂1.2㊀动物分组健康清洁级雄性SD大鼠20只,购自北京中国人民解放军医学科学院实验动物中心,体质量(157.5ʃ7.93)g,购入后在本研究室动物房饲养㊂20只大鼠随机分为对照组(1个月和2个月)和PM2.5暴露组(1个月和2个月),共4组,每组5只㊂在实验过程中,大鼠被安置在对照笼和暴露笼中,分别连续暴露1个月和2个月(每天24h,每周7d)㊂PM2.5组大鼠全身暴露吸入真实环境大气PM2.5;对照组大鼠吸入无颗粒物的空气㊂光照和暗周期均为12h㊂动物可随意获得食物和水㊂每3天更换一次笼内垫料㊂1.3㊀大鼠组织的脏体比与组织病理学分析暴露结束后,大鼠引颈脱臼处死,解剖㊁分离小肠㊁肝脏㊁脾脏㊁胃和肾脏组织㊂分别称量大鼠体质量及肝㊁脾㊁胃和肾的质量,计算其脏器系数(%)(脏器质量/动物体质量ˑ100%)㊂后分别取不同组大鼠5个组织(每个组织约5mm3)固定于福尔马林溶液中,石蜡包埋,制成5μm的切片,脱蜡,苏木精-伊红染色(HE染色),脱水透明,中性树胶固定,封片,显微镜下观察㊂1.4㊀促炎因子和脂质过氧化指标测定取组织样品,剪取0.1g组织㊁加入1mL生理盐水于冰浴上进行组织匀浆,组织匀浆液离心后(3000r㊃min-1,10min,4ħ),收集上清液㊂采用ELISA试剂盒(上海西唐生物科技有限公司)测定大鼠肝㊁肠㊁脾㊁胃和肾组织匀浆液中白细胞介素-6(IL-6)和肿瘤坏死因子-α(TNF-α)水平;采用特定生化试剂盒(南京建成生物工程研究所)测定超氧化物歧化酶(SOD)活性和丙二醛(MDA)含量㊂具体操作步骤按照试剂盒说明书严格进行㊂1.5㊀统计分析实验数据以平均值ʃ标准差(MeanʃSD)表示,应用SPSS19.0软件对数据进行单因素方差分析,用最小显著性差异法(LSD)检验,以P<0.05为差异有统计学意义㊂2㊀结果(Results)2.1㊀PM2.5暴露组和对照组大鼠脏器系数变化与对照相比较,1个月暴露组大鼠肝㊁脾和肾脏器系数增加有显著统计学意义(P<0.01)㊂2个月暴露组和1个月暴露组相比,各组织脏器系数变化的差异均具有统计学意义(P<0.01)㊂如表1所示,PM2.5暴露后动物体质量较相应对照组增加较大;而脏器质量变化较小㊂2.2㊀PM2.5暴露引起各组大鼠组织病理学变化肠HE染色结果表明,对照组大鼠结肠黏膜上皮细胞完好,主要为单层柱状上皮,肠腺有序排列规则整齐,未见损伤及炎性改变(图1(a)和图1(c))㊂与对照组相比,1个月暴露组大鼠结肠黏膜上皮细胞结构紊乱,肠腺排列不整齐(图1(b))㊂2个月暴露组大鼠结肠损伤严重,表现为腺体坏死㊁萎缩,黏液分泌减少,固有层与黏膜下层有部分炎性细胞如中性粒细胞㊁淋巴细胞和单核细胞等浸润(图1(d))㊂肝脏对照组的肝细胞以中央静脉为中心向周围放射状排列,无明显异常(图2(a)和图2(c))㊂与对照组相比,1个月暴露组大鼠肝组织结构紊乱,有少量表1㊀PM2.5暴露对大鼠主要脏器系数的影响Table1㊀Influence of PM2.5exposure on main organ coefficients in rats(%)组别Group肝Liver脾Spleen胃Stomach肾Kidney1个月对照1-month control 2.71ʃ0.120.18ʃ0.03 1.16ʃ0.360.79ʃ0.091个月暴露1-month exposure 3.62ʃ0.97b0.21ʃ0.04b 1.16ʃ0.320.88ʃ0.12b2个月对照2-month control 2.42ʃ0.180.16ʃ0.010.48ʃ0.040.62ʃ0.052个月暴露2-month exposure 2.31ʃ0.20**0.14ʃ0.01**0.54ʃ0.11**0.61ʃ0.03**注:数据表示为xʃSD(n=5);暴露组与对照组比较,b P<0.01;2个月暴露与1个月暴露比较,**P<0.01㊂Note:Data are expressed as xʃSD(n=5);compared with the control group,b P<0.01;comparison of2-month exposure with1-month exposure,**P<0.01.第6期白丽荣等:大气PM2.5全身暴露对大鼠多脏器病理㊁氧化应激指标和炎症因子的影响135㊀出血现象(图2(b))㊂2个月暴露组出现明显脂肪变性,肝细胞水肿,肝细胞浑浊,可见个别细胞坏死和炎细胞浸润(图2(d))㊂脾脏HE染色结果表明,对照组大鼠脾白髓边缘区的界线清晰,均匀致密蓝染且组织排列有序(图3(a)和图3(c));与对照组比较,PM2.5染毒组脾脏有明显的结构异常现象,红髓和白髓的界线不清楚,团块状结构消失且有炎性细胞浸润㊂2个月PM2.5暴图1㊀4组大鼠肠组织病理学改变情况(HE染色,200ˑ)Fig.1㊀Pathological changes of intestine in four groups of rats(HE staining,200ˑ)图2㊀4组大鼠肝组织病理学改变情况(HE染色,200ˑ)Fig.2㊀Pathological changes of liver tissue in four groups of rats(HE staining,200ˑ)136㊀生态毒理学报第15卷露大鼠脾脏损伤比1个月暴露更严重,伴有大量充血现象(图3(b)和图3(d))㊂肾脏对照组无明显异常,各曲管上皮排列有序(图4(a)和图4(c))㊂PM 2.5的1个月暴露组有不同程度肾间质水肿和肾介质炎性细胞浸润,肾小管再生,再生细胞呈现多形,有明显出血现象(图4(b))㊂PM 2.5的2个月暴露组有肾小球固缩,肾小管上皮细胞水肿,有明显出血现象(图4(d))㊂图3㊀4组大鼠脾组织病理学改变情况(HE 染色,200ˑ)Fig.3㊀Pathological changes of spleen in four groups of rats (HE staining,200ˑ)图4㊀4组大鼠肾组织病理学改变情况(HE 染色,200ˑ)Fig.4㊀Pathological changes of kidney in four groups of rats (HE staining,200ˑ)第6期白丽荣等:大气PM2.5全身暴露对大鼠多脏器病理㊁氧化应激指标和炎症因子的影响137㊀㊀㊀HE染色结果表明,对照组大鼠无明显异常(图5(a)和图5(c))㊂1个月暴露组大鼠胃体部浅层黏膜上皮细胞充血(图5(b))㊂2个月暴露组大鼠胃体部浅层黏膜上皮细胞变性坏死,深部细胞排列疏松,间质水肿,有少量上皮细胞核溶解消失,残存少数核碎片(图5(d))㊂2.3㊀PM2.5对大鼠各组织IL-6和TNF-α水平的影响由图6可知,PM2.5暴露组大鼠各组织中的IL-6 (图6(a))和TNF-α(图6(b))表达均上升,与对照组相比有显著性差异(P<0.05,P<0.01)㊂2个月PM2.5暴露组大鼠5个组织中IL-6和TNF-α水平与1个月PM2.5暴露组相比没有显著差异㊂2.4㊀PM2.5对大鼠各组织SOD活性和MDA含量影响由图7(a)可知,大鼠肠㊁脾㊁肾和胃组织中SOD活性均有不同程度的增加,其中,1个月PM2.5暴露组大鼠肠㊁脾和肾组织以及2个月暴露组大鼠胃组织中SOD活性较对照组有显著性差异(P<0.05),2个月PM2.5暴露组大鼠肠㊁脾和肾组织SOD活性与对照组相比具有显著差异(P<0.01)㊂而2个月PM2.5组大鼠肝组织和1个月PM2.5组胃组织中SOD活性较对照组虽有增加但并无显著性差异㊂不同组大鼠5种组织中MDA含量如图7(b)所示㊂与对照组相比,PM2.5暴露后大鼠各组织中MDA含量都显著上升,有显著差异(P<0.05或P<0.01)㊂1个月PM2.5暴露组大鼠肠㊁肝㊁肾和胃组织和2个月PM2.5暴露组大鼠肠㊁脾组织分别较对照组中MDA含量有显著性差异(P<0.05)㊂1个月PM2.5暴露组大鼠脾组织和2个月PM2.5暴露组大鼠肝㊁肾㊁胃组织中MDA含量与对照组相比有显著差异(P<0.01)㊂2个月PM2.5暴露组大鼠肝㊁肾和胃组织中MDA含量与1个月PM2.5暴露组相比,具有显著性差异(P<0.05)㊂3㊀讨论(Discussion)当前我国大气PM2.5污染还未全面解决,许多地区PM2.5污染加重,灰霾现象频繁发生,严重影响空气环境质量和人们的身体健康㊂PM2.5粒径小,比表面积大,含有无机㊁有机和微生物等多种有害物质,进入人体后对健康造成极大危害,PM2.5全身性损伤效应和毒性机制研究是我国大气污染研究关注的焦点㊂早期PM2.5毒理学研究的动物暴露模式多采用气道滴注和口鼻式滴注方式进行颗粒物染毒[11,13,15-16],这种方式的局限性在于动物不是通过自然呼吸空气染毒㊁且滴注过程中对动物有机械损伤,图5㊀4组大鼠胃组织病理学改变情况(HE染色,200ˑ)Fig.5㊀Pathological changes of stomach in four groups of rats(HE staining,200ˑ)138㊀生态毒理学报第15卷图6㊀PM 2.5暴露对大鼠各组织白细胞介素-6(IL-6)和肿瘤坏死因子-α(TNF -α)的影响注:与对照组相比,a P <0.05,b P <0.01㊂Fig.6㊀Effects of PM 2.5exposure on interleukin -6(IL -6)and tumor necrosis factor -α(TNF -α)in rat tissuesNote:Compared with the control,a P <0.05,bP <0.01.图7㊀PM 2.5暴露对大鼠各组织超氧化物歧化酶(SOD )活性和丙二醛(MDA )含量的影响注:与对照组相比,a P <0.05,b P <0.01;与1个月暴露组相比,*P <0.05㊂Fig.7㊀Effects of PM 2.5exposure on superoxide dismutase (SOD)activity and malondialdehyde (MDA)content in different tissues of ratsNote:Compared with the control,aP <0.05,bP <0.01;compared to 1-month of exposure,*P <0.05.不适宜开展较为长期的滴注暴露㊂近年来,全身暴露染毒系统逐渐兴起,通过模拟颗粒物或PM 2.5暴露探讨其对机体的毒性效应及机制,这种方式更为科学,也符合实际颗粒物暴露情境㊂由于城市大气PM 2.5污染受到污染源㊁地形和气象因素影响,PM 2.5污染特征及化学组分存在较大差异,PM 2.5暴露引起的健康危害及毒理学机制较为复杂㊂太原市具有特殊的工业污染㊁盆地地形和冬季辐射逆温特点,污染物不易扩散,大气污染较为严重㊂鉴于此,本研究在太原市这一典型地区,利用全身暴露染毒系统(IVC 通气笼),控制温度㊁湿度㊁气压㊁空气流速和噪音等条件,对大鼠进行PM 2.5实时暴露(暴露期间笼内PM 2.5平均浓度约为61μg ㊃m -3)和亚慢性吸入染毒,研究PM 2.5对不同组织的损伤效应,为PM 2.5与全身性毒性研究提供动物实验基础㊂HE 染色结果是评价动物病理效应的常用方法㊂本研究结果表明,PM 2.5暴露1个月和2个月对大鼠肝㊁肾㊁脾㊁胃和小肠均造成了不同程度的病理损伤,且PM 2.5暴露2个月大鼠5种脏器的病理损伤效应和比暴露1个月严重,提示PM 2.5对多种器官存在病理损伤效应,PM 2.5暴露时间越长,组织损伤越严重,暴露2个月的时间该病理损伤效应未得第6期白丽荣等:大气PM2.5全身暴露对大鼠多脏器病理㊁氧化应激指标和炎症因子的影响139㊀到修复或缓解㊂脏器系数是毒理学研究的常用指标㊂由表1所示,1个月暴露组SD大鼠肝和肾的脏器系数比对照组显著增加,且超出了文献推荐的8~12周SD大鼠正常参考值范围[18],提示肝和肾可能是PM2.5毒作用的靶器官㊂2个月暴露组各组织的脏器系数比1个月暴露组显著下降,这可能与年龄增长㊁动物体质量增加有关㊂MDA是膜脂质过氧化效应和氧化应激的重要标志物[19]㊂PM2.5暴露1个月和2个月后大鼠各组织中MDA含量与对照组相比显著升高,表明PM2.5引起细胞膜脂质过氧化效应,这可能是PM2.5造成组织损伤的原因之一㊂PM2.5暴露2个月后大鼠肝㊁肾和胃组织中MDA含量的增加与1个月暴露组相比有显著统计学意义,提示肝㊁肾和胃组织发生了显著的脂质过氧化效应㊂活性氧(ROS)介导的氧化应激与疾病的发生有关,SOD具有清除超氧阴离子自由基的作用[20]㊂本研究发现,PM2.5暴露引起大鼠5个组织的SOD活性水平增加,在肾㊁脾㊁胃和小肠中的变化与对照组相比有显著性差异㊂我们推测PM2.5暴露导致机体内产生了较多ROS,在太原冬季PM2.5污染的条件下,SOD活性被激发以清除活性氧自由基㊂组织中MDA含量增加,反映出SOD 活性增加没有有效减少ROS堆积,氧化应激效应依然存在㊂文献报道显示,在不同的PM2.5暴露浓度㊁暴露时间和暴露模式下,PM2.5引起了实验动物肝㊁脾㊁肾和肠道的病理损伤和氧化应激效应[10-13],这与本研究结果一致㊂除氧化应激外,炎症反应是决定多器官损伤严重程度和疾病发生的关键因素[6,12,21]㊂由炎症因子测定结果可知,PM2.5引起大鼠5个组织IL-6和TNF-α水平增加,其中,肝和小肠中的水平与对照组相比有显著性差异,这说明PM2.5可使多种组织出现炎症反应,而这种反应有组织差异性,肝和小肠对PM2.5的刺激较为敏感㊂文献报道,PM2.5通过ROS介导线粒体自噬引起肝细胞纤维化,通过内质网应激引起肝细胞凋亡和炎症,通过肝细胞氧化应激和炎症引起的脂质代谢异常[13-15],这些机制可能和肝癌和非酒精性肝病的发病有关㊂PM2.5不仅引起小鼠小肠病理损伤,且PM2.5慢性暴露引起小鼠肠道功能紊乱,并影响了糖代谢异常[22]㊂这说明PM2.5致肝和小肠的损伤机制与代谢和疾病的发生有关,其毒性效应研究应引起重视㊂综上,结合脏器系数㊁HE染色㊁炎症因子和MDA结果,我们发现PM2.5亚慢性暴露造成大鼠多种组织出现不同程度病理损伤㊁脏器系数异常改变㊁脂质过氧化效应和炎症因子水平升高,且PM2.5暴露2个月的损伤效应大于1个月㊂相对而言,肝脏㊁小肠和肾脏受到PM2.5的负面影响较大㊂有研究报道,上海市大气PM2.5暴露对心肌和肺功能有损伤[23]㊂广州㊁深圳㊁东莞和肇庆4个城市的大气PM2.5均可使大鼠肺组织发生氧化应激和炎性反应㊂其中,广州㊁东莞和深圳的PM2.5对指标的影响较肇庆强[24]㊂已有对太原市的研究表明,PM2.5暴露可以通过诱导炎症反应以及炎症因子的表达从而导致心血管疾病发生[25]㊂不同采样点的PM2.5由于成分组成不同而毒性毒理有很大差异,本文结论仅针对所采集具体样本㊂由于目前在全身暴露染毒系统实验条件下,无法分离PM2.5和气态污染物,因此推测雄性SD大鼠暴露器官异常可能是PM2.5和气态污染物共同影响的结果㊂PM2.5对消化系统㊁免疫系统和泌尿系统的健康效应和毒性作用机制还需进一步的实验研究㊂通讯作者简介:李瑞金(1967 ),女,博士,教授,主要研究方向为环境污染与健康㊂参考文献(References):[1]㊀张媛,耿红,张东鹏,等.太原市采暖初期大气PM2.5质量浓度变化分析[J].环境与健康杂志,2016,33(5):385-390,470Zhang Y,Geng H,Zhang D P,et al.Mass concentrationvariation of urban fine particulate matter in initial stage ofheating period in Taiyuan[J].Journal of Environment andHealth,2016,33(5):385-390,470(in Chinese)[2]㊀Li R Y,Zhou R,Zhang J G.Function of PM2.5in the pathogenesis of lung cancer and chronic airway inflamma-tory diseases[J].Oncology Letters,2018,15(5):7506-7514[3]㊀International Agency for Research on Cancer(IARC).IARC Monographs on the Evaluation of CarcinogenicRisks to Humans[M].Lyon:IARC,1990:341[4]㊀Badyda A J,Grellier J,Dᶏbrowiecki P.Ambient PM2.5exposure and mortality due to lung cancer and cardiopul-monary diseases in Polish cities[J].Advances in Experi-mental Medicine and Biology,2017,944:9-17[5]㊀Shou Y K,Huang Y L,Zhu X Z,et al.A review of thepossible associations between ambient PM2.5exposures140㊀生态毒理学报第15卷and the development of Alzheimer s disease[J].Ecotoxi-cology and Environmental Safety,2019,174:344-352 [6]㊀杜鹏瑞,杜睿,任伟珊.城市大气颗粒物毒性效应及机制的研究进展[J].中国环境科学,2016,36(9):2815-2827Du P R,Du R,Ren W S.Research progress on toxicolog-ical characteristics and mechanisms of urban atmosphericparticulate matters[J].China Environmental Science,2016,36(9):2815-2827(in Chinese)[7]㊀Weinmayr G,Pedersen M,Stafoggia M,et al.Particulatematter air pollution components and incidence of cancersof the stomach and the upper aerodigestive tract in theEuropean Study of Cohorts of Air Pollution Effects(ES-CAPE)[J].Environment International,2018,120:163-171[8]㊀Bragg-Gresham J,Morgenstern H,McClellan W,et al.County-level air quality and the prevalence of diagnosedchronic kidney disease in the US Medicare population[J].PLoS One,2018,13(7):e0200612[9]㊀V oPham T,Bertrand K A,Tamimi R M,et al.AmbientPM2.5air pollution exposure and hepatocellular carcinoma incidence in the United States[J].Cancer Causes&Con-trol,2018,29(6):563-572[10]㊀Ge C X,Xu M X,Qin Y T,et al.iRhom2loss alleviatesrenal injury in long-term PM2.5-exposed mice by suppres-sion of inflammation and oxidative stress[J].Redox Biol-ogy,2018,19:147-157[11]㊀Su R J,Jin X T,Lyu L,et al.The potential immunotoxici-ty of fine particulate matter based on SD rat spleen[J].Environmental Science and Pollution Research Interna-tional,2019,26(23):23958-23966[12]㊀Li D C,Zhang R,Cui L H,et al.Multiple organ injury inmale C57BL/6J mice exposed to ambient particulate mat-ter in a real-ambient PM exposure system in Shijiazhuang,China[J].Environmental Pollution,2019,248:874-887 [13]㊀Li R J,Zhang M,Wang Y,et al.Effects of sub-chronicexposure to atmospheric PM2.5on fibrosis,inflammation, endoplasmic reticulum stress and apoptosis in the livers ofrats[J].Toxicology Research,2018,7(2):271-282[14]㊀Qiu Y N,Wang G H,Zhou F,et al.PM2.5induces liver fi-brosis via triggering ROS-mediated mitophagy[J].Eco-toxicology and Environmental Safety,2019,167:178-187[15]㊀Xu M X,Ge C X,Qin Y T,et al.Prolonged PM2.5expo-sure elevates risk of oxidative stress-driven nonalcoholicfatty liver disease by triggering increase of dyslipidemia[J].Free Radical Biology&Medicine,2019,130:542-556 [16]㊀Wang X F,Jiang S F,Liu Y,et prehensive pulmo-nary metabolome responses to intratracheal instillation ofairborne fine particulate matter in rats[J].The Science ofthe Total Environment,2017,592:41-50[17]㊀Xu M X,Zhu Y F,Chang H F,et al.Nanoceria restrainsPM2.5-induced metabolic disorder and hypothalamus in-flammation by inhibition of astrocytes activation relatedNF-κB pathway in Nrf2deficient mice[J].Free RadicalBiology&Medicine,2016,99:259-272[18]㊀董延生,尹纪业,陈长,等.SD大鼠脏器重量及脏器系数正常参考值的确立与应用[J].军事医学,2012,36(5):351-353Dong Y S,Yin J Y,Chen C,et al.Establishment and ap-plication of the normal reference values of organ massesand organ/body coefficients in SD rats[J].Military Medi-cal Sciences,2012,36(5):351-353(in Chinese)[19]㊀Czerska M,Mikołajewska K,Zieliński M,et al.Today soxidative stress markers[J].Medycyna Pracy,2015,66(3):393-405[20]㊀Poprac P,Jomova K,Simunkova M,et al.Targeting freeradicals in oxidative stress-related human diseases[J].Trends in Pharmacological Sciences,2017,38(7):592-607[21]㊀苏瑞军,晋小婷,安全,等.太原市冬季大气PM2.5暴露对小鼠脏器及炎症因子的影响研究[J].环境与健康杂志,2015,32(8):677-679Su R J,Jin X T,An Q,et al.Effects of PM2.5collected in winter of Taiyuan on organs weights and inflammatorycytokines expressions in mice[J].Journal of Environmentand Health,2015,32(8):677-679(in Chinese)[22]㊀Wang W J,Zhou J,Chen M J,et al.Exposure to concen-trated ambient PM2.5alters the composition of gut micro-biota in a murine model[J].Particle and Fibre Toxicolo-gy,2018,15(1):17[23]㊀Xie Y Q,Zhang X,Tian Z Y,et al.Preexposure to PM2.5exacerbates acute viral myocarditis associated with Th17cell[J].International Journal of Cardiology,2013,168(4):3837-3845[24]㊀闫庆倩,赵学彬,杨莉,等.不同地区大气PM2.5致大鼠肺损伤的比较实验研究[J].环境与健康杂志,2012,29(1):7-11,97Yan Q Q,Zhao X B,Yang L,et parative experi-ment of lung injury in rats induced by PM2.5collected from different cities of China[J].Journal of Environmentand Health,2012,29(1):7-11,97(in Chinese)[25]㊀Li R J,Kou X J,Geng H,et al.Mitochondrial damage:An important mechanism of ambient PM2.5exposure-in-duced acute heart injury in rats[J].Journal of HazardousMaterials,2015,287:392-401Ң。
生态毒理学报Asian Journal of Ecotoxicology第16卷第5期2021年10月V ol.16,No.5Oct.2021㊀㊀基金项目:国家自然科学基金资助项目(21677028)㊀㊀第一作者:陈杰(1995 ),男,硕士,研究方向为环境化学,E -mail:***************** ㊀㊀*通讯作者(Corresponding author ),E -mail:**************.cnDOI:10.7524/AJE.1673-5897.20210326002陈杰,马芳芳,郭熙瑞,等.㊃Cl 引发萘的大气氧化机制及动力学[J].生态毒理学报,2021,16(5):14-23Chen J,Ma F F,Guo X R,et al.Atmospheric oxidation mechanism and kinetics of naphthalene initiated by chlorine radicals (㊃Cl)[J].Asian Journal of Ecotoxicology,2021,16(5):14-23(in Chinese)㊃Cl 引发萘的大气氧化机制及动力学陈杰,马芳芳,郭熙瑞,谢宏彬*大连理工大学环境学院,工业生态与环境工程教育部重点实验室,大连116024收稿日期:2021-03-26㊀㊀录用日期:2021-05-23摘要:氯自由基(㊃Cl)内陆来源的新发现增强了其对转化大气有机污染物的贡献,因此,需要更深入地研究㊃Cl 引发有机污染物的转化机制和动力学㊂萘(Nap)是一种重要的化学品,也是城市大气浓度最高的多环芳烃,前人针对羟基自由基(㊃OH)引发Nap 的大气氧化开展了研究㊂然而,目前对于㊃Cl 引发Nap 的大气氧化机制还不清楚㊂本研究通过量子化学计算(ωB97XD/6-311++G(3df,2pd)//ωB97XD/6-31+G(d,p))和动力学模拟相结合的方法研究了㊃Cl 引发Nap 的大气氧化机制与动力学,发现㊃Cl 主要加成到Nap 分子的C 5位置,形成加成中间体㊃C 10H 8Cl(R 1)㊂随后,O 2加成到R 1的C 2和C 6位置生成过氧自由基(RO 2㊃)R 1-2OO -s/a 和R 1-6OO -s/a(s/a =syn/anti ,syn 表示O 2加成方向和㊃Cl 加成方向相同,anti 表示O 2加成方向和㊃Cl 加成方向相反)㊂这4种RO 2㊃的环化㊁氢迁移和氯迁移反应均很难(能垒>20kcal ㊃mol -1)发生㊂因此,在低NO 浓度条件下,RO 2㊃主要和HO 2㊃反应生成氢过氧化合物(QOOH)和烷氧自由基(RO ㊃)R 1-2O -s/a 和R 1-6O -s/a ;在高NO 浓度条件下,RO 2㊃将主要与NO 反应生成RO ㊃(R 1-2O -s/a 和R 1-6O -s/a)和有机硝酸酯(C 10H 8ClNO 3)㊂生成的RO ㊃进一步通过单分子环化反应生成双环产物R 1-21O -s/a 和R 1-61O -s/a ㊂重要的是,生成的有机氢过氧化合物和有机硝酸酯的水生毒性比其母体化合物Nap 更强,表明㊃Cl 引发Nap 反应增加了Nap 释放的环境风险㊂揭示的机制对理解大气Nap 化学及Nap 释放导致的环境风险具有重要意义㊂关键词:萘;氯自由基;量子化学;转化机制;动力学文章编号:1673-5897(2021)5-014-10㊀㊀中图分类号:X171.5㊀㊀文献标识码:AAtmospheric Oxidation Mechanism and Kinetics of Naphthalene Initiated by Chlorine Radicals (㊃Cl )Chen Jie,Ma Fangfang,Guo Xirui,Xie Hongbin *Key Laboratory of Industrial Ecology and Environmental Engineering (Ministry of Education),School of Environmental Science and Technology,Dalian University of Technology,Dalian 116024,ChinaReceived 26March 2021㊀㊀accepted 23May 2021Abstract :The new findings of chlorine radicals (㊃Cl)in mid -continental areas increase the importance of ㊃Cl in transforming atmospheric organic pollutants.Hence,more research should be performed on the atmospheric trans -formation of organic pollutants initiated by ㊃Cl.Naphthalene (Nap)is one type of chemicals and the most abundant polycyclic aromatic hydrocarbons in the urban atmosphere.Atmospheric oxidation of Nap initiated by hydroxyl radicals (㊃OH)have been studied previously.However,the mechanism of ㊃Cl initiated reaction of Nap is not fully understood.Herein,㊃Cl initiated reactions of Nap were investigated by a quantum chemical method (ωB97XD/. All Rights Reserved.第5期陈杰等:㊃Cl 引发萘的大气氧化机制及动力学15㊀311++G(3df,2pd)//ωB97XD/6-31+G(d,p))and kinetics modeling.Results show that ㊃Cl addition to the C 5-position of Nap,forming radicals ㊃C 10H 8Cl(R 1),is the dominant reaction pathway.Subsequently,O 2is mainly added to the C 2and C 6positions of R 1to form peroxy radicals (RO 2㊃)R 1-2OO -s/a and R 1-6OO -s/a depending on the attacking direction of O 2(s/a =syn/anti,syn and anti correspond to the O 2additions from the same and opposite sites of the direction of ㊃Cl addition).In the atmosphere,the isomerization reactions (including cyclization,H -transfer,and Cl -transfer)of the four RO 2㊃proceed very slowly due to the high reaction energy barriers (>20kcal ㊃mol -1).There -fore,the formed four RO 2㊃will mainly react with HO 2㊃to form hydroperoxide (QOOH)and alkoxy radicals (RO ㊃)under the condition of low NO concentration,or react with NO to form organonitrates (C 10H 8ClNO 3)andRO ㊃(R 1-2O -s/a,R 1-6O -s/a)under the condition of high NO concentration.The formed RO ㊃will finally undergo ring -closure reactions forming bi -cyclic intermediates (R 1-21O -s/a,R 1-61O -s/a).More importantly,the predicted toxicity of the formed organonitrates and hydroperoxide is much higher than that of Nap,indicating the ㊃Cl initia -ted transformation of Nap increases the environmental risk caused by Nap emission.The revealed mechanism is ofgreat significance for understanding atmospheric chemistry and environmental risk assessment of Nap.Keywords :naphthalene;chlorine radicals;quantum chemistry;transformation mechanism;kinetics ㊀㊀氯自由基(㊃Cl)具有强的化学反应活性,与大气中大多数污染物的反应速率常数(k Cl )比羟基自由基(㊃OH)对应的反应速率常数(k OH )高10倍~100倍[1-4]㊂过去㊃Cl 被认为主要分布于沿海地区,由海盐经多相化学反应生成,浓度可达㊃OH 的1%~10%[1,5-7]㊂因此,一直认为㊃Cl 对沿海地区污染物的转化起着重要的作用[8-9]㊂然而,最近几年在中国㊁美国和加拿大等国的内陆城郊区域发现较高浓度的㊃Cl 前体物质ClNO 2[10-13]㊂这表明㊃Cl 对有机污染物的转化可能起着比以前更为重要的作用㊂近期关于北京地区大气污染的一项研究指出,㊃Cl 对烷烃类污染物的转化贡献率高于㊃OH [14]㊂此外,值得指出的是,㊃Cl 引发有机污染物的转化可能导致与㊃OH 引发反应不同的大气转化机制,进而可能会导致不同的大气归趋[3,15-17]㊂因此,㊃Cl 引发城市大气有机污染物的转化亟待研究㊂多环芳烃(polycyclic aromatic hydrocarbons,PAHs)是影响城市大气环境质量的重要污染物,主要由化石燃料和生物质燃烧排放,组成成分复杂[18-20]㊂萘(naphthalene,Nap)是分子量最小㊁挥发性最强的PAHs ㊂此外,Nap 还是重要的化学品,被用作化工原料[21]㊂较高的挥发性使Nap 能在生产㊁运输和使用过程中进入大气环境㊂全球尺度的大气模型模拟结果指出,Nap 的含量超过PAHs 总量的2/5[18]㊂特别是在洛杉矶,Reisen 和Arey [19]检测到Nap 的大气浓度高达1600ng ㊃m -3㊂重要的是,前人研究发现,Nap 还是二次有机气溶胶(secondary organic aerosol,SOA)的重要前体物[22-24]㊂如Huang 等[23]研究发现,在北京地区的雾霾事件中,Nap 及Nap 的衍生物转化生成的SOA 最高浓度可达到4.6~6.8μg ㊃m -3,约贡献了SOA 总量的14.9%㊂考虑到Nap 是一种重要化学品及其大气化学的重要性,有必要研究Nap 的大气转化㊂日间被大气中的㊃OH 和㊃Cl 氧化是Nap 的重要大气归趋[22,24-27]㊂前人采用实验和理论计算的方法研究了㊃OH 引发Nap 的大气氧化机制与动力学[26,28-29]㊂研究指出,㊃OH 主要加成到Nap 分子上形成加成中间体Nap -OH ,随后与O 2反应生成萘酚和过氧自由基(RO 2㊃)㊂RO 2㊃既可以经单分子异构化生成肉桂醛,也可与NO 反应生成烷氧自由基(RO ㊃)和有机硝酸酯,RO ㊃经环化后进入下一个O 2/NO 反应,最终生成含羰基的双环中间体[29]㊂然而,针对㊃Cl 引发Nap 的大气氧化,仅有一些实验测量了两者的反应速率常数并鉴别了部分氧化产物的碎片信息[25,30]㊂而对㊃Cl 引发Nap 的具体大气转化机制还不清楚㊂近年来,随着量子化学计算方法计算精度的提高和计算机计算能力的大幅提升,量子化学方法已经成为研究化学品及大气污染物转化机制和动力学的重要工具,在大气化学及化学品风险评价中发挥着重要作用[31-32]㊂如大连理工大学陈景文团队近年来利用量子化学方法成功揭示了多种化学品及大气污染物的转化机制和动力学[32-35]㊂本研究使用量子化学计算结合动力学模拟的方法研究了㊃Cl 引发Nap 的大气氧化机制与动力学,包括㊃Cl 与Nap 的引发反应及引发反应生成的重要. All Rights Reserved.16㊀生态毒理学报第16卷中间体的后续反应(包括加成中间体与O 2/NO 的双分子及RO 2㊃/RO ㊃的解离和异构化)的转化机制与动力学㊂此外,使用化学品水生毒性评估软件ECO -SAR 评估了部分闭壳层产物的水生毒性[36]㊂1㊀计算方法(Computational methods )使用Gaussian 09软件包进行电子结构优化及能量计算[37]㊂在ωB97XD/6-31+G(d,p)计算水平下,对所有反应渠道涉及的反应物(R)㊁反应前络合物(RC)㊁过渡态(TS)㊁反应后络合物(PC)和产物(P)的结构进行优化[38]㊂在相同计算水平下,通过内禀反应坐标计算确定TS 连接的R 和P ㊂使用ωB97XD/6-311++G(3df,2pd)方法计算单点能量㊂通过公式(1)计算每个物种的能量㊂E =E SP +E corr (1)式中:E SP 为单点能;E corr 为ωB97XD/6-31+G(d,p)计算水平获得的零点矫正能㊂本论文所有单分子反应和双分子反应速率常数均采用传统过渡态理论(canonical transition state the -ory,TST)进行计算,TST 计算在MULTIWELL -2014.1软件包中的Thermo 模块进行[39-42],使用Eckart 约近考虑单双分子反应中H 迁移和H 夺取反应的隧道效应[41]㊂2㊀结果与讨论(Results and discussion )2.1㊀㊃Cl 引发反应原则上,㊃Cl 可以夺取Nap 的H 原子或㊃Cl 加成到不饱和键上㊂Nap 分子存在8个H 原子和10个不饱和C 原子(图1),因此,㊃Cl 和Nap 反应有H 夺取和Cl 加成2类,共18条可能的反应渠道㊂由于Nap 具有D 2h 对称性,㊃Cl 与Nap 反应需要考虑的可能渠道可以降低到5条,即:㊃Cl 夺取Nap 的H 12和H 13;㊃Cl 加成到Nap 的C 4㊁C 5和C 6位置(图1)㊂基于计算的热力学数据(表1)可知,㊃Cl 加成到C 5和C 6位置放热显著高于H 夺取反应㊂因此,从热力学上判断,加成反应更容易发生㊂此外,值得指出的是,前人研究也发现㊃Cl 夺取sp 2杂化的C 上的H 很难发生[43]㊂因此,本研究仅考虑㊃Cl 和Nap 的加成反应㊂为确定㊃Cl 和Nap 反应的加成位点,首先,我们进行了二维势能面扫描[44]㊂以Nap 分子C 骨架平行且距离为0.25nm 的平面为扫描平面,以C 3-C 1和C 3-C 4方向分别为X 轴和Y 轴㊂㊃Cl 的移动范围为X =0~0.26nm 和Y =0.07~0.24nm 所包含的矩形,该矩形能覆盖Nap 的C 4㊁C 5和C 6上方区域㊂㊃Cl 在X 和Y 轴上的扫描步长均为0.01nm,共得到约500个扫面点的能量数据,根据此能量数据绘制的二维扫描势能面如图2所示㊂由图2可知,X =0.1~0.26nm 和Y =0.16~0.24nm 构成的区域(C 5㊁C 6附近)的能量显著低于X =0~0.15nm 和Y =0.07~0.10nm 构成的区域(C 4附近),即㊃Cl 在Nap 的C 5和C 6上部区域移动时,体系能量更低㊂结合前文未能优化到C 4加成产物的稳定构型,说明C 4不可能是主要的加成位点,加成反应可能的主要加成位点为C 5和C 6㊂为了进一步确定加成到C 5和C 6哪个位置更可行,对㊃Cl 加成到C 5和C 6位置进行一维势能面扫描㊂C -Cl 垂直距离(r C -Cl)为0.19~0.28nm ,每隔0.005nm 进行扫描㊂㊃Cl 加成到C 5位置的扫描曲线如图3(a)所示,由图可知,当r C 5-Cl 在0.19~0.28nm 范围变化时,随着C 5-Cl 距离的逐渐增加,能量逐渐升高,说明㊃Cl 加成到C 5位置是无能垒的过程㊂㊃Cl图1㊀萘(Nap )的分子结构及原子序号Fig.1㊀Structure and atomic labeling of naphthalene (Nap)表1㊀在ωB97XD/6-311++G(3df,2pd)//ωB97XD/6-31+G(d,p)计算水平下,㊃Cl 和Nap 反应的反应能Table 1㊀Calculated energy changes of ㊃Cl +Nap reaction at the B97XD/6-311++G(3df,2pd)//ωB97XD/6-31+G(d,p)level(kcal ㊃mol -1)产物Products 能变Energy changes(ΔE 0K )自由能变Free energy changes(ΔG 298K )C 4-加成C 4-Addition --C 5-加成C 5-Addition -18.1-10.5C 6-加成C 6-Addition-13.5-6.1H 12-摘氢H 12-Abstraction -0.5-1.8H 13-摘氢H 13-Abstraction -0.1-2.2注:-表示未能优化得到稳定产物,能量信息缺失㊂Note:-represents the missing data for C 4-addition products,and a stable configuration has not been located.. All Rights Reserved.第5期陈杰等:㊃Cl 引发萘的大气氧化机制及动力学17㊀加成到C 6位置的扫描曲线如图3(b)所示,由图可知,在r C 6-Cl =0.275nm 和r C 6-Cl =0.19nm 时,㊃Cl和Nap 已经形成2种不同的局域能量极小点,分别对应可能的C 5和C 6加成㊂因此,即使对㊃Cl 加成到C 6位置进行扫描,也会首先形成加成到C 5位置的中间体,表明C 5位置是最可行的加成位点,与热力学结果一致㊂值得指出的是,在㊃Cl 与Nap 的氧化实验中,㊃Cl 加成到C 5的加成中间体也被推测为主要加成产物,支持计算结果的可靠性[25]㊂接下来仅考虑C 5位置加成中间体㊃C 10H 8Cl(R 1)的后续反应㊂2.2㊀R 1与O 2反应在对流层大气环境中,和其他C 中心自由基反应类似,加成中间体R 1可能进一步与O 2发生双分子反应[29,45]㊂R 1与O 2可以发生直接H 夺取和O 2加成2类反应:(1)O 2直接夺取R 1中 CClH 基团上的H 生成C 10H 7Cl 和HO 2㊃;(2)O 2通过顺式(syns)或反式(anti)2个方向加成到自旋密度较高的C 6(0.59)㊁C 2(0.45)㊁C 10(0.14)㊁C 14(0.16)和C 4(0.13)位置形成过氧自由基中间体R 1-i OO -s/a(i =2,4,6,10,14,s/a =syns/anti ,O 2与㊃Cl 方向相同为syns 加成,相反为anti 加成)㊂表2中列出了R 1与O 2反应的能量信息㊂基于计算的反应能垒(ΔE 0K ʂ)可知,相比于其他反应途径(直接H 夺取和加成到C 10㊁C 14和C 4位点),O 2通过syns 和anti 加成到C 2(ΔE 0K ʂ=6.7kcal ㊃mol -1和5.4kcal ㊃mol -1)和C 6(ΔE 0K ʂ=8.6kcal ㊃mol -1和5.2kcal ㊃mol -1)位点形成R 1-2OO -s/a 和R 1-6OO -s/a是最可行的㊂此外,通过对比加成渠道形成的加成中间体的稳定性(基于ΔE 0K 和ΔG 298K 判断)发现,相比于其他加成位点的加成产物,C 2和C 6加成产物更稳定,这主要是因为这些加成产物中保留着完整的苯环㊂O 2与OH -Nap 的反应中也发现了类似的现象[22]㊂因此,在下文中我们仅考虑O 2加成到C 2和C 6位置形成的加成中间体R 1-2OO -s/a 和R 1-6OO -s/a 的反应㊂图2㊀在ωB97XD/6-31+G(d,p)计算水平下,㊃Cl 加成到Nap 的二维扫描势能面(PES)注:ΔE 为扫描点相对于反应物㊃Cl+Nap 的能量,单位为kcal ㊃mol -1㊂Fig.2㊀Contour plots of the two -dimensional scannedpotential energy surface (PES)for ㊃Cl addition to the Nap atthe ωB97XD/6-31+G(d,p)levelNote:ΔE is the relative energy (kcal ㊃mol -1)withreference to the reactants ㊃Cl+Nap.图3㊀在ωB97XD/6-31+G(d,p)计算水平下,㊃Cl 加成到Nap 的C 5(a)和C 6(b)位点的扫描势能曲线注:r C -Cl 指C -Cl 间距离,单位为10-1nm ;r C -Cl =0.28nm 构型的能量作为ΔE 计算的参考点,单位为kcal ㊃mol -1㊂Fig.3㊀Scanned potential energy curves for ㊃Cl addition to C 5(a)and C 6(b)sites of Nap at the ωB97XD/6-31+G(d,p)levelNote:r C -Cl is the value of Cl -C distance (10-1nm);ΔE is the relative energy (kcal ㊃mol -1)with reference to the configuration (r C -Cl =0.28nm).. All Rights Reserved.18㊀生态毒理学报第16卷2.3㊀R1-2OO-s/a和R1-6OO-s/a的反应氧化加成中间体R1-2OO-s/a和R1-6OO-s/a能发生自身异构化反应或与大气中的NO/HO2㊃发生双分子反应[45-46]㊂R1-2OO-s/a和R1-6OO-s/a所有可能的自身异构反应渠道及对应的ΔE0Kʂ如图4所示,反应类型共可分为3类:(1) OO 基团夺取 CH 或 CClH 的H生成 OOH;(2) OO 基团进攻 CH 的C生成双环中间体;(3) OO 基团夺取 CClH的Cl生成 OOCl(未能找到生成R1-2OOCl-s的TS,参考生成R1-6OOCl-s的能垒很高,认为未找到的TS不影响结论)㊂由图4可知,R1-2OO-s㊁R1-2OO-a和R1-6OO-s㊁R1-6OO-a经最佳自身异构渠道分别生R1-21OO-s㊁R1-2OOH12-a㊁R1-64OO-s和R1-61OO-a,需要克服30.9㊁26.7㊁23.2和24.7kcal㊃mol-1的反应能垒,对应的高压极限速率常数分别是1.7ˑ10-11㊁6.9ˑ10-4㊁9.6ˑ10-6和1.5ˑ10-6 s-1(温度为298K),远低于它们在高浓度NO条件下与NO及HO2㊃的伪一级反应速率常数k NO[NO]= 1.4s-1和k HO2[HO2㊃]=0.02s-1(基于[NO]=1.3ˑ1011 molecules㊃cm-3,[HO2㊃]=1.3ˑ109molecules㊃cm-3,kNO=1.2ˑ10-11cm3㊃molecule-1㊃s-1,k HO2=1.7ˑ10-11 cm3㊃molecule-1㊃s-1,温度为298K)[3,47-49]㊂因此,在高浓度NO条件下,R1-2OO-s/a和R1-6OO-s/a主要与大气中的NO发生双分子反应生成烷氧自由基(RO㊃)(R1-2O-s/a和R1-6O-s/a)和有机硝酸酯(C10H8ClONO2)㊂而在低浓度NO条件下,过氧自由基(R1-2OO-s/a和R1-6OO-s/a)与NO及HO2㊃的伪一级反应速率常数为kNO[NO]=0.001s-1和k HO2[HO2㊃] =0.02s-1(基于[NO]=1.3ˑ108molecules㊃cm-3,[HO2㊃] =1.3ˑ109molecules㊃cm-3,温度为298K)[3]㊂因此,在低浓度NO条件下,R1-2OO-s/a和R1-6OO-s/a主要与HO2㊃反应形成有机氢过氧化物(QOOH)和RO ㊃㊂因此,生成的RO2㊃除转化生成闭壳层产物有机氢过氧化物和有机硝酸酯外,还会生成RO㊃,接下来主要讨论RO㊃中间体R1-2O-s/a和R1-6O-s/a的反应㊂2.4㊀R1-2O-s/a和R1-6O-s/a的反应与其他自由基转化生成的RO㊃的反应类似,R1-2O-s/a和R1-6O-s/a也可以发生单分子自身异构及解离或进一步与大气中的O2反应[29,50]㊂所有考虑的R1-2O-s/a和R1-6O-s/a的单分子自身异构与解离及与O2反应的渠道及对应反应的ΔE0Kʂ如图5所示㊂可以看出,4种RO㊃存在3种单分子异构化反应类型:(1)RO㊃的O原子攻击不同的 CH 的C形成含O双环中间体;(2)C C键断裂生成相应的醛;(3)H迁移和C H键断裂,生成相应的酮㊂对比R1-2O-s/a和R1-6O-s/a所有单分子自身异构化与解离渠道反应ΔE0Kʂ可知,R1-2O-s/a和R1-6O-s/a 经环化生成R1-21O-s/a和R1-61O-s/a双环中间体是最可行的自身异构与解离反应渠道,对应的反应速表2 在ωB97XD/6-311++G(3df,2pd)//ωB97XD/6-31+G(d,p)计算水平下,R1+O2反应的反应能和能垒Table2㊀Reaction energy and barrier heights for reactions of R1with O2at theωB97XD/6-311++G(3df,2pd)//ωB97XD/6-31+G(d,p)level(kcal㊃mol-1)产物Products能垒Barrier heights(ΔE0Kʂ)自由能垒Free energy barriers(ΔG298Kʂ)能变Energy changes(ΔE0K)自由能变Free energy changes(ΔG298K)R1-2OO-syn 6.716.6-4.5 5.7R1-2OO-anti 5.415.7-5.1 5.0R1-4OO-syn26.037.224.635.5R1-4OO-anti23.334.021.632.5R1-6OO-syn8.619.2-5.2 5.8R1-6OO-anti 5.215.8-7.4 3.3R1-10OO-syn26.436.121.231.2R1-10OO-anti26.436.221.731.5R1-14OO-syn24.734.919.829.8R1-14OO-anti23.533.719.329.6C10H7Cl+HO2㊃12.422.4-18.9-12.3 . All Rights Reserved.第5期陈杰等:㊃Cl 引发萘的大气氧化机制及动力学19㊀图4㊀在ωB97XD/6-311++G(3df,2pd)//ωB97XD/6-31+G(d,p)计算水平下,R 1-2OO -s/a (a)和R 1-6OO -s/a (b)自身异构反应渠道和能垒(ΔE 0K ʂ)(单位为kcal ㊃mol -1)Fig.4㊀Barrier heights (ΔE 0K ʂ)for isomerization pathways of R 1-2OO -s/a (a)and R 1-6OO -s/a (b)atthe ωB97XD/6-311++G(3df,2pd)//ωB97XD/6-31+G(d,p)level (all in kcal ㊃mol -1)率常数分别为3.8ˑ108㊁8.3ˑ107㊁6.8ˑ1011和2.7ˑ1011s -1(298K)㊂对于RO ㊃和O 2反应,O 2能分别夺取R 1-2O -s ㊁R 1-2O -a 和R 1-6O -s ㊁R 1-6O -a 的C 2和C 6上的H 原子生成酮,反应的伪一级反应速率常数分别是1.7ˑ103㊁9.1ˑ102㊁6.2ˑ102和3.2ˑ102s -1(双分子反应速率常数分别为3.5ˑ10-16㊁1.9ˑ10-16㊁1.3ˑ10-16和6.5ˑ10-17cm 3㊃molecule -1㊃s -1及O 2浓度为4.9ˑ1018molecule ㊃cm -3,温度为298K),低于对应的最可行的自身异构与解离反应速率常数㊂因此,R 1-2O -s/a 和R 1-6O -s/a 的主要反应渠道是通过单分子环化形成相应的双环中间体R 1-21O -s/a 和R 1-61O -s/a ㊂. All Rights Reserved.20㊀生态毒理学报第16卷图5㊀在ωB97XD/6-311++G(3df,2pd)//ωB97XD/6-31+G(d,p)计算水平下,R 1-2O -s/a和R 1-6O -s/a 单分子反应及与O 2的双分子反应能垒(ΔE 0K ʂ)(单位为kcal ㊃mol -1)Fig.5㊀Barrier heights (ΔE 0K ʂ)for unimolecular reaction of R 1-2O -s/a and R 1-6O -s/a and bimolecular reactions ofR 1-2O -s/a and R 1-6O -s/a with O 2at the ωB97XD/6-311++G(3df,2pd)//ωB97XD/6-31+G(d,p)level (all in kcal ㊃mol -1). All Rights Reserved.第5期陈杰等:㊃Cl引发萘的大气氧化机制及动力学21㊀3㊀环境意义(Implications)本研究基于量子计算化学和过渡态理论,首次揭示了㊃Cl引发Nap的大气转化机制和动力学㊂计算得到的㊃Cl+Nap的主要的反应路径如图6所示㊂可以看出,引发反应主要为㊃Cl加成到Nap的C5位置生成㊃C10H8Cl㊂加成中间体㊃C10H8Cl后续主要与O2通过syns和anti这2个方向加成到C2和C6位置生成4种RO2㊃中间体R1-2OO-s/a和R1-6OO-s/a㊂在低NO浓度条件下,RO2㊃主要和HO2㊃反应生成氢过氧化合物(QOOH)和烷氧自由基(RO㊃);在高NO浓度条件下,RO2㊃将主要与NO反应生成RO㊃(R1-2O-s/ a和R1-6O-s/a)和有机硝酸酯(C10H8ClNO3)㊂重要的是,使用ECOSAR模型预测生成的有机氢过氧化合物㊁有机硝酸酯的水生毒性强于其母体化合物Nap (有机氢过氧化合物㊁有机硝酸酯和Nap对鱼的96h 半致死浓度(LC50)分别为0.92㊁1.54和9.39mg㊃L-1,对绿藻的96h半效应浓度EC50为0.42㊁1.95和6.91mg ㊃L-1),表明㊃Cl引发Nap转化增加了Nap释放的环境风险[36]㊂生成的4种RO㊃主要是通过环化生成相应的双环产物R1-21O-s/a和R1-61O-s/a㊂因此,为进一步完善㊃Cl引发Nap的大气转化机制,将来需要研究双环中间体的后续转化机制和动力学;同时,为准确评估㊃Cl引发Nap转化导致的环境风险,将来有必要对转化过程中生成的闭壳层产物的毒性展开研究㊂本文的研究结果为全面评估㊃Cl引发萘的大气转化及环境风险提供了数据支持㊂图6㊀㊃Cl+Nap的主要反应路径Fig.6㊀Main pathways for the㊃Cl initiated reaction of Nap致谢:感谢大连理工大学超算中心给予的计算资源支持㊂通讯作者简介:谢宏彬(1980 ),男,博士,教授,主要研究方向为环境理论化学㊂参考文献(References):[1]㊀Faxon C B,Allen D T.Chlorine chemistry in urban at-mospheres:A review[J].Environmental Chemistry,2013,10(3):221[2]㊀Young C J,Washenfelder R A,Edwards P M,et al.Chlo-rine as a primary radical:Evaluation of methods to under-stand its role in initiation of oxidative cycles[J].Atmos-pheric Chemistry and Physics,2014,14(7):3427-3440 [3]㊀Liu C,Ma F F,Elm J,et al.Mechanism and predictivemodel development of reaction rate constants for N-centerradicals with O2[J].Chemosphere,2019,237:124411 [4]㊀Atkinson R,Baulch D L,Cox R A,et al.Evaluated kinet-ic and photochemical data for atmospheric chemistry:SupplementⅣ:IUPAC subcommittee on gas kinetic dataevaluation for atmospheric chemistry[J].Atmospheric. All Rights Reserved.22㊀生态毒理学报第16卷Environment Part A General Topics,1992,26(7):1187-1230[5]㊀Keene W C,Khalil M A K,Erickson D JⅢ,et -posite global emissions of reactive chlorine from anthro-pogenic and natural sources:Reactive Chlorine EmissionsInventory[J].Journal of Geophysical Research:Atmos-pheres,1999,104(D7):8429-8440[6]㊀Finlayson-Pitts B J.Chlorine chronicles[J].Nature Chem-istry,2013,5(8):724[7]㊀Knipping E M.Experiments and simulations of ion-en-hanced interfacial chemistry on aqueous NaCl aerosols[J].Science,2000,288(5464):301-306[8]㊀Sommariva R,von Glasow R.Multiphase halogen chem-istry in the tropical Atlantic Ocean[J].EnvironmentalScience&Technology,2012,46(19):10429-10437[9]㊀Lawler M J,Finley B D,Keene W C,et al.Pollution-en-hanced reactive chlorine chemistry in the eastern tropicalAtlantic boundary layer[J].Geophysical Research Letters,2009,36(8):L08810[10]㊀Liu X X,Qu H,Huey L G,et al.High levels of daytimemolecular chlorine and nitryl chloride at a rural site onthe North China plain[J].Environmental Science&Technology,2017,51(17):9588-9595[11]㊀Thornton J A,Kercher J P,Riedel T P,et al.A large a-tomic chlorine source inferred from mid-continental reac-tive nitrogen chemistry[J].Nature,2010,464(7286):271-274[12]㊀Mielke L H,Furgeson A,Osthoff H D.Observation ofClNO2in a mid-continental urban environment[J].Envi-ronmental Science&Technology,2011,45(20):8889-8896[13]㊀Zhou W,Zhao J,Ouyang B,et al.Production of N2O5and ClNO2in summer in urban Beijing,China[J].At-mospheric Chemistry and Physics,2018,18(16):11581-11597[14]㊀Le Breton M,HallquistÅM,Pathak R K,et al.Chlorineoxidation of VOCs at a semi-rural site in Beijing:Signifi-cant chlorine liberation from ClNO2and subsequent gas-and particle-phase Cl-VOC production[J].AtmosphericChemistry and Physics,2018,18(17):13013-13030 [15]㊀Xie H B,Ma F F,Wang Y F,et al.Quantum chemicalstudy on㊃Cl-initiated atmospheric degradation of mono-ethanolamine[J].Environmental Science&Technology,2015,49(22):13246-13255[16]㊀Ma F F,Ding Z Z,Elm J,et al.Atmospheric oxidation ofpiperazine initiated by㊃Cl:Unexpected high nitrosamineyield[J].Environmental Science&Technology,2018,52(17):9801-9809[17]㊀Guo X R,Ma F F,Liu C,et al.Atmospheric oxidationmechanism and kinetics of isoprene initiated by chlorineradicals:A computational study[J].Science of the TotalEnvironment,2020,712:136330[18]㊀Shen H Z,Huang Y,Wang R,et al.Global atmospheric e-missions of polycyclic aromatic hydrocarbons from1960to2008and future predictions[J].Environmental Science&Technology,2013,47(12):6415-6424[19]㊀Reisen F,Arey J.Atmospheric reactions influence season-al PAH and nitro-PAH concentrations in the Los Angelesbasin[J].Environmental Science&Technology,2005,39(1):64-73[20]㊀Lammel G,KlánováJ,Ilic'P,et al.Polycyclic aromatichydrocarbons in air on small spatial and temporalscales II.Mass size distributions and gas-particle parti-tioning[J].Atmospheric Environment,2010,44(38):5022-5027[21]㊀杨春,杨琦,杨素银,等.萘好氧降解菌的筛选及降解特性的初步研究[J].环境科学与技术,2005,28(6):19-21,110Yang C,Yang Q,Yang S Y,et al.Strains isolation andstudy on aerobic biodegradability of naphthalene[J].En-vironmental Science and Technology,2005,28(6):19-21,110(in Chinese)[22]㊀Sasaki J,Aschmann S M,Kwok E S C,et al.Products ofthe gas-phase OH and NO3radical-initiated reactions of naphthalene[J].Environmental Science&Technology,1997,31(11):3173-3179[23]㊀Huang G C,Liu Y,Shao M,et al.Potentially importantcontribution of gas-phase oxidation of naphthalene andmethylnaphthalene to secondary organic aerosol duringhaze events in Beijing[J].Environmental Science&Technology,2019,53(3):1235-1244[24]㊀Bunce N J,Liu L N,Zhu J,et al.Reaction of naphthaleneand its derivatives with hydroxyl radicals in the gas phase[J].Environmental Science&Technology,1997,31(8):2252-2259[25]㊀Riva M,Healy R M,Flaud P M,et al.Gas-and particle-phase products from the chlorine-initiated oxidation ofpolycyclic aromatic hydrocarbons[J].The Journal ofPhysical Chemistry A,2015,119(45):11170-11181 [26]㊀Lee J Y,Lane D A.Unique products from the reaction ofnaphthalene with the hydroxyl radical[J].AtmosphericEnvironment,2009,43(32):4886-4893[27]㊀Kautzman K E,Surratt J D,Chan M N,et al.Chemicalcomposition of gas-and aerosol-phase products from thephotooxidation of naphthalene[J].The Journal of Physi-cal Chemistry A,2010,114(2):913-934. All Rights Reserved.第5期陈杰等:㊃Cl引发萘的大气氧化机制及动力学23㊀[28]㊀Ouyang B,Fang H J,Dong W B,et al.Different mecha-nisms both lead to the production of the naphthalene-OHadduct in the355nm and266nm laser flash photolysisof the mixed aqueous solution of naphthalene and nitrousacid[J].Journal of Photochemistry and Photobiology A:Chemistry,2006,181(2-3):348-356[29]㊀Zhang Z J,Lin L,Wang L M.Atmospheric oxidationmechanism of naphthalene initiated by OH radical.A the-oretical study[J].Physical Chemistry Chemical Physics,2012,14(8):2645-2650[30]㊀Riva M,Healy R M,Flaud P M,et al.Kinetics of the gas-phase reactions of chlorine atoms with naphthalene,ace-naphthene,and acenaphthylene[J].The Journal of Physi-cal Chemistry A,2014,118(20):3535-3540[31]㊀刘聪,马芳芳,付自豪,等.㊃Cl引发3种环状含有NH结构有机化合物的大气转化机制及动力学[J].生态毒理学报,2019,14(4):65-72Liu C,Ma F F,Fu Z H,et al.Atmospheric transformationmechanism and kinetics of three cyclic NH-containingcompounds initiated by㊃Cl[J].Asian Journal of Ecotoxi-cology,2019,14(4):65-72(in Chinese)[32]㊀Li C,Xie H B,Chen J W,et al.Predicting gaseous reac-tion rates of short chain chlorinated paraffins with㊃OH:Overcoming the difficulty in experimental determination[J].Environmental Science&Technology,2014,48(23):13808-13816[33]㊀Xu T,Chen J W,Chen X,et al.Prediction models on pKaand base-catalyzed hydrolysis kinetics of parabens:Ex-perimental and quantum chemical studies[J].Environ-mental Science&Technology,2021,55(9):6022-6031 [34]㊀Xu T,Chen J W,Wang Z Y,et al.Development of predic-tion models on base-catalyzed hydrolysis kinetics ofphthalate esters with density functional theory calculation[J].Environmental Science&Technology,2019,53(10):5828-5837[35]㊀Li C,Chen J W,Xie H B,et al.Effects of atmosphericwater on㊃OH-initiated oxidation of organophosphateflame retardants:A DFT investigation on TCPP[J].Envi-ronmental Science&Technology,2017,51(9):5043-5051 [36]㊀United States Environmental Protection Agency(USEPA).ECOSAR V2.0[CP/OL].[2021-03-26].https:///tsca-screening-tools/ecological-structure-ac-tivity-relationships-ecosar-predictive-model[37]㊀Frisch M J,Trucks G W,Schlegel H B,et al.Gaussian09[CP].Wallingford,CT:Gaussian,Inc,2009[38]㊀Chai J D,Head-Gordon M.Long-range corrected hybriddensity functionals with damped atom-atom dispersioncorrections[J].Physical Chemistry Chemical Physics,2008,10(44):6615[39]㊀Barker J R,Nguyen T L,Stanton J F,et al.MultiWellprogram suite[CP].Ann Arbor,MI:University of Michi-gan,2014[40]㊀Barker J R.Multiple-well,multiple-path unimolecular re-action systems.Ⅰ.MultiWell computer program suite[J].International Journal of Chemical Kinetics,2001,33(4):232-245[41]㊀Carl E.The penetration of a potential barrier by electrons[J].Physical Review,1930,35(11):1303-1309[42]㊀Raman S,Ashcraft R W,Vial M,et al.Oxidation of hy-droxylamine by nitrous and nitric acids.Model develop-ment from first principle SCRF calculations[J].The Jour-nal of Physical Chemistry A,2005,109(38):8526-8536 [43]㊀Braña P,Sordo J A.Mechanistic aspects of the abstractionof an allylic hydrogen in the chlorine atom reaction with2-methyl-1,3-butadiene(isoprene)[J].Journal of the A-merican Chemical Society,2001,123(42):10348-10353 [44]㊀Sun C H,Xu B E,Zhang S W.Atmospheric reaction ofCl+methacrolein:A theoretical study on the mechanism,and pressure-and temperature-dependent rate constants[J].The Journal of Physical Chemistry A,2014,118(20):3541-3551[45]㊀Wang L Y,Wang L M.Atmospheric oxidation mechanismof acenaphthene initiated by OH radicals[J].AtmosphericEnvironment,2020,243:117870[46]㊀Dang J,He M X.Mechanisms and kinetic parameters forthe gas-phase reactions of anthracene and pyrene with Clatoms in the presence of NOx[J].RSC Advances,2016,6(21):17345-17353[47]㊀Wennberg P O,Bates K H,Crounse J D,et al.Gas-phasereactions of isoprene and its major oxidation products[J].Chemical Reviews,2018,118(7):3337-3390[48]㊀Patchen A K,Pennino M J,Elrod M J.Overall rate con-stant measurements of the reaction of chloroalkylperoxyradicals with nitric oxide[J].The Journal of PhysicalChemistry A,2005,109(26):5865-5871[49]㊀Fu Z H,Xie H B,Elm J,et al.Formation of low-volatileproducts and unexpected high formaldehyde yield fromthe atmospheric oxidation of methylsiloxanes[J].Envi-ronmental Science&Technology,2020,54(12):7136-7145[50]㊀Atkinson R.Atmospheric reactions of alkoxy and-hydroxyalkoxy radicals[J].International Journal ofChemical Kinetics,1997,29(2):99-111Ң. All Rights Reserved.。
atmospheric environment审稿三天Atmospheric environment refers to the combination of air, water, soil, climate and other environmental elements that surround us. It is an important factor that affects human health, the environment and economic development.Environmental protection is the key to maintain a good atmospheric environment. To start with, advanced technology and effective management should be employed in industrial and agricultural production to reduce pollutants. Besides, the public should be educated on environmental awareness, reduce the waste of energy and natural resource, and persuade people to abandon outdated and polluting equipment. Moreover, advanced ecological engineering should be employed to control environmental degradation, such as water and air pollution caused by industrial and residential activities.Apart from protecting the existing environment, remember to restore the degraded environment. Reforestation and grassland restoration should be conducted to increase the vegetation cover and improve the local environment. Local governments should also take necessary actions to protect the environment from human activities.In conclusion, it is vital that we take effective steps to protect the atmospheric environment. With the joint efforts of the governments, enterprises, and individuals, hopefully, a better environment can be achieved.。
Submesoscale Coastal Ocean Flows Detected By Very High Frequency Radar and Autonomous Underwater VehiclesLynn K.Shay1,Thomas M.Cook1,P.Edgar An21Division of Meteorology and Physical OceanographyRosenstiel School of Marine and Atmospheric ScienceUniversity of MiamiMiami,FL331492Department of Ocean EngineeringFlorida Atlantic UniversityDania Beach,FL33004Accepted in Journal of Atmospheric and Oceanic TechnologyCorresponding Author Address:Dr.Lynn K.Shay4600Rickenbacker CausewayDivision of Meteorology and Physical OceanographyRosenstiel School of Marine and Atmospheric ScienceMiami,FL33149Email:nick@April10,20031AbstractOver a29-d time series in July1999,an Ocean Surface Current Radar(OSCR)in very high frequency (VHF)mode mapped the surface velocityfield at250m resolution at700cells offuderdale,Florida. During the experiment,autonomous underwater vehicles(AUV),equipped with upward and downward-looking1.2MHz Acoustic Doppler Current Profilers(ADCP),measured subsurface current structure over four to six radar cells during two mixed layer patterns on9and27July1999.As these AUV sampling patterns were conducted over500m x500m,and500m x750m areas,these missions required about80 to90min(four radar sample intervals)to form four and seven synoptic snapshots,respectively.Based on autocorrelation analyses of the profiler data,along-AUV track subsurface profiles were aver-aged at10-sec intervals,mapped to a surface from1.5to6.5m,and compared to surface currents at more than500points for each parisons between the surface and subsurface currents from the AUV revealed spatially averaged differences ranging from4to26cm s−1during these two experiments. The largest differences occurred when the surface and subsurface current vectors were orthogonal,other-wise differences were O(10cm s−1).Scatterplots between2-m and radar-derived surface currents indicated a consistent relationship with mooring data.From the seven spatial snapshots acquired during the second experiment,current profiles suggested a time-dependent oscillation that was corroborated by radar and moored ADCP data.Least-squarefits of these profiles from sequential AUV snapshots to a simple model isolated an∼9.2±1h oscillation where the along-shelf current was O(50cm s−1).Spatially averaged current profiles from four and seven snapshots were subsequently time-averaged to form a mean profile from each experiment.In the down-wind directions,these mean profiles were compared to a wind-driven,logarithmic layer profile in the upper6.5m based on a10-m surface winds.Regression analyses suggest a slope of≈1.16between the theoretical and observed mean profiles with a bias of about 3cm s−1.In this context,the averaged winds played a role in driving the coastal ocean circulation. These results further suggest that the spatial averaging by the radar is consistent when subsurface current variations are averaged over similar space and time scales.21IntroductionGiven escalating interest in establishing coastal ocean observatories,there has been an increased em-phasis on connecting observatories using High Frequency(HF)radar techniques to map submesoscale to mesoscale surface velocityfields.These measurements have provided spatial context for the more conven-tional observational approaches using moorings,drifters,and ships during recent coastal experiments.An important measurement issue emerging from previous studies has been the range of the differences(4to 25cm s−1)between the HF radar current measurements and those observed by Eulerian and Lagrangian methods(Paduan and Graber1997).Comparisons to conventional techniques,however,do not necessarily provide sufficent spatial coverage at high temporal resolution to understand areal averaging by a Doppler radar and to assess its relationship to upper ocean subsurface structure over similar spatial scales.For example,RMS differences between surface currents from an Ocean Surface Current Radar(OSCR)and subsurface currents(4-m)from a Vector Measuring Current Meters(VMCM)were7cm s−1over a1m s−1range from observations acquired during the Duck94experiment(Shay et al.1998a).Since measurement accuracies of a VMCM are2cm s−1(Weller and Davis1980),these differences are within the cited accuracy of4to5cm s−1for OSCR measurements.Teague et al.(2001)recently described a set of measurements acquired during the third Chesapeake Bay Outflow Plume Experiment(COPE-3)from the University of Michigan’s MultiChannel Radar system(MCR),which measures bulk surface currents over differing parisons atfive Acoustic Doppler Current Profiler(ADCP)moorings in the radar domain revealed RMS differences ranging from5to20cm s−1depending on the frequency and location within the radar domain.During this same period,subsurface current comparisons to OSCR-derived surface currents over a14-d time series ranged between7to11cm s−1RMS where the range of current variability exceeded1m s−1(Shay et al.2003).Previous results point to a reasonably good estimate of a radar-derived surface current compared to a subsurface current observation.However,a moored measurement is a single point observation whereas a radar-derived surface current represents a spatial average over areas with dimensions of0.6to4km2 depending on the horizontal resolution of the radar in either Very High Frequency(VHF)or HF modes,3respectively.A key scientific issue in using HF-radars has been the relationship between radar-derived sur-face current estimates and those determined from subsurface measurements over similar spatial footprints. The approach described herein attempts to examine this radar measurement issue.Recent surface current observations from OSCR in VHF mode revealed complex surface current patterns in the South Florida Ocean Measurement Center(SFOMC)area located offuderdale,Florida (Shay et al.2000).These high-resolution surface current observations provided spatial context for repeated velocity profiles acquired by an Autonomous Underwater Vehicle(AUV)over small-scale grids at a speed of1.2m s−1at9-m depth(An et al.2001;Dhanak et al.2001).Downward and upward-looking ADCP moorings and ship-based ADCP and conductivity-temperature-depth(CTD)measurements were also acquired during these AUV-based sampling schemes to examine submesoscale coastal ocean processes forced by energetic Florida Current(FC)intrusions and vortices across the shelf break(Peters et al.2002; Shay et al.2002).An AUV provides a novel approach to acquire oceanographic data to observe and understand sub-mesoscale variability in the current and density structure(An et al.2001).AUV-based measurements allow for a diversity of spatial-temporal resolution of oceanic variables over kilometer scales or less.This subsurface sampling approach benefits from improved platform stability and maneuverability of the AUV with enhanced positional capabilities in oceanic environments.To enable effective operations,an adaptive sampling strategy has to be chosen carefully to maximize the utility of AUV sampling strategy relative to the imposed surface boundary conditions derived from radar-derived surface measurements.A fundamen-tal step in developing AUV sampling algorithms is to characterize the spatially evolving environmental structure at the appropriate scales(Curtin et al.1993).To improve our understanding of this spatial averaging by VHF(and HF)radars,subsurface current profiles from repeated ADCP transects from an AUV underneath4to6cells are compared to radar-derived surface currents.As these subsurface mapping patterns required≈80min(4∆t,where∆t is the radar sample at a20-min interval)to complete,spatial variability in coastal oceanicflows was acquired at high-spatial resolution inshore of the FC as it intruded across the shelf break.Accordingly,the experimental4design is described in Section2.In Section3,spatially averaged current snapshots from two AUV exper-iments are used to isolate a10-h oscillation found in the radar,mooring and AUV data,and to compare the radar-derived surface current to the AUV-derived subsurface current.In section4,time-mean profiles from each experiment are compared to theoretical log-layer estimates in the down-wind directions.Results are summarized in Section5with concluding remarks.2Surface and Subsurface MeasurementsThe experiment was conducted in July1999from Port Everglades to Hollywood Beach,Florida which encompassed the SFOMC domain.In this section,VHF radar and AUV sensor suite are described where the emphasis is on observations to resolve submesoscale coastal ocean current variability(Table1).Given the north-south alignment of the south Florida coastline,currents will be referred to as cross-shelf(U)and along-shelf(V)to represent east-west and north-south directions,respectively.2.1Surface Wind ObservationsSurface atmospheric conditions observed from NOAA Coastal Monitoring Automated Network(CMAN) stations at Fowey Rocks and Lake Worth,and at the Northeast(NE)Mooring(Table1)were used to adjust surface winds to the10-m level.While the CMAN stations recorded data at hourly intervals,the NE surface mooring acquired data at15-min intervals beginning on15July1999.During thefirst mixed layer experiment(ML1),surface winds at Fowey Rocks ranged between5to9m s−1towards the west between yeardays190to191(9to10July)(Fig.2a).This onshore wind averaged about6.5m s−1.Based on Fairall et al.(1996),surface friction velocities ranged from0.15to0.27m s−1at Fowey Rocks(Fig. 2c).By contrast,northward surface winds were weaker during the ML2experiment(yeardays208to209; 27to28July)as shown in Figs.2a,b.Winds were directed towards the north at speeds of4to5m s−1 where u over this12-h period ranged from0.02m s−1to as high as0.16m s−1at Fowey Rocks.Surface frictional velocities at the NE mooring(shown below)ranged between0.05to0.2m s−1based on adjusted 10-m records.Notice the marked similarities in the adjusted surface winds between Fowey Rocks and the NE mooring site separated by a distance of65km(An et al.2001).Surface winds observed at the Lake Worth CMAN station(not shown)indicated similar trends,suggestive of coherent wind patterns.52.2VHF RadarThe OSCR radar system was deployed from25June to10August1999.During this period,a29-d continuous time series of vector surface currents were acquired at20-min intervals starting on9July and ending7August1999(Shay et al.2002).The system consisted of two VHF radar transmit and receive stations operating at49.945MHz that sensed electromagnetic signals scattered from surface gravity waves with a Bragg wavelength of2.95m(radar wavelength is5.9m).Note that the Bragg frequency is0.72Hz, and the offset of thefirst-order spectral peak from this Bragg frequency is proportional to the radial current for a Bragg wave advancing(positive)or receding(negative)from the radar station(Stewart and Joy1974). The radar system mapped coastal ocean currents over a7.5km×8km domain with a horizontal resolution of250m at700grid points(Fig.1).Radar sites were located in John U.Lloyd State Park(Master)and at an oceanfront site in Hollywood Beach FL(Slave),yielding a baseline distance of6.7km.Each site consisted of a four-element transmit and thirty-element receiving array oriented at an angle of37◦(SW-NE at Master)and160◦(SE-NW at Slave).A total of2078snapshots of the2-dimensional surface current vector were acquired from0320GMT9July1999until2340GMT6August yielding a29-d time series.Of the2078samples,only69samples were missing from this time series,equating to a3.3%data loss in the snapshots.Previous experiments have yielded similar results ranging from93to97%data return(Haus et al.1998).The geometric dilution of precision(GDOP)provides an estimate of the spatial dependence of observed surface current differences based on geometrical constraints(Chapman et al.1997).Using the radar’s mean look direction(α),and the half angle(φ)between intersecting beams,expressions for the currentcomponents areσv= 2 sin2(α)sin2(φ)+cos2(α)cos2(φ)sin(2φ) 1σ,(1)andσu= 2 cos2(α)sin2(φ)+sin2(α)cos2(φ)sin2(2φ) 12σ,(2)whereσrepresents rms current differences.The GDOP is a nondimensional number defined by the ratiosofσvσandσuσin the along-shelf and cross-shelf components,respectively.As shown in Fig.3a,the GDOP6ranged from0.75to2for each component.In the core of the domain where a large fraction of the subsurface measurements were acquired,the GDOP for each radial current component was O(1).Close to the coast, however,the GDOP increased to a value of2over a1.5to2km distance as intersection angles(Fig.3b) approached the outer limits of the envelope required to construct a vector currentfield from radial current data.In this VHF domain,optimal intersections angles are defined here as between30◦and150◦.These angles encompassed most of the domain except for the grid points closest to the shore,and those just beyond the40◦limits in the farfield(Shay et al.2002).2.3Surface Velocity FieldsAs shown in Fig.4,an energetic surface current regime encompassed the second AUV experiment (ML2)from27to28July(yeardays208to209)(Table2).The AUV measurement domain will be referred to hereafter as the Intensive Observational Domain(IOD).In the inner region of the IOD,surface currents of30cm s−1were towards the south whereas offshore currents associated with the FC exceeded1.4m s−1(Fig.4a).This surface current pattern indicated large spatial gradients over small scales of<2km particularly along the shelf break where vorticities were large compared to f(Peters et al.2002).Four hours later,surface currents decreased over the IOD,however,the FC remained in the radar footprint where currents were1.4m s−1shown in Fig.4b.A divergence region formed between the inner and outer parts of the IOD as currents were towards the northwest and east directions,respectively(Fig.4c). This surface current divergence was associated with a cyclonically rotating ocean feature evident in the central portion of the radar domain.As shown in Fig.4d,surfaceflows indicated northwestwardflows over the inner and outer portions of the IOD.By0000UTC28July(yearday208),the surfaceflow pattern suggested a cyclonically rotating feature with60to70cm s−1tangential currents(Fig.4e)and had marked asymmetric structure similar to the previously observed vortices.In the inner portion of the IOD,surface currents were about30cm s−1directed towards the northeast,but were weaker(≈10cm s−1)than those over the outer part of the IOD.Forty minutes later(Fig.4f),a cyclonically rotating feature influenced the surfaceflows over the IOD.This asymmetric surface velocity structure supported wave-like structures and vortices that appeared to be propagating along the inner edge of the FC(Peters et al.2002).These features may be more similar to the submesoscale vortices(Shay et al.1998b)than to large spin-offeddies7(Lee and Mayer1977).2.4StratificationTo examine the upper-ocean stratification,the NE mooring data were used to determine the mean temperature and density conditions observed from a series of micro-cats at1,5,10,15,20,30and40m sampling at5-min intervals(Fig.5).Temperatures between5and10m indicated similar time variations during ML2when weakly stratified conditions persisted.Since the AUV measuredfields at9-m depth, mean profiles(defined here as the average conditions over the12-h period)revealed oceanic conditions between1and10m,with small vertical density differences of less than5x10−4gm cm−3.Ship-board CTD data indicated agreement in the temperature profiles,but with slightly different density profiles due to salinity differences.Stronger stratification was evident beneath10m depth in the CTD and mooring data(Peters et al.2002).This implies that the ML2measurements beginning at1.5m and extending to6.5m were contained within a weakly stratified layer as opposed to a well mixed surface layer.The episodic nature of the FC intrusions across the shelf break makes it difficult to define a well-mixed surface layer.For the purposes of this manuscript,the AUV missions will be referred to as thefirst(ML1)and second(ML2)mixed layer experiments.2.5AUV-Based MeasurementsThe Ocean Explorer(OEX)is a compact AUV that was used to map subsurface coastal ocean structure (An et al.2001;Dhanak et al.2001).Upgrades to the AUV sensor suite may be done in a plug-and-play mode,resulting in scalable and reconfigurable systems.At a speed of1.2m s−1at9m±5cm beneath the surface,the AUV mapped submesoscale,3-dimensional current variability using a pair of upward-and downward-looking1.2MHz ADCPs for nearly continuous periods ranging from6to24h within the IOD (Table2).The upward-looking ADCP was programmed to have0.5m bins for the Mixed Layer(ML)and Tidal(T)Missions and1.5m bins in the downward-looking direction.For the Bottom Boundary Layer (BBL)Mission,higher resolution bins of0.5m were set in the downward-looking ADCP.The downward-looking ADCP also performed bottom-tracking,which provided underwater navigation.Sampling rates for these ADCPs were set to1Hz for the1.5and0.5m bins.A Falmouth Scientific Instrument Conductivity8Temperature and Depth(CTD)was used to map the temperature,salinity and densityfields at9m. During all missions,scientific and vehicle data were acquired by the OEX including time,temperature, salinity,current,AUV position(latitude,longitude,depth,altitude),and motion parameters(velocity, acceleration,angular rate and attitude).The OEX was acoustically tracked from a chase boat using its onboard Trackpoint II system where data were integrated into the boat’s Global Positioning System(GPS) andflux-gate compass.Within the IOD(see Fig.1),the OEX sampled subsurface current velocities in water depths ranging from20to35m.For both ADCPs,only raw data were acquired(no ensemble averaging)from the AUV. The ML1mission was conducted over a period of6h and15min,and ML2was completed over a12-h period,that also included a turbulence package as part of the payload(Dhanak et al.2001).The AUV pre-programmed,lawn-mower pattern was repeatedly executed until its battery energy was depleted,resulting in four and seven snapshots for the ML1(500m x500m)and ML2(500m x750m)missions,respectively. During these missions,the OEX surfaced to acquire differential GPSfixes that bounded positional errors of≈25m(typically10%of radar grid resolution).The time to complete one complete pattern was approximately80-min or4radar samples(∆t=20min).To facilitate a comparison to the radar-derived surface current measurements,AUV-derived current profiles were smoothed in the along AUV-track directions.As shown in Fig.6,differing time scales of10, 20and30s,corresponding to horizontal length scales of12,24and36m for a steadily moving vehicle, were tested using autocorrelation methods.For an AUV speed of1.2m s−1,the10-s averaging period was optimal since autocorrelation levels were above0.9,yielding a12-m averaging interval for each profile. Based upon a1Hz sample interval for the0.5m bin data,measurement uncertainty over a10-s interval was approximately4cm s−1whereas it was less than2cm s−1for the downward-looking ADCP averaged over a1.5m bin.The smallest resolvable horizontal wavelength in the measurements was24m.For each synoptic snapshot,these10-s averaged data were gridded to a surface at the various depths beginning at 1.5m depth(kinetic energy differences of the mapped currentfields were typically1to7%less than the observedfield).Contamination by side-lobe,surface reflections may induce additional uncertainty into near-surface ADCP measurements.For example,the range is given by H(1-cos(θ))where H is∼9m and9the tranducer angle(θ)is30◦(Teague et al.2001).Thus,thefirst bin of usable AUV’s ADCP data is estimated to be1.5to2m.Radar-derived surface velocityfields were mapped to the same surface for comparison purposes over each80-min synoptic snapshot at more than500points as described below.3Surface and Subsurface Comparisons3.1Vertical StructureAlong one section during the ML2mission(BB in Fig.1),velocity profiler data from1.5to6.5m were spatially averaged along an east-west transect from each snapshot,and compared to the southwest (SW)and northeast(NE)ADCP mooring data(averaged over80-min time intervals)as shown in Fig.7. Cross-shelfflows were40to50%less energetic than along-shelfflows associated with the FC.The upper ocean data suggests a depth-independentflow,consistent with ship transect data.Although AUV profiling was conducted over a12-h period,these spatially averaged profiles indicated a current oscillation in the time-dependent variability from snapshot to snapshot.This energetic oscillation was apparent in the along-shelfflows with an≈10-h period,which can be distinguished from weaker semidiurnal tidal currents of 4to6cm s−1(An et al.2001).These motions were also observed at the NE and SW moorings located along the50-m and20-m isobaths,respectively(see Fig.1).The AUV snapshots captured this wave-like variability where the wavelength was found to be about30km for a wave period of≈10-h period based on spectral analysis(Peters et al.2002).At the SW mooring,this oscillation was similar but with a shifted phase and a smaller amplitude presumably due to the importance of bottom friction in the shallower water. These time-dependent current transects from the AUV are similar to the oscillatory motions observed at the moorings.To isolate this observed wave-like variability,observed currents from the moorings and AUV were layer averaged(as per Fig.7)from3to6m and1.5to6.5m,respectively over the≈80-min time scale.These layer-averaged data(from each sequential snapshot)werefit to a model(Rossby and Sanford1976)to determine the carrier frequency of the wave form.For each trial frequency starting at0.8σo and ending at1.2σo,layer-averaged,cross-shelf currents werefit to an expression of the formU i(t)={A u cos(σoδt+θu)+B u sin(σoδt+θu)}+u r(t),(3)10and,with a similar expression for the along-shelf componentV i(t)={A v sin(σoδt+θv)+B v cos(σoδt+θv)}+v r(t),(4) where A u,v and B u,v represent Fourier amplitudes for the cross-shelf and along-shelf components,respec-tively,δt≈80-min represents the time interval to complete one snapshot,θu,v is the phase angle,and u r and v r represents unresolved currents(Marquardt1963).The carrier frequency is defined as the frequency minimizing the covariance between the unresolved components(<u r,v r>).Sensitivity tests have been reported elsewhere using a wave of known frequency(period)with an amplitude of1and0.2m s−1and a random noise component of10%of these input signals.The uncertainty in resolving the period was found to be±1h.Results of the least-squaresfit are listed in Table3and shown in Fig.8.In the cross-shelf direction, amplitudes from the mooring data were∼16cm s−1at the NE mooring(50-m)whereas the amplitude at the SW mooring(20-m)was6cm s−1.The cross-shelf amplitude(A u)derived from the AUV was 14cm s−1,which agreed with the A y determined from the50-m mooring as phase angles(θu)increased shoreward.By contrast,along-shelf current amplitudes ranged from32cm s−1at the SW mooring to44 cm s−1at the NE mooring.The amplitude derived from the seven AUV snapshots was52cm s−1.In this component,the phase angle decreased towards the coast.The periods of the carrier signals at the mooring sites ranged from10to11h compared to about9.2h from the12-h of AUV ing a29-d time series of surface current data over the NE mooring,Peters et al.(2002)found a10.2h period oscillation in variance-conserving spectra from the surface current data where the horizontal wavelength was30km. Thus,this9.2h period is within the uncertainty of the approach of≈1h in resolving the wave period. RMS amplitudes defined here as (<u r,v r>)ranged from4to14cm s−1.These larger differences atthe mooring sites were due to the lack of time dependence in model amplitudes(3,4).That is,amplitudes were assumed to have the same value over a3-d time series used in thefit.The4cm s−1difference suggestsa goodfit between the model and the AUV current data.3.2Radar and AUV-Derived Current SnapshotsThe mapped velocityfields from the radar and the AUV are compared at each of these points to assess11the relationship between subsurface structure and radar-derived surface velocityfield.During thefirst snapshot of the ML1experiment(Fig.9a),currents revealed a predominant northwardflow and evidence of an anticylonic current veering with depth.The anticyclonic current veering during the second snapshot (Fig.9b)became more pronounced as the angle between the surface and subsurface current vectors increased to more than45◦.As shown in Fig.9c,weaker surface currents were directed towards the coast compared to the northwardflows in thefirst snapshot.The2-m subsurface currents decreased from50 cm s−1to less than15cm s−1.Current veering with depth was erratic at times due to weaker cross-shelf flows.By the fourth snapshot(Fig.9d),currents reversed direction from a northward to a southwardflow, but did not indicate any dominate spatial trends.Over the six hours,these snapshots indicated a gradual cyclonic turning of the currentfield as observedflows were towards the north,northwest,west,and south to southwest directions.During the ML2experiment,surface and subsurface currents were aligned in the same direction with weak currents at both levels(Fig.10a).As shown in Fig.10b,currentsflowed towards the south as the 2-m currents were more energetic than those observed in the surface layer.Surface currents subsequently increased and were directed towards the southeast whereas subsurface currents remained directed towards the south(Fig.10c).Surfaceflows were less energetic than those observed at2-m during the fourth snapshot(Fig.10d).During the third and fourth snapshots,subsurface currents were nearly orthogonal to the surface currents.However,over the next averaging cycle,surface velocities of30cm s−1were towards the northeast direction whereas subsurface current directions were less than30◦to the right of the surface current(Fig.10e).Cross-shelf current differences of5to10cm s−1occurred during this period when a surface current convergence zone developed within the IOD.In the northern part of the IOD,these directional differences exceeded45◦between surface and subsurface currents(Fig.10f).By contrast,the alignment between the surface and subsurface current vector was less than20◦in the central and southern portions of the IOD.In the last snapshot(Fig.10g),weaker currents were directed towards the west with magnitudes of20cm s−1in the core of the IOD.During this experiment,northward surface winds prevailed as shown in Fig.2.123.3Areal AveragingAs per thesefields(Figs.9,10),the areal sampling of radar signals over the cells is estimated by differencing surface and subsurface currents and spatially averaging them over500m x500m and500m x 750m for the ML1and ML2experiments,respectively.For the areal averaged currents,differences ranged between7to8cm s−1during thefirst snapshot of ML1experiment(Table4)as currents wereflowing in the same direction.A difference of22cm s−1occurred during the second snapshot due to weaker cross-shelf surfaceflows directed towards the coast compared to the more energetic northward subsurfaceflows. By the third snapshot,these differences decreased to16and6cm s−1for the cross-shelf and along-shelf currents,respectively.For the fourth snapshot,averaged differences decreased to7to9cm s−1.In the ML2experiment,similar results were found as suggested by Fig.10.The largest differences(20to26cm s−1)occurred during snapshots3and4when directional differences of up to90◦were observed.During the other snapshots,current differences,based on more than500spatial points,were in the10cm s−1 range.Based on the seven snapshots,averaged current differences(in an arithmetic sense)ranged from8 to14cm s−1in a regime where FC surface velocities as high as1.4m s−1were observed.This equates to an uncertainty of approximately5to20%,which is similar to previous observations.These results suggest a linkage between subsurface ocean structure and radar-derived current signals from80-min(≈1.3hr) sample intervals over the same spatial scales.Given the10-s averaging period of the AUV and the20-min sample interval from the VHF radar,the scatter between the two platforms from the11snapshots is shown in Fig.11.In the cross-shelf direction, the scatter ranged between±30cm s−1(Fig.11a).Based on30samples from the NE-mooring data comparisons during ML2,the regression curve had a slope of1.5,centered approximately in the scatter. The histogram revealed that95%of the current differences were from-30cm s−1to20cm s−1.Similarly, along-shelf currents indicated multiple subsurface values(at≈35cm s−1)during one radar sample(Fig. 11b).Regression analysis of the NE-mooring data revealed a slope of O(1)with a bias of9cm s−1.This slope generally followed with the scatter data,but was skewed towards more positive values.Notice that 95%of the data ranged between-20to25cm s−1,but with more current differences located between -15and10cm s−1.Even though the10-s averaged data were subsampled at30-s intervals,the analysis13。