Application of ray-traced tropospheric slant delays to geodetic VLBI analysis

The correction of tropospheric influences via so-called path delays is critical for the analysis of observations from space geodetic techniques like the very long baseline interferometry (VLBI). In standard VLBI analysis, the a priori slant path delays are determined using the concept of zenith delays, mapping functions and gradients. The a priori use of ray-traced delays, i.e., tropospheric slant path delays determined with the technique of ray-tracing through the meteorological data of numerical weather models (NWM), serves as an alternative way of correcting the influences of the troposphere on the VLBI observations within the analysis. In the presented research, the application of ray-traced delays to the VLBI analysis of sessions in a time span of 16.5 years is investigated. Ray-traced delays have been determined with program RADIATE (see Hofmeister in Ph.D. thesis, Department of Geodesy and Geophysics, Faculty of Mathematics and Geoinformation, Technische Universität Wien. http://resolver.obvsg.at/urn:nbn:at:at-ubtuw:1-3444, 2016) utilizing meteorological data provided by NWM of the European Centre for Medium-Range Weather Forecasts (ECMWF). In comparison with a standard VLBI analysis, which includes the tropospheric gradient estimation, the application of the ray-traced delays to an analysis, which uses the same parameterization except for the a priori slant path delay handling and the used wet mapping factors for the zenith wet delay (ZWD) estimation, improves the baseline length repeatability (BLR) at 55.9% of the baselines at sub-mm level. If no tropospheric gradients are estimated within the compared analyses, 90.6% of all baselines benefit from the application of the ray-traced delays, which leads to an average improvement of the BLR of 1 mm. The effects of the ray-traced delays on the terrestrial reference frame are also investigated. A separate assessment of the RADIATE ray-traced delays is carried out by comparison to the ray-traced delays from the National Aeronautics and Space Administration Goddard Space Flight Center (NASA GSFC) (Eriksson and MacMillan in http://lacerta.gsfc.nasa.gov/tropodelays, 2016) with respect to the analysis performances in terms of BLR results. If tropospheric gradient estimation is included in the analysis, 51.3% of the baselines benefit from the RADIATE ray-traced delays at sub-mm difference level. If no tropospheric gradients are estimated within the analysis, the RADIATE ray-traced delays deliver a better BLR at 63% of the baselines compared to the NASA GSFC ray-traced delays.     

Nocturnal Near-Surface Temperature,but not Flow Dynamics,can be Predicted by Microtopography in a Mid-Range Mountain Valley

We investigate nocturnal flow dynamics and temperature behaviour near the surface of a 170-m long gentle slope in a mid-range mountain valley. In contrast to many existing studies focusing on locations with significant topographic variations, gentle slopes cover a greater spatial extent of the Earth’s surface. Air temperatures were measured using the high-resolution distributed-temperature-sensing method within a two-dimensional fibre-optic array in the lowest metre above the surface. The main objectives are to characterize the spatio-temporal patterns in the near-surface temperature and flow dynamics, and quantify their responses to the microtopography and land cover. For the duration of the experiment, including even clear-sky nights with weak winds and strong radiative forcing, the classical cold-air drainage predicted by theory could not be detected. In contrast, we show that the airflow for the two dominant flow modes originates non-locally. The most abundant flow mode is characterized by vertically-decoupled layers featuring a near-surface flow perpendicular to the slope and strong stable stratification, which contradicts the expectation of a gravity-driven downslope flow of locally produced cold air. Differences in microtopography and land cover clearly affect spatio-temporal temperature perturbations. The second most abundant flow mode is characterized by strong mixing, leading to vertical coupling with airflow directed down the local slope. Here variations of microtopography and land cover lead to negligible near-surface temperature perturbations. We conclude that spatio-temporal temperature perturbations, but not flow dynamics, can be predicted by microtopography, which complicates the prediction of advective-heat components and the existence and dynamics of cold-air pools in gently sloped terrain in the absence of observations.     

Impact of in-consistency between the climate model and its initial conditions on climate prediction

We investigated the influence of dynamical in-consistency of initial conditions on the predictive skill of decadal climate predictions. The investigation builds on the fully coupled global model “Coupled GCM for Earth Simulator” (CFES). In two separate experiments, the ocean component of the coupled model is full-field initialized with two different initial fields from either the same coupled model CFES or the GECCO2 Ocean Synthesis while the atmosphere is initialized from CFES in both cases. Differences between both experiments show that higher SST forecast skill is obtained when initializing with coupled data assimilation initial conditions (CIH) instead of those from GECCO2 (GIH), with the most significant difference in skill obtained over the tropical Pacific at lead year one. High predictive skill of SST over the tropical Pacific seen in CIH reflects the good reproduction of El Niño events at lead year one. In contrast, GIH produces additional erroneous El Niño events. The tropical Pacific skill differences between both runs can be rationalized in terms of the zonal momentum balance between the wind stress and pressure gradient force, which characterizes the upper equatorial Pacific. In GIH, the differences between the oceanic and atmospheric state at initial time leads to imbalance between the zonal wind stress and pressure gradient force over the equatorial Pacific, which leads to the additional pseudo El Niño events and explains reduced predictive skill. The balance can be reestablished if anomaly initialization strategy is applied with GECCO2 initial conditions and improved predictive skill in the tropical Pacific is observed at lead year one. However, initializing the coupled model with self-consistent initial conditions leads to the highest skill of climate prediction in the tropical Pacific by preserving the momentum balance between zonal wind stress and pressure gradient force along the equatorial Pacific.

     

Efficiency of reduced-order, time-dependent adjoint data assimilation approaches

Applications of adjoint data assimilation, which is designed to bring an ocean circulation model into consistency with ocean observations, are computationally demanding. To improve the convergence rate of an optimization, reduced-order optimization methods that reduce the size of the control vector by projecting it onto a limited number of basis functions were suggested. In this paper, we show that such order reduction can indeed speed up the initial convergence rate of an assimilation effort in the eastern subtropical North Atlantic using in situ and satellite data as constraints. However, an improved performance of the optimization was only obtained with a hybrid approach where the optimization is started in a reduced subspace but is continued subsequently using the full control space. In such an experiment about 50% of the computational cost can be saved as compared to the optimization in the full control space. Although several order-reduction approaches seem feasible, the best result was obtained by projecting the control vector onto Empirical Orthogonal Functions (EOFs) computed from a set of adjusted control vectors estimated previously from an optimization using the same model configuration.     

Discharge experiments simulating chemical evolution on the surface of Titan

Saturn's moon Titan has been considered as one of the few places in our Solar System, where atmospheric and surface conditions could have produced organic compounds essential as precursors for an evolution of life. The Cassini-Huygens mission has provided new data on Titan's atmosphere and surface, which enabled us to simulate the chemical processes occurring under these conditions. Possible lightning events on Titan cannot only produce higher hydrocarbons, but also allow surface water ice to participate in the reaction scenario, resulting in CHO, CHN, and CHON compounds including several molecules relevant for the formation of amino acids and nucleic acids.     

FORSAT: a 3D forest monitoring system for cover mapping and volumetric 3D change detection

ABSTRACT

A 3D forest monitoring system, called FORSAT (a satellite very high resolution image processing platform for forest assessment), was developed for the extraction of 3D geometric forest information from very high resolution (VHR) satellite imagery and the automatic 3D change detection. FORSAT is composed of two complementary tasks: (1) the geometric and radiometric processing of satellite optical imagery and digital surface model (DSM) reconstruction by using a precise and robust image matching approach specially designed for VHR satellite imagery, (2) 3D surface comparison for change detection. It allows the users to import DSMs, align them using an advanced 3D surface matching approach and calculate the 3D differences and volume changes (together with precision values) between epochs. FORSAT is a single source and flexible forest information solution, allowing expert and non-expert remote sensing users to monitor forests in three and four (time) dimensions. The geometric resolution and thematic content of VHR optical imagery are sufficient for many forest information needs such as deforestation, clear-cut and fire severity mapping. The capacity and benefits of FORSAT, as a forest information system contributing to the sustainable forest management, have been tested and validated in case studies located in Austria, Switzerland and Spain.     

A study to explain the emission of nitric oxide from a marsh soil

In the period 18–21 September 1989, soil NO emission was studied at Halvergate Marshes, Norfolk (U.K.) within the framework of the BIATEX-LOVENOX joint field experiment. Using a dynamic chamber technique, 186 measurements at four plots were performed showing a net NO flux of 7.2–14.6×10^–12 kgN m^–2 s^–1. Soil samples from a soil profile (1.0 m) at a representative site and from the uppermost layer (0.1 m) of each of the four plots were sent to the laboratory for (a) detailed physical and chemical soil analysis, (b) determination of NO production rates, NO uptake rate constants, and NO compensation mixing ratios, and (c) characterization of the microbial processes involved. A diffusive model (Galbally and Johansson, 1989) was applied to the laboratory results to infer NO fluxes of the individual soil samples. When we compared these fluxes with those measured in the field, we found agreement within a factor 2–4. Furthermore, laboratory studies showed, that NO was produced and consumed only in the upper soil layer (0–0.1 m depth) and that the NO production and consumption activities observed in the Halvergate marsh soil were most probably due to the anaerobic metabolism of denitrifying bacteria operating in anaerobic microniches within the generally aerobic soil.     

Different age response of zircon and monazite during the tectono-metamorphic evolution of a high grade paragneiss from the Ruhla Crystalline Complex,central Germany

Results of TIMS, SIMS and SEM analyses show that zircon and monazite in a high-grade paragneiss of the Ruhla Crystalline Complex, central Germany, were formed and/or altered during different stages of a tectono-metamorphic history between Early Devonian and Permian times. Detrital zircon cores of >460 Ma place an older limit on the age of anatexis, and show that the paragneiss sequence contains rocks at least as young as early Cambrian. Metamorphic zircon growth commenced at ~365 Ma, peaking at ~360–355 Ma at the same time that granite dykes were emplaced. In contrast, monazite in the paragneiss preserves little record of the metamorphic peak. Most monazite grains grew or were recrystallised in the Lower Carboniferous at ~339 Ma, contemporaneous with the emplacement of voluminous diorite and granite bodies. These intrusions and related tectonics caused some of the high-U zircon overgrowths to undergo moderate to severe Pb loss. A second Pb loss event, between 300 and 280 Ma, can be related to Late Carboniferous/Early Permian large-scale block faulting.Editorial responsibility: J. Hoefs     

The formation of high-grade ignimbrites, I: Experiments on high- and low-concentration transport systems containing sticky particles

 High-grade ignimbrites are thought to be deposited by pyroclastic flows at temperatures exceeding minimum welding temperature or even solidus temperature. Corresponding pyroclastic-flow particles range from plastic to partially liquid and are able to aggregate or coalesce. This contrasts with particles in pyroclastic flows producing unwelded ignimbrite, which are capable of elastic grain interactions. The low aspect ratio and great areal extent of high-grade ignimbrites requires transport in a particulate state either by (a) high-concentration mass flow facilitated by fluidizing gas reducing internal friction, or by (b) expanded turbulent flow of low but downward increasing concentration. This paper presents experiments designed to investigate the effects of plastic to liquid particles on these two contrasting transport mechanisms. Gas fluidization experiments using polyethyleneglycole (PEG) powders heated above minimum sintering (T_ms) and melting (T_m) temperatures cover a wide range of fluidization velocities (U_mf>U_a>0.6·U_t) but are always in the bubbly fluidization regime similar to fluidized ignimbrite ash, where particle volume concentration outside the bubbles is high (≈10^–1). When the powders reach a critical temperature T_m≥T≥T_ms, defluidization by catastrophic particle aggregation immediately commences in both stationary and laterally moving fluidized beds as well as in experiments using mixtures of high- and low-T_m (≥30 wt.%) PEG powders, when T≥T_ms of the lower-T_m powder. This indicates that extended particulate transport at T≥T_ms is not possible at such high particle concentrations. In the turbulent flow experiments, liquid sprays of molten PEG or water, vertically injected into a high-Re (>10⁴) horizontal air flow, form a low-concentration (10^–5 to 10^–4) turbulent suspension current. Proximal formation of partially coalesced aggregates, which settle faster than individual particles, causes the measured downstream decay of sedimentation rate to be steeper than predicted by theory of single solid-particle sedimentation from turbulent suspensions. As particles become finer downstream and coalescence efficiency decreases in response to cooling, more distally formed aggregates become too small and rare to modify sedimentation-rate decay from that of suspension flows containing solid particles. The key difference between the two transport systems is particle concentration, C. Since particle collision rate R_coll∝C², collision rates in fluidized beds are so high that all particles immediately aggregate when coalescence efficiency (1≥E_coal≥0) is larger than 10^-3. Low-concentration suspensions, on the other hand, require much higher values of E_coal for significant aggregation to occur. Dilute pyroclastic flows will have higher particle volume fractions (≈10^–3) than the experimental currents, but then viscous pyroclasts should have lower coalescence efficiencies than PEG droplets. Experimental results thus support an expanded turbulent transport mechanism of pyroclastic flows generating extensive high-grade ignimbrite sheets. Received: 28 August 1996 / Accepted: 3 December 1997