A generalized transformation approach for simulating steady-state variably-saturated subsurface flow

A numerical method is proposed to accurately and efficiently compute a direct steady-state solution of the nonlinear Richards equation. In the proposed method, the Kirchhoff integral transformation and a complementary transformation are applied to the governing equation in order to separate the nonlinear hyperbolic characteristic from the linear parabolic part. The separation allows the transformed governing equation to be applied to partially- to fully-saturated systems with arbitrary constitutive relations between primary (pressure head) and secondary variables (relative permeability). The transformed governing equation is then discretized with control volume finite difference/finite element approximations, followed by inverse transformation. The approach is compared to analytical and other numerical approaches for variably-saturated flow in 1-D and 3-D domains. The results clearly demonstrate that the approach is not only more computationally efficient but also more accurate than traditional numerical solutions. The approach is also applied to an example flow problem involving a regional-scale variably-saturated heterogeneous system, where the vadose zone is up to 1 km thick. The performance, stability, and effectiveness of the transform approach is exemplified for this complex heterogeneous example, which is typical of many problems encountered in the field. It is shown that computational performance can be enhanced by several orders of magnitude with the described integral transformation approach.     

The K-band luminosity function of nearby field galaxies

We present a measurement of the K -band luminosity function (LF) of field galaxies obtained from near-infrared imaging of a sample of 345 galaxies selected from the Stromlo-APM Redshift Survey. The LF is reasonably well fitted over the 10-mag range −26 M _K −16 by a Schechter function with parameters α =−1.16±0.19, M *=−23.58±0.42 and φ *=0.012±0.008 Mpc⁻³, assuming a Hubble constant of H ₀=100 km s⁻¹ Mpc⁻¹. We have also estimated the LF for two subsets of galaxies subdivided by the equivalent width of the H α emission line at EW(H α )=10 Å. There is no significant difference in LF shape between the two samples, although there is a hint (∼1 σ significance) that emission-line galaxies (ELGs) have M * roughly 1 mag fainter than non-ELGs. Contrary to the optical LF, there is no difference in faint-end slope α between the two samples.     

GCM-simulated hydrology in the Arctic duringn the past 21,000 years

The Paleoclimates from Arctic Lakes and Estuaries (PALE) project has chosen to conduct high resolution data-model comparisons for the Arctic region at 21 and 10 (calendar) ka BP. The model simulations for 21, 10, and 0 ka BP were conducted with the GENESIS 2.0 GCM. The 10 ka BP simulation was coupled to the EVE vegetation model. The primary boundary conditions differing from present at 21 ka BP were the northern hemisphere ice sheets and lower CO2, and at 10 ka BP were the orbital insolation and smaller northern hemisphere ice sheets. The purpose of this article is to discuss the hydrological consequences of these simulations.At the Last Glacial Maximum (21 ka BP) the large ice sheets over North America and Eurasia and the lower CO2 levels produced a colder climate than present, with less precipitation throughout the Arctic, except where circulation was altered by the ice sheets. At 10 ka BP greater summer insolation resulted in a warmer and wetter Beringia, but conditions remained cold and dry in the north Atlantic sector, in the vicinity of the remnant ice sheets. Less winter insolation at 10 ka BP resulted in colder and drier conditions throughout the Arctic. Precipitation - evaporation generally correlated with precipitation except where changes in the surface type (ice sheets, vegetation at 10 ka BP, or sea level at 21 ka BP) caused large changes in the evaporation rate. The primary hydrological differences (from present) at 21 and 10 ka BP correlated with the temperature differences, which were a direct result of the large-scale boundary condition changes.     

Macrocrystal phlogopite Rb–Sr dates for the Ekati property kimberlites, Slave Province, Canada: evidence for multiple intrusive episodes in the Paleocene and Eocene

New Rb–Sr age determinations using macrocrystal phlogopite are presented for 27 kimberlites from the Ekati property of the Lac de Gras region, Slave Province, Canada. These new data show that kimberlite magmatism at Ekati ranges in age from at least Late Paleocene (61 Ma) to Middle Eocene time (45 Ma). Older, perovskite-bearing kimberlites from Ekati extend this age range to Late Cretaceous time (74 Ma). Within this age range, emplacement episodes at 48, 51–53, 55–56 and 59–61 Ma can be recognized. Middle Eocene kimberlite magmatism of the previously dated Mark kimberlite (47.5 Ma) is shown to include four other pipes from the east-central Ekati property. A single kimberlite (Aaron) may be younger than the 47.5 Ma Mark kimberlite. The economically important Panda kimberlite is precisely dated in this study to be 53.3±0.6 Ma using the phlogopite isochron method, and up to six additional kimberlites from the central Ekati property have Early Eocene ages indistinguishable from that of Panda, including the Koala and Koala North occurrences. Late Paleocene 55–56 Ma kimberlite magmatism, represented by the Diavik kimberlite pipes adjacent to the southeastern Ekati property, is shown to extend onto the southeastern Ekati property and includes three, and possibly four, kimberlites. A precise eight-point phlogopite isochron for the Cobra South kimberlite yields an emplacement age of 59.7±0.4 Ma; eight other kimberlites from across the Ekati property have similar Late Paleocene Rb–Sr model ages. The addition of 27 new emplacement ages for kimberlites from the Ekati property confirms that kimberlite magmatism from the central Slave Province is geologically young, despite ages ranging back to Cambrian time from elsewhere in the Slave Province. With the available geochronologic database, Lac de Gras kimberlites with the highest diamond potential are currently restricted to the 51–53 and 55–56 Ma periods of kimberlite magmatism.     

The geology of kimberlite pipes of the Ekati property, Northwest Territories, Canada

This paper reviews key characteristics of kimberlites on the Ekati property, NWT, Canada. To date 150 kimberlites have been discovered on the property, five of which are mined for diamonds. The kimberlites intrude Archean basement of the central Slave craton. Numerous Proterozoic diabase dykes intrude the area. The Precambrian rocks are overlain by Quaternary glacial sediments. No Phanerozoic rocks are present. However, mudstone xenoliths and disaggregated sediment within the kimberlites indicate that late-Cretaceous and Tertiary cover (likely <200 m) was present at the time of emplacement. The Ekati kimberlites range in age from 45 to 75 Ma. They are mostly small pipe-like bodies (surface area mostly <3 ha but up to 20 ha) that typically extend to projected depths of 400–600 m below current surface. Pipe morphologies are strongly controlled by joints and faults. The kimberlites consist primarily of variably bedded volcaniclastic kimberlite (VK). This is dominated by juvenile constituents (olivine and lesser kimberlitic ash) and variable amounts of exotic sediment (primarily mud), with minor amounts of xenolithic wall-rock material (generally <5%). Kimberlite types include: mud-rich resedimented VK (mRVK); olivine-rich VK (oVK); sedimentary kimberlite; primary VK (PVK); tuffisitic kimberlite (TK) and magmatic kimberlite (MK). The presence and arrangement of these rock types varies widely. The majority of bodies are dominated by oVK and mRVK, but PVK is prominent in the lower portions of certain kimberlites. TK is rare. MK occurs primarily as precursor dykes but, in a few cases, forms pipe-filling intrusions. The internal geology of the kimberlites ranges from simple single-phase pipes (RVK or MK), to complex bodies with multiple, distinct units of VK. The latter include pipes infilled with steep, irregular VK blocks/wedges and at least one case in which the pipe is occupied by well-defined sub-horizontal VK phases, including a unique, 100-m-thick graded sequence. The whole-rock compositions of VK samples suggest significant loss of kimberlitic fines during eruption followed by variable dilution by surface sediment and concurrent incorporation of kimberlitic ash. Diamond distribution within the kimberlites reflects the amount and nature of mantle material sampled by individual kimberlite phases, but is modified considerably by eruption and depositional processes. The characteristics of the Ekati kimberlites are consistent with a two-stage emplacement process: (1) explosive eruption/s causing vent clearing followed by formation of a significant tephra rim/cone of highly fragmented, olivine-enriched juvenile material with varying amounts of kimberlitic ash and surface sediments (predominantly mud); and (2) infilling of the vent by direct deposition from the eruption column and/or resedimentation of crater rim materials. The presence of less fragmented, juvenile-rich PVK in the lower portions of certain pipes and the intrusion of large volumes of MK to shallow levels in some bodies suggest emplacement of relatively volatile-depleted, less explosive kimberlite in the later stages of pipe formation and/or filling. Explosive devolatilisation of CO₂-rich kimberlite magma is interpreted to have been the dominant eruption mechanism, but phreatomagmatism is thought to have played a role and, in certain cases, may have been dominant.     

Summer temperature variability across four urban neighborhoods in Knoxville,Tennessee, USA

The urban heat island (UHI) is a well-documented effect of urbanization on local climate, identified by higher temperatures compared to surrounding areas, especially at night and during the warm season. The details of a UHI are city-specific, and microclimates may even exist within a given city. Thus, investigating the spatiotemporal variability of a city’s UHI is an ongoing and critical research need. We deploy ten weather stations across Knoxville, Tennessee, to analyze the city’s UHI and its differential impacts across urban neighborhoods: two each in four neighborhoods, one in more dense tree cover and one in less dense tree cover, and one each in downtown Knoxville and Ijams Nature Center that serve as control locations. Three months of temperature data (beginning 2 July 2014) are analyzed using paired-sample t tests and a three-way analysis of variance. Major findings include the following: (1) Within a given neighborhood, tree cover helps negate daytime heat (resulting in up to 1.19 ^°C lower maximum temperature), but does not have as large of an influence on minimum temperature; (2) largest temperature differences between neighborhoods occur during the day (0.38–1.16 ^°C difference), but larger differences between neighborhoods and the downtown control occur at night (1.04–1.88 ^°C difference); (3) presiding weather (i.e., air mass type) has a significant, consistent impact on the temperature in a given city, and lacks the differential impacts found at a larger-scale in previous studies; (4) distance from city center does not impact temperature as much as land use factors. This is a preliminary step towards informing local planning with a scientific understanding of how mitigation strategies may help minimize the UHI and reduce the effects of extreme weather on public health and well-being.     

GGR Critical Review of Analytical Developments in 2003

The articles that comprise this critical review serve to draw attention to research papers published in specific fields of interest during 2003, provide critical comment on the relevance and importance of individual publications in these fields, and offer an overview of the comparative importance of advances in particular areas. In this way, these articles aim to assist experts in the field by keeping them informed of relevant recent publications, as well as providing an important resource for students or early career researchers who are embarking on studies in an area new to them. This year, five papers provide summaries of developments in bulk sample determinations employing (1) ICP-AES and ICP-MS (trace elements), (2) XRF and atomic absorption spectrometry and INAA, (3) isotope ratio measurements (TIMS, MC-ICP-MS, ICP-MS, ToF), as well as in situ measurements conducted using (4) secondary ion mass spectrometry and (5) laser ablation ICP-MS (trace element and isotope ratio determinations).     

The last glaciation of Germania Land and adjacent areas,northeast Greenland

A profile across the unglaciated coast of northeast Greenland at 77°N was studied with regard to the Quaternary stratigraphy and glacial history. The Germania Land peninsula is characterised by till-covered lower ground which contrasts sharply with the blockfields and extensive gelifluction deposits of its higher altitudes. Two glaciations are distinguished. The older one extended over the entire area and had its margin on the continental shelf. The younger one, of Late Weichselian age, reached the present coastline and several mountains and high plateaus on western Germania Land formed nunataks. The Late Weichselian glaciation was more extensive and occurred later on the Germania Land peninsula than on the coast further south. Radiocarbon dates suggest that the glacier margin rested to the east of the present coastline until ca. 10 000 yr BP. This correlates with the Late Weichselian Milne Land Stage, which is found as a late glacial readvance along the coast of East Greenland. A series of recessional moraines formed during the deglaciation were probably caused by glacier dynamics, as opposed to being of climatic origin.     

Precise measurement of seismic travel times from an air gun source

Summary. Seismic travel times for extrema, zero-crossings, or entire body waves need to be determined precisely to one part in 10³ or better in several varieties of seismic studies employing an impulsive artificial source. Examples are crosshole surveys which delineate rock crack distribution separating the holes and monitoring of crustal seismic travel times in earthquake precursor studies. A timing resolution of one part in 10³ has been achieved previously using digitally recorded seismic data. These methods, however, do not use interpolation between digitized data points as a method to increase the timing resolution. We report travel-time determinations based on interpolation between digitized points which achieve a precision of two parts in 10⁴, a five-fold improvement over the existing methods. In addition, the effects of seismic noise on travel-time measurement have been compared for the extremum location, the unnormalized correlation, and the normalized correlation method. The following conclusions are drawn from this comparison: (1) the normalized correlation method provides an 18–55 per cent improvement in the standard deviation of the mean over the extremum location method, and (2) results as accurate as those by the normalized correlation can be obtained by the unnormalized correlation if a complete up-and-down swing of the waveform is used as the master trace and if the master trace is close to being sinusoidal. The advantage of the unnormalized correlation over the normalized correlation is speed; the unnormalized correlation is faster by a factor of 28 in computing time.     

Storylines: an alternative approach to representing uncertainty in physical aspects of climate change

Investigating the relationships between climate extremes and crop yield can help us understand how unfavourable climatic conditions affect crop production. In this study, two statistical models, multiple linear regression and random forest, were used to identify rainfall extremes indices affecting wheat yield in three different regions of the New South Wales wheat belt. The results show that the random forest model explained 41–67% of the year-to-year yield variation, whereas the multiple linear regression model explained 34–58%. In the two models, 3-month timescale standardized precipitation index of Jun.–Aug. (SPI_JJA), Sep.–Nov. (SPI_SON), and consecutive dry days (CDDs) were identified as the three most important indices which can explain yield variability for most of the wheat belt. Our results indicated that the inter-annual variability of rainfall in winter and spring was largely responsible for wheat yield variation, and pre-growing season rainfall played a secondary role. Frequent shortages of rainfall posed a greater threat to crop growth than excessive rainfall in eastern Australia. We concluded that the comparison between multiple linear regression and machine learning algorithm proposed in the present study would be useful to provide robust prediction of yields and new insights of the effects of various rainfall extremes, when suitable climate and yield datasets are available.