Fast static correction methods for high-frequency multichannel marine seismic reflection data: A high-resolution seismic study of channel-levee systems on the Bengal Fan

Very high-frequency marine multichannel seismic reflection data generated by small-volume air- or waterguns allow detailed, high-resolution studies of sedimentary structures of the order of one to few metres wavelength. The high-frequency content, however, requires (1) a very exact knowledge of the source and receiver positions, and (2) the development of data processing methods which take this exact geometry into account. Static corrections are crucial for the quality of very high-frequency stacked data because static shifts caused by variations of the source and streamer depths are of the order of half to one dominant wavelength, so that they can lead to destructive interference during stacking of CDP sorted traces. As common surface-consistent residual static correction methods developed for land seismic data require fixed shot and receiver locations two simple and fast techniques have been developed for marine seismic data with moving sources and receivers to correct such static shifts. The first method – called CDP static correction method – is based on a simultaneous recording of Parasound sediment echosounder and multichannel seismic reflection data. It compares the depth information derived from the first arrivals of both data sets to calculate static correction time shifts for each seismic channel relative to the Parasound water depths. The second method – called average static correction method – utilises the fact that the streamer depth is mainly controlled by bird units, which keep the streamer in a predefined depth at certain increments but do not prevent the streamer from being slightly buoyant in-between. In case of calm weather conditions these streamer bendings mainly contribute to the overall static time shifts, whereas depth variations of the source are negligible. Hence, mean static correction time shifts are calculated for each channel by averaging the depth values determined at each geophone group position for several subsequent shots. Application of both methods to data of a high-resolution seismic survey of channel-levee systems on the Bengal Fan shows that the quality of the stacked section can be improved significantly compared to stacking results achieved without preceding static corrections. The optimised records show sedimentary features in great detail, that are not visible without static corrections. Limitations only result from the sea floor topography. The CDP static correction method generally provides more coherent reflections than the average static correction method but can only be applied in areas with rather flat sea floor, where no diffraction hyperbolae occur. In contrast, the average static correction method can also be used in regions with rough morphology, but the coherency of reflections is slightly reduced compared to the results of the CDP static correction method.     

The growth, collapse and quiescence of Teno volcano, Tenerife: new constraints from paleomagnetic data

Tenerife basically consists of three Miocene shield volcanoes, the Anaga, the Teno and Central shield, as well as the Pliocene Cañadas volcano. The temporal evolution and structural significance of each volcano with respect to the history of Tenerife is still a matter of debate. We present paleomagnetic results in order to enhance the view of the volcanic history of the Teno volcano by means of magnetostratigraphy. It is found that the initial subaerial phase shows reverse magnetizations throughout. After two major sector collapses, dominantly normally magnetized lavas extruded. Comparisons of observed magnetic polarities with the geomagnetic polarity timescale show that these volcanic activities occurred within 0.4 Myr between 6.3 and 5.9 Ma. Significantly younger flows, ～ 5.3 Myr old according to their radiometric age, revealed again normal polarity throughout. The absence of inversely magnetized lavas in-between the two normal periods indicates a volcanic hiatus or erosional phase. The evolutionary sequence and the estimated high production rates for the initial building phase are similar as would be expected for a hotspot volcano. The average geomagnetic field for 6.0 ± 0.2 Ma is close to an axial dipole field showing a slight far-sided/right-handed effect. The field strength, determined by Thellier-type intensity determinations, corresponds to a virtual axial dipole moment of 4.9 × 10²² A m². This value is approximately half of the present day field strength, but similar to values obtained for the mid-Miocene. It also corresponds to the proposed tertiary low-field level of the geomagnetic dipole moment.     

地震静校正全局最优化问题的求解

针对地震静校正存在的非线性和多参数性等问题，本文提出了一种全局优化的剩余静校正方法.在模型更新规则上，把模拟退火方法和均匀优化设计方法结合在一起，使模型的更新更为合理，从而加快了搜寻全局最优解的速度.同时提出了模型的不同分量的温度参数和退火过程的选取方式，把温度参数的选取与地震剖面的能量联系起来，使温度参数的选取具有自适应的特点.该方法克服了常规模拟退火方法所具有的寻优空间不均匀以及退火参数需通过多次试验选取的缺陷.通过实际地震资料的计算证明，本文提出的方法对地震静校正问题合理而有效.     

GPS精密定位中的海潮位移改正

根据海洋负荷潮理论，利用NAO99b全球海潮模型，计算了中国部分IGS站的海潮位移改正，并将海潮位移改正应用到GPS数据处理当中。在GAMIT软件的解算过程中，分别按加入和不加入海潮位移改正，对GPS基线分量和测站坐标分别进行了计算和比较分析。结果表明，海潮位移改正无论是对GPS基线分量还是对测站坐标，都有一定的影响。     

Absolute velocity of North America during the Mesozoic from paleomagnetic data

We use paleomagnetic data to map Mesozoic absolute motion of North America, using paleomagnetic Euler poles (PEP). First, we address two important questions: (1) How much clockwise rotation has been experienced by crustal blocks within and adjacent to the Colorado Plateau? (2) Why is there disagreement between the apparent polar wander (APW) path constructed using poles from southwestern North America and the alternative path based on poles from eastern North America? Regarding (1), a 10.5° clockwise rotation of the Colorado Plateau about a pole located near 35°N, 102°W seems to fit the evidence best. Regarding (2), it appears that some rock units from the Appalachian region retain a hard overprint acquired during the mid-Cretaceous, when the geomagnetic field had constant normal polarity and APW was negligible.We found three well-defined small-circle APW tracks: 245–200 Ma (PEP at 39.2°N, 245.2°E, R=81.1°, root mean square error (RMS)=1.82°), 200–160 Ma (38.5°N, 270.1°E, R=80.4°, RMS=1.06°), 160 to 125 Ma (45.1°N, 48.5°E, R=60.7°, RMS=1.84°). Intersections of these tracks (the “cusps” of Gordon et al. [Tectonics 3 (1984) 499]) are located at 59.6°N, 69.5°E (the 200 Ma or “J1” cusp) and 48.9°N, 144.0°E (the 160 Ma or “J2” cusp). At these times, the absolute velocity of North America appears to have changed abruptly.North America absolute motion also changed abruptly at the beginning and end of the Cretaceous APW stillstand, currently dated at about 125 and 88 Ma (J. Geophys. Res. 97 (1992b) 19651). During this interval, the APW path degenerates into a single point, implying rotation about an Euler pole coincident with the spin axis.Using our PEP and cusp locations, we calculate the absolute motion of seven points on the North American continent. Our intention is to provide a chronological framework for the analysis of Mesozoic tectonics. Clearly, if APW is caused by plate motion, abrupt changes in absolute motion should correlate with major tectonic events. This follows because large accelerations reflect important changes in the balance of forces acting on the plate, the most important of which are edge effects (subduction, terrane accretion, etc.). Some tectonic interpretations: (1) The J1 cusp may be associated with the inception of rifting of North America away from land masses to the east; the J2 cusp seems to mark the beginning of rapid spreading in the North Atlantic. (2) The J2 cusp signals the beginning of a period of rapid northwestward absolute motion of western North America; motion of tectonostratigraphic terranes in the westernmost Cordillera seems likely to have been directed toward the south during this interval. (3) The interval 88 to 80 Ma saw a rapid decrease in the paleolatitude of North America; unless this represents a period of true polar wander, terrane motion during this time should have been relatively northward.     

长距离单历元单频非差北斗网络差分方法

针对北斗卫星导航系统现有的网络差分方法多是基于GPS系统并采用双差定位模型,而双差定位模型的采用使得网络差分方法存在一些缺陷的问题,该文提出一种长距离单历元单频非差北斗网络差分方法;首先长距离北斗基准站网向覆盖区域内的北斗用户提供非差误差改正数;然后用户内插计算出自己所需的非差误差改正数,并改正北斗单频伪距观测值的误差;最后进行单历元单频网络差分定位。实验表明,该方法能够使用单历元单频伪距观测数据实现长距离北斗网络差分定位。     

Accuracy assessment of NLCD 2011 percent impervious cover for selected USA metropolitan areas

The emergence of high-resolution land cover data has created the opportunity to assess the accuracy of impervious cover (IC) provided by the National Land Cover Database (NLCD). We assessed the accuracy of the 900 m² NLCD2011 %IC for 18 metropolitan areas throughout the conterminous United States using reference data from 1 m² land cover data developed as part of the United States Environmental Protection Agency’s EnviroAtlas project. Agreement was assessed from two perspectives: 1) sensitivity to the size of the assessment unit used for the comparison, and 2) utility of NLCD %IC to serve as a proxy for high-resolution IC. The former perspective was considered because statistical relationships can be sensitive to assessment unit size and shape, and the latter perspective was considered because high resolution (reference) %IC data are not available nationwide. The utility of NLCD %IC as a proxy for the high resolution data was assessed for seven lattice (square) cell sizes ranging from 1 ha to 200 ha using four EnviroAtlas IC indicators: 1) %IC per 100 ha (1 km²); 2) %IC by Census block group; 3) %IC within a 15 m (radius) of the riparian zone, and; 4) %IC within a 50 m (radius) of the riparian zone. Agreement was quantified as per assessment unit deviation (NLCD %IC – reference %IC) and summarized as Mean Absolute Deviation (MAD) and Mean Deviation (MD) both within and across the 18 metropolitan areas. Ordinary least squares (OLS) regression (y = reference %IC and x = NLCD %IC) was also used to evaluate the quality of the NLCD %IC data. MAD was ≤ 5% for six of the seven lattice cell sizes. MAD was also ≤ 5% for Census block groups > 100 ha and for both riparian units. These results suggest that uncertainty attributable to the measurement of %IC was no greater than the uncertainty related to the effect of IC on aquatic resources that have been derived from studies of aquatic condition (e.g., benthic fauna) over a range of %IC. Overall, agreement was variable from one metropolitan area to the next. Agreement improved as assessment unit size increased and declined as the level of urbanization (NLCD %IC) increased. NLCD %IC tended to underestimate reference %IC overall, but NLCD %IC was sometimes greater than reference %IC in urbanized settings.     

联合BDS/GPS的北斗广域差分实时电离层延迟格网改正方法研究

研究并实现中国区域北斗广域差分实时电离层延迟格网改正算法,针对北斗单系统实时建立电离层延迟格网的不足,提出联合BDS/GPS的北斗广域差分实时电离层延迟改正方法。比较BDS、GPS以及BDS/GPS观测的实时电离层延迟格网,并利用北斗单频单点伪距定位进行精度验证。结果表明,在有效格网点覆盖充足的地区,BDS与GPS实时电离层延迟格网定位效果相当;而联合BDS/GPS观测的实时电离层延迟格网极大提高了偏远地区的单频定位精度与可定位历元数,也使原本单系统格网覆盖充足的地区定位效果得到进一步增强和稳定。     

青岛陆上观测点的气溶胶光学厚度特征分析及与渤海气溶胶光学厚度的比较

阐述了气溶胶散射对太阳光衰减的影响和对观测数据的统计分析.作者使用手持式五波段太阳光度计和臭氧监测仪对青岛崂山区陆域和渤海海域的大气衰减进行了较长时间的观测.分析表明,对可见光和近红外波段电磁渡的大气衰减贡献最大的是气溶胶散射,气溶胶散射的光学厚度随波长的增大而减小,但是它在大气总光学厚度中占的比例随波长的增大而增加.在青岛崂山区冬季,对应于440nm波长的气溶胶散射光学厚度主要分布在0.4～0.6和1.5～1.6区间,浑浊度系数α主要分布在0.15～0.45区间,埃斯特朗指数β主要分布在1.0～1.4区间.从时间上看,α和β在1d内的变化不大,但逐日变化非常显著.与晴天相比,青岛崂山区冬季雾天的气溶胶散射光学厚度增加了大约一倍,较大的浑浊度系数和较小的埃斯特朗指数表明,雾天的大气具有气溶胶散射的光学厚度大、粒子浓度大和粒径尺度大的特点.通过与海上观测数据的比较可见,青岛崂山区陆地冬季在雾天时的气溶胶分布与渤海夏季时相似,青岛崂山区陆地冬季在晴天无云时的气溶胶分布与渤海秋季时相似.     

Two different methods of geoidal height determinations using a spherical harmonic representation of the geopotential, topographic corrections and the height anomaly–geoidal height difference

 It is suggested that a spherical harmonic representation of the geoidal heights using global Earth gravity models (EGM) might be accurate enough for many applications, although we know that some short-wavelength signals are missing in a potential coefficient model. A `direct' method of geoidal height determination from a global Earth gravity model coefficient alone and an `indirect' approach of geoidal height determination through height anomaly computed from a global gravity model are investigated. In both methods, suitable correction terms are applied. The results of computations in two test areas show that the direct and indirect approaches of geoid height determination yield good agreement with the classical gravimetric geoidal heights which are determined from Stokes' formula. Surprisingly, the results of the indirect method of geoidal height determination yield better agreement with the global positioning system (GPS)-levelling derived geoid heights, which are used to demonstrate such improvements, than the results of gravimetric geoid heights at to the same GPS stations. It has been demonstrated that the application of correction terms in both methods improves the agreement of geoidal heights at GPS-levelling stations. It is also found that the correction terms in the direct method of geoidal height determination are mostly similar to the correction terms used for the indirect determination of geoidal heights from height anomalies. Received: 26 July 2001 / Accepted: 21 February 2002