Investigating a Hierarchy of Eulerian Closure Models for Scalar Transfer Inside Forested Canopies

Modelling the transfer of heat, water vapour, and CO₂ between the biosphere and the atmosphere is made difficult by the complex two-way interaction between leaves and their immediate microclimate. When simulating scalar sources and sinks inside canopies on seasonal, inter-annual, or forest development time scales, the so-called well-mixed assumption (WMA) of mean concentration (i.e. vertically constant inside the canopy but dynamically evolving in time) is often employed. The WMA eliminates the need to model how vegetation alters its immediate microclimate, which necessitates formulations that utilize turbulent transport theories. Here, two inter-related questions pertinent to the WMA for modelling scalar sources, sinks, and fluxes at seasonal to inter-annual time scales are explored: (1) if the WMA is to be replaced so as to resolve this two-way interaction, how detailed must the turbulent transport model be? And (2) what are the added predictive skills gained by resolving the two-way interaction vis-à-vis other uncertainties such as seasonal variations in physiological parameters. These two questions are addressed by simulating multi-year mean scalar concentration and eddy-covariance scalar flux measurements collected in a Loblolly pine (P. taeda L.) plantation near Durham, North Carolina, U.S.A. using turbulent transport models ranging from K-theory (or first-order closure) to third-order closure schemes. The multi-layer model calculations with these closure schemes were contrasted with model calculations employing the WMA. These comparisons suggested that (i) among the three scalars, sensible heat flux predictions are most biased with respect to eddy-covariance measurements when using the WMA, (ii) first-order closure schemes are sufficient to reproduce the seasonal to inter-annual variations in scalar fluxes provided the canonical length scale of turbulence is properly specified, (iii) second-order closure models best agree with measured mean scalar concentration (and temperature) profiles inside the canopy as well as scalar fluxes above the canopy, (iv) there are no clear gains in predictive skills when using third-order closure schemes over their second-order closure counterparts. At inter-annual time scales, biases in modelled scalar fluxes incurred by using the WMA exceed those incurred when correcting for the seasonal amplitude in the maximum carboxylation capacity (V _{cmax, 25}) provided its mean value is unbiased. The role of local thermal stratification inside the canopy and possible computational simplifications in decoupling scalar transfer from the generation of the flow statistics are also discussed.

“The tree, tilting its leaves to capture bullets of light; inhaling, exhaling; its many thousand stomata breathing, creating the air”. Ruth Stone, 2002, In the Next Galaxy

     

Post‐breakup burial and exhumation of the southern margin of Africa

Despite many years of study, the processes involved in the development of the continental margin of southern Africa and the distinctive topography of the hinterland remain poorly understood. Previous thermochronological studies carried out within a monotonic cooling framework have failed to take into account constraints provided by Mesozoic sedimentary basins along the southern margin. We report apatite fission track analysis and vitrinite reflectance data in outcrop samples from the Late Jurassic to Early Cretaceous sedimentary fill of the Oudtshoorn, Gamtoos and Algoa Basins (Uitenhage Group), as well as isolated sedimentary remnants further west, plus underlying Paleozoic rocks (Cape Supergroup) and Permian‐Triassic sandstones from the Karoo Supergroup around the Great Escarpment. Results define a series of major regional cooling episodes. Latest Triassic to Early Jurassic cooling which began between 205 and 180 Ma is seen dominantly in basement flanks to the Algoa and Gamtoos Basins. This episode may have affected a wider region but in most places any effects have been overprinted by later events. The effects of Early Cretaceous (beginning between 145 and 130 Ma) and Early to mid‐Cretaceous (120–100 Ma) cooling are both delimited by major structures, while Late Cretaceous (85–75 Ma) cooling appears to have affected the whole region. These cooling events are all interpreted as dominantly reflecting exhumation. Higher Late Cretaceous paleotemperatures in samples from the core of the Swartberg Range, coupled with evidence for localised Cenozoic cooling, are interpreted as representing Cenozoic differential exhumation of the mountain range. Late Cretaceous paleotemperatures between 60°C and 90°C in outcropping Uitenhage Group sediments from the Oudtshoorn, Gamtoos and Algoa Basins require burial by between 1.2 and 2.2 km prior to Late Cretaceous exhumation. Because these sediments lie in depositional contact with underlying Paleozoic rocks in many places, relatively uniform Late Cretaceous paleotemperatures across most of the region, in samples of both basin fill and underlying basement, suggest the whole region may have been buried prior to Late Cretaceous exhumation. Cenozoic cooling (beginning between 30 and 20 Ma) is focussed mainly in mountainous regions and is interpreted as representing denudation which produced the modern‐day relief. Features such as the Great Escarpment are not related to continental break up, as is often supposed, but are much younger (post‐30 Ma). This history of post‐breakup burial and subsequent episodic exhumation is very different from conventional ideas of passive margin evolution, and requires a radical re‐think of models for development of continental margins.     

Discrete element modelling of material non‐coaxiality in simple shear flows

We investigate the quasi‐static simple shear flow of a two‐dimensional assembly of cohesionless particles using discrete element method (DEM) simulations. We focus on the unsteady flow regime where the solid would experience significant evolution of stresses, mobilised shear strength and dilation. We construct the DEM model using a discretised‐wall confined granular cell where the apparent boundary is allowed to dilate or contract synchronously with the confined solid. A rather uniform simple shear field is achieved across the whole assembly, which benefits rheological studies in generalising constitutive laws for continuum methods. We examine two aspects of the simple shear behaviour: macroscopic stress and strain rate evolution, particularly the non‐coaxiality between the principal directions of the two; and micromechanics such as evolution of fabric. For an initially anisotropic specimen sheared under constant normal pressure, the direction of principal stress rotates towards that of the principal strain rate, gradually reducing the degree of non‐coaxiality from about 45° to fluctuating around 0°. The rate in approaching coaxiality is slower in samples with larger initial porosity, stress ratio and mean stress. Generally, a faster rate in approaching coaxiality in simple shear is observed in a more dilatant sample, which often shows a larger degree of mobilised fabric anisotropy, suggesting the possible important role of instantaneous internal friction angle. The evolution of principal fabric direction resembles that of the principal stress direction. © 2013 The Authors. International Journal for Numerical and Analytical Methods in Geomechanics published by John Wiley & Sons, Ltd.     

Microwave Remote Sensing of the Temperature and Distribution of Sulfur Compounds in the Lower Atmosphere of Venus

A multi-wavelength radio frequency observation of Venus was performed on April 5, 1996, with the Very Large Array to investigate potential variations in the vertical and horizontal distribution of temperature and the sulfur compounds sulfur dioxide (SO₂) and sulfuric acid vapor (H₂SO₄(g)) in the atmosphere of the planet. Brightness temperature maps were produced which feature significantly darkened polar regions compared to the brighter low-latitude regions at both observed frequencies. This is the first time such polar features have been seen unambiguously in radio wavelength observations of Venus. The limb-darkening displayed in the maps helps to constrain the vertical profile of H₂SO₄(g), temperature, and to some degree SO₂. The maps were interpreted by applying a retrieval algorithm to produce vertical profiles of temperature and abundance of H₂SO₄(g) given an assumed sub-cloud abundance of SO₂. The results indicate a substantially higher abundance of H₂SO₄(g) at high latitudes (above 45°) than in the low-latitude regions. The retrieved temperature profiles are up to 25 K warmer than the profile obtained by the Pioneer Venus sounder probe at altitudes below 40 km (depending on location and assumed SO₂ abundance). For 150 ppm of SO₂, it is more consistent with the temperature profile obtained by Mariner 5, extrapolated to the surface via a dry adiabat. The profiles obtained for H₂SO₄(g) at high latitudes are consistent with those derived from the Magellan radio occultation experiments, peaking at around 8 ppm at an altitude of 46 km and decaying rapidly away from that altitude. At low latitudes, no significant H₂SO₄(g) is observed, regardless of the assumed SO₂ content. This is well below that measured by Mariner 10 (Lipa and Tyler 1979, Icarus39, 192-208), which peaked at ∼14 ppm near 47 km. Our results favor ≤100 ppm of SO₂ at low latitudes and ≤50 ppm in polar regions. The low-latitude value is statistically consistent with the results of Bézard et al. (1983, Geophs. Res. Lett.20, 1587-1590), who found that a sub-cloud SO₂ abundance of 130±40 ppm best matched their observations in the near-IR. The retrieved temperature profile and higher abundance of H₂SO₄(g) in polar regions are consistent with a strong equatorial-to-polar, cloud-level flow due to a Hadley cell in the atmosphere of Venus.     

A geocomputational analysis of Twitter activity around different world cities

The penetration and use of social media services differs from city to city. This paper is aimed to provide a comparison of the use of Twitter between different cities of the world. We present a temporal analysis of activity on Twitter in 15 cities. Our study consists of two parts: First, we created temporal graphs of the activity in the 15 cities, through which hours of high and low activity could be identified. Second, we created heat map visualizations of the Twitter activities during the period of 19 September 2012–25 September 2013. The heat map visualizations make the periods of intense and sparse activity apparent and provide a snapshot of the activity during the whole year.     

An experimental study of anorthosite dissolution in lunar picritic magmas: Implications for crustal assimilation processes

Experiments characterizing the kinetics of anorthosite dissolution in lunar picritic magmas (very low-Ti, low-Ti, and high-Ti picritic glasses) were conducted at 0.6 GPa and 1250-1400 °C using the dissolution couple method. Reaction between the anorthosite and lunar picritic magmas at 1250-1300 °C produced a spinel + melt layer. Reaction between the anorthosite and an olivine-saturated low-Ti magma at 1250-1300 °C produced a crystal-free region between the spinel + melt layer and the olivine-saturated magma. The anorthosite dissolution experiments conducted at 1400 °C simply dissolved anorthosite and did not result in a crystal-bearing region. The rate of anorthosite dissolution strongly depends on temperature and composition of the reacting melt. Concentration profiles that develop during anorthosite dissolution are nonlinear and extend from the picritic glass compositions to anorthite. These profiles feature a large and continuous variation in melt density and viscosity from the anorthosite-melt interface to the initial picritic magmas. In both the low-Ti and high-Ti magmas the diffusive fluxes of TiO₂, Al₂O₃, and SiO₂ are strongly coupled to the concentration gradients of CaO and FeO. Anorthosite dissolution may play an important role in producing the chemical variability of the lunar picritic magmas, the origin of spinel in the lunar basalts and picritic glasses, and the petrogenesis of the high-Al basalts.     

Trends in northern midlatitude atmospheric wave power from 1950 to 2099

In a warming climate, atmospheric wave activity and associated weather patterns may change, although conflicting results have been reported on this topic. Additionally, atmospheric wave changes in a future climate have mainly focused on waves of a specified spatial scale, rather than a particular spatiotemporal scale. Here, changes in the variability of Rossby waves of multiple spatiotemporal scales are analyzed using the wavenumber-frequency power spectrum, a tool commonly applied to analyze atmospheric equatorial waves. Daily 500 hPa geopotential height data over 40°–60°N from historical (1950–2005) and future (2006–2099) simulations from 20 models in the Coupled Model Intercomparison Project Phase 5 (CMIP5) under the RCP8.5 scenario were analyzed. When compared to the historical period, the late 21st century climate projections showed a decline in spectral power for both eastward and westward propagating waves with wavenumbers greater than 8 that spanned over all frequencies in all seasons, but an increase in mean power for eastward propagating waves with wavenumbers 1–7 over all frequencies was shown in winter and spring. This increase in power was accompanied by increased variance, i.e., an increased meridional extent of 500 hPa ridges and troughs, and was the result of increases in the mean number of high amplitude events and duration of activity within this wave band. These results indicate that large-scale (~ 104 km) eastward propagating weather systems may intensify with higher amplitudes for ridges and troughs, while short-scale (102–103 km) weather systems may decrease in their intensity due to reduced variability in the late 21st century under the high emissions scenario. Potential mechanisms for these changes are discussed, including enhanced Arctic warming and midlatitude-tropical interactions.     

Models of corundum origin from alkali basaltic terrains: a reappraisal

Corundums from basalt fields, particularly in Australia and Asia, include a dominant blue-green-yellow zoned “magmatic” suite (BGY suite) and subsidiary vari-coloured “metamorphic” suites. The BGY corundums have distinctive trace element contents (up to 0.04 wt% Ga₂O₃ and low Cr/Ga and Ti/Ga ratios <1). Different melt origins for BGY corundums are considered here from their inclusion and intergrowth mineralogy, petrologic associations and tectonic setting. Analysed primary inclusion minerals (over 100 inclusions) cover typical feldspars, zircon and Nb-Ta oxides and also include hercynite-magnetite, gahnospinel, rutile-ilmenite solid solution, calcic plagioclase, Ni-rich pyrrhotite, thorite and low-Si and Fe-rich glassy inclusions. This widens a previous inclusion survey; New England, East Australia corundums contain the most diverse inclusion suite known from basalt fields (20 phases). Zircon inclusion, intergrowth and megacryst rare earth element data show similar patterns, except for Eu which shows variable depletion. Temperature estimates from magnetite exsolution, feldspar compositions and fluid inclusion homogenization suggest that some corundums crystallized between 685–900 °C. Overlap of inclusion Nb, Ta oxide compositions with new comparative data from niobium-yttrium-fluorine enriched granitic pegmatites favour a silicate melt origin for the corundums. The feasibility of crystallizing corundum from low-volume initial melting of amphibole-bearing mantle assemblages was tested using the MELTS program on amphibole-pyroxenite xenolith chemistry from basalts. Corundum appears in the calculations at 720–880 °C and 0.7–1.1 GPa with residual feldspathic assemblages that match mineral compositions found in corundums and their related xenoliths. A model that generates melts from amphibole-bearing lithospheric mantle during magmatic plume activity is proposed for BGY corundum formation. Received: 3 January 1997 / Accepted: 8 July 1998     

Thermal maturity and organic composition of Pennsylvanian coals and carbonaceous shales, north-central Texas: Implications for coalbed gas potential

Thermal maturity was determined for about 120 core, cuttings, and outcrop samples to investigate the potential for coalbed gas resources in Pennsylvanian strata of north-central Texas. Shallow (< 600 m; 2000 ft) coal and carbonaceous shale cuttings samples from the Middle-Upper Pennsylvanian Strawn, Canyon, and Cisco Groups in Archer and Young Counties on the Eastern Shelf of the Midland basin (northwest and downdip from the outcrop) yielded mean random vitrinite reflectance (R_o) values between about 0.4 and 0.8%. This range of R_o values indicates rank from subbituminous C to high volatile A bituminous in the shallow subsurface, which may be sufficient for early thermogenic gas generation. Near-surface (< 100 m; 300 ft) core and outcrop samples of coal from areas of historical underground coal mining in the region yielded similar R_o values of 0.5 to 0.8%. Carbonaceous shale core samples of Lower Pennsylvanian strata (lower Atoka Group) from two deeper wells (samples from ~ 1650 m; 5400 ft) in Jack and western Wise Counties in the western part of the Fort Worth basin yielded higher R_o values of about 1.0%. Pyrolysis and petrographic data for the lower Atoka samples indicate mixed Type II/Type III organic matter, suggesting generated hydrocarbons may be both gas- and oil-prone. In all other samples, organic material is dominated by Type III organic matter (vitrinite), indicating that generated hydrocarbons should be gas-prone. Individual coal beds are thin at outcrop (< 1 m; 3.3 ft), laterally discontinuous, and moderately high in ash yield and sulfur content. A possible analog for coalbed gas potential in the Pennsylvanian section of north-central Texas occurs on the northeast Oklahoma shelf and in the Cherokee basin of southeastern Kansas, where contemporaneous gas-producing coal beds are similar in thickness, quality, and rank.