Large-Eddy Simulation of Flows over Random Urban-like Obstacles

Further to our previous large-eddy simulation (LES) of flow over a staggered array of uniform cubes, a simulation of flow over random urban-like obstacles is presented. To gain a deeper insight into the effects of randomness in the obstacle topology, the current results, e.g. spatially-averaged mean velocity, Reynolds stresses, turbulence kinetic energy and dispersive stresses, are compared with our previous LES data and direct numerical simulation data of flow over uniform cubes. Significantly different features in the turbulence statistics are observed within and immediately above the canopy, although there are some similarities in the spatially-averaged statistics. It is also found that the relatively high pressures on the tallest buildings generate contributions to the total surface drag that are far in excess of their proportionate frontal area within the array. Details of the turbulence characteristics (like the stress anisotropy) are compared with those in regular roughness arrays and attempts to find some generality in the turbulence statistics within the canopy region are discussed.     

Unique isoprene oxidation products demonstrate chlorine atom chemistry occurs in the Houston,Texas urban area

As part of the 2000 Texas Air Quality Study (TexAQS), we studied the isoprene oxidation process under ambient conditions to discern the presence of chlorine atom (Cl) chemistry in the Houston, Texas urban area. By measuring chloromethylbutenone (CMBO) and an isomer of chloromethylbutenal (CMBA), we clearly observed sixteen episodes of active Cl chemistry during the 24-day experiment. Estimated median Cl concentration during each of these episodes was between the detection limit of ~10² atoms cm⁻³ and 50 _{- 30} ^{+ 70} ×10⁴ 50_{ - 30}^{ + 70} \times {10^4} atoms cm⁻³. Cl concentration during all the episodes averaged 7.6 _{- 2.0} ^{+ 4.7} ×10⁴ 7.6_{ - 2.0}^{ + 4.7} \times {10^4} atoms cm⁻³ and thus amounted to less than 3% of the OH concentration during the same periods. During the episodes, the fraction of oxidation chemistry initiated by Cl ranged from 3–43% and was strongly dependent on the quantity and type of hydrocarbons present in the atmosphere. Because of its intermittent presence and low concentration, Cl is not a broadly influential oxidant in the Houston, Texas urban area.     

Current status of ENSO prediction skill in coupled ocean–atmosphere models

The overall skill of ENSO prediction in retrospective forecasts made with ten different coupled GCMs is investigated. The coupled GCM datasets of the APCC/CliPAS and DEMETER projects are used for four seasons in the common 22 years from 1980 to 2001. As a baseline, a dynamic-statistical SST forecast and persistence are compared. Our study focuses on the tropical Pacific SST, especially by analyzing the NINO34 index. In coupled models, the accuracy of the simulated variability is related to the accuracy of the simulated mean state. Almost all models have problems in simulating the mean and mean annual cycle of SST, in spite of the positive influence of realistic initial conditions. As a result, the simulation of the interannual SST variability is also far from perfect in most coupled models. With increasing lead time, this discrepancy gets worse. As one measure of forecast skill, the tier-1 multi-model ensemble (MME) forecasts of NINO3.4 SST have an anomaly correlation coefficient of 0.86 at the month 6. This is higher than that of any individual model as well as both forecasts based on persistence and those made with the dynamic-statistical model. The forecast skill of individual models and the MME depends strongly on season, ENSO phase, and ENSO intensity. A stronger El Niño is better predicted. The growth phases of both the warm and cold events are better predicted than the corresponding decaying phases. ENSO-neutral periods are far worse predicted than warm or cold events. The skill of forecasts that start in February or May drops faster than that of forecasts that start in August or November. This behavior, often termed the spring predictability barrier, is in part because predictions starting from February or May contain more events in the decaying phase of ENSO.     

Modeled seasonality of glacial abrupt climate events

Greenland ice cores, as well as many other paleo-archives from the northern hemisphere, recorded a series of 25 warm interstadial events, the so-called Dansgaard-Oeschger (D-O) events, during the last glacial period. We use the three-dimensional coupled global ocean–atmosphere–sea ice model ECBILT-CLIO and force it with freshwater input into the North Atlantic to simulate abrupt glacial climate events, which we use as analogues for D-O events. We focus our analysis on the Northern Hemisphere. The simulated events show large differences in the regional and seasonal distribution of the temperature and precipitation changes. While the temperature changes in high northern latitudes and in the North Atlantic region are dominated by winter changes, the largest temperature increases in most other land regions are seen in spring. Smallest changes over land are found during the summer months. Our model simulations also demonstrate that the temperature and precipitation change patterns for different intensifications of the Atlantic meridional overturning circulation are not linear. The extent of the transitions varies, and local non-linearities influence the amplitude of the annual mean response as well as the response in different seasons. Implications for the interpretation of paleo-records are discussed.     

气候变化对中国水资源影响的适应性评估与管理框架

 通过论述气候变化对中国水资源影响的适应性评估与管理框架,提出一个气候变化影响决策评估工具,它包括：未来气候变化对中国水资源潜在影响的定性描述分析、半定量与定量分析以及适应性对策评估。由于不同气候区域所面临的水资源问题不同,选择中国4个典型案例区域,并确定不同的目标进行气候变化适应性管理综合研究,提出了甄别气候变化影响和适应性管理的新的思路、框架与方法论。该项研究为应对未来气候变化影响的水资源规划与风险管理提供了途径与方法。     

Dynamical passage to approximate equilibrium shapes for spinning, gravitating rubble asteroids

Many asteroids are thought to be particle aggregates held together principally by self-gravity. Here we study — for static and dynamical situations — the equilibrium shapes of spinning asteroids that are permitted for rubble piles. As in the case of spinning fluid masses, not all shapes are compatible with a granular rheology. We take the asteroid to always be an ellipsoid with an interior modeled as a rigid-plastic, cohesion-less material with a Drucker-Prager yield criterion. Using an approximate volume-averaged procedure, based on the classical method of moments, we investigate the dynamical process by which such objects may achieve equilibrium. We first collapse our dynamical approach to its statical limit to derive regions in spin-shape parameter space that allow equilibrium solutions to exist. At present, only a graphical illustration of these solutions for a prolate ellipsoid following the Drucker-Prager failure law is available [Sharma, I., Jenkins, J.T., Burns, J.A., 2005a. Bull. Am. Astron. Soc. 37, 643; Sharma, I., Jenkins, J.T., Burns, J.A., 2005b. Equilibrium shapes of ellipsoidal soil asteroids. In: García-Rojo, R., Hermann, H.J., McNamara, S. (Eds.), Proceedings of the 5th International Conference on Micromechanics of Granular Media, vol. 1. A.A. Balkema, UK; Holsapple, K.A., 2007. Icarus 187, 500-509]. Here, we obtain the equilibrium landscapes for general triaxial ellipsoids, as well as provide the requisite governing formulae. In addition, we demonstrate that it may be possible to better interpret the results of Richardson et al. [Richardson, D.C., Elankumaran, P., Sanderson, R.E., 2005. Icarus 173, 349-361] within the context of a Drucker-Prager material. The graphical result for prolate ellipsoids in the static limit is the same as those of Holsapple [Holsapple, K.A., 2007. Icarus 187, 500-509] because, when worked out, his final equations will match ours. This is because, though the formalisms to reach these expressions differ, in statics, at the lowest level of approximation, volume-averaging and the approach of Holsapple [Holsapple, K.A., 2007. Icarus 187, 500-509] coincide. We note that the approach applied here was obtained independently [Sharma, I., Jenkins, J.T., Burns, J.A., 2003. Bull. Am. Astron. Soc. 35, 1034; Sharma, I., 2004. Rotational Dynamics of Deformable Ellipsoids with Applications to Asteroids. Ph.D. thesis, Cornell University] and it provides a general, though approximate, framework that is amenable to systematic improvements and is flexible enough to incorporate the dynamical effects of a changing shape, different rheologies and complex rotational histories. To demonstrate our technique, we investigate the non-equilibrium dynamics of rigid-plastic, spinning, prolate asteroids to examine the simultaneous histories of shape and spin rate for rubble piles. We have succeeded in recovering most results of Richardson et al. [Richardson, D.C., Elankumaran, P., Sanderson, R.E., 2005. Icarus 173, 349-361], who obtained equilibrium shapes by studying numerically the passage into equilibrium of aggregates containing discrete, interacting, frictionless, spherical particles. Our mainly analytical approach aids in understanding and quantifying previous numerical simulations.     

Cosmography with cluster strong lensing

By stacking an ensemble of strong lensing clusters, we demonstrate the feasibility of placing constraints on the dark energy equation of state. This is achieved by using multiple images of sources at two or more distinct redshift planes. The sample of smooth clusters in our simulations is based on observations of massive clusters and the distribution of background galaxies is constructed using the Hubble Deep Field . Our source distribution reproduces the observed redshift distribution of multiply imaged sources in Abell 1689. The cosmology recovery depends on the number of image families with known spectroscopic redshifts and the number of stacked clusters. Our simulations suggest that constraints comparable to those derived from other competing established techniques on a constant dark energy equation of state can be obtained using 10–40 clusters with five or more families of multiple images. We have also studied the observational errors in the image redshifts and positions. We find that spectroscopic redshifts and high-resolution Hubble Space Telescope ( HST ) images are required to eliminate confidence contour relaxation relative to the ideal case in our simulations. This suggests that the dark energy equation of state, and other cosmological parameters, can be constrained with existing HST images of lensing clusters coupled with dedicated ground-based arc spectroscopy.     

Chemical evolution with radial mixing

Models of the chemical evolution of our Galaxy are extended to include radial migration of stars and flow of gas through the disc. The models track the production of both iron and α-elements. A model is chosen that provides an excellent fit to the metallicity distribution of stars in the Geneva–Copenhagen survey (GCS) of the solar neighbourhood and a good fit to the local Hess diagram. The model provides a good fit to the distribution of GCS stars in the age–metallicity plane, although this plane was not used in the fitting process. Although this model's star formation rate is monotonically declining, its disc naturally splits into an α-enhanced thick disc and a normal thin disc. In particular, the model's distribution of stars in the ([O/Fe], [Fe/H]) plane resembles that of Galactic stars in displaying a ridge line for each disc. The thin-disc's ridge line is entirely due to stellar migration, and there is the characteristic variation of stellar angular momentum along it that has been noted by Haywood in survey data. Radial mixing of stellar populations with high  σ _z   from inner regions of the disc to the solar neighbourhood provides a natural explanation of why measurements yield a steeper increase of  σ _z   with age than predicted by theory. The metallicity gradient in the interstellar medium is predicted to be steeper than in earlier models, but appears to be in good agreement with data for both our Galaxy and external galaxies. The models are inconsistent with a cut-off in the star formation rate at low gas surface densities. The absolute magnitude of the disc is given as a function of time in several photometric bands, and radial colour profiles are plotted for representative times.     

Annual cycles in the interstellar scintillation time-scales of PKS B1519–273 and PKS B1622–253

We have used the University of Tasmania's 30-m radio telescope at Ceduna in South Australia to regularly monitor the flux density of a number of southern blazars. We report the detection of an annual cycle in the variability time-scale of the centimetre radio emission of PKS B1622−253. Observations of PKS B1519−273 over a period of nearly 2 yr confirm the presence of an annual cycle in the variability time-scale in that source. These observations prove that interstellar scintillation is the principal cause of inter-day variability at radio wavelengths in these sources. The best-fitting annual cycle model for both sources implies a high degree of anisotropy in the scattering screen and that it has a large velocity offset with respect to the local standard of rest. This is consistent with a greater screen distance for these 'slow' intra-day variability (IDV) sources than for rapid scintillators such as PKS B0405−385 or J1819+3845.