期刊文献+
共找到17篇文章
< 1 >
每页显示 20 50 100
卫星遥感海表温度资料和高度计资料的变分同化 被引量:16
1
作者 肖贤俊 何娜 +2 位作者 张祖强 刘怀明 王东晓 《热带海洋学报》 CAS CSCD 北大核心 2011年第3期1-8,共8页
利用国家气候中心正在发展的第二代全球海洋资料同化系统(BCC_GODAS2.0),针对多变量同化的协调性问题,发展了一种基于三维变分框架(3DVAR)下的高度计和海表温度(SST)相互约束的同化方法。该方法使海面高度和SST资料在同一个动力约束关... 利用国家气候中心正在发展的第二代全球海洋资料同化系统(BCC_GODAS2.0),针对多变量同化的协调性问题,发展了一种基于三维变分框架(3DVAR)下的高度计和海表温度(SST)相互约束的同化方法。该方法使海面高度和SST资料在同一个动力约束关系下进行同化。在一般方法中,海面高度和SST观测项是代价函数中2个独立的观测项,海面高度项引入动力高度计算公式,海表温度项用统计关系进行垂向投影。在代价函数的实际求解的计算过程中,虽然其总体积分效应受海面高度观测的约束,但整个水柱中各层温盐分析变量的调整是无序的。针对这个问题,文章提出一种新的同化方案。该方案将SST的观测项并入海面高度观测项中,海面高度的一部分,确切说是上层海洋部分,由SST决定,因此至少在SST的统计关系能影响到深度的上层海洋,在代价函数的求解过程中,温盐的调整是受较强的统计关系约束的,而这种统计关系的有效性已经在很多SST的同化试验中被其他学者广泛应用并证明。利用该方法,对1993—1997年的AVHRR卫星遥感海表温度资料进行变分同化试验,用TAO、OISST和SODA数据集进行检验证明,通过对卫星遥感资料的同化能够有效改进对海洋温度和盐度的估计。海洋表层的月平均温、盐度的总均方根误差相对同化前分别降低了0.67℃和0.2‰。在混合层中,同化效果较好。 展开更多
关键词 海洋资料同化 三维变分 海表温度 高度计 遥感资料 垂向投影
下载PDF
2002年3月19日沙尘暴爆发条件分析 被引量:17
2
作者 肖贤俊 刘还珠 +1 位作者 宋振鑫 何金海 《应用气象学报》 CSCD 北大核心 2004年第1期1-9,共9页
使用沙尘暴预报模式的输出资料 ,对 2 0 0 2年 3月 1 9~ 2 1日的一次特强沙尘暴过程作动力诊断分析。结果表明 ,该次沙尘暴由蒙古气旋后部冷锋锋生产生的偏西北大风引发。近地面风速的垂直切变和地面热通量的加大 ,都可使边界层湍流加... 使用沙尘暴预报模式的输出资料 ,对 2 0 0 2年 3月 1 9~ 2 1日的一次特强沙尘暴过程作动力诊断分析。结果表明 ,该次沙尘暴由蒙古气旋后部冷锋锋生产生的偏西北大风引发。近地面风速的垂直切变和地面热通量的加大 ,都可使边界层湍流加强扬起地面沙尘。地面锋区附近风场的强水平切变 ,锋面垂直环流及锋后斜压转换的作用 ,将地面卷起的沙尘带到高空 。 展开更多
关键词 沙尘暴模式 边界层湍流 斜压转换 冷锋 天气过程分析
下载PDF
南海三维变分海洋同化模式及其验证 被引量:3
3
作者 肖贤俊 王东晓 +1 位作者 闫长香 朱江 《自然科学进展》 北大核心 2007年第3期353-361,共9页
基于Princeton Ocena Model(POM)建立了一个南海三维变分海洋同化模式,该模式主要用于卫星遥感高度计沿轨资料同化,考虑了每条轨道上的误差相关性,并且引入递归滤波优化运算.采用1998年春季和夏季、2000年夏季3个航次温盐资料对这个系... 基于Princeton Ocena Model(POM)建立了一个南海三维变分海洋同化模式,该模式主要用于卫星遥感高度计沿轨资料同化,考虑了每条轨道上的误差相关性,并且引入递归滤波优化运算.采用1998年春季和夏季、2000年夏季3个航次温盐资料对这个系统进行评估表明,相对于未同化的模式模拟结果,无论对温盐异常年份还是正常年份,该系统的温盐均方根误差和偏差都显著减小,误差结构更符合正态分布.特别是在盐度场上能重现南海由于降水和径流所造成的淡水盖,该系统也能将海表动力高度信息有效传递进模式次表层;温盐资料验证各航次评估结果差别不大,表明该同化系统性能是稳定的. 展开更多
关键词 南海 三维变分 海洋同化 递归滤波 卫星高度计
下载PDF
GTS的温盐资料在BCC_GODAS中的同化结果分析 被引量:4
4
作者 刘向文 李维京 +1 位作者 吴统文 肖贤俊 《应用气象学报》 CSCD 北大核心 2010年第5期558-569,共12页
分析了从GTS(全球无线通讯系统)获得的2002-2007年海洋温盐观测资料在国家气候中心第2代全球海洋资料同化系统(BCC_GODAS 2.0)中的同化结果。与SODA(简易海洋同化数据)资料的比较表明:GTS中的海洋温盐资料同化对模式温盐场的改进之处主... 分析了从GTS(全球无线通讯系统)获得的2002-2007年海洋温盐观测资料在国家气候中心第2代全球海洋资料同化系统(BCC_GODAS 2.0)中的同化结果。与SODA(简易海洋同化数据)资料的比较表明:GTS中的海洋温盐资料同化对模式温盐场的改进之处主要表现在混合层暖区的范围和中心强度、温跃层中温度槽脊的深度、温跃层附近的温度梯度以及盐度高、低值区的范围和中心强度等方面,同化后全球温盐场的均方根误差得到一定程度的降低。挑选位于不同海区的单点温盐廓线与ARGO(地转海洋学实时观测阵)观测作了进一步比较,结果表明:大多数情况下,同化后温盐廓线的均方根误差得到明显降低,模拟的温盐垂向分布特征也更为准确。与TAO(热带大气海洋观测网)资料的比较也同样表明:同化后的温盐场特征会得到一定程度改善。 展开更多
关键词 GTS BCC_GODAS 温度 盐度 剖面 廓线
下载PDF
基于三层嵌套网格的珠江口冬季盐度层化的数值模拟 被引量:4
5
作者 周巍 王东晓 +2 位作者 俎婷婷 李海洋 肖贤俊 《海洋与湖沼》 CAS CSCD 北大核心 2013年第3期545-556,共12页
运用三层嵌套网格的ROMS模式较好地模拟了2009年冬季珠江口的主要水动力过程和盐度分布。结果表明,珠江口盐度层化具有明显的潮周期变化特征,涨潮时表底层盐度差较小,层化较弱;落潮时层化较强。利用势能异常变化平衡方程分析影响层化的... 运用三层嵌套网格的ROMS模式较好地模拟了2009年冬季珠江口的主要水动力过程和盐度分布。结果表明,珠江口盐度层化具有明显的潮周期变化特征,涨潮时表底层盐度差较小,层化较弱;落潮时层化较强。利用势能异常变化平衡方程分析影响层化的贡献项,结果表明势能异常的平流项和应变项是影响珠江口势能异常变化的主要因素。 展开更多
关键词 ROMS模式 盐度层化 势能异常 珠江口
下载PDF
海洋环流模式中卫星遥感资料同化的应用进展 被引量:3
6
作者 何娜 张祖强 +1 位作者 肖贤俊 郭品文 《海洋科学》 CAS CSCD 北大核心 2008年第12期92-96,共5页
关键词 海洋环流模式 卫星遥感技术 资料同化技术 卫星观测数据 数值模拟分析 观测资料 观测技术
下载PDF
从GTS获得的海洋温、盐资料在BCC海洋同化系统中的质量控制及同化结果初步分析 被引量:2
7
作者 刘向文 李维京 +1 位作者 吴统文 肖贤俊 《气象学报》 CAS CSCD 北大核心 2011年第4期672-681,共10页
介绍了从全球电信系统(GTS)上获得的海洋温度、盐度观测资料在中国国家气候中心(BCC)新一代海洋同化系统中的应用情况。通过资料的质量控制判断温、盐观测的重复记录、观测深度、地形、极端值、气候变率、层结、空间差异,有效地过滤了... 介绍了从全球电信系统(GTS)上获得的海洋温度、盐度观测资料在中国国家气候中心(BCC)新一代海洋同化系统中的应用情况。通过资料的质量控制判断温、盐观测的重复记录、观测深度、地形、极端值、气候变率、层结、空间差异,有效地过滤了错误的或不可靠的观测信息。质量控制后,将温、盐观测资料加入同化系统,有效地改进了模块化全球海洋环流模式MOM4中的全球热带、副热带海洋,尤其是太平洋地区的多年平均海表温度、盐度场分布特征;此外,同化温、盐资料对南北半球中纬度地区的海表温度分布特征也有明显的改进。对比同化前后的均方根误差(RMSE)发现,同化后大部分海区,尤其是热带海洋的海表温度/盐度的均方根误差明显降低,降幅通常在0.1—1.0℃/psu,模拟与观测的海表温、盐分布特征也更为接近。进一步分析指出,同化明显地改善了模式对Nino3、Nino4区海温时间演变特征的模拟,同化后的Nino3海温与最优插值海表温度的差异减小,但其通常在上半年改进较多(差值绝对值多在0.5℃左右),而在下半年则改进较少(差值绝对值常达1℃左右);Nino4区的海温特征则改进明显,其与最优插值海表温度的差值绝对值通常都控制在0.5℃以下。 展开更多
关键词 BCC_GODAS 质量控制 海表温度 海表盐度
下载PDF
南海北部集合卡曼滤波同化SST试验 被引量:1
8
作者 舒业强 隋丹丹 +1 位作者 王伟文 肖贤俊 《热带海洋学报》 CAS CSCD 北大核心 2010年第5期10-16,共7页
基于POM(Princeton Ocean Model)建立一个南海北部集合卡曼滤波的同化模式,主要用于卫星海表面温度的同化。模式的平均水平分辨率为5km,垂向分层为20层;侧边界条件嵌套到一个大范围的南海海洋模式,在同化方案上采用一个均方根集合卡曼... 基于POM(Princeton Ocean Model)建立一个南海北部集合卡曼滤波的同化模式,主要用于卫星海表面温度的同化。模式的平均水平分辨率为5km,垂向分层为20层;侧边界条件嵌套到一个大范围的南海海洋模式,在同化方案上采用一个均方根集合卡曼滤波算法,避免观测的扰动;适当引入局地化算子,消除样本在空间上的虚假相关,同时增加集合样本的自由度。该同化试验同化了2008年夏季6月到7月的GHRSST(Global High-Resolution Sea Surface Temperature),然后采用2008年夏季SCOPE(Northern South China Sea Coastal Oceanographic Process Ex-periment)航次的温、盐数据对同化结果进行评估。结果表明,相对于未同化模式模拟结果,同化模式温度的改善比较明显,表现在加强了南海北部的上升流,校正了海表温度的偏差,改善了温度的垂向分布。由于集合卡曼滤波是一种多变量调整的同化方法,同化SST不仅能改善表层与次表层的温度分布,而且对流场和盐度的调整也比较明显。 展开更多
关键词 集合卡曼滤波 同化 南海北部 SST
下载PDF
背景误差协方差矩阵不同求逆方案在高度计资料同化试验中的应用比较
9
作者 刘向文 李维京 +1 位作者 吴统文 肖贤俊 《海洋学报》 CAS CSCD 北大核心 2012年第5期65-73,共9页
针对海表高度计资料的同化,考查了背景误差协方差矩阵的不同求逆方案对同化效果的影响。所使用的求逆方案包括避免求逆的经验正交函数方案(EOF/EOF_var)、递归滤波方案(RF/RF_var)以及采用初等变换法直接求逆的方案(Inv)。基于上述方案... 针对海表高度计资料的同化,考查了背景误差协方差矩阵的不同求逆方案对同化效果的影响。所使用的求逆方案包括避免求逆的经验正交函数方案(EOF/EOF_var)、递归滤波方案(RF/RF_var)以及采用初等变换法直接求逆的方案(Inv)。基于上述方案开展了热带太平洋地区2002年1-7月的TOPEX/Poseidon高度计资料同化试验,并利用SODA再分析资料和TAO观测资料评估了各方案对温度场的同化效果,主要得到如下结果:与SODA相比较,Inv方案对模式温度场改进甚微,其余四种方案在100~300m深度之间对温度场改进较多,在其它深度范围内则改进较少;与TAO观测相比较,EOF_var、RF_var方案对模式温度场改进最多,EOF和RF方案次之,Inv方案则对温度场改进甚少。 展开更多
关键词 海表高度 经验正交函数 递归滤波 初等变换法直接求逆 温度 均方根误差
下载PDF
GTSPP资料在BCC海洋资料同化系统中的同化效果分析
10
作者 武丰民 张祖强 +2 位作者 肖贤俊 何金海 刘怀明 《大气科学学报》 CSCD 北大核心 2014年第1期38-46,共9页
介绍了国家气候中心第二代海洋资料同化系统(BCC_GODAS2.0)所用的模式和同化方案,并设计了两个试验对同化系统的效果和稳定性进行检验:一个是加入GTSPP资料的1990—2008年19 a的同化试验,另一个是同时段内不加任何资料的模拟试验。利用O... 介绍了国家气候中心第二代海洋资料同化系统(BCC_GODAS2.0)所用的模式和同化方案,并设计了两个试验对同化系统的效果和稳定性进行检验:一个是加入GTSPP资料的1990—2008年19 a的同化试验,另一个是同时段内不加任何资料的模拟试验。利用OISST和SODA(Simple Ocean Data Assimilation)数据集对试验结果进行评估。评估结果表明,加入GTSPP资料的同化结果可以有效改进原模式对温度和盐度的模拟。在表层,整个海洋的温度和盐度场误差都有所减小,其中热带太平洋的改进效果最为明显;此外,同化结果能更好地反映Nino3和Nino4区海表温度的时间演变。垂直层次上的修正效果也很明显,在200 m左右的混合层效果改进最好,温度的均方根误差可减小1.5℃,盐度的均方根误差可减小0.6 psu。 展开更多
关键词 BCC—GODAS2 0 资料同化 海表温度 海表盐度 评估
下载PDF
Argo温盐资料在BCC_GODAS2.0中的同化试验 被引量:1
11
作者 刘怀明 张祖强 +1 位作者 肖贤俊 何金海 《安徽农业科学》 CAS 北大核心 2010年第26期14521-14526,共6页
基于国家气候中心新一代的海洋资料同化系统(BCCGODAS2.0),对20042008年Argo温盐资料进行了同化试验,并利用TAO(International Tropical Ocean and Global Atmosphere/Tropical Atmosphere Ocean)、OISST(OptimallyZnterPolated Se... 基于国家气候中心新一代的海洋资料同化系统(BCCGODAS2.0),对20042008年Argo温盐资料进行了同化试验,并利用TAO(International Tropical Ocean and Global Atmosphere/Tropical Atmosphere Ocean)、OISST(OptimallyZnterPolated Sea Surface Temperature)和SODA(Simple Ocean Data Assimilation)数据集进行了检验和评估。同化结果对01000m的整层海洋的模拟都有改善,尤其在混合层更为明显,能够使温度的均方根误差(RMSE,Root Mean Squared Error)减小3℃,盐度的均方根误差减小0.2psu。尤其是减小了北半球中高纬地区、北大西洋温度模拟的显著误差。模拟结果中在中西太平洋、印度洋地区明显的暖异常也在同化结果中得到了校正。通过TAO时间序列资料的验证表明,同化后能很好地模拟出热带附近上层海洋的季节变化震荡幅度。此外,同化之后,模拟的效果能得到较为稳定的改善,证明系统是稳定的,这对在业务上的应用打好了良好的基础。 展开更多
关键词 资料同化 BCCGODAS2.0 ARGO资料
下载PDF
短中期降水温度天气过程区域分布的研究 被引量:5
12
作者 刘还珠 肖贤俊 《干旱气象》 2003年第3期14-20,共7页
运用旋转主分量方法对分布于全国621个站点的候降水量和2214个站的逐日最高最低气温进行分析,得到4个季节的降水气温主特征模态及其相对应的时间变率。分析结果表明,该方法所分解得到的特征模态较好地反映了全国不同区域降水温度演变的... 运用旋转主分量方法对分布于全国621个站点的候降水量和2214个站的逐日最高最低气温进行分析,得到4个季节的降水气温主特征模态及其相对应的时间变率。分析结果表明,该方法所分解得到的特征模态较好地反映了全国不同区域降水温度演变的差别,依据各模态的相关系数将全国划分为不同的降水温度天气区。最后,获得表征全国不同天气特点的252个代表站,这些站的资料可作为提高要素客观预报、扩展服务领域的基本资料。 展开更多
关键词 候降水量 旋转主分量 特征模态 天气区 代表站 降水天气过程 预报
下载PDF
基于一维变分算法的红外高光谱(IASI)卫星遥感大气温湿廓线研究 被引量:8
13
作者 官元红 任杰 +3 位作者 鲍艳松 陆其峰 刘辉 肖贤俊 《大气科学学报》 CSCD 北大核心 2019年第4期602-611,共10页
大气温湿度廓线是大气重要参数,在数值天气预报及天气预警中具有重要的应用价值。为获得高精度的大气温度与水汽混合比廓线数据,研究了基于Metop-A/IASI红外高光谱资料的大气温度与水汽混合比廓线变分反演方法。利用IASI高光谱传感器温... 大气温湿度廓线是大气重要参数,在数值天气预报及天气预警中具有重要的应用价值。为获得高精度的大气温度与水汽混合比廓线数据,研究了基于Metop-A/IASI红外高光谱资料的大气温度与水汽混合比廓线变分反演方法。利用IASI高光谱传感器温度和水汽探测通道资料,结合CRTM模式和WRF模式预报技术,使用一维变分方法,研究了卫星资料质量控制、背景误差协方差本地化、观测误差协方差计算等方法,构建了大气温度及水汽混合比廓线变分反演系统,并在北京、青岛、沈阳3个地区开展了反演试验。以探空为标准的反演结果对比显示,使用WRF模式预报值作为背景场,温度的平均误差绝对值小于0.6K,均方根误差为0.89K;水汽混合比的平均误差绝对值小于0.021g/kg,均方根误差为0.02g/kg。试验结果表明:基于一维变分方法,可以利用Metop-A/IASI红外高光谱资料进行大气温度与水汽混合比廓线高精度探测。 展开更多
关键词 Metop-A/IASI 温湿廓线 一维变分 反演
下载PDF
参加2007年AOGS会议总结
14
作者 张秀芝 朱艳峰 肖贤俊 《气象科技合作动态》 2007年第6期15-16,共2页
亚洲一大洋洲地球科学学会(AOGS)是由亚洲、大洋洲以及部分北美和欧洲的地球科学科技工作者组成的国际学术组织,会员超过1000人。其宗旨是为从事地球科学研究工作的组织和个人提供合作和交流的平台。亚洲一大洋洲地球科学学会第4届... 亚洲一大洋洲地球科学学会(AOGS)是由亚洲、大洋洲以及部分北美和欧洲的地球科学科技工作者组成的国际学术组织,会员超过1000人。其宗旨是为从事地球科学研究工作的组织和个人提供合作和交流的平台。亚洲一大洋洲地球科学学会第4届年会(AOGS 2007 Asia Oceania 展开更多
关键词 会议总结 地球科学 科学研究工作 科学学会 学术组织 科技工作者 大洋洲 ASIA
下载PDF
2000年夏季南海西南部环流同化试验 被引量:2
15
作者 肖贤俊 王东晓 徐建军 《科学通报》 EI CAS CSCD 北大核心 2006年第B11期27-32,共6页
对2000年夏季南海南部进行同化模拟,对比航次资料的分析表明,通过在模式中同化卫星遥感海面高度信息,改善了南海南部流场的模拟结果,减少了由于风应力强迫场误差带来的模式模拟的不确定性,得到了较为真实的2000年夏季南海南部温盐场.... 对2000年夏季南海南部进行同化模拟,对比航次资料的分析表明,通过在模式中同化卫星遥感海面高度信息,改善了南海南部流场的模拟结果,减少了由于风应力强迫场误差带来的模式模拟的不确定性,得到了较为真实的2000年夏季南海南部温盐场.为改善海洋环流模拟提供了一条新思路. 展开更多
关键词 南海 环流同化试验 边缘海
原文传递
与南海西边界流有关的区域海洋学进展 被引量:10
16
作者 王东晓 刘钦燕 +7 位作者 谢强 贺志刚 庄伟 舒业强 肖贤俊 宏波 吴湘玉 隋丹丹 《科学通报》 EI CAS CSCD 北大核心 2013年第14期1277-1288,共12页
从动力学角度,回顾了与南海西部主流系及其涡旋研究有关的南海物理海洋学的进展.主要内容涉及南海西边界流漂流浮标观测、西边界流源区涡旋遥感观测、西边界流源区跨陆架交换、涡流相互作用、西边界流年际变化、西边界流区海气相互作用... 从动力学角度,回顾了与南海西部主流系及其涡旋研究有关的南海物理海洋学的进展.主要内容涉及南海西边界流漂流浮标观测、西边界流源区涡旋遥感观测、西边界流源区跨陆架交换、涡流相互作用、西边界流年际变化、西边界流区海气相互作用、南海贯穿流等方面的进展.西边界流是一个强流系,其与季节和年际变化相关的结构特征、变率及跟某些动力过程的关联有待研究.在西边界流变率、涡旋环流相互作用、海气过程以及南海贯穿流取得了以下成果:(1)利用漂流浮标观测样本对南海西边界环流进行分析,指出南海西边界表层环流在到达越南中部沿岸后伴随着流幅变窄的同时流速加强;探讨了南海北部环流变化机理,定量诊断南海西边界流北支冬季逆风流产生的动力机制;利用航次数据对18°N断面的经向地转流进行诊断,表明南海西边界流的经向输送年际变化明显;(2)结合航次观测数据,对2003/2004年冬季南海北部2个反气旋涡旋和2007年夏季18°N附近的3个反气旋涡旋进行研究,指出冬季2个涡旋产生后以罗斯贝(Rossby)波速度(约0.1m/s)沿北部陆坡向西南方向传播,并初步揭示了南海西边界环流与夏季3个涡旋的相互作用;南海北部陆架区涡旋西南向传播最大(最大为0.09m/s),而越南以东海域涡动能(EKE)最大,这都说明涡旋活动与南海西边界流存在强的相互作用;(3)发现南海西边界流附近海表面温度(SST)强的季节内振荡特征,进一步研究表明此区域冬季SST季节内变化会使潜热季节内信号减弱20%;(4)探讨了南海贯穿流的长期变化特征以及与整个太平洋环流系统的相互关联. 展开更多
关键词 南海 西边界流 涡流相互作用 季节内变化 南海贯穿流
原文传递
一套全球海洋再分析资料初步结果分析 被引量:4
17
作者 王东晓 秦英豪 +2 位作者 肖贤俊 张祖强 武丰民 《科学通报》 CAS CSCD 北大核心 2012年第22期2102-2110,共9页
基于全球海洋资料同化系统(BCC_GODAS2.0)的20年(1990~2009年)再分析资料,利用OISST和SODA数据集从整体上进行检验,分析同化产品误差特征表明,同化系统能够有效改进海洋温度和盐度的估计,尤其在赤道东太平洋地区最为显著.垂向温盐总体R... 基于全球海洋资料同化系统(BCC_GODAS2.0)的20年(1990~2009年)再分析资料,利用OISST和SODA数据集从整体上进行检验,分析同化产品误差特征表明,同化系统能够有效改进海洋温度和盐度的估计,尤其在赤道东太平洋地区最为显著.垂向温盐总体RMSE相对同化前分别减小0.53℃和0.28psu,相对第一代同化系统BCC_GODAS1.0,温盐各层RMSE均有不同程度的减小.然后利用HadISST数据对比分析了热带太平洋海温异常变化特征.同化结果的Nio3指数与HadISST的海表温度异常有很好的一致性,相关系数高达93.6%;同化结果很好地模拟出了观测HadISST赤道太平洋季节变化特征,而赤道东太平洋地区逐月SSTA标准差季节变化的幅度变化较美国GODAS更接近HadISST;此外,同化在热带太平洋SSTA年际变率的EOF第一主模态所表现的ElNio空间结构,改善了冷舌地区的正异常年际变率强度,第二模态所表现的El Nio Modoki空间结构,改善了赤道东太平洋负异常局限于狭窄赤道地区的问题,并且赤道中太平洋海温正异常向南北高纬地区延伸的幅度和强度都得到加强,与观测HadISST比较吻合. 展开更多
关键词 BCC_GODAS2.0 温盐 EOF EL Nio EL Nio Modoki
原文传递
上一页 1 下一页 到第
使用帮助 返回顶部