Uncertainties in simulating regional climate of Southern Africa: sensitivity to physical parameterizations using WRF

This study aims at quantifying seasonal biases of regional climate model outputs during southern African summer, against a dense in situ measurement network (daily rain-gauge and surface air temperature records, and 12?h UTC radiosondes), and uncertainties associated with some physical parameterizations. Using the non-hydrostatic Advanced Research Weather Forecast (WRF) laterally forced by ERA40 reanalysis, twenty-seven experiments configured with three schemes of cumulus (CU), planetary boundary layer (PBL) and microphysics (MP), are performed at 35?km horizontal resolution during the core of a summer rainy season (December 1993 to February 1994 season) representative of the South African rainfall climatology. WRF simulates accurately seasonal large-scale rainfall patterns, as well as seasonal gradients of South African rainfall and 2-m temperature, and seasonal vertical profiles of the air temperature and humidity. However seasonal biases fluctuate strongly from an experiment to another, denoting considerable uncertainties generated by the physical package. Rainfall amounts are the most sensitive parameter to the tested schemes. Their geography, intensity, and intraseasonal characteristics are predominantly sensitive to CU schemes, and much less to PBL and MP schemes. Some CU-PBL combinations produce additive effects, which can dramatically either reduce or increase biases. Satisfactory configurations are found for South African climate, which would not have been possible without testing numerous physical parameterizations.     

High Resolution DEM Generation in High‐Alpine Terrain Using Airborne Remote Sensing Techniques

Up‐to‐date and accurate digital elevation models (DEMs) are essential for many applications such as numerical modeling of mass movements or mapping of terrain changes. Today the Federal Department of Topography, swisstopo, provides Digital Terrain Models (DTMs) and Digital Surface Models (DSMs) derived from airborne LiDAR data with a high spatial resolution of 2 m covering the entire area of Switzerland below an elevation of 2000 m a.s.l.. However, above an elevation of 2000 m a.s.l., which is typical for high‐alpine terrain, the best product available is the a DTM with a spatial resolution of 25 m. This spatial resolution is insufficient for many applications in complex terrain. In this study, we investigate the quality of DSMs derived from opto‐electronic scanner data (ADS80; acquired in autumn 2010) using photogrammetric image correlation techniques based on the multispectral nadir and backward looking sensor data. As reference, we take a high precision airborne LiDAR data set with a spatial resolution of ca. 0.5 m, acquired in late summer 2010, covering the Grabengufer/Dorfbach catchment near Randa, VS. We find the deviations between the two datasets are surprisingly low. In terrain with inclination angles of less than 30° the RMSE is below 0.5 m. In extremely steep terrain of more than 50° the RMSE goes up to 2 m and outliers increase significantly. We also find dependencies of the deviations on illumination conditions and ground cover classes. Finally we discuss advantages and disadvantages of the different data acquisition methods.     

Thermo‐tectonic evolution of the south‐central Pyrenees from rifting to orogeny: insights from detrital zircon U/Pb and (U‐Th)/He thermochronometry

Constraining the thermal and denudational evolution of continental margins from extensional episodes to early orogenic stages is critical in the objective to better understand the sediment routing during the growth of orogenic topography. Here, we report 160 detrital zircon U/Pb ages and 73 (U‐Th)/He ages from Albian, Upper Cretaceous and Eocene sandstones from the south‐central Pyrenees. All samples show dominant zircon U/Pb age peaks at 310–320 Ma, indicating a primary contribution from Variscan granites of the central Pyrenean Axial Zone. A secondary population at 450–600 Ma documents zircon grains sourced from the eastern Pyrenees. Zircon (U‐Th)/He ages recovered from older samples document, a Triassic age peak at ca. 241 Ma, corresponding to denudation coeval with the initiation of Atlantic rifting. An Early Cretaceous cooling event at ca. 133 Ma appears consistent with rift‐related exhumation and thermal overprint on the Iberian margin. The (U‐Th)/He age peaks from ca. 80 Ma to ca. 68 Ma with decreasing depositional ages are interpreted to reflect the southward‐migrating thrust‐related exhumation on the pro‐wedge side of the Pyrenean orogen. The increase in lag times, from ca. 15 Ma in the Tremp Formation (ca. 65 Ma) to 28 Ma in the Escanilla Formation (ca. 40 Ma), suggests decreasing exhumation rates from 0.4 km Myr^–1 to 0.2 km Myr^–1. The apparent inconsistency with convergence rates is used to infer that rocks cooled at 68 Ma may have resided in the crust before final exhumation to the surface. Finally, the cooling event observed at 68 Ma provides support to the inferred acceleration of convergence, shortening and exhumation during Late Cretaceous times.     

Steerable wavelet analysis of CMB structures alignment

This paper reviews the application of a novel methodology for analysing the isotropy of the universe by probing the alignment of local structures in the CMB. The strength of the proposed methodology relies on the steerable wavelet filtering of the CMB signal. One the one hand, the filter steerability renders the computation of the local orientation of the CMB features affordable in terms of computation time. On the other hand, the scale-space nature of the wavelet filtering allows to explore the alignment of the local structures at different scales, probing possible different phenomena. We present the WMAP first-year data analysis recently performed by the same authors (Wiaux et al.), where an extremely significant anisotropy was found. In particular, a preferred plane was detected, having a normal direction with a northern end position at (θ, ) = (34°, 331°), close to the northern end of the CMB dipole axis. In addition, a most preferred direction was found in that plane, with a northern end direction at (θ, ) = (71°, 91°), very close to the north ecliptic pole. This result synthesised for the first time previously reported anomalies identified in the direction of the dipole and the ecliptic poles axes. In a forthcoming paper (Vielva et al.), we have extended our analysis to the study of individual frequency maps finding first indications for discarding foregrounds as the origin of the anomaly. We have also tested that the preferred orientations are defined by structures homogeneously distributed in the sky, rather than from localised regions. We have also analysed the WMAP 3-year data, finding the same anomaly pattern, although at a slightly lower significance level.     

Water temperature modelling in a controlled environment: comparative study of heat budget equations

One of the challenges when modelling a complex variable such as water temperature in rivers is that it can be difficult to determine the sources of error and to ensure that the simulations are truly representative of the reality. Therefore, a heat budget study was completed in a controlled environment, which excluded advection and bottom fluxes but enabled observation of all the other fluxes. A 21.42 m³ pool was installed and insulated to limit heat exchange through the sides and bottom. All the major energy fluxes were monitored for a 50‐day period. Different equations for individual heat budget terms were tested to determine their ability to reproduce the observations. This experiment also permitted to assess the relative importance of each component of the heat budget. Performance of each semi‐empirical equation was determined by comparing predictions and measured values. It was thus possible to choose the formulae that best represented the measured heat exchange processes, while understanding the limits of some of the semi‐empirical representations of heat exchange processes. The results highlight the importance of radiative terms into the heat budget because they controlled the major sources and sinks. The study also showed the importance of the wind function determination into the calculation of latent heat flux. The resulting water temperature model returned simulated hourly water temperature with an overall root mean square error of 0.71 °C/h and a modified Nash–Sutcliffe coefficient of 0.97. Copyright © 2012 John Wiley & Sons, Ltd.     

Removing attitude-related variations in the line-of-sight for kinematic GPS positioning

We are presenting a method for removing attitude-related variations (ATTRVs) in dynamic 1 Hz GPS positioning. The ATTRVs are separated into vertical and horizontal components. These result from the translational and rotational motions of the vehicle, which is a marine research vessel in our case. We have developed new observation equations that use corrected pseudoranges and carrier phases to account for ATTRVs. In the present contribution, we are only focusing on the vertical signals. These modeled ATTRVs are included as corrections in the line-of-sight (LoS) to each GPS satellite in the RINEX data sets using a new software called RNXATTCOR. Precise IGS sp3-orbits are used as inputs together with precise lever arm coordinates of the onboard GPS antennas, observations from the marine inertial navigation system and a priori 3D position of the vessel. The corrected RINEX data sets are then processed using kinematic processing or sequential processing in Precise Point Positioning (PPP) mode. The reduction of the standard deviation from a running mean in the final processed GPS time series is about 95%. The current method is being proposed for marine geodesy science applications such as along-track calibration/validation of altimetry products, and also for applications related to offshore kinematic high precision GPS positioning such as drilling, offshore platforms stability, marine pipeline positioning, earthquake monitoring and tsunami detection.     

CBM and CO2-ECBM related sorption processes in coal: A review

This article reviews the state of research on sorption of gases (CO₂, CH₄) and water on coal for primary recovery of coalbed methane (CBM), secondary recovery by an enhancement with carbon dioxide injection (CO₂-ECBM), and for permanent storage of CO₂ in coal seams.Especially in the last decade a large amount of data has been published characterizing coals from various coal basins world-wide for their gas sorption capacity. This research was either related to commercial CBM production or to the usage of coal seams as a permanent sink for anthropogenic CO₂ emissions. Presently, producing methane from coal beds is an attractive option and operations are under way or planned in many coal basins around the globe. Gas-in-place determinations using canister desorption tests and CH₄ isotherms are performed routinely and have provided large datasets for correlating gas transport and sorption properties with coal characteristic parameters.Publicly funded research projects have produced large datasets on the interaction of CO₂ with coals. The determination of sorption isotherms, sorption capacities and rates has meanwhile become a standard approach.In this study we discuss and compare the manometric, volumetric and gravimetric methods for recording sorption isotherms and provide an uncertainty analysis. Using published datasets and theoretical considerations, water sorption is discussed in detail as an important mechanisms controlling gas sorption on coal. Most sorption isotherms are still recorded for dry coals, which usually do not represent in-seam conditions, and water present in the coal has a significant control on CBM gas contents and CO₂ storage potential. This section is followed by considerations of the interdependence of sorption capacity and coal properties like coal rank, maceral composition or ash content. For assessment of the most suitable coal rank for CO₂ storage data on the CO₂/CH₄ sorption ratio data have been collected and compared with coal rank.Finally, we discuss sorption rates and gas diffusion in the coal matrix as well as the different unipore or bidisperse models used for describing these processes.This review does not include information on low-pressure sorption measurements (BET approach) to characterize pore sizes or pore volume since this would be a review of its own. We also do not consider sorption of gas mixtures since the data base is still limited and measurement techniques are associated with large uncertainties.     

Solubility of CO2 in a Ca-rich leucitite: effects of pressure,temperature, and oxygen fugacity

The solubility of carbon dioxide in a Ca-rich leucitite has been investigated as a function of pressure (0.1–2.0 GPa), temperature (1200–1600°C), and oxygen fugacity. The experiments were done in a rapid-quench internally-heated pressure vessel (0.1 GPa) and a piston cylinder (0.5–2.0 GPa). The leucitite glass, previously equilibrated at NNO, and silver oxalate were loaded in Fe-doped Pt capsules (oxidized conditions) and graphitelined Pt capsules (reduced conditions). Secondary Ion Mass Spectrometry and bulk carbon analyses were used to determine the amount of dissolved carbon. Speciation of carbon was characterized by Fourier transform microinfrared spectroscopy. At oxidized conditions, only CO₃ ^2- is observed as a dissolved species. The solubility is high with CO₂ contents in the melt attaining 6.2 wt% at 2.0 GPa and 1350°C. The solubility increases with pressure and shows a significant negative temperature dependence. An excellent correlation is obtained when the data are fit to a model, based on the simplified solubility reaction CO₂ ^(vapor)+O^2-(melt)CO₃ ^2-(melt), which describes the solubility of CO₂ as a function of pressure, temperature and fCO₂. At reduced conditions, the amount of carbon dissolved is significantly lower, and CO₃ ^2- is still the only species present in the melt. If the solubility model established at oxidized conditions is applied, the carbon dissolved appears to be essentially a function of fCO₂ alone although divergence increases in a consistent manner with pressure and temperature. This could suggest a low but significant solubility of CO with a positive temperature dependence or a departure of the calculated fluid compositions determined by the equation of state from the actual ones. The strong preferential solubility of carbon in its oxidized C⁴⁺ form, even at reduced conditions, implies that ascending melts with high CO₂ solubility can experience significant oxidation through degassing. This could reconcile the oxidized nature of some Ca-rich alkaline magmas with more reduced mantle source regions.     

Distribution of222Rn over the north Pacific: Implications for continental influences

The atmospheric distribution of²²²Rn over the north Pacific is simulated with a three-dimensional chemical tracer model using meteorological input from the NASA-GISS general circulation model (4°×5° resolution). Radon-222 (half-life 3.8 days) is a tracer of continental air. Model results are in good agreement with measurements from ships and aircraft. Strong Asian influence is found throughout the tropospheric column over the north Pacific in spring, reflecting a combination of frequent convection over the continent, strong westerly winds at altitude, and subsidence over the ocean. In summer, the upper troposphere over the north Pacific is heavily affected by deep convection over China; however, Asian influences at the surface are then at their yearly minimum. In winter, strong Asian influence is found near the surface but not at high altitudes. Transport of American air over the Pacific is important only at low latitudes. American sources account for 11% of total²²²Rn in the model at Midway, 30% at Mauna Loa and 59% at Oahu. Results for Hawaii indicate two seasonal peaks of American influence, one in summer and one in winter. The tropical western Pacific is particularly remote from continental influences year round.