Cloud detection and analysis: A review of recent progress

The major types of cloud retrieval algorithms are reviewed with special emphasis being placed upon recent (i.e., post 1981/1982) developments and novel techniques. Satellite-based retrieval algorithms can be grouped into three classes: threshold methods, statistical procedures, and radiative transfer techniques, although each algorithm depends upon implicit, if not overt, inversion of the radiative transfer equation. The fourth type of retrieval differs very considerably from the satellite-based techniques as it depends upon surface-based measurements which are generally, but not always, human rather than computer based. There is a tendency to assume that surface-based observations, especially of total cloud amount, are “correct” but that they may differ from satellite-based retrievals because of the differences in viewing geometry. Actually the literature reveals surprisingly few intercomparison studies. None of the satellite-based techniques have yet been well-validated in a variety of situations, while surface-based observations are made in all terrain and climate regimes. Overall there seems to be good reason to believe that surface-based observations of cloud amount and especially of low cloud amount, character, and base height, can add significantly to satellite-based global nephanalyses now operational or currently being planned.     

A three-dimensional,iterative mapping procedure for the implementation of an ionosphere-magnetosphere anisotropic Ohm’s law boundary condition in global magnetohydrodynamic simulations

The mathematical formulation of an iterative procedure for the numerical implementation of an ionosphere-magnetosphere (IM) anisotropic Ohm’s law boundary condition is presented. The procedure may be used in global magnetohydrodynamic (MHD) simulations of the magnetosphere. The basic form of the boundary condition is well known, but a well-defined, simple, explicit method for implementing it in an MHD code has not been presented previously. The boundary condition relates the ionospheric electric field to the magnetic field-aligned current density driven through the ionosphere by the magnetospheric convection electric field, which is orthogonal to the magnetic field B, and maps down into the ionosphere along equipotential magnetic field lines. The source of this electric field is the flow of the solar wind orthogonal to B. The electric field and current density in the ionosphere are connected through an anisotropic conductivity tensor which involves the Hall, Pedersen, and parallel conductivities. Only the height-integrated Hall and Pedersen conductivities (conductances) appear in the final form of the boundary condition, and are assumed to be known functions of position on the spherical surface R=R₁ representing the boundary between the ionosphere and magnetosphere. The implementation presented consists of an iterative mapping of the electrostatic potential , the gradient of which gives the electric field, and the field-aligned current density between the IM boundary at R=R₁ and the inner boundary of an MHD code which is taken to be at R₂>R₁. Given the field-aligned current density on R=R₂, as computed by the MHD simulation, it is mapped down to R=R₁ where it is used to compute by solving the equation that is the IM Ohm’s law boundary condition. Then is mapped out to R=R₂, where it is used to update the electric field and the component of velocity perpendicular to B. The updated electric field and perpendicular velocity serve as new boundary conditions for the MHD simulation which is then used to compute a new field-aligned current density. This process is iterated at each time step. The required Hall and Pedersen conductances may be determined by any method of choice, and may be specified anew at each time step. In this sense the coupling between the ionosphere and magnetosphere may be taken into account in a self-consistent manner.     

Numerical analyses of deformability tests in jointed rock — “Joint Perturbation” and “No Tension” Finite Element solutions

SummaryNumerical Analyses of Deformability Tests in Jointed Rock — Joint Perturbation and No Tension Finite Element Solutions Sound design of structures in or upon rock requires a thorough knowledge of the rock mass deformability. Field test data interpretation has generally relied upon analytical or Finite Element linear elastic solutions. However, when testing in jointed rock, these can no longer be readily used since the medium cannot resist the tension induced upon loading. Accordingly, other solutions are required. If the discontinuities of the rock mass can all be mapped and their surface properties determined, the Finite Element model will use a Joint Perturbation solution where the joint elements have variable stiffnesses to account for movements along the fractures. If the discontinuities cannot be satisfactorily mapped, a global approach is provided by a No Tension analysis of the stress transfer type. The two techniques are compared in the case of bore hole jack deformability tests and found to agree remarkably well, indicating an appreciable difference from results of linear elasticity.

---

ZusammenfassungNumerische Analyse von Verformungsversuchen in geklüfteten Fels — Lösungen mittels Kluft-Perturbation und No Tension Analysis auf Grund der Methode der Finiten Elemente Der Entwurf von Bauwerken in oder auf Fels erfordert eine genaue Kenntnis des Formänderungsverhaltens der Felsmasse. Im allgemeinen basiert die Interpretation der Ergebnisse von Feldversuchen auf den Lösungen des zugehörigen Randwertproblems der linearen Elastostatik mittels klassischer analytischer Methoden oder der Methode der Finiten Elemente. Wenn die Versuche aber in geklüftetem Fels durchgeführt werden, müssen die theoretischen Methoden modifiziert werden, da das Medium keine Zugspannungen aufnehmen kann. Zwei Fälle können unterschieden werden. Erstens: die Kluftflächen der Felsmasse können geometrisch erfaßt und ihre Oberflächeneigenschaften können bestimmt werden; in diesem Falle kann die Finite-Elemente-Lösung eine Kluft-Perturbations-Methode verwenden, in welcher spezielle Kluft-Elemente verschiedener Steifigkeiten normal und tangential zu den Kluftflächen die Bewegungen entlang derselben bestimmen. Zweitens: die Diskontinuitäten der Felsmasse können nicht hinreichend beschrieben werden; für diesen Fall ist vorgekehrt durch die Finite-Element-Analyse vom Typ der Spannungsumlagerung unter Ausschaltung der Zugspannungen (No Tension Analysis), um das globale Verhalten der geklüfteten Felsmasse zu erfassen. Die Lösungen dieser beiden Methoden wurden für den Fall der Bohrloch-Verformungs-Tests verglichen; der Vergleich zeigt ausgezeichnete Übereinstimmung, im Gegensatz zu den Lösungen als lineares Elastizitätsproblem.

---

RésuméLa méthode des Eléments Finis appliquée aux résultats d'essais de vérins cylindriques dans les roches fissurées — solutions dites Joint Pertubé et Sans Traction La construction de structures sur ou dans des roches nécessite une connaissance approfendie de la déformabilité du massif rocheux. L'analyse des résultats d'essais en place se fait généralement à l'aide de solutions analytiques ou par Eléments Finis, qui supposent une élasticité linéaire du milieu. De fait, dans les roches fissurées, ces solutions sont inadéquates puisque le milieu ne peut résister aux sollicitations de traction. Il s'agit donc d'utiliser d'autres méthodes d'analyse. Lors qu'il est possible de déterminer la géométrie et les propriétés de toutes les discontinuités, l'on pourra utiliser une solution par Eléments Finis, dite Joint Perturbé dans laquelle la rigidité de surface du joint (ou fissure) peut varier lorsque des movements prennent place au niveau du joint. Si cela est impossible, il faudra utiliser une approche globale dite Sans Traction qui redistribue les contraintes de traction et les dissipe. On compare ces deux types de solutions dans le cas des essais de vérins cylindriques unidirectionnels. Leurs résultats sont identiques et ont nettement différents de ceux obtenus avec une solution élastique linéaire.

---

---

Research supported by the National Aeronautics and Space Administration — Contract NAS 8-21432.

---

With 8 Figures     

中国当代土地利用对区域气候影响的数值模拟

使用RegCM3区域气候模式，嵌套欧洲数值预报中心(ECMWF)ERA40再分析资料，分别进行了中国区域在实际植被和理想植被分布情况下各15年时间长度（1987-2001）的积分试验，以研究我国土地利用状况对气候的影响。通过两个试验结果的对比，研究了我国土地利用状况对气候的影响。分析主要集中于气温、降水等的变化上，并对结果进行了统计显著性检验。结果表明，当代土地利用/植被覆盖变化加强了中国地区冬、夏季的季风环流，同时改变了地表能量平衡状况，从而对各气候要素产生重要影响。冬季，植被改变引起长江以南降水减少、气温降低，长江以北降水增加。夏季，植被改变显著影响了南方地区的气候，使得这里降水增多，黄淮、江淮气温降低，华南气温上升；同时引起中国北方降水减少，气温在西北部分植被退化地区升高。植被变化对日最低、最高气温的影响更大。总体来说，土地利用引起了年平均降水在南方增加、北方减少，年平均气温在南方显著降低。     

Sea level,climate change and coastal evolution in Morar,northwest Scotland

Analyses of geomorphologically contrasting sites in Morar, NW Scotland, describe the forcing mechanisms of coastal change. Isolation basins (i.e. basins behind rock sills and now isolated from the sea following isostatic uplift) accumulated continuous marine and freshwater sediments from c.12 to 2 ka BP. Raised dune, marsh and wetland sites register breaching, migration and stability of dunes from c. 9 to 2 ka BP. High-resolution methods designed to address issues of macroscale and microscale sea-level changes and patterns of storminess include 1-mm sampling for pollen, dinocyst and diatom analyses, infra-red photography, X-ray photography and thin-section analysis. The data enhance the record of relative sea-level change for the area. Major phases of landward migration of the coast occurred during the period of low sea-level rise in the mid-Holocene as the rate of rise decreased from c. 3 to < 1 mm/year. Relative sea-level change controls the broad pattern of coastal evolution at each site; local site-specific factors contribute to short-term process change. There is no record of extreme events such as tsunami. Within a system of dynamic metastable equilibrium, the Holocene records show that site-specific factors determine the exact timing of system breakdown, e.g. dune breaching, superimposed on regional sea-level rise. The global average sea-level rise of 3 to 6 mm/yr by AD 2050 predicted by IPCC would only partly be offset in the Morar area by isostatic uplift of about 1 mm/yr. A change from relative sea-level fall to sea-level rise, in areas where the regional rate of uplift no longer offsets global processes, is a critical factor in the management of coastal resources.     

Fault and system stiffnesses and stick-slip phenomena

Summary This paper discusses the influence of system stiffness on the dynamic instability of fault surfaces under laboratory conditions for a number of test modes. In conjunction with shear load stiffness, the normal load stiffness, often neglected, is shown to have a considerable effect on the stick-slip process —its presence or absence and its characteristics. Also appropriate stiffnesses are suggested for an earthquake sequence modeled as a growing dislocation.     

Jet Deflection by a Quasi-Steady-State Side Wind in the Laboratory

We present experimental data on the steady state deflection of a highly supersonic jet by a side-wind in the laboratory. The use of a long interaction region enables internal shocks to fully cross the jet, leading to the development of significantly more structure in the jet than in previous work with a similar setup (Lebedev et al., 2004). The ability to control the length of the interaction region in the laboratory allows the switch between a regime representing a clumpy jet or wind and a regime similar to a slowly varying mass loss rate. The results indicate that multiple internal oblique shocks develop in the jet and the possible formation of a second working surface as the jet attempts to tunnel through the ambient medium.