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.     

Equation of state under nuclear statistical equilibrium conditions

We investigate a three-parameter equation of state for stellar matter under nuclear statistical equilibrium conditions in the ranges of temperatures 3×10⁹–10¹¹ K and densities 10⁴–10¹³ g cm^?3 and for various ratios of the total number of neutrons to the total number of protons within the range 1–1.5. These conditions correspond to the initial stages of the gravitational collapse of iron stellar cores that are accompanied by nonequilibrium matter neutronization. We analyze the effect of the excited levels of atomic nuclei on the thermodynamic properties of the matter. We show that this effect is insignificant at low densities, ρ?10¹⁰ g cm^?3, but it leads to an expansion of the instability region, γ<4/3, at higher densities. The incorporated effects of the Fermi degeneracy of free nucleons prove to be insignificant, because their concentrations are low at low temperatures. In the future, we plan to investigate the effects of Coulomb interactions and neutron-rich nuclei on the thermodynamic properties of the matter.     

Primitive island arc and oceanic lavas from the hunter ridge-hunter fracture zone. Evidence from glass,olivine and spinel compositions

Summary At its southernmost end, the main spreading centre of the North Fiji Basin is propagating into arc crust of the poorly-known Hunter Ridge. We define nine magmatic groups from major element glass chemistry and olivine and spinel compositions in samples dredged from twenty six sites in this area by the R/V Academician A. Nesmeyanov in 1990. These include groups of boninites, island arc tholeiites (IAT), mid-ocean ridge basalts (MORB), enriched mid-ocean ridge basalts (E-MORB), olivine porphyritic andesites and basaltic andesite and Na-rhyolites. Primitive lavas containing highly forsteritic olivine phenocrysts are common in all the groups, except for the rhyolites.We report over 100 glass analyses for dredged rocks from this region, and about 300 olivine-spinel pairs for representatives of all the magmatic groups identified, except the Na-rhyolites.The MORB in this region are probably produced at the propagating spreading centre in the southern part of the North Fiji Basin. Juxtaposition of shallow, hot MORB-source diapirs supplying the MORB in this area, and the sub-arc damp, refractory upper mantle beneath the Hunter Ridge, provides suitable petrogenetic conditions to produce a range of magma types, from island arc tholeiites through to high-Ca boninites. The latter were recovered in eleven dredges.The E-MORB lavas recovered from the extreme southern margin of the North Fiji Basin are shown to be essentially identical to those dredged from adjacent older South Fiji Basin crust. It is hypothesized that the former were either scraped off the South Fiji Basin crust during an episode of oblique subduction that may have generated the Hunter Ridge during the last 5 Myr, or alternatively, that slices of the South Fiji Basin crust were trapped and incorporated into the North Fiji Basin as the subduction zones fronting the Vanuatu arc stepped or propagated southward.

---

Primitive Inselbogen- und ozeanische Laven von Hunter Rücken und der Hunter BruchZone: die Bedeutung der Zusammensetzung von Glas, Olivin und Spinell

Zusammenfassung Das spreading centre des Nord-Fidschi-Beckens setzt sich an seinem südlichsten Ende in die Inselbogen-Kruste des noch wenig bekannten Hunter-Rückens fort.Wir definieren 9 magmatische Gruppen auf der Basis der Hauptelement-Zusammensetzung von Gläsern und der Zusammensetzung von Olivin und Spinell in Proben die das Forschungsschiff R/V Akademiker A. Nesmeyanov von 26 Stellen im Jahr 1990 aufgesammelt hat. Zu diesen gehören Gruppen von Boniniten, Inselbogentholeiiten (IAT), Basalten zentralozeanischer Rücken (MORB), angereicherte zentralozeanische Rücken (E-MORB), Olivin-porphyritische Andesite und basaltische Andesite, sowie Na-Rhyolite. Primitive Laven mit Olivinkristallen, die reich an Forsterit-Komponenten sind, kommen in allen diesen Gruppen, mit Ausnahme der Rhyolite, vor.Wir legen über 100 Glas-Analysen von Gesteinen aus diesem Gebiet vor und über 300 Olivin-Spinell-Paare für Vertreter aller der hier identifizierten magmatischen Gruppen, mit Ausnahme der Natriumrhyolite.Die MORB in diesem Gebiet sind wahrscheinlich an dem aktiven Spreading Center im Südteil des Nord-Fidschi-Beckens entstanden. Das nebeneinander Vorkommen von seichten heißen MORB-Quellen Diapiren, die MORB in diesem Gebiet erzeugen, und der Sub-Inselbogen, wasserhaltige, refraktäre obere Mantel unter dem Hunter Rücken führen zu geeigneten petrogenetischen Bedingungen für die Entstehung einer Vielfalt von Magmatypen, von Inselbogentholeiiten bis zu sehr Kalziumreichen Boniniten. Die letzteren wurden in 11 Dredge-Proben gefunden.Die E-MORS Laven, die aus dem extremen Südteil des Nord-Fidschi-Beckens entnommen wurden, sind im wesentlichen mit jenen identisch, die aus der benachbarten älteren Kruste des Süd-Fidschi-Beckens stammen. Es wird vermutet, daß die letzteren entweder aus dem Süd-Fidschi-Becken während einer Episode obliger Subduktion, die in den letzten 5 Millionen Jahren zur Entstehung des Hunter Ridge geführt hat, abgeschert wurden, oder auch daß Teile der Kruste des Süd-Fidschi-Beckens in das Nord-Fidschi-Becken inkorporiert wurden als Subduktionszonen gegenüber dem Vanuato-Bogen sich nach Süden fortbewegten.

     

High-pressure partial melting of gabbro and its role in the Hawaiian magma source

We have conducted high-pressure experiments on a natural oceanic gabbro composition (Gb108). Our aim was to test recent proposals that Sr-enrichment in rare primitive melt inclusions from Mauna Loa, Hawaii, may have resulted from melting of garnet pyroxenite formed in the magma source regions by reaction of peridotite with siliceous, Sr-enriched partial melts of eclogite of gabbroic composition. Gb108 is a natural, Sr-enriched olivine gabbro, which has a strong positive Sr anomaly superimposed on an overall depleted incompatible trace element pattern, reflecting its origin as a plagioclase-rich cumulate. At high pressures it crystallises as a coesite eclogite assemblage, with the solidus between 1,300 and 1,350°C at 3.5 GPa and 1,450 and 1,500°C at 4.5 GPa. Clinopyroxenes contain 4–9% Ca-eskolaite component, which varies systematically with pressure and temperature. Garnets are almandine and grossular-rich. Low degree partial melts are highly siliceous in composition, resembling dacites. Coesite is eliminated between 50 and 100°C above the solidus. The whole-rock Sr-enrichment is primarily hosted by clinopyroxene. This phase dominates the mode (>75 wt%) at all investigated PT conditions, and is the major contributor to partial melts of this eclogite composition. Hence the partial melts have trace element patterns sub-parallel to those of clinopyroxene with ≈10× greater overall abundances and with strong positive Sr anomalies. Recent studies of primitive Hawaiian volcanics have suggested the incorporation into their source regions of eclogite, formerly gabbroic material recycled through the mantle at subduction zones. The models suggest that formerly gabbroic material, present as eclogite in the Hawaiian plume, partially melted earlier than surrounding peridotite (i.e. at higher pressure) because of the lower solidus temperature of eclogite compared with peridotite. This produced highly siliceous melts which reacted with surrounding peridotite producing hybrid pyroxene + garnet lithologies. The Sr-enriched nature of the formerly plagioclase-rich gabbro was present in the siliceous partial melts, as demonstrated by these experiments, and was transferred to the reactive pyroxenite. These in turn partially melted, producing Sr-enriched picritic liquids which mixed with normal picritic partial melts of peridotite before eruption. On rare occasions these mixed, relatively Sr-rich melts were trapped as melt inclusions in primitive olivine phenocrysts. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

MODELING HEADCUT DEVELOPMENT AND MIGRATION IN UPLAND CONCENTRATED FLOWS

On hillslopes and agricultural fields, discrete areas of intense, localized soil erosion commonly take place in the form of migrating headcuts. These erosional features significantly increase soil loss and landscape degradation, yet the unsteady, transient, and migratory habits of headcuts complicate their phenomenological and erosional characterization. Here a unique experimental facility was constructed to examine actively migrating headcuts typical of upland concentrated flows. Essential components of the facility include a deep soil cavity with external drainage, rainfall simulator, capacity for overland flow, and a video recording technique for data collection. Results from these experiments show that： （1） after a short period of adjustment, headcut migration attained a steady-state condition, where the rate of migration, scour hole geometry, and sediment discharge remain constant with time; （2） boundary conditions of higher rates of overland flow, steeper bed slopes, and larger initial headcut heights produced systematically larger scour holes with higher rates of soil erosion; and （3） during migration, the turbulent flow structure within the scour hole remained unchanged, consisting of an overfall nappe at the brink transitioning into a reattached wall jet with two recirculation eddies within the plunge pool. The systematic behavior of headcut development and migration enabled the application of modified jet impingement theory to predict with good success the characteristics of the impinging jet, the depth of maximum scour, the rate of headcut migration, and the rate of sediment erosion. These laboratory data and the analytical formulation can be used in conjunction with soil erosion prediction technology to improve the management of agricultural areas impacted by headcut development and ephemeral gully erosion.     

On the circulation in the East Frisian Wadden Sea: numerical modeling and data analysis

 In this paper we use a combination of numerical modeling and data analysis to gain a better understanding of the major characteristics of the circulation in the East Frisian Wadden Sea. In particular, we concentrate on the asymmetry of the tidal wave and its modulation in the coastal area, which results in a complex pattern of responses to the sea-level forcing from the North Sea. The numerical simulations are based on the 3-D primitive equation General Estuarine Transport Model (GETM) with a horizontal resolution of 200 m and terrain-following vertical coordinates. The model is forced at its open boundaries with sea-level data from an operational model for the German Bight (German Hydrographic Office). The validation data for our model simulations include time series of tidal gauge data and surface currents measured at a pile in the back-barrier basin of the Island Langeoog, as well as several ADCP transects in the Accumer Ee tidal inlet. Circulation and turbulence characteristics are investigated for typical situations driven by spring and neap tides, and the analysis is focused on dominating temporal and spatial patterns. By investigating the response of five back-barrier basins with rather different morphologies to external forcing, an attempt is made to elucidate the dominating physical balances controlling the circulation in the individual sub-basins. It is demonstrated that the friction at the seabed tends to slow down the tidal signal in the shallow water. This leads to the establishment of flood dominance in the shallow sea north of the barrier islands. South of the islands, where the water volume of the channels at low tide is smaller than the tidal prism, the asymmetry of the tidal signal is shifted towards ebb dominance, a feature which is particularly pronounced at spring tide. At the northern open boundary, the tidal wave propagating from west to east generates a sea-level difference of ∼1 m along the boundary, and thereby triggers vigorous alongshore currents. The frictional control in the model is located in the inlets, as well as along the northern boundary. The correlation between velocity and turbulent kinetic energy tends to the establishment of a net southward transport, giving theoretical support to the observed accumulation of sediments on the intertidal flats. Weak turbulence along the northern shores of the barrier islands and the small magnitude of the residual currents there promote accumulation of suspended matter in these areas, although wave action will generally counteract this effect. Received: 29 May 2002 / Accepted: 26 September 2002 Responsible Editor: Jean-Marie Beckers Acknowledgements We are indebted to S. Dick for providing the data from the operational model of BSH and to B. Flemming for the useful discussions. The topography data and Fig. 1 have been prepared in cooperation with F. Meyer. Figure 2 has been prepared by G. Brink-Spalink. We also thank for the comments from an anonymous reviewer which helped to improve our paper.     

Geochemistry of sedimentary rocks from Permian–Triassic boundary sections of Tethys Himalaya: implications for paleo-weathering,provenance, and tectonic setting

The geochemical characteristics of two sections—the Permian–Triassic boundary (PTB) Guryul Ravine section, Kashmir Valley, Jammu and Kashmir, India; and the Attargoo section, Spiti Valley, Himachal Pradesh, India—have been studied in the context of provenance, paleo-weathering, and plate tectonic setting. These sections represent the siliciclastic sedimentary sequence from the Tethys Himalaya. The PTB siliciclastic sedimentary sequence in these regions primarily consists of sandstones and shales with variable thickness. Present studied sandstones and shales of both sections had chemical index of alteration values between 65 and 74; such values reveal low-to-moderate degree of chemical weathering. The chemical index of weathering in studied samples ranged from 71 to 94, suggesting a minor K-metasomatism effect on these samples. Plagioclase index of alteration in studied sections ranged from 68 to 92, indicating a moderate degree of weathering of plagioclase feldspars. The provenance discriminant function diagram suggests that the detritus involved in the formation of present studied siliciclastic sedimentary rocks fall in quartzose sedimentary and felsic igneous provenances. These sediments were deposited in a passive continental margin plate tectonic setting according to their location on a Si₂O versus K₂O/Na₂O tectonic setting diagram.     

Computer processing of LANDSAT-1 data and imagery generation

LANDSAT-1 (ERTS-1) data of the complete scene, 1218-04320, covering an area 185 Km x 185 Km in the Godavari Valley, India was computer processed in the Computer Division of the Oil and Natural Gas Commissíon to produce an imagery of readily recognizable land forms. Comparison with the original maps supplied by NASA showed an amazing agreement in quality. Computer oriented data enhancement techniques were developed and the data was reprocessed with the help of these techniques. The resulting imagery reveal an improvement over the originals. The computer software package developed to-date comprises nine different programmes to process the LANDSAT data on the IBM 360, 370 series.     

Performance analysis of the spaceborne laser ranging system

The Spaceborne Laser Ranging System is a proposed short pulse laser on board an orbiting spacecraft.^1,2,3,4 It measures the distance between the spacecraft and many laser retroreflectors (targets) deployed on the Earth’s surface. The precision of these range measurements is assumed to be about ±2 cm (M. W. Fitzmaurice, private communication). These measurements are then used together with the orbital dynamics of the spacecraft, to derive the relative position of the laser ground targets. Assuming a six day observing period with 50% cloud cover, uncertainties in the baseline for target separations of 50 km to 1200 km were estimated to be on the order of 1 to 3 cm and the corresponding values in the vertical direction, ranged from 1 cm to 12 cm. By redetermining the measurements of the relative target positions, the estimated precision in the baseline for a target separation of 50 km is less than 0.3 cm and for a separation of 1200 km is less than 1 cm. In the vertical direction, the estimated precision ranged from 0.4 cm to 4.0 cm respectively. As a result of the repeated estimation of the relative laser target positions, most of the non-temporal effects of error sources as exemplified by the errors in geopotential are reduced. The Spaceborne Laser Ranging System’s capability of determining baselines to a high degree of precision provides a measure of strain and strain rate as shown byCohen, 1979. These quantities are essential for crustal dynamic studies which include determination and monitoring of strain near seismic zones, land subsidence, and edifice building preceding volcanic eruptions. It is evident that such a system can also be used for geodetic surveys where such precisions are more than adquate.     

The source of sulfur in sulfide deposits in the Siberian Platform traps (from isotope data)

The source of sulfur in giant Norilsk-type sulfide deposits is discussed. A review of the state of the problem and a critical analysis of existing hypotheses are made. The distribution of δ³⁴S in sulfides of ore occurrences and small and large deposits and in normal sedimentary, metamorphogenic, and hypogene sulfates is considered. A large number of new δ³⁴S data for sulfides and sulfates in various deposits, volcanic and terrigenous rocks, coals, graphites, and metasomatites are presented. The main attention is focused on the objects of the Norilsk and Kureika ore districts. The δ³⁴S value varies from -14 to + 22.5‰ in sulfides of rocks and ores and from 15.3 to 33‰ in anhydrites. In sulfide-sulfate intergrowths and assemblages, δ³⁴S is within 4.2-14.6‰ in sulfides and within 15.3-21.3‰ in anhydrites. The most isotopically heavy sulfur was found in pyrrhotite veins in basalts (δ³⁴S = 21.6‰), in sulfate veins cutting dolomites (δ³⁴S = 33‰), and in subsidence caldera sulfates in basalts (δ³⁴S = 23.2-25.2‰). Sulfide ores of the Tsentral’naya Shilki intrusion have a heavy sulfur isotope composition (δ³⁴S = + 17.7‰ (n = 15)). Thermobarogeochemical studies of anhydrites have revealed inclusions of different types with homogenization temperatures ranging from 685 °C to 80 °C. Metamorphogenic and hypogene anhydrites are associated with a carbonaceous substance, and hypogene anhydrites have inclusions of chloride-containing salt melts. We assume that sulfur in the trap sulfide deposits was introduced with sulfates of sedimentary rocks (δ³⁴S = 22-24‰). No assimilation of sulfates by basaltic melt took place. The sedimentary anhydrites were “steamed” by hydrocarbons, which led to sulfate reduction and δ³⁴S fractionation. As a result, isotopically light sulfur accumulated in sulfides and hydrogen sulfide, isotopically heavy sulfur was removed by aqueous calcium sulfate solution, and “residual” metamorphogenic anhydrite acquired a lighter sulfur isotope composition as compared with the sedimentary one. The wide variations in δ³⁴S in sulfides and sulfates are due to changes in the physicochemical parameters of the ore-forming system (first of all, temperature and P_ch4) during the sulfate reduction. The regional hydrocarbon resources were sufficient for large-scale ore formation.