不同湿度环境对岩石风化产生的影响
Influence of Different Humidity Environments on the Rock Weathering
-
Keywords:
- humidity change /
- rock weathering /
- relative humidity /
- water vapor /
- numerical simulation /
- diffusion
-
类受体激酶(receptor-like kinase,RLK)是植物受体最主要的组成成分,在细胞间的跨膜转导信号中起重要作用,控制植物的繁殖、生长、发育及对各种环境的适应[1]。RLK家族一般包含胞外结构域(extracellular domain,ED)、跨膜区(transmembrane region,TM)和胞内激酶结构域(Ser/Thr kinase domain,STK)[2]。在植物病原菌互作的信号转导通路中,RLK作为模式识别受体 (pattern-recognition receptors,PRRs)可以识别病原菌分子模式 (pathogen-associated molecular patterns,PAMPs),从而在细胞膜上响应植物的天然免疫系统(pattern-triggered immunity,PTI),激活植物的抗病能力[3-5]。当外界信号PAMPs分子被受体类激酶的ED识别后,通过STK将信号传导给下游的信号通路,引发一系列免疫反应,包括ROS爆发、MAPK激活、胼胝质沉淀、钙离子流入、水杨酸积累等[3-5]。如油菜素内酯信号转导通路中,共受体BAK1和受体FLS2共同识别病原菌的鞭毛蛋白,BRI1受体类激酶和BAK1相互作用并经交互磷酸化形成1个活跃的受体复合物,从而激活油菜素内酯抵御病原菌。越来越多的RLKs已被证明是PRRs[3-5]。
细胞壁受体激酶(wall-associated kinases,WAK)及其类基因WAKLs (WAK like genes,WAKLs)属于RLK家族的亚家族[6],合称为WAK/WAKL家族。WAKs和WAKLs的结构非常相似,主要是胞外表皮生长因子结构域(epidermal growth factor domain,EGF)、跨膜结构域和胞内STK,但部分WAKLs缺少一些结构域,如WAKL7缺少跨膜域和激酶结构域[7]。在植物防御反应中,拥有胞外结构域的WAKs/WAKLs基因可以作为细胞外界信号接收员,而胞内激酶部位参与胞质信号级联反应[7],并且在生长发育和对非生物压力的耐受性中起作用。WAKs/WAKLs基因在植物防御病原体方面发挥着关键作用,如:AtWAK1编码WAKs蛋白,与病原菌互作,在拟南芥的防御反应中起重要作用[8];从玉米(Zea mays L.)中克隆的主效数量性状基因ZmWAK可抗玉米丝黑穗病[9];在棉花(Gossypium hirsutum L.)中鉴定到的GhWAK7A参与调节对镰刀菌和黄萎病的反应[10];在山新杨(Populus davidiana×P. bolleana)中过表达克隆到PdPapWAK基因后,其抵抗非生物和生物胁迫的能力进一步增强[11]。
植物WAK/WAKL家族是响应病原菌信号的重要成员。长久以来,基因对基因假说提供了植物抗病模式PTI与ETI系统、植物PRRs与抗病基因(resistance gene,R)的明确界限,但这一界限面临着大量研究的挑战。WAK/WAKL家族作为信号转导途径中的重要受体,可以激活植物的PTI免疫反应模式,检测到病原菌的效应蛋白,可能原因之一是其在植物R基因与病原菌无毒基因互作过程中发挥“监测”作用[12-13]。在小麦对禾谷丝核菌的防御反应中,鉴定到WAK7的转录丰度升高[14],积极响应了小麦禾谷丝核菌的防御反应。小麦基因组中定位到1个WAKL4/stb6基因,该基因识别病原菌效应物Avr Stb6,且绕过过敏性反应赋予小麦对枯叶病的抗性,符合基因对基因假说[15]。当过表达OsWAK91时,水稻激活ETI系统,产生H2O2反应并增强相关R基因的表达[16]。SlWAK1通过番茄中的FLS2/FLS3复合物参与调节PRRs介导的免疫应答[17]。总之,WAK/WAKL家族在PTI和ETI系统中的调控机制很复杂,具体情况尚不清楚。
本研究从课题组前期测定的小麦响应条锈病侵染的转录组数据中筛选出4个在抗感材料中差异表达的WAKs/WAKLs基因,对其进行基因结构、蛋白基本理化性质、保守结构域、系统进化等生物信息学分析,结合实时荧光定量 PCR (RT-qPCR)在抗病和感病材料中分析其响应条锈菌的表达模式,为进一步研究小麦WAKs/WAKLs基因响应条锈菌的侵染作用机制及功能提供参考。
1. 材料与方法
1.1 试验材料
供试小麦品种为抗条锈病品种云337和感条锈病品种云402;小麦条锈菌(Puccinia striiformis f.sp.tritici)为CYR32。将小麦播种在人工培养室中,温度控制在昼16~18 ℃、夜12~15 ℃,每天光照16 h;当小麦生长到二叶期时,将CYR32与滑石粉按体积比1∶20混合后进行撒粉接种;分别采集条锈菌侵染后0、24、48、72、120和168 h的叶片样本,于−80 ℃保存备用。
1.2 转录组测序及差异基因表达分析
转录组测序由北京百迈客生物科技有限公司完成。文库质检合格的测序数据采用DESeq2进行差异基因分析,以差异倍数(fold change,FC)≥1.5且错误发现率(false discovery rate,FDR)<0.01作为差异基因筛选标准,从中挖掘在抗感材料中差异表达的小麦WAKs/WAKLs基因,并对其进行COG分类、GO功能富集、KEGG注释和KEGG通路富集分析。
1.3 小麦WAK/WAKL家族基因鉴定
下载拟南芥、水稻、大麦和小麦基因组数据库文件,采用TBtools v1.108 从中比对出差异表达的小麦WAKs/WAKLs基因的同源基因组和氨基酸序列,并筛选同源性较高的基因;利用NCBI的Batch CD-Search在线软件(https//www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi)对筛选得到的候选基因编码蛋白进行蛋白结构域分析。
1.4 小麦WAK/WAKL家族成员生物信息学分析
利用GSDS在线软件(http//gsds.cbi.pku.edu.cn/)分析和绘制小麦WAKs/WAKLs基因结构;利用ExPASyProtParam Tool (http//web.expasy.org/protparam)预测蛋白质的等电点和分子质量;利用ExPASy (http//web.expasy.org/protscale)预测蛋白质的亲水性和疏水性;利用MEME 在线软件(https://meme-suite.org/meme/tools/meme)预测蛋白基序。为了研究小麦WAKs/WAKLs蛋白与其他物种间的系统进化关系,对小麦、大麦、拟南芥和水稻4个物种中同源性较高的WAKs/WAKLs蛋白全长序列进行比对,采用DNAMAN V6比对蛋白质多序列,采用MEGA 7.0的邻接法构建系统进化树,采用 TBtools 绘制系统进化树与基序预测图。
1.5 小麦WAKs/WAKLs基因的表达模式分析
根据从转录组筛选出的小麦WAKs/WAKLs基因序列,使用Beacon Designer 8设计RT-qPCR引物序列(表1),引物由北京擎科生物股份有限公司合成。采用RNA Easy Fast 植物组织 RNA快速提取试剂盒(天根生化科技有限公司)提取小麦叶片的总RNA,再采用Aid lab公司反转录试剂盒(TURE script 1st Stand cDNA SYNTHESIS Kit)进行反转录。以小麦GAPDH基因为内参,使用2×SYBR® Green预混液(中国DF)试剂盒在Analytikjena-qTOWER2.2型荧光定量PCR仪(德国)上进行RT-qPCR反应。反应体系总体积10.0 μL,包括cDNA 模板1.0 μL,SYBR Green 1 Master 5.0 μL,上、下游引物各 0.5 μL,ddH2O 3.0 μL。反应程序为:95 ℃ 3 min;95 ℃ 10 s,60 ℃ 30 s,39 个循环 ,溶解曲线(60~95 ℃,每个循环增加1 ℃,每次循环时间为4 s)设置2个生物学重复,并使用2-ΔΔCt[18]计算各个样品及各组中的基因相对表达量。采用Excel 2010 处理数据;采用GraphPad Prism 8.4.2 作图。
表 1 供试引物Table 1. Primers for test引物名称
primer names引物序列 (5′→3′)
primer sequencesGAPDH-F GCCTTGCTCCTCTTGCTA GAPDH-R TCCACCTCTCCAGTCCTT TraesCS1B02G004100-F CAATCAGCGGACGGTCTA TraesCS1B02G004100-R ATGTCGCCGATGACTTGT TraesCS6B02G459100-F TTCATCGGAATGTGGTCAAG TraesCS6B02G459100-R GCTGAGTGGAGATAAGAAAGG TraesCSU02G079300-F GAGATGTCGCTGAGAGGTT TraesCSU02G079300-R CCATGCTGTCCAGTGCTA TraesCS3A02G122300-F AGCTCGACACAATGAAGG TraesCS3A02G122300-R CCCAGCCTCAATGTTGAT 2. 结果与分析
2.1 4个WAKs/WAKLs基因的结构
经转录组筛选出4个小麦WAKs/WAKLs基因,分别为TraesCS1B02G004100、TraesCS6B02G459100、TraesCSU02G079300和TraesCS3A02G122300基因(图1)。其中,基因组序列最长的是TraesCS6B02G459100,长度为5 117 kb,包含4个外显子和3个内含子;序列最短的是TraesCS3A02G122300,长度为2 934 kb,包含3个外显子和2个内含子;TraesCSU02G079300和TraesCS1B02G004100的长度相近,分别为3208和3194 bp,均包含3个外显子和2个内含子。
2.2 WAKs/WAKLs蛋白的理化性质和亚细胞定位
由表2可知:TraesCS1B02G004100、TraesCS6B02G459100、TraesCSU02G079300和TraesCS3A02G122300编码的氨基酸序列长度分别为 714、960、737和652 aa;蛋白分子质量分别为78.68、106.29、81.86和69.67 ku;理论等电点分别为8.50、5.98、5.65和6.25;脂肪族氨基酸指数分别为89.59、85.33、81.38和90.49;4个基因的蛋白亲疏水性值均为负值,为亲水性蛋白。亚细胞定位预测显示:4个基因均定位于细胞膜上,为膜蛋白。
表 2 4个WAKs/WAKLs 基因理化性质Table 2. Physicochemical properties of four WAKs/WAKLs genes基因
gene蛋白质序列号
accession number of protein多肽长度/aa
peptide length蛋白分子质量/ku
protein molecular weight等电点
pI氨基酸指数
amino acid
index蛋白亲水性
protein
hydrophilicityTraesCS1B02G004100 XP_044361219.1 714 78.68 8.50 89.59 −0.119 TraesCS6B02G459100 XP_044410226.1 960 106.29 5.98 85.33 −0.120 TraesCSU02G079300 XP_044407301.1 737 81.86 5.65 81.38 −0.214 TraesCS3A02G122300 XP_044339162.1 652 69.67 6.25 90.49 −0.024 2.3 WAKs/WAKLs的系统进化以及基序和功能结构域
小麦、大麦、拟南芥和水稻中同源性较高的43个WAKs/WAKLs蛋白质序列按照亲缘关系被分为4个分支,TraesCS6B02G459100和来自水稻、大麦、拟南芥的成员聚为一类,处于分支Ⅰ,且与水稻Os05g04460亲缘关系最近;TraesCSU02G079300和TraesCS3A02G122300处于分支Ⅱ,分别与水稻Os02g02120、拟南芥AtWAKL15亲缘关系最近;TraesCS1B02G004100处于分支Ⅲ,与大麦HORVU.MOREX.r3.5HG0503880亲缘关系最近;分支Ⅳ仅包括小麦成员(图2a)。
图 2 小麦与拟南芥、水稻、大麦 WAK/WAKLs基因中保守蛋白质的系统进化关系(a)、蛋白结构域(b)和基序结构(c)注:红色框为从转录组中筛选的4个WAKs/WAKLs基因;下同。Figure 2. Phylogenetic relationship (a), protein domain(b) and motifs structure (c) based on the amino acid sequences of WAKs/WAKLs from Triticum aestivum and Arabidopsis thaliana, Oryza sativa, and Hordeum vulgareNote: In the red box are four WAKs/WAKLs genes screened from the transcriptome; the same as below.蛋白结构域预测结果显示:小麦、大麦、拟南芥和水稻WAKs/WAKLs蛋白序列均包括激酶部位(图2b)。在进化分支Ⅰ中,小麦WAKs/WAKLs蛋白都包括2个GUB_WAK_bind结构域,1个EGF和1个STK,除水稻Os01g26300和拟南芥AtWAKL16编码的蛋白质不具有胞外结构域外,其他成员均具有胞外结构域。在进化分支Ⅱ中,大部分水稻成员不具备胞外结构域,而其他成员均具备;部分小麦成员还具有跨膜TM_FGFR-like结构域。在进化分支Ⅲ中,除小麦TraesCS1B02G004100和大麦HORVU.MOREX.r2.1HG0004690、HORVU.MOREX.r3.5HG0503880编码的蛋白质具有完整的胞内和胞外结构域外,其他成员均不具备。在进化分支Ⅳ中,只有TraesCS3A02G049500编码的蛋白质具有胞外和胞内结构域。说明胞内STK为所有成员共有,预测WAKs/WAKLs家族基因发挥功能与激酶部位有较大关联。
基序预测结果显示:小麦、大麦、拟南芥和水稻的WAKs/WAKLs蛋白序列共包括12种基序(图2c)。进化分支Ⅰ包括11种基序,TraesCS6B02G459100与大麦HORVU.MOREX.r3.3HG0499530编码的蛋白质基序相似;进化分支Ⅱ包括10种基序,TraesCSU02G079300与水稻Os02g02120、大麦HORVU.MOREX.r3.7HG0649430编码的蛋白基序高度相似,TraesCS3A02G122300与拟南芥AtWAKL15编码的蛋白基序高度相似;进化分支Ⅲ包括10种基序,TraesCS1B02G004100与大麦HORVU.MOREX.r3.5HG0503880编码的蛋白基序高度相似;进化分支Ⅳ包括9种基序,均是来自小麦WAKs/WAKLs家族的成员,且编码的蛋白基序基本相似。此外,基序1、2、3、4、5、7、8、11为所有成员共有,均包括STK,预测胞内激酶部位保守性较高,是决定蛋白质功能的关键基序;基序9、10、12包括胞外EGF和半乳糖醛酸结合结构域(galacturonic acid binding domain,GUB),说明胞外结构域多样,但具有一定的保守性。
2.4 WAKs/WAKLs蛋白激酶部位的同源性比对
为了进一步研究激酶部位的保守性,预测其发挥功能的部位,对涉及的成员进行激酶部位蛋白序列同源性比对分析。该部位通常有2个与激酶发挥功能有关的重要区域,包括蛋白激酶ATP结合位点、丝氨酸/苏氨酸蛋白激酶活性位点和酪氨酸蛋白激酶活性位点(图3),分析结果显示:胞内结构域同源性达73.34%,蛋白激酶ATP结合位点、丝氨酸/苏氨酸蛋白激酶活性位点和酪氨酸蛋白激酶活性位点保守性强,均符合这2个结合位点序列构成特征。根据蛋白激酶ATP结合位点中不变的天冬氨酸(D)附近是否存在保守精氨酸(R),将激酶分为RD型激酶和非RD型激酶,由图4可知:TraesCS6B02G459100和TraesCS3A02G122300为RD型激酶,TraesCS1B02G004100和TraesCSU02G079300为非RD型激酶。
2.5 小麦 WAKs/WAKLs基因响应条锈病的表达
由图5可知:4个基因在小麦接种条锈菌亲和互作和非亲和互作过程中差异表达。侵染0 h,TraesCS1B02G004100和TraesCS6B02G459100在小麦与条锈菌非亲和互作中被短暂诱导表达后下调,随后在侵染早期(24~48 h)上调表达,而TraesCS3A02G122300则表达相反;TraesCSU02G079300在早期表达量无显著差异。接种48 h后,TraesCS1B02G004100和TraesCSU02G079300表达受到抑制而下调直至接种168 h;TraesCS6B02G459100表达则呈上调直至接种120 h。综上所述,亲和互作过程中,TraesCS1B02G004100和TraesCS6B02G459100几乎不表达,TraesCS3A02G122300早期表达趋势和非亲和互作一致,TraesCSU02G079300在早期表达量低于非亲和互作。
图 5 4个WAKs/WAKLs基因在抗病品种云337 (Y337)和感病品种云402 (Y402)接种CYR32后的转录组表达量注:“**”表示与上一个时间相比,在0.01水平上差异极显著;下同。Figure 5. Transcriptome expression level of four WAKs/WAKLs genes in resistant variety Yun337 (Y337) and susceptible variety Yun402 (Y402) post inoculation with CYR32Note: “**” indicate extremely significant differences at the 0.01 level compared to the previous time; the same as below.由图6可知:非亲和互作过中,在条锈菌侵染小麦的24~48 h,TraesCS1B02G004100上调表达,TraesCS6B02G459100显著上调表达,TraesCS3A02G122300和TraesCSU02G079300表达受到抑制而下调;在接种条锈菌48 h后,TraesCS1B02G004100和TraesCS6B02G459100先下调后上调表达;TraesCS6B02G459100、TraesCS3A02G122300和TraesCSU02G079300的表达量分别在接种条锈菌120、24和168 h后达到峰值。亲和互作过程中,TraesCS1B02G004100几乎不表达;转录组中几乎未表达的TraesCS6B02G459100,在接种0~24 h被抑制表达,24 h后上调表达,168 h表达量达到峰值,但由于该基因转录丰度较低,可能导致RT-qPCR检测结果不能有效地反映其实际表达情况;TraesCS3A02G122300和TraesCSU02G079300在接种早期下调表达,接种120 h后表达上调。综上所述,4个基因在小麦和条锈菌互作过程中利用转录组和RT-qPCR分析的表达模式基本一致。
3. 讨论
WAKs/WAKLs基因在植物中属于大家族,目前从拟南芥中鉴定了26个WAKs/WAKLs基因,水稻中125个,大麦中91个,小麦中320个[19-21],作为病原菌的模式识别受体,它在植物信号转导过程中发挥着重要作用。本课题组从转录组中筛选的4个WAKs/WAKLs基因,其结构域存在一定差异,暗示这4个基因功能上存在差异;通过RT-qPCR验证,阐明TraesCS1B02G004100基因参与调控小麦对条锈病的抗性,研究结果为揭示WAKs/WAKLs基因调控抗病性的功能和机制提供了参考。
本研究的4个WAKs/WAKLs基因具有典型的胞外富含半胱氨酸的GUB_WAK_bind和胞内丝氨酸/苏氨酸激酶结构域。从进化分支上看,小麦TraesCS6B02G459100与水稻OsWAK91的亲缘关系最近,OsWAK91正向调控水稻抵御稻瘟病菌侵染[16]。已有研究报道:2个胞外结构域变异的OsWAK10基因以不同的亲和力结合果胶寡糖,以控制水稻次生壁的合成[22],由此可以说明:GUB作为首先感知到外界信号的受体,可以与果胶结合激活WAKs/WAKLs基因,以调节植物细胞形态。TraesCS6B02G459100与上述2个水稻WAKs基因同为RD激酶,由此可以推测,该基因可能编码与调节植物生长发育、抗病相关的激酶,相较于其他3个基因,其具有的2个GUB结构域在信号转导过程中的作用机制值得探索。TraesCS3A02G122300与未知功能的激酶聚为同一分支,在抗感材料中都被抑制表达,但在后期感病材料的表达量高于抗病材料,尽管属于RD激酶,但已有报道中也出现能够调控植物抗病反应的RD激酶,如OsWAK91、OsWAK92、OsWAK94和OsWAK14基因是水稻对稻瘟病菌数量抗性的调节因子[16],因此TraesCS3A02G122300的具体功能机制值得进一步研究。
WAK/WAKL家族胞内的ATP结合位点及丝氨酸/苏氨酸结合位点是调节植物先天免疫的重要结构[23],尤其是丝氨酸/苏氨酸蛋白激酶活性位点中保守的天冬氨酸附近是否存在精氨酸很大程度上决定了激酶的功能。已知在先天免疫中作为PRRs的15个激酶中有12个是非RD型[24],表明所有非RD 激酶可能都参与先天免疫。TraesCS1B02G004100和TraesCSU02G079300为非RD型激酶。其中,TraesCS1B02G004100在亲和互作体系中几乎不表达;而在非亲和体系中,该基因早期被诱导,24 h以后被抑制,这与已报道的OsWAK112d基因表达模式一致,且过表达OsWAK112d基因及其突变体分别增加了水稻的易感性和抗性,表明其在植物感染期间是抑制基础抗性的负调控因子[16],故推测TraesCS1B02G004100可能是小麦响应条锈菌侵染的负调控基因,并且早期可能被PAMPs诱导。TraesCSU02G079300与OsWAK11处于同一进化分支,OsWAK11通过GUB结合果胶并监测果胶甲基化的变化以解毒过量的铜,并且可以调节细胞的伸长率以适应外界环境的变化[25],其转录丰度的趋势与TraesCS1B02G004100相似,并且在2种互作体系中差异表达,该基因可能也在防御反应中发挥作用。
本研究筛选的小麦WAK1基因TraesCS1B02G004100和已经在小麦中鉴定到的Snn1基因同源性最高,达99.86%。在腐生病原菌(Parastagonospora nodorum)和小麦互作过程中,SnTox1作为病原菌分泌的效应蛋白与Snn1基因相互作用,引起小麦细胞程序性死亡并赋予其对叶斑病的易感性,Snn1-SnTox1是小麦和叶斑病互作的多个重要模式之一[26-27]。本研究发现:在小麦与条锈菌非亲和互作过程中,TraesCS1B02G004100基因响应条锈菌侵染且表达受到抑制,而在亲和互作过程中不受条锈菌诱导表达,说明该基因极大可能参与调控小麦对条锈病的抗性。
4. 结论
4个小麦WAKs/WAKLs基因中,TraesCS3A02G122300和TraesCSU02G079300在小麦与条锈菌亲和及非亲和互作体系中差异表达;TraesCS6B02G459100和TraesCS1B02G004100在小麦非亲和互作中特异表达。非RD型激酶TraesCS1B02G004100在表达过程中受到抑制,作为病原菌效应蛋白的受体,其在参与调控小麦对条锈病的抗性功能值得进一步研究;其他3个WAKs/WAKLs基因在植物生长发育过程中的作用也值得探索。
-
图 6 玄武岩风化速率(v)—相对湿度(RH)与混凝土水分扩散系数(D)—内部相对湿度关系对比图[24]
Figure 6. Comparison between basalt weathering rate (v)-relative humidity (RH) and concrete moisture diffusion coefficient (D)-internal relative humidity
表 1 材料模型参数
Table 1 Parameters of material model
材料名称
material name孔隙率/%
porosity密度/(kg·m–3)
density压强/Pa
pressure扩散系数/(m·s)
diffusion coefficient玄武岩 basalt 0.24 2.943 4×103 101 325 8.6×10–15 表 2 不同时间段岩石的风化速率
Table 2 The weathering rate of rock in different time periods
时间/a
time位置/m
position位移/m
displacement速率/(μm·a–1)
rate0~1 000 0.015 2 0.015 2 15.2 1 000~2 000 0.023 0 0.007 8 7.8 2 000~3 000 0.029 5 0.006 5 6.5 3 000~4 000 0.035 7 0.006 2 6.2 4 000~5 000 0.043 7 0.008 0 8.0 5 000~5 213 0.050 0 0.006 3 29.6 -
[1] 孙家齐. 工程地质[M]. 武汉: 武汉工业大学出版社, 2000. [2] 刘作程. 岩石学[M]. 北京: 冶金工业出版社, 1992. [3] MOON V, JAYAWARDANE J. Geochemical changes during early stages of weathering of Karamu Basalt, New Zealand[J]. Engineering Geology, 2004, 74(1): 57. DOI: 10.1016/j.enggeo.2004.02.002.
[4] 尚彦军, 吴宏伟, 曲永新. 花岗岩风化程度的化学指标及微观特征对比——以香港九龙地区为例[J]. 地质科学, 2001, 36(3): 279. DOI: 10.3321/j.issn:0563-5020.2001.03.003. [5] 巫锡勇, 罗健, 魏有仪. 岩石风化与岩石化学成分的变化研究[J]. 地质与勘探, 2004, 40(4): 85. DOI: 10.3969/j.issn.0495-5331.2004.04.019. [6] 徐则民, 黄润秋, 唐正光. 硅酸盐矿物溶解动力学及其对滑坡研究的意义[J]. 岩石力学与工程学报, 2005, 24(9): 1479. DOI: 10.3321/j.issn:1000-6915.2005.09.004. [7] 江世贵. 结晶学及矿物学[M]. 北京: 冶金工业出版社, 1993. [8] 汤连生, 王思敬. 水—岩化学作用对岩体变形破坏力学效应研究进展[J]. 地球科学进展, 1999, 14(5): 433. . DOI: 10.3321/j.issn:1001-8166.1999.05.002. [9] 袁和忠, 王运生, 房冬恒, 等. 西南某水电站深风化现象及其机理浅析[J]. 水土保持研究, 2006, 13(1): 103. DOI: 10.3969/j.issn.1005-3409.2006.01.034. [10] 徐则民, 唐正光. 石灰岩腐岩的基本特征及其形成机制[J]. 地质论评, 2012, 53(3): 421. DOI: 10.3321/j.issn:0371-5736.2007.03.017. [11] MATSUKURA Y, HASHIZUME K, OGUCHI C T. Effect of microstructure and weathering on the strength anisotropy of porous rhyolite[J]. Engineering Geology, 2002, 63(1): 39. DOI: 10.1016/S0013-7952(01)00067-9.
[12] YANG J Q, XU Z M, ZHANG R, et al. Formation and evolution of Emeishan basalt saprolite in vadose zones of Touzhai landslide source rockmass[J]. Journal of Mountain Science, 2017, 14(6): 1174. DOI: 10.1007/s11629-016-4169-y.
[13] BUSENBERG E, CLEMENCY C V. The dissolution kinetics of feldspars at 25 °C and 1 atm CO2 partial pressure[J]. Geochimica et Cosmochimica Acta, 1976, 40(1): 41. DOI: 10.1016/0016-7037(76)90192-7.
[14] TWIDALE C R, BOURNE J A. International science "down under": the British association meeting in Australia, August 1914, with special reference to related activities in Adelaide[J]. Earth Sciences History: Journal of the History of the Earth Sciences Society, 2002, 21(21): 166.
[15] FRIEDMAN I, LONG W. Hydration rate of obsidian[J]. Science, 1976, 191(4225): 347. DOI: 10.1126/science.191.4225.347.
[16] 王承, 陶瑛. 硅酸盐玻璃的风化[J]. 硅酸盐学报, 2003, 31(1): 78. DOI: 10.3321/j.issn:0454-5648.2003.01.017. [17] CARMELIET J, DESCAMPS F, HOUVENAGHEL G. A multiscale network model for simulating moisture transfer properties of porous media[J]. Transport in Porous Media, 1999, 35(1): 67. DOI: 10.1023/A:1006500716417.
[18] 梅宁, 尹凤, 陆虹涛. 湿度变化对气体污染物扩散影响的研究[J]. 中国海洋大学学报(自然科学版), 2006, 36(6): 987. DOI: 10.3969/j.issn.1672-5174.2006.06.027. [19] 张君, 侯东伟. 基于内部湿度试验的早龄期混凝土水分扩散系数求解[J]. 清华大学学报(自然科学版), 2008(12): 2033. [20] 唐春安, 唐世斌. 岩体中的湿度扩散与流变效应分析[J]. 采矿与安全工程学报, 2010, 27(3): 292. DOI: 10.3969/j.issn.1673-3363.2010.03.002. [21] 唐世斌, 唐春安, 林皋, 等. 水泥基复合材料湿度扩散特性细观数值模拟研究[J]. 大连理工大学学报, 2011(6): 868. DOI: doi:10.7511/dllgxb201106015. [22] 杨佳玲. 湿度梯度作用下天然气泄漏扩散浓度场研究[D]. 重庆: 重庆交通大学, 2014. [23] AUGUSTITHIS S S, VGENOPOULOS A. On the hawleyite-sphalerite-wurtzite-galena paragenesis from Ragada, Komotini, (Rhodope) North Greece[M]// AMSTUTZ G C, FRENZEL G, KLUTH C, et al. Ore genesis: special publication of the society for geology applied to mineral deposits, Vol 2. Berlin Heidelberg: Springer-Verlag, 1982.
[24] OZAKA Y, FUJIWARA T, AKITA H. A practical procedure for the analysis of moisture transfer within concrete due to drying[J]. Magazine of Concrete Research, 1997, 49(179): 129. DOI: 10.1680/macr.1997.49.179.129.