Earth Rotation and Gravity Field Parameters from Satellite Laser Ranging

We present the results from a simultaneous estimation of the gravity field, Earth rotation parameters, and station coordinates from combined SLR solutions incorporating up to nine geodetic satellites: LAGEOS-1/2, Starlette, Stella, AJISAI, Beacon-C, Lares, Blits and LARES. These solutions cover all three pillars of satellite geodesy and ensure full consistency between the Earth rotation parameters, gravity field coefficients, and geometry-related parameters. We address benefits emerging from such an approach and discuss particular aspects and limitations of the gravity field recovery using SLR data. The current accuracy of SLR-derived polar motion, by the means of WRMS w.r.t. IERS-08-C04 series, is at a level of 118-149 μas, which corresponds to 4 to 5 mm on the Earth’s surface. The WRMS of SLR-derived Length-of-Day, when the gravity field parameters are simultaneously estimated, is 56 μs/day, corresponding to about 26 mm on the ground, and the mean bias of SLR-derived Length-of-Day is 6.3 μs/day, corresponding to 3 mm.

Time variable Earth’s gravity field from SLR satellites

The time variable Earth’s gravity field contains information about the mass transport within the system Earth, i.e., the relationship between mass variations in the atmosphere, oceans, land hydrology, and ice sheets. For many years, satellite laser ranging (SLR) observations to geodetic satellites have provided valuable information of the low-degree coefficients of the Earth’s gravity field. Today, the Gravity Recovery and Climate Experiment (GRACE) mission is the major source of information for the time variable field of a high spatial resolution. We recover the low-degree coefficients of the time variable Earth’s gravity field using SLR observations up to nine geodetic satellites: LAGEOS-1, LAGEOS-2, Starlette, Stella, AJISAI, LARES, Larets, BLITS, and Beacon-C. We estimate monthly gravity field coefficients up to degree and order 10/10 for the time span 2003–2013 and we compare the results with the GRACE-derived gravity field coefficients. We show that not only degree-2 gravity field coefficients can be well determined from SLR, but also other coefficients up to degree 10 using the combination of short 1-day arcs for low orbiting satellites and 10-day arcs for LAGEOS-1/2. In this way, LAGEOS-1/2 allow recovering zonal terms, which are associated with long-term satellite orbit perturbations, whereas the tesseral and sectorial terms benefit most from low orbiting satellites, whose orbit modeling deficiencies are minimized due to short 1-day arcs. The amplitudes of the annual signal in the low-degree gravity field coefficients derived from SLR agree with GRACE K-band results at a level of 77 %. This implies that SLR has a great potential to fill the gap between the current GRACE and the future GRACE Follow-On mission for recovering of the seasonal variations and secular trends of the longest wavelengths in gravity field, which are associated with the large-scale mass transport in the system Earth.

CODE’s new solar radiation pressure model for GNSS orbit determination

The Empirical CODE Orbit Model (ECOM) of the Center for Orbit Determination in Europe (CODE), which was developed in the early 1990s, is widely used in the International GNSS Service (IGS) community. For a rather long time, spurious spectral lines are known to exist in geophysical parameters, in particular in the Earth Rotation Parameters (ERPs) and in the estimated geocenter coordinates, which could recently be attributed to the ECOM. These effects grew creepingly with the increasing influence of the GLONASS system in recent years in the CODE analysis, which is based on a rigorous combination of GPS and GLONASS since May 2003. In a first step we show that the problems associated with the ECOM are to the largest extent caused by the GLONASS, which was reaching full deployment by the end of 2011. GPS-only, GLONASS-only, and combined GPS/GLONASS solutions using the observations in the years 2009–2011 of a global network of 92 combined GPS/GLONASS receivers were analyzed for this purpose. In a second step we review direct solar radiation pressure (SRP) models for GNSS satellites. We demonstrate that only even-order short-period harmonic perturbations acting along the direction Sun-satellite occur for GPS and GLONASS satellites, and only odd-order perturbations acting along the direction perpendicular to both, the vector Sun-satellite and the spacecraft’s solar panel axis. Based on this insight we assess in the third step the performance of four candidate orbit models for the future ECOM. The geocenter coordinates, the ERP differences w. r. t. the IERS 08 C04 series of ERPs, the misclosures for the midnight epochs of the daily orbital arcs, and scale parameters of Helmert transformations for station coordinates serve as quality criteria. The old and updated ECOM are validated in addition with satellite laser ranging (SLR) observations and by comparing the orbits to those of the IGS and other analysis centers. Based on all tests, we present a new extended ECOM which substantially reduces the spurious signals in the geocenter coordinate z (by about a factor of 2–6), reduces the orbit misclosures at the day boundaries by about 10 %, slightly improves the consistency of the estimated ERPs with those of the IERS 08 C04 Earth rotation series, and substantially reduces the systematics in the SLR validation of the GNSS orbits.

The impact of common versus separate estimation of orbit parameters on GRACE gravity field solutions

Gravity field parameters are usually determined from observations of the GRACE satellite mission together with arc-specific parameters in a generalized orbit determination process. When separating the estimation of gravity field parameters from the determination of the satellites’ orbits, correlations between orbit parameters and gravity field coefficients are ignored and the latter parameters are biased towards the a priori force model. We are thus confronted with a kind of hidden regularization. To decipher the underlying mechanisms, the Celestial Mechanics Approach is complemented by tools to modify the impact of the pseudo-stochastic arc-specific parameters on the normal equations level and to efficiently generate ensembles of solutions. By introducing a time variable a priori model and solving for hourly pseudo-stochastic accelerations, a significant reduction of noisy striping in the monthly solutions can be achieved. Setting up more frequent pseudo-stochastic parameters results in a further reduction of the noise, but also in a notable damping of the observed geophysical signals. To quantify the effect of the a priori model on the monthly solutions, the process of fixing the orbit parameters is replaced by an equivalent introduction of special pseudo-observations, i.e., by explicit regularization. The contribution of the thereby introduced a priori information is determined by a contribution analysis. The presented mechanism is valid universally. It may be used to separate any subset of parameters by pseudo-observations of a special design and to quantify the damage imposed on the solution.

Time Variable Gravity: Contributions of GOCE Satellite Data to Monthly and Bi-monthly GRACE Gravity Estimates

A feasibility study by Pail et al. (Can GOCE help to improve temporal gravity field estimates? In: Ouwehand L (ed) Proceedings of the 4th International GOCE User Workshop, ESA Publication SP-696, 2011b) shows that GOCE (‘Gravity field and steady-state Ocean Circulation Explorer’) satellite gravity gradiometer (SGG) data in combination with GPS derived orbit data (satellite-to-satellite tracking: SST-hl) can be used to stabilize and reduce the striping pattern of a bi-monthly GRACE (‘Gravity Recovery and Climate Experiment’) gravity field estimate. In this study several monthly (and bi-monthly) combinations of GRACE with GOCE SGG and GOCE SST-hl data on the basis of normal equations are investigated. Our aim is to assess the role of the gradients (solely) in the combination and whether already one month of GOCE observations provides sufficient data for having an impact in the combination. The estimation of clean and stable monthly GOCE SGG normal equations at high resolution ( >  d/o 150) is found to be difficult, and the SGG component, solely, does not show significant added value to monthly and bi-monthly GRACE gravity fields. Comparisons of GRACE-only and combined monthly and bi-monthly solutions show that the striping pattern can only be reduced when using both GOCE observation types (SGG, SST-hl), and mainly between d/o 45 and 60.

Satellite laser ranging to GPS and GLONASS

Satellite laser ranging (SLR) to the satellites of the global navigation satellite systems (GNSS) provides substantial and valuable information about the accuracy and quality of GNSS orbits and allows for the SLR-GNSS co-location in space. In the framework of the NAVSTAR-SLR experiment two GPS satellites of Block-IIA were equipped with laser retroreflector arrays (LRAs), whereas all satellites of the GLONASS system are equipped with LRAs in an operational mode. We summarize the outcome of the NAVSTAR-SLR experiment by processing 20 years of SLR observations to GPS and 12 years of SLR observations to GLONASS satellites using the reprocessed microwave orbits provided by the center for orbit determination in Europe (CODE). The dependency of the SLR residuals on the size, shape, and number of corner cubes in LRAs is studied. We show that the mean SLR residuals and the RMS of residuals depend on the coating of the LRAs and the block or type of GNSS satellites. The SLR mean residuals are also a function of the equipment used at SLR stations including the single-photon and multi-photon detection modes. We also show that the SLR observations to GNSS satellites are important to validate GNSS orbits and to assess deficiencies in the solar radiation pressure models. We found that the satellite signature effect, which is defined as a spread of optical pulse signals due to reflection from multiple reflectors, causes the variations of mean SLR residuals of up to 15 mm between the observations at nadir angles of 0∘ and 14∘. in case of multi-photon SLR stations. For single-photon SLR stations this effect does not exceed 1 mm. When using the new empirical CODE orbit model (ECOM), the SLR mean residual falls into the range 0.1–1.8 mm for high-performing single-photon SLR stations observing GLONASS-M satellites with uncoated corner cubes. For best-performing multi-photon stations the mean SLR residuals are between −12.2 and −25.6 mm due to the satellite signature effect.

Evaluating gillnetting protocols to characterize lacustrine fish communities

Ecological research and monitoring of lacustrine ecosystems often requires a whole-lake assessment of fish communities. Gillnet sampling offers an efficient means of estimating abundance, biomass and fish community composition. However the choice of gillnet sampling protocol may influence lake characterization via physical properties of the nets and allocation of sampling effort between littoral, benthic and pelagic habitats. This paper compares two commonly used, whole-lake sampling protocols applied across 17 prealpine, subalpine and alpine European lakes ranging widely in size, depth and altitude to determine their relative strength for research and management applications. Effort-corrected estimates of abundance, biomass and species richness were correlated between the protocols and both distinguished the trout-dominated alpine communities from subalpine and prealpine lakes dominated by whitefish and perch. A considerable amount of variance remained unexplained between the two protocols however, which seemed to correspond with differences in the proportion of effort among benthic and pelagic habitats. We suggest that both the European standard (CEN) and vertical (VERT) netting protocols are suitable for assessing ecological status and monitoring changes in lake fish communities through time. However the details of each protocol should be kept in mind when comparing fish communities between lakes. Mesh sizes used in CEN nets produce a more even size frequency distribution, suggesting that this protocol is most appropriate for assessing size structure of fish assemblages. The high proportion of netting effort in benthic habitats shallower than 70 m depth under the CEN protocol means that, particularly in larger lakes, outcomes will be disproportionately influenced by the ecological condition of this habitat. The VERT protocol presumably provides a more accurate estimate of whole-lake CPUE and community composition because effort, in terms of net area, is more evenly distributed across the entire volume of the lake. This is particularly important in large and deep lakes where pelagic habitats occupy a high proportion of the lake volume.

Active eosinophilic esophagitis is characterized by epithelial barrier defects and eosinophil extracellular trap formation.

BACKGROUND Eosinophilic esophagitis (EoE) exhibits esophageal dysfunction owing to an eosinophil-predominant inflammation. Activated eosinophils generate eosinophil extracellular traps (EETs) able to kill bacteria. There is evidence of an impaired barrier function in EoE that might allow pathogens to invade the esophagus. This study aimed to investigate the presence and distribution of EETs in esophageal tissues from EoE patients and their association with possible epithelial barrier defects. METHODS Anonymized tissue samples from 18 patients with active EoE were analyzed. The presence of DNA nets associated with eosinophil granule proteins forming EETs and the expression of filaggrin, the protease inhibitor lympho-epithelial Kazal-type-related inhibitor (LEKTI), antimicrobial peptides, and cytokines were evaluated by confocal microscopy following immune fluorescence staining techniques. RESULTS Eosinophil extracellular trap formation occurred frequently and was detected in all EoE samples correlating with the numbers of infiltrating eosinophils. While the expression of both filaggrin and LEKTI was reduced, epithelial antimicrobial peptides (human beta-defensin-2, human beta-defensin-3, cathelicidin LL-37, psoriasin) and cytokines (TSLP, IL-25, IL-32, IL-33) were elevated in EoE as compared to normal esophageal tissues. There was a significant correlation between EET formation and TSLP expression (P = 0.02) as well as psoriasin expression (P = 0.016). On the other hand, a significant negative correlation was found between EET formation and LEKTI expression (P = 0.016). CONCLUSION Active EoE exhibits the presence of EETs. Indications of epithelial barrier defects in association with epithelial cytokines are also present which may have contributed to the activation of eosinophils. The formation of EETs could serve as a firewall against the invasion of pathogens.

A consistent combination of GNSS and SLR with minimum constraints

In this article, the realization of a global terrestrial reference system (TRS) based on a consistent combination of Global Navigation Satellite System (GNSS) and Satellite Laser Ranging (SLR) is studied. Our input data consists of normal equation systems from 17 years (1994– 2010) of homogeneously reprocessed GPS, GLONASS and SLR data. This effort used common state of the art reduction models and the same processing software (Bernese GNSS Software) to ensure the highest consistency when combining GNSS and SLR. Residual surface load deformations are modeled with a spherical harmonic approach. The estimated degree-1 surface load coefficients have a strong annual signal for which the GNSS- and SLR-only solutions show very similar results. A combination including these coefficients reduces systematic uncertainties in comparison to the singletechnique solution. In particular, uncertainties due to solar radiation pressure modeling in the coefficient time series can be reduced up to 50 % in the GNSS+SLR solution compared to the GNSS-only solution. In contrast to the ITRF2008 realization, no local ties are used to combine the different geodetic techniques.We combine the pole coordinates as global ties and apply minimum constraints to define the geodetic datum. We show that a common origin, scale and orientation can be reliably realized from our combination strategy in comparison to the ITRF2008.

Estimating whole-lake fish catch per unit effort

The European standard for gillnetsampling to characterize lake fish communities stratifies sampling effort (i.e., number of nets) within depth strata. Nets to sample benthic habitats are randomly distributed throughout the lake within each depth strata. Pelagic nets are also stratified by depth, but are set only at the deepest point of the lake. Multiple authors have suggested that this design under-represents pelagic habitats, resulting in estimates of whole-lake CPUE and community composition which are disproportionately influenced by ecological conditions of littoral and benthic habitats. To address this issue, researchers have proposed estimating whole-lake CPUE by weighting the catch rate in each depth-compartment by the proportion of the volume of the lake contributed by the compartment. Our study aimed to assess the effectiveness of volume-weighting by applying it to fish communities sampled according to the European standard (CEN), and by a second whole-lake gillnetting protocol (VERT), which prescribes additional fishing effort in pelagic habitats. We assume that convergence between the protocols indicates that volume-weighting provides a more accurate estimate of whole-lake catch rate and community composition. Our results indicate that volume-weighting improves agreement between the protocols for whole-lake total CPUE, estimated proportion of perch and roach and the overall fish community composition. Discrepancies between the protocols remaining after volume-weighting maybe because sampling under the CEN protocol overlooks horizontal variation in pelagic fish communities. Analyses based on multiple pelagic-set VERT nets identified gradients in the density and biomass of pelagic fish communities in almost half the lakes that corresponded with the depth of water at net-setting location and distance along the length of a lake. Additional CEN pelagic sampling effort allocated across water depths and distributed throughout the lake would therefore help to reconcile differences between the sampling protocols and, in combination with volume-weighting, converge on a more accurate estimate of whole-lake fish communities.

Pathology: Classification and Immunoprofile

The classification of neuroendocrine neoplasms (NENs) has been evolving steadily over the last decades. Important prognostic factors of NENs are their proliferative activity and presence/absence of necrosis. These factors are reported in NENs of all body sites; however, the terminology as well as the exact rules of classification differ according to the location of the primary tumor. Only in gastroenteropancreatic (GEP) NENs a formal grading is performed. This grading is based on proliferation assessed by the mitotic count and/or Ki-67 proliferation index. In the lung, NEN grading is an intrinsic part of the tumor designation with typical carcinoids corresponding to neuroendocrine tumor (NET) G1 and atypical carcinoids to NET G2; however, the presence or absence of necrotic foci is as important as proliferation for the differentiation between typical and atypical carcinoids. Immunohistochemical markers can be used to demonstrate neuroendocrine differentiation. Synaptophysin and chromogranin A are, to date, the most reliable and most commonly used for this purpose. Beyond this, other markers can be helpful, for example in the situation of a NET metastasis of unknown primary, where a hormonal profile or a panel of transcription factors can give hints to the primary site. Many immunohistochemical markers have been shown to correlate with prognosis but are not used in clinical practice, for example cytokeratin 19 and KIT expression in pancreatic NETs. There is no predictive biomarker in use, with the exception of somatostatin receptor (SSTR) 2 expression for predicting the amenability of a tumor to in vivo SSTR targeting for imaging or therapy.

CODE ultra-rapid product series for the IGS

CODE, the Center for Orbit Determination in Europe, is a joint venture of the following four institutions: Astronomical Institute, University of Bern (AIUB), Bern, Switzerland; Federal Office of Topography swisstopo, Wabern, Switzerland; Federal Agency of Cartography and Geodesy (BKG), Frankfurt a. M., Germany; Institut für Astronomische und Physikalische Geodäsie, Technische Universität München (IAPG, TUM), Munich, Germany. It acts as a global analysis center of the International GNSS Service (IGS). The operational computations are performed at AIUB using the latest development version of the Bernese GNSS Software. In this context an ultra-rapid solution series is generated considering GPS and GLONASS satellites. It is updated several times per day and contains 24 hours of observed and 24 hours of predicted orbit interval. More details are available in: Lutz, S., G. Beutler, S. Schaer, R. Dach, A. Jäggi; 2014: CODE's new ultra-rapid orbit and ERP products for the IGS. GPS Solutions. DOI 10.1007/s10291-014-0432-2

CODE repro2 product series for the IGS

CODE, the Center for Orbit Determination in Europe, is a joint venture of the following four institutions: • Astronomical Institute, University of Bern (AIUB), Bern, Switzerland • Federal Office of Topography swisstopo, Wabern, Switzerland • Federal Agency of Cartography and Geodesy (BKG), Frankfurt a. M., Germany • Institut für Astronomische und Physikalische Geodäsie, Technische Universität München (IAPG, TUM), Munich, Germany It acts as a global analysis center of the International GNSS Service (IGS, Dow et al, 2009). The operational computations are performed at AIUB using the latest development version of the Bernese GNSS Software (Dach et al., 2015). In this context the contribution to the IGS repro02 effort is generated considering only the GPS satellites between 1994 and 2001 as well as the GPS and GLONASS satellites from 2002 to the end of 2013.