Single frequency processing of Ørsted GPS radio occultation measurements

The Global Positioning System (GPS) radio occultation measurements obtained using the TurboRogue GPS receiver on the Danish satellite Ørsted have been processed using the single frequency method. Atmospheric profiles of refractivity and temperature are derived and validated against numerical weather prediction data from the European Centre for Medium-Range Weather Forecast (ECMWF). Results from the Ørsted GPS measurement campaign in February 2000 indicate that the single frequency method can provide retrievals with accuracy comparable to that of using two frequencies. From comparisons between measured dry temperature profiles and corresponding dry temperature profiles derived from ECMWF analysis fields, we find a mean difference of less than 0.5 K and a standard deviation of 2–4 K between 500 and 30 hPa in height. Above 30 hPa the impact of the ionosphere becomes more dominant and more difficult to eliminate using the single frequency method, and the results show degraded accuracy when compared to previous analysis results of occultation data from other missions using the dual frequency method. At latitudes less than 40° (denoted low latitudes), the standard deviation is generally smaller than at latitudes higher than 40° (denoted high latitudes). A small temperature bias is observed centered at 200 hPa for low latitudes and at 300 hPa for high latitudes. This indicates that the ECMWF analyses do not adequately resolve the tropopause temperature minimum. In the lowest part of the troposphere an observed warm bias is thought to be due to erroneous tracking of the GPS signal in cases of atmospheric multipath propagation.     

反演寿县TMI微波频率的地表比辐射率

利用一个微波辐射传输模式，晴天下的大气探测廓线，地面的部分观测资料以及卫星观测的亮温，计算了淮河流域能量和水循环实验(HUBEx)中，寿县地区的热带降雨观测卫星微波成像仪(TRMM／TMI)微波频率的地表比辐射率。通过比较所计算的微波地表比辐射率随地表状况的改变，发现地表微波比辐射率随地表状况的变化存在敏感性，并且其变化是合理的。本中，地表肤温设定等于地表空气温度，并且仅计算了寿县的地表比辐射率，在HUBEx区域上的地表比辐射率的计算需要更多的观测资料。     

The onset of Pleistocene glaciation in the Barents Sea: implications for glacial isostatic adjustment

In this paper the effect of a delayed onset of glaciation in the Barents Sea on glacial isostatic adjustment is investigated. The model calculations solve the sea-level equation governing the total mass redistributions associated with the last glaciation cycle on a spherically symmetric, linear, Maxwell viscoelastic earth for two different scenarios for the growth phase of the Barents Sea ice sheet. In the first ice model a linear growing history is used for the Barents Sea ice sheet, which closely relates its development to the build-up of other major Late Pleistocene ice sheets. In the second ice model the accumulation of the Barents Sea ice sheet is restricted to the last 6 ka prior to the last glacial maximum.
The calculations predict relative sea levels, present-day radial velocities, and gravity anomalies for the area formerly covered by the Weichselian ice sheet. The results show that observed relative sea levels in the Barents Sea are appropriate for distinguishing between the different glaciation histories. In particular, present-day observables such as the free-air gravity anomaly over the Barents Sea, and the present-day radial velocities are sensitive to changes in the glaciation history on this scale.
A palaeobathymetry derived from relative sea-level predictions before the last glacial maximum based on the second ice model essentially agrees with a palaeobathymetry derived by Lambeck (1995). The additional emerged areas provide centres for the build-up of an ice sheet and thus support the theory of Hald, Danielsen & Lorentzen (1990) and Mangerud et al. (1992) that the Barents Sea was an essentially marine environment shortly before the last glacial maximum.     

Life cycle of Agapetus fuscipes (Trichoptera, Glossosomatidae) in a first-order upland stream in central Germany

The number of immature stages and the seasonal patterns of development are basic life history features of a stream dwelling species and knowledge about these important components are essential for understanding its adaptations to its dynamic environment. The life cycle of Agapetus fuscipes (Trichoptera, Glossosomatidae), one of the dominant scrapers in the upper and middle reaches of the Breitenbach, a first-order upland stream in central Germany, was analysed. The pronotum length and the relationship between pronotum length, larval biomass and case length showed seven distinct larval instars, contrary to earlier findings from the Breitenbach. In addition to a few trichopteran species from other functional feeding groups, A. fuscipes is the only scraping caddis fly reported to have more than five larval instars. The moult increments of pronotum length and larval biomass were distinctly lower than in glossosomatid species with five larval instars. A. fuscipes is clearly univoltine in the Breitenbach. First-instar larvae were found from July to the beginning of December, and second-instar larvae from July to January. At the beginning of December the population consisted of the instars I to V, and development did not cease during winter. The sixth-instar larvae occurred mostly in January, and the seventh-instar larvae were never present before January. The prepupae and pupae occurred in April. The last pupae were found at the beginning of September, although most of the emergence took place in June and July. At least five different immature stages with different ecological demands were present at any time throughout the year. The ecological advantage having two additional larval instars compared to other glossosomatid species may be to compensate for the high rate of mouthpart wear that occurs while the larvae feed on the rough Bunter Sandstone substratum. A further advantage may be to spread the risk of high mortality under unfavourable environmental conditions.     

Discrete element model for crack propagation in brittle materials

We propose a discrete element model for brittle rupture. The material consists of a bidimensional set of closed‐packed particles in contact. We explore the isotropic elastic behavior of this regular structure to derive a rupture criterion compatible to continuum mechanics. We introduce a classical criterion of mixed mode crack propagation based on the value of the stress intensity factors, obtained by the analysis of two adjacent contacts near a crack tip. Hence, the toughness becomes a direct parameter of the model, without any calibration procedure. We verify the consistency of the formulation as well as its convergence by comparison with theoretical solutions of tensile cracks, a pre‐cracked beam, and an inclined crack under biaxial stress. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical modeling of fracking fluid migration through fault zones and fractures in the North German Basin

Gas production from shale formations by hydraulic fracturing has raised concerns about the effects on the quality of fresh groundwater. The migration of injected fracking fluids towards the surface was investigated in the North German Basin, based on the known standard lithology. This included cases with natural preferential pathways such as permeable fault zones and fracture networks. Conservative assumptions were applied in the simulation of flow and mass transport triggered by a high pressure boundary of up to 50 MPa excess pressure. The results show no significant fluid migration for a case with undisturbed cap rocks and a maximum of 41 m vertical transport within a permeable fault zone during the pressurization. Open fractures, if present, strongly control the flow field and migration; here vertical transport of fracking fluids reaches up to 200 m during hydraulic fracturing simulation. Long-term transport of the injected water was simulated for 300 years. The fracking fluid rises vertically within the fault zone up to 485 m due to buoyancy. Progressively, it is transported horizontally into sandstone layers, following the natural groundwater flow direction. In the long-term, the injected fluids are diluted to minor concentrations. Despite the presence of permeable pathways, the injected fracking fluids in the reported model did not reach near-surface aquifers, either during the hydraulic fracturing or in the long term. Therefore, the probability of impacts on shallow groundwater by the rise of fracking fluids from a deep shale-gas formation through the geological underground to the surface is small.     

Regional geopotential model improvement for the Iranian geoid determination

Spherical harmonic expansions of the geopotential are frequently used for modelling the earth’s gravity field. Degree and order of recently available models go up to 360, corresponding to a resolution of about50 km. Thus, the high degree potential coefficients can be verified nowadays even by locally distributed sets of terrestrial gravity anomalies. These verifications are important when combining the short wavelength model impact, e.g. for regional geoid determinations by means of collocation solutions. A method based on integral formulae is presented, enabling the improvement of geopotential models with respect to non-global distributed gravity anomalies. To illustrate the foregoing, geoid computations are carried out for the area of Iran, introducing theGPM2 geopotential model in combination with available regional gravity data. The accuracy of the geoid determination is estimated from a comparison with Doppler and levelling data to ±1.4m.     

Sr-gorceixite,a weathering product in rich iron ores from the Córrego do Feijão Mine,Minas Gerais,Brazil

A sulfate-bearing Sr-Ba phosphate (Sr-gorceixite) formed by supergene alteration in rich iron ores from Minas Gerais, Brazil, has been studied by means of microprobe techniques. Analytical data and the structural formula are presented.