The 2007 dry spell in Luzon (Philippines): its cause, impact and corresponding response measures

Climate change, involving both human-induced global warming and natural climate variability, has been called upon to explain the occurrences of weather extremes with their associated natural hazards. The Philippines experienced a dry spell in 2007 specifically in parts of Luzon which occurred during the rainy season. On the other hand, areas in Mindanao, southern Philippines which were supposed to be dry, were wet due to the non-migration of the inter-tropical convergence zone northward. The 2007 dry spell in Luzon, northern Philippines affected the agricultural, power, water and health sectors. The local effects of this weather extreme have to be documented and studied to ensure that the appropriate response measures are adopted should there be a recurrence. The building up of the database on this weather extreme and related natural hazards will definitely help the country cope with future similar events.     

Modeling heavy precipitation in complex terrain

Summary A numerical prediction model is described which uses the full set of prognostic equations on a domain roughly the size of the United States with a 96 km horizontal grid resolution and six sigma-coordinate levels. Within this grid resides a nested domain of approximately 1000×1000 km with 24 km horizontal resolution. In this nested grid only modifications to the wind field by the better resolved terrain are considered on the lowest two sigma levels. The terrain effects necessitate adjustments in the location of these two sigma levels. Adjusted wind fields cause modifications in the mass and moisture divergence fields, hence in precipitation. These modifications are averaged into the appropriate meteorological fields on the larger grid.The algorithms used by our model allow continuous interaction between both grids with high computational efficiency.The relative importance of synoptic forcing and terrain is demonstrated for the cases of the Big Thompson, Colorado, flood of 1976 and the Cheyenne, Wyoming, flood of 1985.With 15 Figures     

The atmospheric tide as a continuous spectrum: Lunar semidiurnal tide in surface pressure

Summary The standard equations for the theory of atmospheric tides are solved here by an integral representation on the continuous spectrum of free oscillations. The model profile of back-ground temperature is that of the U.S. Standard Atmosphere in the lower and middle atmosphere, and in the lower thermosphere, above which an isothermal top extends to arbitrarily great heights. The top is warm enough to bring both the Lamb and the Pekeris modes into the continuous spectrum.Computations are made for semidiurnal lunar tidal pressure at sea level at the equator, and the contributions are partitioned according to vertical as well as horizontal structure. Almost all the response is taken up by the Lamb and Pekeris modes of the slowest westward-propagating gravity wave. At sea level, the Lamb-mode response is direct and is relatively insensitive to details of the temperature profile. The Pekeris mode at sea level has an indirect response-in competition with the Lamb mode-and, as has been known since the time of its discovery, it is quite sensitive to the temperature profile, in particular to stratopause temperature. In the standard atmosphere the Lamb mode contributes about +0.078 mb to tidal surface pressure at the equator and the Pekeris mode about –0.048 mb.The aim of this investigation is to illustrate some consequences of representing the tide in terms of the structures of free oscillations. To simplify that task as much as possible, all modifying influences were omitted, such as background wind and ocean or earth tide. Perhaps the main defect of this paper's implementation of the free-oscillation spectrum is that, in contrast to the conventional expansion in the structures of forced oscillations, it does not include dissipation, either implicity or explicity, and thus does not satisfy causality. Dissipation could be added implicity by means of an impedance condition, for example, which would cause up-going energy flux to exceed downgoing flux at the base of the isothermal top layer. To achieve complete causality, however, the dissipation must be modeled explicity. Nevertheless, since the Lamb and Pekeris modes are strongly trapped in the lower and middle atmosphere, where dissipation is rather weak (except possibly in the surface boundary layer), more realistic modeling is not likely to change the broad features of the present results.Symbols a earth's mean radius; expansion coefficient in (5.3) - b recursion variable in (7.4); proximity to resonance in (9.2) - c sound speed in (2.2); specific heatc _p in (2.2) - f Coriolis parameter 2sin in (2.2) - g standard surface gravity - h equivalent depth - i ; discretization index in (7.3) - j index for horizontal structure - k index for horizontal structure; upward unit vectork in (2.2) - m wave number in longitude - n spherical-harmonic degree; number of grid layers in a model layer - p tidal pressure perturbation; background pressurep ₀ - q heating function (energy per mass per time) - r tidal state vector in (2.1) - s tidal entropy perturbation; background entropys ₀ - t time - u tidal horizontal velocityu - w tidal vertical component of velocity - x excitation vector defined in (2.3); vertical coordinate lnp _*/p ₀ [except in (3.8), where it is lnp /p ₀] - y vertical-structure function in (7.1) - z geopotential height - A constant defined in (6.2) - C spherical-harmonic expansion coefficient in (3.6) - D vertical cross section defined in (5.6) and (5.9) - E eigenstate vector - F vertical-structure function for eigenstate pressure in (3.2) [re-defined with WKB scaling in (7.2)] - G vertical-structure function for eigenstate vertical velocity in (3.2) [re-defined with WKB scaling in (7.2)] - H pressure-scale height - I mode intensity defined in (8.1) - K quadratic form defined in (4.4) - L quadratic form defined in (4.4); horizontal-structure magnification factor defined in (5.11) - M vertical-structure magnification factor defined in (4.6) - P eigenstate pressure in (3.2); tidal pressure in (6.2) - R tidal state vector in (5.1) - S eigenstate entropy in (3.2); spherical surface area, in differential dS - T background molecular-scale (NOAA, 1976) absolute temperatureT ₀ - U eigenstate horizontal velocityU in (3.2); coefficient in (7.3) - V horizontal-structure functionV for eigenstate horizontal velocity in (3.2); recursion variable in (7.3) - W eigenstate vertical velocity in (3.2) - X excitation vector in (5.1) - Y surface spherical harmonic in (3.7) - Z Hough function defined in (3.6) - +dH/dz - (1––)/2 - Kronecker delta; Dirac delta; correction operator in (7.6) - equilibrium tide elevation - (square-root of Hough-function eigenvalue) - ratio of specific gas constant to specific heat for air=2/7 - longitude - - - background density ₀ - eigenstate frequency in (3.1) - proxy for heating functionq =c _P/t - latitude - tide frequency - operator for the limitz - horizontal-structure function for eigenstate pressure in (3.2) - Hough function defined in (6.2) - earth's rotation speed - horizontal gradient operator - ()₀ background variable - ()_* surface value of background variable - () value at base of isothermal top layer - Õ state vector with zerow-component - , energy product defined in (2.4) - | | energy norm - ()^* complex conjugate With 10 Figures     

Atmosphere-ocean interaction in mid-latitude storms

Summary Surface fluxes of heat, latent heat, and momentum, and entrainment fluxes and vertical motion at the top of the boundary layer have been calculated for limited regions of several mid-latitude ocean storms. Results have been combined to describe distributions of boundary layer processes which are characteristic of such storms. Surface heat fluxes have important effects in the region west of cold or occluded fronts and are relatively unimportant within a band of about 200 km width east of fronts. Entrainment in pre-frontal regions is driven largely by vertical shear at the top of the boundary layer, while in post-frontal regions it is driven largely by surface heat flux. Boundary layers are well defined in regions more than roughly 200 km east or west of fronts; but closer to fronts boundary layers are not well defined due to the combined effects of entrainment, condensation, and vertical motion associated with the distribution of surface stress.With 12 Figures     

CO2 feedbacks and the 100 K year cycle in climate

Summary In an earlier paper (Lindzen, 1986), it was shown that allowing CO₂ to vary with snow/sea ice position could lead to a greatly enhanced response in glaciation to 100 K year orbital forcing—even when 20 K year forcing was much stronger. In that model, snow/sea ice position (SSIP) and glaciation were different: the former was the forcing for the latter. However, SSIP and glaciation were not decorrelated. Observations (Berner et al., 1979; Lorius et al., 1985; Neftel et al., 1982) suggest that CO₂ may be independently related to both SSIP and glaciation. In the present paper, we allow (in a highly simplified manner) such independent dependence, and show how it alters the earlier results. Briefly, the dependence of CO₂ on glaciation can contribute to and even cause a highly enhanced response to the 100 K year component of the forcing. However, the CO₂ dependence on SSIP is, on the whole, more effective in this regard. Thus, we expect time series of CO₂ to show variation on the faster time scales than does glaciation.With 5 Figures     

Taxonomic sufficiency for soft-bottom sublittoral mollusks assemblages in a tropical estuary, Guanabara Bay, Southeast Brazil

Guanabara Bay (GB) is considered to be one of the most polluted environments of the southern Brazilian coastline. This typical estuarine system is impacted by the heavy discharge of both industrial and domestic waste from the Rio de Janeiro metropolitan area. The mollusc community structure and distribution was investigated between 2000 and 2001, using a three month sampling design of 38 stations, according to austral seasons. Species abundance was aggregated into progressively higher taxa matrices (genus, family, order) and were analysed using multivariate techniques. Mollusc distribution in GB varied significantly in space and time and was probably ruled by the organic enrichment effects of hypoxia and altered redox conditions coupled with prevailing patterns of circulation. Within the sectors of GB an increasing gradient in mollusc diversity and occurrence was observed, ranging from the azoic and impoverished stations in the inner sector to a well-structured community in terms of species composition and abundance inhabiting the outer sector. The non-metric multidimensional scaling (nMDS) and cluster analysis showed similar results when species were aggregated into genera and families, while greater difference occurred at coarser taxonomic identification (order). The literature about taxonomic sufficiency has demonstrated that faunal patterns at different taxonomic levels tend to become similar with increased pollution. In Guanabara Bay, an analysis carried out solely at family level is perfectly adequate to describe the ecophysiological stress. Further aggregation to order level changed the perceived patterns of differences. However, a different taxonomic resolution can be chosen depending on the type of ecological patterns investigated.     

Gravity variations along the Southeast Bohol Ophiolite Complex (SEBOC), Central Philippines: Implications on Ophiolite Emplacement

Abstract The basement complex of Bohol Island consists of the Southeast Bohol Ophiolite Complex (SEBOC), Cansiwang Melange and Alicia Schist. The SEBOC is a complete, but dismembered ophiolite with outcrops generally trending northeast– southwest and dipping north-west. The harzburgite units of the SEBOC are almost always observed to be thrusted onto the Cansiwang Melange, which in turn is thrusted onto the Alicia Schist. Bouguer gravity values on Bohol range from about +60 mGal in the west to +120 mGal in the east, in the region to the north-east of the SEBOC outcrops. Based on the present distribution of the SEBOC units and their thrust fault relationship with the Cansiwang Melange and Alicia Schist, it is proposed that the SEBOC was emplaced by onramping towards the south-eastward direction. However, the orientation of the Bouguer highs suggests that the thrusting direction of the ophiolite units is towards the south-west and not towards the south-east.     

Crustal thickness and adakite occurrence in the Philippines: Is there a relationship?

Adakites are increasingly being recognized worldwide in a variety of tectonic settings. Models on the formation of this geochemically distinct class of volcanic rocks have evolved from partial melting of subducted young, hot oceanic slabs to magmatism resulting from oblique subduction, low‐angle or flat subduction, or even slab‐tearing. Some workers have also pointed to the partial melting of thickened crust to explain the generation of adakitic melts. Rare earth element ratios from adakites and adakitic rocks in the Philippines were used in this study to obtain approximations of the levels where they were generated. These were tied to available geophysical data that defines the crustal thickness of the areas where the samples were collected. High Sm/Yb and La/Yb ratios denote the involvement of amphiboles, and in some cases garnet, in the generation of adakites and adakitic magmas. The presence of amphibole and garnet as residual phases suggests high pressures corresponding to thicker crust (～30 to 45 km). Adakites and adakitic rocks formed through processes other than melting of subducted young oceanic crust would need ≥30 km to account for the heavy rare earth element signatures. If mantle fractionation is not the process involved, crustal thickness is critical to generate adakites and adakitic rocks.     

Surface complexation modeling of Cu(II) adsorption on mixtures of hydrous ferric oxide and kaolinite

Background  

The application of surface complexation models (SCMs) to natural sediments and soils is hindered by a lack of consistent models and data for large suites of metals and minerals of interest. Furthermore, the surface complexation approach has mostly been developed and tested for single solid systems. Few studies have extended the SCM approach to systems containing multiple solids.