A modified mix-differenced approach for estimating multi-GNSS real-time satellite clock offsets

A modified mixed-differenced approach for estimating multi-GNSS real-time clock offsets is presented. This approach, as compared to the earlier presented mixed-differenced approach which uses epoch-differenced and undifferenced observations, further adds a satellite-differenced process. The proposed approach, based on real-time orbit products and a mix of epoch-differenced and satellite-differenced observations to estimate only satellite clock offsets and tropospheric zenith wet delays, has fewer estimated parameters than other approaches, and thus its implementing procedure is efficient and can be performed and extended easily. To obtain high accuracy, the approach involves three steps. First, the high-accuracy tropospheric zenith wet delay of each station is estimated using mixed-differenced carrier phase observations. Second, satellite clock offset changes between adjacent epochs are estimated using also mixed-differenced carrier phase observations. Third, the satellite clock offsets at the initial epoch are estimated using satellite-differenced pseudorange observations. Finally, the initial epoch clock results and clock offset changes are concatenated to obtain the clock results of the current epoch. To validate the real-time satellite clock results, multi-GNSS post-processing clock products from IGS ACs were selected for comparison. From the comparison, the standard deviations of the GPS, GLONASS, BeiDou and Galileo systems clock results are approximately 0.1–0.4 ns, except for the BeiDou GEO satellites. The root mean squares are about 0.4–2.3 ns, which are similar to those of other international real-time products. When the clock estimates were assessed based on a pseudo-kinematic PPP procedure, the positioning accuracies in the East, North and Up components reach 5.6, 5.5 and 7.6 cm, respectively, which meet the centimeter level and are comparable to the application of other products.     

Fe-Ni and Mn-Cr isotopic systems in sulfides from unequilibrated enstatite chondrites

In situ measurements of ⁶⁰Fe-⁶⁰Ni and ⁵³Mn-⁵³Cr isotopic systems with an ion microprobe have been carried out for sulfide assemblages from unequilibrated enstatite chondrites (UECs). Evidence for the initial presence of ⁶⁰Fe has been observed in nine sulfide inclusions from three UECs: ALHA77295, MAC88136, and Qingzhen. The inferred initial (⁶⁰Fe/⁵⁶Fe) [(⁶⁰Fe/⁵⁶Fe)₀] ratios show a large variation range, from ∼2 × 10⁻⁷ to ∼2 × 10⁻⁶. The sulfide inclusions with high Fe/Ni ratios yield (⁶⁰Fe/⁵⁶Fe)₀ ratios of ∼(2-7) × 10⁻⁷, similar to most of the (⁶⁰Fe/⁵⁶Fe)₀ values of troilite and pyroxene observed in unequilibrated ordinary chondrites (UOCs). Inclusions with high inferred (⁶⁰Fe/⁵⁶Fe)₀ ratios (∼1-2 × 10⁻⁶) have low Fe/Ni ratios and the magnitude of the ⁶⁰Ni excesses is similar in two MAC88136 assemblages in spite of a difference of a factor of two in their Fe/Ni ratios. The inferred high (⁶⁰Fe/⁵⁶Fe)₀ ratios were probably the result of Fe-Ni re-distribution in the sulfides during later alteration processes.The ⁵³Mn-⁵³Cr system was measured in five of the sulfide assemblages that were examined for their ⁶⁰Fe-⁶⁰Ni systematics. The ⁵³Mn-⁵³Cr isochrons yielded variable initial (⁵³Mn/⁵⁵Mn) [(⁵³Mn/⁵⁵Mn)₀] ratios from ∼(2-7) × 10⁻⁷. There is no obvious correlation between the (⁶⁰Fe/⁵⁶Fe)₀ and (⁵³Mn/⁵⁵Mn)₀ ratios. The variable ⁵³Mn-⁵³Cr isochrons probably also indicate later disturbance to the isotopic systems in these sulfides. Even though no chronological information can be extracted from the ⁶⁰Fe-⁶⁰Ni and ⁵³Mn-⁵³Cr systems in these UEC sulfides, our results indicate that ⁶⁰Fe was present in the enstatite chondrite formation region of the early Solar System.     

一种小波与谱分析结合的GPS精密卫星钟差拟合方法研究

根据卫星钟固有的特性,即频移、频漂和频漂率以及周期性,提出一种新的拟合方法。新方法首先构造一个合理的钟差模型,该模型包含一个二次多项式和多个周期项。然后利用IGS提供的钟差产品数据,先粗略补充缺失数据,用二次多项式拟合;并对拟合后的残差进行小波分析,作降噪处理;最后通过谱分析,确定其主要的周期项,从而构建出适当的拟合模型,实现精密GPS卫星钟差拟合和预报。多天数据的实验结果表明,采用本文提出的新方法能够有效地对卫星钟差进行拟合和预报,满足不同目的的需求。     

北斗三号全球导航卫星系统全球广播电离层延迟修正模型(BDGIM)应用性能评估

2020年6月23日,我国北斗三号全球导航卫星系统正式完成星座全球组网。北斗三号全球导航卫星系统采用新一代全球广播电离层延迟修正模型(BDGIM),为用户提供电离层延迟改正服务。本文利用高精度全球电离层格网(GIM)以及实测BDS/GPS数据提供的电离层TEC作为参考,从延迟改正精度及北斗单频伪距单点定位应用、模型系数性能等方面,对北斗三号系统组网前后(2020年5月1日至2020年7月20日)BDGIM模型的改正精度等应用性能进行了分析与研究,并将其与美国GPS播发的Klobuchar模型和北斗二号卫星导航系统播发的BDS Klobuchar模型进行对比。研究表明,BDGIM模型在对北斗三号系统组网完成前后电离层延迟修正精度没有发生显著变化。上述时段内,以国际GNSS服务(IGS)发布的最终GIM产品为参考,BDGIM模型在中国区域、亚太地区和全球范围内的电离层修正百分比分别达到84.45%、74.74%和64.57%;以选取的全球83个GNSS检测站BDS、GPS双频数据实测电离层TEC为参考,BDGIM在中国区域、亚太地区和全球范围内的电离层修正百分比分别为73.12%、70.18%及68.06%;当BDGIM模型应用于北斗单频伪距单点定位时,在中国区域、亚太地区和全球范围内分别实现了2.22、2.66和2.96 m的三维定位精度。     

多频多模GNSS接收机差分相位偏差的短期时变特性

随着中国北斗三号全球导航卫星系统(BeiDou-3 Navigation Satellite System,BDS-3)的建成、欧盟伽利略系统(Galileo)及日本准天顶卫星系统(quasi-zenith satellite system,QZSS)的发展,越来越多的卫星可用于反演大气电离层.通常,接收机差分码偏差(differential code biases,DCB)的短时变化被认为是利用全球导航卫星系统(Global Navigation Satellite System,GNSS)反演电离层的重要误差来源,然而,有研究表明,接收机差分相位偏差(differential phase biases,DPB)的短时变化也有可能影响电离层反演的精度和可靠性.为此,本文提出了基于站间单差模型并采用不变换参考星策略来估计接收机DPB的方法,可实现接收机DPB的连续估计.基于几台可跟踪BDS-3信号的多频多模接收机采集的数据,对BDS-3、Galileo、GPS和QZSS重叠频率组合的DPB进行了分析.结果表明,四系统的接收机DPB日变化都是很明显的,并且和温度有很强的相关性;基于不同系统重叠频率组合的DPB之间存在强相关;基于相同类型接收机的DPB的变化也存在明显的相关性.     

Petrogenesis of the Northwest Africa 4898 high‐Al mare basalt

Northwest Africa (NWA) 4898 is the only low‐Ti, high‐Al basaltic lunar meteorite yet recognized. It predominantly consists of pyroxene (53.8 vol%) and plagioclase (38.6 vol%). Pyroxene has a wide range of compositions (En_12–62Fs_25–62Wo_11–36), which display a continuous trend from Mg‐rich cores toward Ca‐rich mantles and then to Fe‐rich rims. Plagioclase has relatively restricted compositions (An_87–96Or_0–1Ab_4–13), and was transformed to maskelynite. The REE zoning of all silicate minerals was not significantly modified by shock metamorphism and weathering. Relatively large (up to 1 mm) olivine phenocrysts have homogenous inner parts with Fo ~74 and sharply decrease to 64 within the thin out rims (~30 μm in width). Four types of inclusions with a variety of textures and modal mineralogy were identified in olivine phenocrysts. The contrasting morphologies of these inclusions and the chemical zoning of olivine phenocrysts suggest NWA 4898 underwent at least two stages of crystallization. The aluminous chromite in NWA 4898 reveals that its high alumina character was inherited from the parental magma, rather than by fractional crystallization. The mineral chemistry and major element compositions of NWA 4898 are different from those of 12038 and Luna 16 basalts, but resemble those of Apollo 14 high‐Al basalts. However, the trace element compositions demonstrate that NWA 4898 and Apollo 14 high‐Al basalts could not have been derived from the same mantle source. REE compositions of its parental magma indicate that NWA 4898 probably originated from a unique depleted mantle source that has not been sampled yet. Unlike Apollo 14 high‐Al basalts, which assimilated KREEPy materials during their formation, NWA 4898 could have formed by closed‐system fractional crystallization.     

Origin of a metamorphosed lithic clast in CM chondrite Grove Mountains 021536

Abstract– A metamorphosed lithic clast was discovered in the CM chondrite Grove Mountains 021536, which was collected in the Antarctica by the Chinese Antarctic Research Exploration team. The lithic clast is composed mainly of Fe‐rich olivine (Fo₆₂) with minor diopside (Fs_9.7–11.1Wo_48.3–51.6), plagioclase (An_43–46.5), nepheline, merrillite, Al‐rich chromite (21.8 wt% Al₂O₃; 4.43 wt% TiO₂), and pentlandite. Δ¹⁷O values of olivine in the lithic clast vary from ?3.9‰ to ?0.8‰. Mineral compositions and oxygen isotopic compositions of olivine suggest that the lithic clast has an exotic source different from the CM chondrite parent body. The clast could be derived from strong thermal metamorphism of pre‐existing chondrule that has experienced low‐temperature anhydrous alteration. The lithic clast is similar in mineral assemblage and chemistry to a few clasts observed in oxidized CV3 chondrites (Mokoia and Yamato‐86009) and might have been derived from the interior of the primitive CV asteroid. The apparent lack of hydration in the lithic clast indicates that the clast accreted into the CM chondrite after hydration of the CM components.     

基于精选基准消秩亏的GNSS参考网数据处理方法

全球导航卫星系统（GNSS）参考网多用于估计卫星轨道/钟差、监测地表形变和速度场、确定精密地球自转参数等方面。相关数据处理模式包括：双差基线解（DD）和非差精密单点定位（PPP）等。本文首先从GNSS基本观测方程出发,通过选取两组基准参数,导出了上述两模式下的列满秩观测方程,然后分析了它们的不足,例如：相位偏差在DD模式中吸收了钟差,丧失了时不变特性;模糊度在PPP模式中吸收了相位偏差,失去了整数性。基于上述分析,本文提出了一种新的参考网数据处理方案,以充分融合DD和PPP模式的优势。它的关键策略是精选基准参数,以达到消秩亏的目的,具体优点体现在：相位偏差独立可估,若合理约束为时不变参数,可充分减少参数个数,提高网解精度;待估模糊度具备整周特性,经由模糊度固定,可改善网解可靠性。     

Petrogenesis of the new lherzolitic shergottite Grove Mountains 99027: Constraints of petrography,mineral chemistry,and rare earth elements

Abstract— We report petrography, mineral chemistry, and microdistribution of rare earth elements (REE) in a new lherzolitic shergottite, Grove Mountains (GRV) 99027. The textural relationship and REE patterns of minerals suggest precipitation of cumulus olivine and chromite, followed by equilibrium crystallization of a closed system with a bulk composition of the inferred intercumulus melt. Subsolidus equilibrium temperatures of pyroxenes and olivine range from 1100 to 1210 °C, based on a two‐pyroxene thermometry and Ca partitioning between augite and olivine. Oxygen fugacity of the parent magma is 1.5–2.5 (av. 2.0 ± 0.4) log units below the quartz‐fayalite‐magnetite (QFM) buffer at 960–1360 °C, according to the olivine‐orthopyroxene‐chromite barometer. The ilmenite‐chromite barometer and thermometer show much wider ranges of oxygen fugacity (1.0–7.0 log unit below QFM) and temperature (1130–480 °C), suggesting subsolidus equilibration of the oxides at low temperatures, probably due to deep burial of GRV 99027 on Mars. The low oxygen fugacity and LREE depletion of the parent magma of GRV 99027 suggest low contamination by martian crust. Characteristics of GRV 99027 demonstrate similarity of lherzolitic shergottites, suggesting a high possibility of launch pairing or a homogeneous upper mantle of Mars if they were ejected by individual impact events. However, GRV 99027 probably experienced severe post‐shock thermal metamorphism in comparison with other lherzolitic shergottites, based on the re‐crystallization of maskelynite, the homogeneity of minerals, and the low subsolidus equilibrium temperatures between chromite and ilmenite.     

A grid-based tropospheric product for China using a GNSS network

Tropospheric delay accounts for one source of error in global navigation satellite systems (GNSS). To better characterize the tropospheric delays in the temporal and spatial domain and facilitate the safety-critical use of GNSS across China, a method is proposed to generate a grid-based tropospheric product (GTP) using the GNSS network with an empirical tropospheric model, known as IGGtrop. The prototype system generates the GTPs in post-processing and real-time modes and is based on the undifferenced and uncombined precise point positioning (UU-PPP) technique. GTPs are constructed for a grid form (\(2.0{^{\circ }}\times 2.5{^{\circ }}\) latitude–longitude) over China with a time resolution of 5 min. The real-time GTP messages are encoded in a self-defined RTCM3 format and broadcast to users using NTRIP (networked transport of RTCM via internet protocol), which enables efficient and safe transmission to real-time users. Our approach for GTP generation consists of three sequential steps. In the first step, GNSS-derived zenith tropospheric delays (ZTDs) for a network of GNSS stations are estimated using UU-PPP. In the second step, vertical adjustments for the GNSS-derived ZTDs are applied to address the height differences between the GNSS stations and grid points. The ZTD height corrections are provided by the IGGtrop model. Finally, an inverse distance weighting method is used to interpolate the GNSS-derived ZTDs from the surrounding GNSS stations to the location of the grid point. A total of 210 global positioning system (GPS) stations from the crustal movement observation network of China are used to generate the GTPs in both post-processing and real-time modes. The accuracies of the GTPs are assessed against with ERA-Interim-derived ZTDs and the GPS-derived ZTDs at 12 test GPS stations, respectively. The results show that the post-processing and real-time GTPs can provide the ZTDs with accuracies of 1.4 and 1.8 cm, respectively. We also apply the GTPs in real-time kinematic GPS PPP, and the results show that the convergence time of the PPP solutions is shortened. These results confirm that the GTPs can act as an efficient information source to augment GNSS positioning over China.