Grid transformation for incorporating the Arctic in a global ocean model

A grid transformation is described which isolates the Arctic and North Atlantic, rotates the spherical grid to pass an equator up the Atlantic through the north pole and remaps the Coriolis parameter. Boundary condition information is exchanged along the equatorial Atlantic so that the Arctic-Atlantic model is dynamically coupled to a model of the rest of the global ocean (which remains on the geographic spherical grid). The transform produces a more regular grid over the Arctic and eliminates the need for filtering or special treatment at the pole. The transform has been implemented in the GFDL Modular Ocean Model. After testing with idealized geometry, a 300 y global integration is compared to an integration using the geographic spherical grid and Fourier filtering. Results are similar, with differences in the Arctic and western North Atlantic regions leading to smaller air-sea heat flux near the Gulf Stream separation latitude for the transform case. Use of the transform also leads to a reduction in computation time.     

Bistability of the thermohaline circulation identified through comprehensive 2-parameter sweeps of an efficient climate model

The effect of changes in zonal and meridional atmospheric moisture transports on Atlantic overturning is investigated. Zonal transports are considered in terms of net moisture export from the Atlantic sector. Meridional transports are related to the vigour of the global hydrological cycle. The equilibrium thermohaline circulation (THC) simulated with an efficient climate model is strongly dependent on two key parameters that control these transports: an anomaly in the specified Atlantic–Pacific moisture flux (F_a) and atmospheric moisture diffusivity (K_q). In a large ensemble of spinup experiments, the values of F_a and K_q are varied by small increments across wide ranges, to identify sharp transitions of equilibrium THC strength in a 2-parameter space (between Conveyor On and Off states). Final states from this ensemble of simulations are then used as the initial states for further such ensembles. Large differences in THC strength between ensembles, for identical combinations of F_a and K_q, reveal the co-existence of two stable THC states (Conveyor On and Off)—i.e. a bistable regime. In further sensitivity experiments, the model is forced with small, temporary freshwater perturbations to the mid-latitude North Atlantic, to establish the minimum perturbation necessary for irreversible THC collapse in this bistable regime. A threshold is identified in terms of the forcing duration required. The model THC, in a Conveyor On state, irreversibly collapses to a Conveyor Off state under additional freshwater forcing of just 0.1 Sv applied for around 100 years. The irreversible collapse is primarily due to a positive feedback associated with suppressed convection and reduced surface heat loss in the sinking region. Increased atmosphere-to-ocean freshwater flux, under a collapsed Conveyor, plays a secondary role.     

Propagation of the sea-surface temperature anomalies in the Tropical Atlantic

Summary The evolving modes of the sea-surface temperature (SST) in the Tropical Atlantic on the short interannual (IA) timescale were obtained by performing the extended empirical orthogonal function (EEOF) analyses on this variable separately for the 106-year (1871–1976) and 20-year (1881–1900; 1901–1920; 1921–1940; 1941–1960) periods. The equatorial and inter-hemispheric patterns manifest in the first EEOF mode of each analysis as part of the short IA evolution of the SST anomalies in the Tropical Atlantic. Another outstanding feature of the first EEOF mode of each analysis concerns the propagations of the SST anomalies in the meridional direction within the 20°N–20°S band and in the zonal direction in the sector 40°W–20°W. For all analyses, the SST anomalies propagate northward from the equator to 15°N and southward from 20°N to 15°N, with the same sign anomalies merging approximately at 15°N. On the other hand, the SST anomalies propagate westward in the sector 40°W–20°W with a propagation rate close to that of the phase speed of the fastest baroclinic Rossby wave in the ocean. So, the observed propagations of the SST anomalies in the 20°N–20°S band might result from the combined effect of the surface oceanic currents in this band and the baroclinic Rossby waves in the ocean.     

A study of the sensitivity of bare soil evaporation schemes to soil surface wetness,using the coupled soil moisture and surface temperature prediction model,BARESOIL

Summary The performance of evaporation schemes with and approach and their combination within resistance representation of evaporation from bare soil surface is discussed. For this purpose nine schemes, based on different functions of or , on the ratio of the volumetric soil moisture content and its saturated value are used.The quality of the chosen schemes has been evaluated using the results of time integration by the coupled soil moisture and surface temperature prediction model, BARESOIL, using in situ data. A sensitivity analysis was made using two sets of data derived from the volumetric soil moisture content of the top soil layer. One with values below the wilting point (0.17 m³m^–3) and the second with values above 0.20m³m^–3. Data sets were obtained at the experimental site Rimski anevi, Yugoslavia, from the bare surface of a chernozem soil.With 4 Figures     

Seasonal variations of anthropogenic sulfate aerosol and direct radiative forcing over China

Summary Regional climate model (RegCM2) and sulfur transport model (NJUADMS) were combined to simulate the distribution of anthropogenic sulfate aerosol burden over China, where a look up table method was applied to illustrate sulfate formation from SO₂-oxidation. Direct radiative forcing of sulfate aerosol was further estimated using the scheme suggested by Charlson et al (1991). Investigations show that the annual average total sulfate column over mainland China is 2.01mg/m² with high value in East and Central areas (more than 7mg/m²). The annual average direct radiative forcing of China is about –0.85W/m². The forcing can reach –7W/m² in Central and East China during the winter season. Total sulfate column shows significant seasonal variations with winter maximum-summer minimum in the Southern part of China and spring maximum-autumn minimum in the northern part of China. Strong seasonal cycles of direct radiative forcing are also found due to the influence of total sulfate column, cloud, relative humidity and the reflectivity of underlying surfaceReceived May 16, 2001; accepted August 5, 2002 Published online: May 8, 2003     

Re-analysis of the nuclear winter phenomenon

Summary An analysis of the report of the (U.S.) National Academy of Sciences (NAS) on atmospheric effects of a nuclear exchange leads to conclusions that differ from those of the NAS and of the earlier TTAPS and AMBIO studies. Any cooling of the earth's surface is likely to beshort-lived because of rapid removal of the smoke clouds originating from nuclear burst-initiated fires, andminor because of appreciable green-house effects due to several distinct physical causes. (One of these, neglected in prior analyses, is the infrared absorption from cirrus clouds produced directly by the nuclear bursts.) Taken together, these effects may even induce slight surface warming (nuclear summer) instead of cooling (nuclear winter). The consequences to atmospheric ozone are similarly ambiguous; depending on the detailed nuclear scenario, the net ozone content may increase-rather than decrease as argued by TTAPS. Experiments could settle some uncertainties.With 2 Figures     

Trace gas exchange at the air/water interface: Measurements of mass accommodation coefficients

A liquid jet of 90 m diameter and variable length has been utilized to determine absorption rates and, hence, mass accommodation coefficients , of atmospheric trace gases. The compounds investigated are HCl (0.01), HNO₃ (0.01), N₂O₅ (0.005), peroxyacetyl nitrate (>0.001), and HONO (0.005). It is concluded that the absorption of these trace gases by liquid atmospheric water is not significantly retarded by interfacial mass transport. The strengths and limitations of the liquid jet technique for measuring mass accommodation coefficients are explored.     

Application of a “Big-Tree” model to regional climate modeling: a sensitivity study

Summary ¶In order to better understand land-atmosphere interactions and increase the predictability of climate models, it is important to investigate the role of forest representation in climate modeling. Corresponding to the big-leaf model commonly employed in land surface schemes to represent the effects of a forest, a so called big-tree model, which uses multi-layer vegetation to represent the vertical canopy heterogeneity, was introduced and incorporated into the National Center for Atmospheric Research (NCAR) regional climate model RegCM2, to make the vegetation model more physically based. Using this augmented RegCM2 and station data for China during 1991 Meiyu season, we performed 10 experiments to investigate the effects of the application of the big-tree model on the summer monsoon climate.With the big-tree model incorporated into the regional climate model, some climate characteristics, e.g. the 3-month-mean surface temperature, circulation, and precipitation, are significantly and systematically changed over the model domain, and the change of the characteristics differs depending on the area. Due to the better representation of the shading effect in the big-tree model, the temperature of the lower layer atmosphere above the plant canopy is increased, which further influences the 850hPa temperature. In addition, there are significant decreases in the mean latent heat fluxes (within 20–30W/m²) in the three areas of the model domain.The application of the big-tree model influences not only the simulated climate of the forested area, but also that of the whole model domain, and its impact is greater on the lower atmosphere than on the upper atmosphere. The simulated rainfall and surface temperature deviate from the originally simulated result and are (or seem to be) closer to the observations, which implies that an appropriate representation of the big-tree model may improve the simulation of the summer monsoon climate.We also find that the simulated climate is sensitive to some big-tree parameter values and schemes, such as the shape, height, zero-plane displacement height and mixing-length scheme. The simulated local/grid differences may be very large although the simulated areal-average differences may be much lower. The area-average differences in the monthly-mean surface temperature and heat fluxes can amount to 0.5°C and 4W/m², respectively, which correspond to maximum local/grid differences of 3.0°C and 40W/m² respectively. It seems that the simulated climate is most sensitive to the parameter of the zero-plane displacement among the parameters studied.     

Arctic hazes in summer over Greenland and the North American Arctic: II. Nature and concentrations of accumulation-mode and giant particles

Airborne measurements made during August 1985 over Greenland and its environs show that both accumulation-mode (0.1 m D2.0 m) and giant (D2 m) particles were present in relatively high concentrations in arctic haze layers and that the accumulation-mode particles dominated light scattering. Particles with diameters (D) between 1 and 4 m consisted predominately of mixed materials, small and dense inclusions, and probably organic compounds containing sulfur. Many of the particles from 0.1 to 1 m in diameter were also of mixed composition, with sulfuric acid, ammonium sulfate and organics probably the dominant constituents.     

Numerical modeling of the nocturnal urban boundary layer

A numerical case study with a second-order turbulence closure model is proposed to study the role of urban canopy layer (UCL) for the formation of the nocturnal urban boundary layer (UBL). The turbulent diffusion coefficient was determined from an algebraic stress model. The concept of urban building surface area density is proposed to represent the UCL. Calculated results were also compared with field observation data. The height of the elevated inversion above an urban center was simulated and found to be approximately twice the average building height. The turbulent kinetic energy k, energy dissipation rate , and turbulence intensities u ² and w ² increase rapidly at the upwind edge of the urban area. The Reynolds stress uw displayed a nearly uniform profile inside the UBL, and the vertical sensible heat flux w had a negative value at the inversion base height. This indicates that the downward transport of sensible heat from the inversion base may play an important role in the formation of the nocturnal UBL.     

The validation of various schemes for parameterizing evaporation from bare soil for use in meteorological models: A numerical study using in situ data

Parameterization of evaporation from a non-plant-covered surface is very important in the hierarchy strategy of modelling land surface processes. One of the representations frequently used in its computation is the resistance formulation. The performance of the evaporation schemes using the , , and their combination resistance approaches to parameterize evaporation from bare soil surfaces is discussed. For that purpose, the nine schemes, based on a different dependence of and on volumetric soil moisture content and its saturated value, are used.The tests of performances of the considered schemes are based on time integrations by the land surface module (BARESOIL) using observed data. The 23 data sets at a bare surface experimental site in Rimski anevi, Yugoslavia on chernozem soil, were used for the resistance algorithm evaluation. The quality of the schemes was compared with the observed values of the latent heat flux using several statistical parameters.     

ALPEX-Simulation

Summary In a project ALPEX-Simulation, sponsored by the Österreichischer Fond zur Förderung der wissenschaftlichen Forschung (FWF), all eight cases of ALPEX-SOP cyclones were numerically simulated with a fine mesh isentropic model of the atmosphere. These numerical simulations in six-hourly intervals allow a deeper insight into the synoptics and dynamics of the cyclogeneses in the Western Mediterranean, especially into the genesis of the two basic types of cyclones: the so-called Überströmungs-type and Vorderseiten-type. In the first phase of cyclogenesis of the Überströmungs-type, the blocking and flow splitting of the cold air due to the Alps and the canalization between the Alps and the Massif Central are important. Cold air flows cyclonically around the western part of the Alps, creating a vorticity maximum at the south western edge of the Alpine, bow and leads also to an enhanced PV. In connection with warm air in the Mediterranean, a strong baroclinic zone is generated. The interaction between the arriving PV maximum in the upper troposphere and the enhanced PV at the bottom leads to cyclogenesis in the Western Mediterranean. In the case of the Vorderseiten-type warm air advection dominates with the exception of a shallow layer of cold air in the inner Po-Valley, which is shielded by the Alpine ridge. A well-pronounced PV maximum builds up and couples with the PV maximum arriving at upper levels, even before the cold air, coming from the north-west, has surrounded the Alps. The cold air only intensifies the development by raising the baroclinity. Therefore, the Vorderseiten-cyclogenesis is an orographically modified cyclogenesis, in the course of which the cyclonic development is triggered by the Alps, whereas the Überströmungs-cyclogenesis is an orographically induced cyclogenesis i.e. a true lee cyclogenesis.With 14 FiguresDied in a tragic traffic accident on June 6, 1993.     

Bioaccumulation of Cobalt in Parmelia sulcata

Six locations across mainland Portugal were selected for exposing Parmelia sulcata, for a one-year period (8 months for one site), with simultaneous measurement of total (dry + wet) deposition (one-month periods). The exposed lichens and the total (dry + wet) deposition were analysed for cobalt contents by INAA (instrumental neutron activation analysis) and ICP-MS (inductively coupled plasma mass spectroscopy), respectively. The designated wet deposition was evaluated through the collected water volume; the designated dry deposition was assessed after the (dried) residual mass of the wet deposition. An excellent agreement between Co contents in exposed lichens and the cumulative (1) Co contents in the dry deposition, (2) dry deposition, and (3) wet deposition has been found for the locations with alternate drought and precipitation months, high dry deposition, and high Co contents in the latter. Continuous rainfall was found to hinder the Co accumulation in the lichen due to its release from the lichen and/or lower Co contents in the dry deposition. At three locations, P. sulcata Co contents, after subtraction of the background (before exposure), equalled or exceeded the Co contents in the cumulative dry deposition at the end of the exposure time. The optimal exposure period for this species likely depends on the exposure conditions.     

Orbital forcing of Arctic climate: mechanisms of climate response and implications for continental glaciation

Progress in understanding how terrestrial ice volume is linked to Earths orbital configuration has been impeded by the cost of simulating climate system processes relevant to glaciation over orbital time scales (10³–10⁵ years). A compromise is usually made to represent the climate system by models that are averaged over one or more spatial dimensions or by three-dimensional models that are limited to simulating particular snapshots in time. We take advantage of the short equilibration time (10 years) of a climate model consisting of a three-dimensional atmosphere coupled to a simple slab ocean to derive the equilibrium climate response to accelerated variations in Earths orbital configuration over the past 165,000 years. Prominent decreases in ice melt and increases in snowfall are simulated during three time intervals near 26, 73, and 117 thousand years ago (ka) when aphelion was in late spring and obliquity was low. There were also significant decreases in ice melt and increases in snowfall near 97 and 142 ka when eccentricity was relatively large, aphelion was in late spring, and obliquity was high or near its long term mean. These glaciation-friendly time intervals correspond to prominent and secondary phases of terrestrial ice growth seen within the marine ¹⁸O record. Both dynamical and thermal effects contribute to the increases in snowfall during these periods, through increases in storm activity and the fraction of precipitation falling as snow. The majority of the mid- to high latitude response to orbital forcing is organized by the properties of sea ice, through its influence on radiative feedbacks that nearly double the size of the orbital forcing as well as its influence on the seasonal evolution of the latitudinal temperature gradient.     

Hybrid modeling in meteorological applications Part I: Concepts and approaches

Summary Hybrid modeling entails the combination of a numerical weather prediction model and a symbolic model. The latter uses symbolic objects, their characterizing attributes and sets of behavioral constrains which prescribe changes in the states of these objects as functions of time, space, and other imposed quantitative or heuristic conditions. Integration of these two modeling components for an on-line, real-time, operational system is feasible only if both the numerical and the symbolic model can be executed in a distributed mode, i.e. at a user's location rather than in a central weather service office. This condition entails the design of a numerical model that can run on relatively inexpensive desktop workstations or high-end PCs. Given such a capability, the output from the numerical model can be used to satisfy a number of behavioral constraints of objects (such as thunderstorm, blizzard, etc.), defined in the symbolic model. These constraints can be embedded invisibly as functions of time and pixel location on the computer screen, to be called upon as soon as the respective object is activated, e.g. by placing an icon on the screen.To make such a hybrid weather prediction model responsive to details in topography, it will have to be able to interface with a geographic information system (GIS) database. Since such databases can be very voluminous, management procedures for indexing and rapid information retrieval have to be instituted. The approach discussed here involves restructuring of given GIS data into B+ trees.The hybrid prediction model whose design is described in this paper, executes very quickly on a PC (e.g. a 33 megahertz Intel 80486 chip based machine). It allows assimilation of locally generated observational data to improve forecast quality, and can respond to queries of a highly specialized nature in support of tactical decisions within the time frame between nowcasting (3 hours) and 24 hours.With 6 Figures     

The role of the Atlantic freshwater balance in the hysteresis of the meridional overturning circulation

We have studied the response of the Atlantic meridional overturning circulation to surface freshwater forcing using an ocean GCM coupled to an energy-moisture-balance atmosphere model. The overturning collapses rapidly when a slowly increasing forcing applied to the North Atlantic passes a positive threshold, and spins up equally quickly when the forcing falls below a negative threshold. This well-known behaviour is referred to as hysteresis because the thresholds in forcing are different for the transitions in opposite directions. However, we argue that the behaviour of the Atlantic salinity is more fundamental than the forcing. Hysteresis as a function of freshwater forcing occurs because the states with North Atlantic overturning on and off each tend to reinforce their associated salinity distributions and inhibit the transition to the other state. During the collapse, the Atlantic becomes less saline because of the import of 80 Sv year of freshwater by ocean transports across 30°S; during the spin-up this freshwater is exported again. We show that qualitatively similar hysteresis behaviour can be produced by perturbing the system without any net freshwater forcing. The salinity flip-flop is associated with the appearance and disappearance of a shallow reverse overturning circulation south of the Equator, which is present while the northern overturning is absent, and may provide the mechanism for the ocean freshwater influx during collapse.     

On the turbulence structure in the stable boundary layer over the Greenland ice sheet

Turbulence measurements from a 30 m tower in the stably stratifiedboundary layer over the Greenland ice sheet are analyzed. The observationsinclude profile and eddy-correlation measurements at various levels. Atfirst, the analysis of the turbulence data from the lowest level (2 m aboveground) shows that the linear form of the non-dimensional wind profile(m) is in good agreement with the observations for z/L <0.4, whereL represents the Obukhov length. A linear regression yieldsm=1+5.8z/L. The non-dimensional temperature profile (h) at the2m level shows no tendency to increase with increasing stability. The datafrom the upper levels of the tower are analyzed in terms of both localscaling and surface-layer scaling. The m and the h values show atendency to level off at large stability (z/>0.4) where represents the local Obukhov length. Hence, the linear form of the functions is no longer appropriate under such conditions. The bestcorrespondence to the data can be achieved when using the expression ofBeljaars and Holtslag for m and h. The vertical profiles of theturbulent fluxes, the wind velocity variances and temperature variance arealso determined. The momentum flux profile and the profiles of wind speedvariances are in general agreement with other observations if a welldeveloped low-level wind maximum occurs, and the height of this maximum isused as a height scale.     

Effect of land management on ecosystem carbon fluxes at a subalpine grassland site in the Swiss Alps

Summary The influence of agricultural management on the CO₂ budget of a typical subalpine grassland was investigated at the Swiss CARBOMONT site at Rigi-Seebodenalp (1025m a.s.l.) in Central Switzerland. Eddy covariance flux measurements obtained during the first growing season from the mid of spring until the first snow fall (17 Mai to 25 September 2002) are reported. With respect to the 10-year average 1992–2001, we found that this growing season had started 10 days earlier than normal, but was close to average temperature with above-normal precipitation (100–255% depending on month). Using a footprint model we found that a simple approach using wind direction sectors was adequate to classify our CO₂ fluxes as being controlled by either meadow or pasture. Two significantly different light response curves could be determined: one for periods with external interventions (grass cutting, cattle grazing) and the other for periods without external interventions. Other than this, meadow and pasture were similar, with a net carbon gain of –128±17g Cm^–2 on the undisturbed meadow, and a net carbon loss of 79±17g Cm^–2 on the managed meadow, and 270±24g Cm^–2 on the pasture during 131 days of the growing season, respectively. The grass cut in June reduced the gross CO₂ uptake of the meadow by 50±2% until regrowth of the vegetation. Cattle grazing reduced gross uptake over the whole vegetation period (37±2%), but left respiration at a similar level as observed in the meadow.     

Transport-Dissipation Analytical Solutions to the E-∈Turbulence Model and their Role in Predictions of the Neutral ABL

E- turbulence model predictions of the neutralatmospheric boundary layer (NABL) are reinvestigated to determine thecause for turbulence overpredictions found in previous applications. Analytical solutions to the coupled E and equations for the case of steady balance between transport and dissipation terms, the dominant balance just below the NABL top, are derived. It is found that analytical turbulence profiles laminarizeat a finite height only for values of closure parameter ratio c ₂ /_e equal toor slightly greater than one, with laminarization as z for greater . The point = 2 is additionally foundthat where analytical turbulent length scale (l) profilesmade a transition from ones ofdecreasing ( < 2) to increasing ( > 2)values with height. Numerically predicted profiles near the NABL topare consistent with analytical findings. The height-increasingvalues of l predicted throughout the NABL with standard values ofclosure parameters thus appear a consequence of 2.5(> 2), implied by these values (c ₂ = 1.92, _{= 1.3}, _{e = 1}). Comparison of numericalpredictions with DNS data shows that turbulence overpredictions obtained with standard-valued parameters are rectifiedby resetting and _e to 1.1 and 1.6, respectively, giving, with c ₂ = 1.92, 1.3, and laminarization of the NABL's cappingtransport-dissipation region at a finite height.     

Nitrate in the atmospheric boundary layer of the tropical South Pacific: Implications regarding sources and transport

Weekly bulk aerosol samples collected at Funafuti, Tuvalu (8°30S, 179°12E), American Samoa (14°15S, 170°35W), and Rarotonga (21°15S, 159°45W), from 1983 through most of 1987 have been analyzed for nitrate and other constituents. The mean nitrate concentration is about 0.11 g m^–3 at each of these stations: 0.107±0.011 g m^–3 at Funafuti; 0.116±0.008 at American Samoa; and 0.117±0.010 at Rarotonga. Previous measurements of mineral aerosol and trace metal concentrations at American Samoa are among the lowest ever recorded for the near-surface troposphere and indicate that this region is minimally affected by transport of soil material and pollutants from the continents. Consequently, the nitrate concentration of 0.11 g m^–3 can be regarded as the natural level for the remote marine boundary layer of the tropical South Pacific Ocean. In contrast, over the tropical North Pacific which is significantly impacted by the transport of material from Asia and North America, the mean nitrate concentrations are about three times higher, 0.29 and 0.36 g m^–3 at Midway and Oahu, respectively. The major sources of the nitrate over the tropical South Pacific are still very uncertain. A very significant correlation between the nitrate concentrations at American Samoa and the concentrations of ²¹⁰Pb suggests that transport from continental sources might be important. This continental source could be lightning, which occurs most frequently over the tropical continents. A near-zero correlation with ⁷Be indicates that the stratosphere and upper troposphere are probably not the major sources. A significant biogenic source would be consistent with the higher mean nitrate concentrations, 0.16 to 0.17 g m^–3, found over the equatorial Pacific at Fanning Island (3°55N, 159°20W) and Nauru (0°32S, 166°57E). The lack of correlation between nitrate and nss sulfate at American Samoa does not necessarily preclude an important role for marine biogenic sources.