Evaluation of the Interpretation of Ceilometer Data with RASS and Radiosonde Data

Since 2006 different remote monitoring methods for determining mixing-layer height have been operated in parallel in Augsburg (Germany). One method is based on the operation of eye-safe commercial mini-lidar systems (ceilometers). The optical backscatter intensities recorded with ceilometers provide information about the range-dependent aerosol concentration; gradient minima within this profile mark the tops of mixed layers. Special software for these ceilometers provides routine retrievals of lower atmospheric layering. A second method, based on sodar observations, detects the height of a turbulent layer characterized by high acoustic backscatter intensities due to thermal fluctuations and a high variance of the vertical velocity component. This information is extended by measurements with a radio-acoustic sounding system (RASS) that directly provides the vertical temperature profile from the detection of acoustic signal propagation and thus temperature inversions that mark atmospheric layers. Ceilometer backscatter information is evaluated by comparison with parallel measurements. Data are presented from 2 years of combined ceilometer and RASS measurements at the same site and from comparison with a nearby (60 km) radiosonde for larger-scale humidity information. This evaluation is designed to ensure mixing-layer height monitoring from ceilometer data more reliable.     

A gridded monthly upper-air data set from 1918 to 1957

Significant efforts have been devoted in recent years towards extending observation-based three-dimensional atmospheric data sets back in time. Such data sets form an important basis for a better understanding of the climate system. Here we present a new monthly three-dimensional global data set that is based on historical upper-air data and surface data. We use statistical reconstruction techniques, calibrated using ERA-40 data, to obtain gridded data from the numerous, but scattered and heterogeneous historical upper-air observations. In contrast to previous work, in which we used hemispheric principal components on both the predictor and the predictand side to reconstruct spatially complete fields back to 1880, we restrict the procedure to places and times where upper-air observations are available. Each grid column (consisting of four variables at six levels) is then reconstructed independently using only predictor variables in the vicinity (i.e., only local stationarity is required rather than stationary large-scale patterns). The product, termed REC2, is a gridded, global monthly data set of geopotential height, temperature, and u and v wind from 850 to 100?hPa back to 1918. The data set is sparse (i.e., many grid cells are empty), but provides an alternative to large-scale reconstructions as it allows for non-stationary teleconnections. We show the results of several validation experiments, compare our new data set with a number of existing data sets, and demonstrate that it is suitable for analysing large-scale climate variability on interannual time-scales.     

The effects of boreal forest expansion on the summer Arctic frontal zone

Over the last 100?years, Arctic warming has resulted in a longer growing season in boreal and tundra ecosystems. This has contributed to a slow northward expansion of the boreal forest and a decrease in the surface albedo. Corresponding changes to the surface and atmospheric energy budgets have contributed to a broad region of warming over areas of boreal forest expansion. In addition, mesoscale and synoptic scale patterns have changed as a result of the excess energy at and near the surface. Previous studies have identified a relationship between the positioning of the boreal forest-tundra ecotone and the Arctic frontal zone in summer. This study examines the climate response to hypothetical boreal forest expansion and its influence on the summer Arctic frontal zone. Using the Weather Research and Forecasting model over the Northern Hemisphere, an experiment was performed to evaluate the atmospheric response to expansion of evergreen and deciduous boreal needleleaf forests into open shrubland along the northern boundary of the existing forest. Results show that the lower surface albedo with forest expansion leads to a local increase in net radiation and an average hemispheric warming of 0.6°C at and near the surface during June with some locations warming by 1–2°C. This warming contributes to changes in the meridional temperature gradient that enhances the Arctic frontal zone and strengthens the summertime jet. This experiment suggests that continued Northern Hemisphere high-latitude warming and boreal forest expansion might contribute to additional climate changes during the summer.     

The corridor correction concept

Atmospheric delays are contributors to the GNSS error budget in precise GNSS positioning that can reduce positioning accuracy considerably if not compensated appropriately. Both ionospheric and tropospheric delay corrections can be determined with help of reference stations in active GNSS networks. One approach to interpolate these error terms to the user’s location that is employed in Germany’s SAPOS network is the determination of area correction parameters (ACP, German: “Fl?chenkorrekturparameter—FKP”). A 2D interpolation scheme using data from at least 3 reference stations surrounding the rover is employed. A modification of this method was developed which only makes use of as few as 2 reference stations and provides 1D linear correction parameters along a “corridor” in which the user’s rover is moving. We present the results of a feasibility study portraying results from use of corridor correction parameters for precise RTK-like positioning. The differences to the reference coordinates (3D) attained in average for 1 h of data employing selected network nodes in Germany are between 0.8 and 2.0 cm, which compares well with the traditional area correction method that yields an error of 0.7 up to 1.1 cm.     

Primitive andesites from the Taupo Volcanic Zone formed by magma mixing

Andesites with Mg# >45 erupted at subduction zones form either by partial melting of metasomatized mantle or by mixing and assimilation processes during melt ascent. Primitive whole rock basaltic andesites from the Pukeonake vent in the Tongariro Volcanic Centre in New Zealand’s Taupo Volcanic Zone contain olivine, clino- and orthopyroxene, and plagioclase xeno- and antecrysts in a partly glassy matrix. Glass pools interstitial between minerals and glass inclusions in clinopyroxene, orthopyroxene and plagioclase as well as matrix glasses are rhyolitic to dacitic indicating that the melts were more evolved than their andesitic bulk host rock analyses indicate. Olivine xenocrysts have high Fo contents up to 94%, δ¹⁸O_(SMOW) of +5.1‰, and contain Cr-spinel inclusions, all of which imply an origin in equilibrium with primitive mantle-derived melts. Mineral zoning in olivine, clinopyroxene and plagioclase suggest that fractional crystallization occurred. Elevated O isotope ratios in clinopyroxene and glass indicate that the lavas assimilated sedimentary rocks during stagnation in the crust. Thus, the Pukeonake andesites formed by a combination of fractional crystallization, assimilation of crustal rocks, and mixing of dacite liquid with mantle-derived minerals in a complex crustal magma system. The disequilibrium textures and O isotope compositions of the minerals indicate mixing processes on timescales of less than a year prior to eruption. Similar processes may occur in other subduction zones and require careful study of the lavas to determine the origin of andesite magmas in arc volcanoes situated on continental crust.     

Turbidite stacking patterns in salt‐controlled minibasins: Insights from integrated analogue models and numerical fluid flow simulations

The sea floor of intraslope minibasins on passive continental margins plays a significant role in controlling turbidity current pathways and the resulting sediment distribution. To address this, laboratory analogue modelling of intraslope minibasin formation is combined with numerical flow simulations of multi‐event turbidity currents. This approach permits an improved understanding of evolving flow–bathymetry–deposit interactions and the resulting internal stacking patterns of the infills of such minibasins. The bathymetry includes a shelf to slope channel followed by an upper minibasin, which are separated by a confining ridge from two lower minibasins that compares well with analogous bathymetries reported from natural settings. From a wider range of numerical flow experiments, a series of 100 consecutive flows is reported in detail. The turbidity currents are released into the channel and upon reaching the upper minibasin follow a series of stages from short initial ponding, ‘filling and spilling’ and an extended transition to long retrogradational ponding. Upon reaching the upper minibasin floor, the currents undergo a hydraulic jump and therefore much sediment is deposited in the central part of the minibasin and the counterslope. This modifies the bathymetry such that in the fill and spill stage, flow stripping and grain‐size partitioning cause some finer sediment to be transported across the confining ridge into the lower minibasins. Throughout the basin infill process, the sequences retrograde upstream, accompanied by lateral switching into locally formed depressions in the upper minibasin. After the fill and spill stage, significant deposition occurs in the channel where retrograding cyclic steps with wavelengths of 1 to 2 km develop as a function of pulsating flow criticality. These results are at variance with conventional schemes that emphasize sequential downstream minibasin filling through ponding dominated by vertical aggradation. Comparison of these results with published field and experimental examples provides support for the main conclusions.     

Continental rift-setting and evolution of Neoproterozoic Sindreth Basin in NW-India

The Neoproterozoic Sindreth Basin, NW India, and its surrounding area represent a half graben structure situated between the undeformed Malani Igneous Suite (MIS) in the west and a corridor of coeval Cryogenian ductile deformation, anatexis and granite intrusion in the east. The main lithologies observed in the basin are conglomerate, fanglomerate, debris flow and lake deposits derived from a nearby continental provenance, intercalated with concurrent mafic and felsic lava flows. Based on geological traverses across the strike of the basin, we propose a three-fold classification comprising Lower Clastic Unit and an Upper Clastic Unit and a Bimodal (basalt–rhyolite) Volcanic Unit separating the two. Tilting due to basin inversion and faulting has been observed; however, the rocks are unmetamorphosed and show undisturbed primary sedimentary features. The stratigraphic record of the basin is characteristic for deposition and magmatism in a fault-related continental setting. Implications of the findings have been discussed in the context of Neoproterozoic crustal dynamics in NW India. This study provides conclusive evidence for a continental setting for Sindreth Basin evolution and contests the recent models of active subduction setting (either back-arc basin or accretionary sediments over a subduction zone).     

Novel basin modelling concept for simulating deformation from mechanical compaction using level sets

As sedimentation progresses in the formation and evolution of a depositional geologic basin, the rock strata are subject to various stresses. With increasing lithostatic pressure, compressional forces act to compact the porous rock matrix, leading to overpressure buildup, changes in the fluid pore pressure and fluid flow. In the context of petroleum systems modelling, the present study concerns the geometry changes that a compacting basin experiences subject to deposition. The purpose is to track the positions of the rock layer interfaces as compaction occurs. To handle the challenge of potentially large geometry deformations, a new modelling concept is proposed that couples the pore pressure equation with a level set method to determine the movement of lithostratigraphic interfaces. The level set method propagates an interface according to a prescribed speed. The coupling term for the pore pressure and level-set equations consists of this speed function, which is dependent on the compaction law. The two primary features of this approach are the simplicity of the grid and the flexibility of the speed function. A first evaluation of the model concept is presented based on an implementation for one spatial dimension accounting for vertical effective stress. Isothermal conditions with a constant fluid density and viscosity were assumed. The accuracy of the implemented numerical solution for the case of a single stratigraphic unit with a linear compaction law was compared to the available analytical solution [38]. The multi-layer setup and the nonlinear case were tested for plausibility.     

A comparison of two physics-based numerical models for simulating surface water–groundwater interactions

Problems in hydrology and water management that involve both surface water and groundwater are best addressed with simulation models that can represent the interactions between these two flow regimes. In the current generation of coupled models, a variety of approaches is used to resolve surface–subsurface interactions and other key processes such as surface flow propagation. In this study we compare two physics-based numerical models that use a 3D Richards equation representation of subsurface flow. In one model, surface flow is represented by a fully 2D kinematic approximation to the Saint–Venant equations with a sheet flow conceptualization. In the second model, surface routing is performed via a quasi-2D diffusive formulation and surface runoff follows a rill flow conceptualization. The coupling between the land surface and the subsurface is handled via an explicit exchange term resolved by continuity principles in the first model (a fully-coupled approach) and by special treatment of atmospheric boundary conditions in the second (a sequential approach). Despite the significant differences in formulation between the two models, we found them to be in good agreement for the simulation experiments conducted. In these numerical tests, on a sloping plane and a tilted V-catchment, we examined saturation excess and infiltration excess runoff production under homogeneous and heterogeneous conditions, the dynamics of the return flow process, the differences in hydrologic response under rill flow and sheet flow parameterizations, and the effects of factors such as grid discretization, time step size, and slope angle. Low sensitivity to vertical discretization and time step size was found for the two models under saturation excess and homogeneous conditions. Larger sensitivity and differences in response were observed under infiltration excess and heterogeneous conditions, due to the different coupling approaches and spatial discretization schemes used in the two models. For these cases, the sensitivity to vertical and temporal resolution was greatest for processes such as reinfiltration and ponding, although the differences between the hydrographs of the two models decreased as mesh and step size were progressively refined. In return flow behavior, the models are in general agreement, with the largest discrepancies, during the recession phase, attributable to the different parameterizations of diffusion in the surface water propagation schemes. Our results also show that under equivalent parameterizations, the rill and sheet flow conceptualizations used in the two models produce very similar responses in terms of hydrograph shape and flow depth distribution.