Performance assessment of three convective parameterization schemes in WRF for downscaling summer rainfall over South Africa

Austral summer rainfall over the period 1991/1992 to 2010/2011 was dynamically downscaled by the weather research and forecasting (WRF) model at 9 km resolution for South Africa. Lateral boundary conditions for WRF were provided from the European Centre for medium-range weather (ECMWF) reanalysis (ERA) interim data. The model biases for the rainfall were evaluated over the South Africa as a whole and its nine provinces separately by employing three different convective parameterization schemes, namely the (1) Kain–Fritsch (KF), (2) Betts–Miller–Janjic (BMJ) and (3) Grell–Devenyi ensemble (GDE) schemes. All three schemes have generated positive rainfall biases over South Africa, with the KF scheme producing the largest biases and mean absolute errors. Only the BMJ scheme could reproduce the intensity of rainfall anomalies, and also exhibited the highest correlation with observed interannual summer rainfall variability. In the KF scheme, a significantly high amount of moisture was transported from the tropics into South Africa. The vertical thermodynamic profiles show that the KF scheme has caused low level moisture convergence, due to the highly unstable atmosphere, and hence contributed to the widespread positive biases of rainfall. The negative bias in moisture, along with a stable atmosphere and negative biases of vertical velocity simulated by the GDE scheme resulted in negative rainfall biases, especially over the Limpopo Province. In terms of rain rate, the KF scheme generated the lowest number of low rain rates and the maximum number of moderate to high rain rates associated with more convective unstable environment. KF and GDE schemes overestimated the convective rain and underestimated the stratiform rain. However, the simulated convective and stratiform rain with BMJ scheme is in more agreement with the observations. This study also documents the performance of regional model in downscaling the large scale climate mode such as El Niño Southern Oscillation (ENSO) and subtropical dipole modes. The correlations between the simulated area averaged rainfalls over South Africa and Nino3.4 index were ?0.66, ?0.69 and ?0.49 with KF, BMJ and GDE scheme respectively as compared to the observed correlation of ?0.57. The model could reproduce the observed ENSO-South Africa rainfall relationship and could successfully simulate three wet (dry) years that are associated with La Niña (El Niño) and the BMJ scheme is closest to the observed variability. Also, the model showed good skill in simulating the excess rainfall over South Africa that is associated with positive subtropical Indian Ocean Dipole for the DJF season 2005/2006.     

A progress report on the design of the 8 m aperture Optical/Infrared Telescope

Progress of the design study carried out y the Instituto de Astrophysica de Canarias and Royal Greenwich Observatory joint project office for an 8 m aperture Optical/Infrared Telescope. This includes the development of an active pneumatic support system for both the radial and axial supports of an 8 m diameter thin meniscus mirror and the investigations of new technology encoding systems.Paper presented at the Symposium on the JNLT and Related Engineering Developments, Tokyo, November 29–December 2, 1988.     

Boreal temperature variability inferred from maximum latewood density and tree-ring width data, Wrangell Mountain region, Alaska

Variations in both width and density of annual rings from a network of tree chronologies were used to develop high-resolution proxies to extend the climate record in the Wrangell Mountain region of Alaska. We developed a warm-season (July–September) temperature reconstruction that spans A.D. 1593–1992 based on the first eigenvector from principal component analysis of six maximum latewood density (MXD) chronologies. The climate/tree-growth model accounts for 51% of the temperature variance from 1958 to 1992 and shows cold in the late 1600s–early 1700s followed by a warmer period, cooling in the late 1700s–early 1800s, and warming in the 20th century. The 20th century is the warmest of the past four centuries. Several severely cold warm-seasons coincide with major volcanic eruptions. The first eigenvector from a ring-width (RW) network, based on nine chronologies from the Wrangell Mountain region (A.D. 1550–1970), is correlated positively with both reconstructed and recorded Northern Hemisphere temperatures. RW shows a temporal history similar to that of MXD by increased growth (warmer) and decreased growth (cooler) intervals and trends. After around 1970 the RW series show a decrease in growth, while station data show continued warming, which may be related to increasing moisture stress or other factors. Both the temperature history based on MXD and the growth trends from the RW series are consistent with well-dated glacier fluctuations in the Wrangell Mountains and some of the temperature variations also correspond to variations in solar activity.     

Geological Relationships Between the Intrusions,Country Rocks,and Ni-Cu-PGE Sulfides of the Kharaelakh Intrusion,Noril′sk Region:Implications for the Roles of Sulfide Differentiation and Metasomatism in Their Genesis

The Ni-Cu-platinum group element sulfide ore deposits of the Kharaelakh Intrusion,Noril′sk Region,Siberia,represent a large concentration of sulfides associated with a small differentiated intrusion formed at the edge of the Siberian Craton in the roots of the Siberian Trap flood basalt.The deposit is associated with an intrusion that occupies a flanking periclinal structure adjacent to the Noril′sk-Kharaelakh Fault.The intrusion is strongly differentiated and comprises taxitic gabbrodolerites,picritic gabbrodolerites,and gabbrodolerites within the main body which in turn forms a chonolith within a sheet-like intrusion that extends laterally to form extensive undifferentiated sills of gabbrodolerite.The intrusion substantially replaces the stratigraphy of the country rocks,and although it appears to have exploited the axis of structures developed in response to transtension,the intrusion has created space by both mechanical dilation of stratigraphy and magmatic replacement of pre-existing sedimentary rocks.The frontal lobes of the main intrusion have complex apophyses of gabrodolerite on a range of scales that demonstrate replacement of the sedimentary rocks and link to the development of an extensive metamorphic halo in the country rocks.This halo is much narrower over the main body of the intrusion,and these observations have implications for the thermal history of the intrusion.Mg-skarns and breccias are developed in the roof of the main body of the intrusion.Within the intrusion,the taxitic rocks contain vesicles and the blebby sulfides developed in the picritic and taxitic gabbrodolerites appear to have a linkage to volatile phases.Cuprous sulfide mineralization developed at the roof of the Kharaelakh Intrusion is associated with metamorphosed and skarn-bearing country rocks,and appears to have been generated by a combination of sulfide fractionation and associated metasomatism.The geological relationships appear consistent with a chonolith model for the development of the differentiated intrusion and mineralization,but the extent of metasmorphism of the country rocks appears to be related to the unusual thickness of gabbrodolerite apophyses at the flanks of the intrusion rather than metamorphism produced by the passage of mafic magma through the intrusion.Variations in disseminated sulfide compositions and metasomatic textures in the skarns are described,and a model is proposed which balances traditional views on the evolution of the magma conduits with the impact of magmatic fluids transported through the magma column(i.e.transmagmatic fluids).The importance of structures in controlling the nature of the conduit,and the resultant small intrusions with excess sulfide is a feature of many other Ni-Cu sulfide deposits including Voisey′s Bay,and it is suggested that the sulfides are more likely to have beentransported from depth into their final resting place rather than developed by in-situ equilibration of sulfide with fresh magma in the chonolith.     

Geochemical processes controlling arsenic mobility in groundwater: A case study of arsenic mobilization and natural attenuation

The behavior of As in the subsurface environment was examined along a transect of groundwater monitoring wells at a Superfund site, where enhanced reductive dechlorination (ERD) is being used for the remediation of groundwater contaminated with chlorinated solvents. The transect was installed parallel to the groundwater flow direction through the treatment area. The ERD technology involves the injection of organic C (OC) to stimulate in situ microbial dechlorination processes. A secondary effect of the ERD treatment at this site, however, is the mobilization of As, as well as Fe and Mn. The concentrations of these elements are low in groundwater collected upgradient of the ERD treatment area, indicating that, in the absence of the injected OC, the As that occurs naturally in the sediment is relatively immobile. Batch experiments conducted using sediments from the site inoculated with an Fe(III)- and As(V)-reducing bacterium and amended with lactate resulted in mobilization of As, Fe and Mn, suggesting that As mobilization in the field is due to microbial processes.     

Formation of chlorite rims and the impact of pore-lining chlorite on reservoir quality: a case study from Shiqianfeng sandstones in upper Permian of Dongpu Depression,Bohai Bay Basin,eastern China

The source material, precursor and formation processes of chlorite rims, and impact of pore-lining chlorite on reservoir quality of the Shiqianfeng sandstones, Dongpu Depression, Bohai Bay Basin, China, are studied using an integrated approach, including core observation, point-count analysis of thin-sections, scanning electron microscopy, electron microprobe analysis and cathode luminescence. The petrographic analysis shows that chlorite rims consist of grain-coating chlorite, poorly crystalline pore-lining chlorite and euhedral-crystallised pore-lining chlorite. The chemical composition shows that pore-lining chlorite is mainly Fe-rich with an average of 0.785 for Fe/(Fe+Mg) ratio. Petrographic analysis shows a large amount of volcanic dust (3.0～16.0 vol%, average of 7.93 vol%) in Shiqianfeng sandstones, which determines the formation of poorly crystalline pore-lining chlorite. Transformation of volcanic dust to smectite rims started with shallow burial depth at an early diagenetic stage, followed by in situ alteration of the smectite rims to poorly crystalline chlorite rims. Euhedral-crystallised chlorite mainly develops in sandstones with high porosity, high permeability and open flow systems. Pore-lining chlorite can inhibit quartz overgrowth but cannot effectively prevent pore-filling of authigenic quartz, carbonate and kaolinite cements, and therefore cannot prevent porosity destruction. However, the occurrence of pore-lining chlorite is a barometer of good reservoir quality and intense hydrodynamic conditions.     

Centennial to millennial geomagnetic secular variation

A time-varying spherical harmonic model of the palaeomagnetic field for 0–7 ka is used to investigate large-scale global geomagnetic secular variation on centennial to millennial scales. We study dipole moment evolution over the past 7 kyr, and estimate its rate of change using the Gauss coefficients of degree 1 (dipole coefficients) from the CALS7K.2 field model and by two alternative methods that confirm the robustness of the predicted variations. All methods show substantial dipole moment variation on timescales ranging from centennial to millennial. The dipole moment from CALS7K.2 has the best resolution and is able to resolve the general decrease in dipole moment seen in historical observations since about 1830. The currently observed rate of dipole decay is underestimated by CALS7K.2, but is still not extraordinarily strong in comparison to the rates of change shown by the model over the whole 7 kyr interval. Truly continuous phases of dipole decrease or increase are decadal to centennial in length rather than longer-term features. The general large-scale secular variation shows substantial changes in power in higher spherical harmonic degrees on similar timescales to the dipole. Comparisons are made between statistical variations calculated directly from CALS7K.2 and longer-term palaeosecular variation models: CALS7K.2 has lower overall variance in the dipole and quadrupole terms, but exhibits an imbalance between dispersion in   g ¹₂  and   h ¹₂  , suggestive of long-term non-zonal structure in the secular variations.     

Receiver function method in reflection seismology

The receiver function method was originally developed to analyse earthquake data recorded by multicomponent (3C) sensors and consists in deconvolving the horizontal component by the vertical component. The deconvolution process removes travel path effects from the source to the base of the target as well as the earthquake source signature. In addition, it provides the possibility of separating the emergent P and PS waves based on adaptive subtraction between recorded components if plane waves of constant ray parameters are considered. The resulting receiver function signal is the local PS wave's impulse response generated at impedance contrasts below the 3C receiver.We propose to adapt this technique to the wide‐angle multi‐component reflection acquisition geometry. We focus on the simplest case of land data reflection acquisition. Our adapted version of the receiver function approach consists in a multi‐step procedure that first removes the P wavefield recorded on the horizontal component and next removes the source signature. The separation step is performed in the τ?p domain while the source designature can be achieved in either the τ?p or the t?x domain. Our technique does not require any a priori knowledge of the subsurface. The resulting receiver function is a pure PS‐wave reflectivity response, which can be used for amplitude versus slowness or offset analysis. Stack of the receiver function leads to a high‐quality S wave image.