A numerical ocean circulation model of the Norwegian and Greenland Seas

The dynamics and thermodynamics of the Norwegian and Greenland Seas are investigated using a three-dimensional primitive equation ocean circulation model. The horizontal resolution of the model is 1° in the zonal direction and 0.5° in the meridional direction. The vertical structure is described by 15 levels. The model is driven by both annual mean and seasonally varying wind and thermohaline forcing. The connections of the Norwegian and Greenland Seas with the North Atlantic and Arctic Ocean are modelled with an open boundary condition. The simulated currents are in reasonable agreement with the observed circulation

Polynomial-smoothing and derivative-estimating formulas for functions of one or two independent variables

Application of polynomial-smoothing formulas and related derivative-estimating formulas simplifies certain linear least-squares problems. These problems can then be solved by a computer method

Magic, Mystery, Illusion and Magic

This document vaguely describes what happened to me and my paintings during my time spent at VCU

The ACPI Project, Element 1: Initializing a Coupled Climate Model from Observed Conditions

A problem for climate change studies with coupled ocean-atmosphere models has been how to incorporate observed initial conditions into the ocean, which holds most of the ‘memory’ of anthropogenic forcing effects. The first difficulty is the lack of comprehensive three-dimensional observations of the current ocean temperature (T) and salinity (S) fields to initialize to. The second problem is that directly imposing observed T and S fields into the model results in rapid drift back to the model climatology, with the corresponding loss of the observed information. Anthropogenic forcing scenarios therefore typically initialize future runs by starting with pre-industrial conditions. However, if the future climate depends on the details of the present climate, then initializing the model to observations may provide more accurate forecasts. Also, this ∼130 yr spin up imposes substantial overhead if only a few decades of predictions are desired. A new technique to address these problems is presented. In lieu of observed T and S, assimilated ocean data were used. To reduce model drift, an anomaly coupling scheme was devised. This consists of letting the model’s climatological (pre-industrial) oceanic and atmospheric heat contents and transports balance each other, while adding on the (much smaller) changes in heat content since the pre-industrial era as anomalies. The result is model drift of no more than 0.2 K over 50 years, significantly smaller than the forced response of 1.0 K. An ensemble of runs with these assimilated initial conditions is then compared to a set spun up from pre-industrial conditions. No systematic differences were found, i.e., the model simulation of the ocean temperature structure in the late 1990s is statistically indistinguishable from the assimilated observations. However, a model with a worse representation of the late 20th century climate might show significant differences if initialized in this way.This work was supported by the Department of Energy under grant DE-FG03– 98ER62505

Equatorial currents and transports in the upper central Indian Ocean: Annual cycle and interannual variability

The zonal circulation south of Sri Lanka is an important link for the exchange of water between the Bay of Bengal and the Arabian Sea. Results from a first array of three moorings along 80 degrees 30'E north of 4 degrees 10'N from January .1991 to March 1992 were used to investigate the Monsoon Current regime [Schott et al., 1994]. Measurements from a second array of six current meter moorings are presented here. This array was deployed along 80 degrees 30'E between 45'S and 5 degrees N from July 1993 to September 1994 to investigate the annual cycle and interannual variability of the equatorial currents at this longitude. Both sets of moorings contribute to the Indian Ocean current meter array ICM8 of the World Ocean Circulation Experiment. The semiannual equatorial jet (EJ) was showing a large seasonal asymmetry, reaching a monthly mean eastward transport of 35 Sv (1 Sv = 1 x 10(6) m(3) s(-1)) in November 1993, but just 5 Sv in May 1994. The Equatorial Undercurrent (EUC) had a maximum transport of 17 Sv in March to April 1994. Unexpectedly, compared to previous observations and model studies, the EUC was reappearing again in August 1994 at more than 10 Sv transport and was still flowing when the moorings were recovered. In addition, monthly mean ship drifts near the equator are evaluated to support the interpretation of the moored observations. Interannual variability of the EJ in our measurements and ship drift data appears to be related to the variability of the zonal winds and Southern Oscillation Index. The output of a global numerical model (Parallel Ocean Climate Model) driven by the winds for 1993/1994 is used to connect our observations to the larger scale. The model reproduces the EJ asymmetry and shows the existence of the EUC and its reappearance during summer 1994

Possible Mechanisms of the Exclusion of Johnson Grass by Tall Grass Prairies

Historically, plant distribution typically has been studied with the purpose of learning why a species grows and survives where it does; but why a species does not survive in a particular habitat has rarely been studied, although it may be just as important. According to the US Department of Agriculture, Johnsongrass [Sorghum halepense (L.) Pers.; formerly Johnson grass] is listed as an agricultural pest in most states south of the 42nd parallel. Control of Johnsongrass in agricultural fields involves various labor intensive cultural, mechanical, and chemical means. Release of a bio-control agent has not been suitable for intensively cropped areas. An agriculturally important weed and prominent member of early stage secondary succession, Johnsongrass is not present in later stages of prairie succession. Various environmental factors (biotic and abiotic) that might be involved in restricting Johnsongrass survival were examined in this research. In two sites in Oklahoma, soil conditions were found to be more favorable for survival and growth of Johnsongrass in undisturbed prairie than in the disturbed areas in which Johnsongrass was found vigorously growing. However, even when its rhizomes were introduced into mature prairie, Johnsongrass did not thrive. In laboratory and field trials, presence of the living dominant prairie grasses or leachate from living or dead leaf blades seemed to influence growth and survival of Johnsongrass rhizomes. The prairie grasses, little bluestem [Schizachyrium scoparium (Michx.) Nash] and Indian grass [Sorghastrum nutans (L.) Nash], seem to play a similar allelopathic role in restricting the growth of Johnsongrass to outside of the prairies. Looking at this past study might lead to new methods for the future. (Semtner 2012

Deep currents and the eastward salinity tongue in the equatorial Atlantic: Results from an eddy-resolving, primitive equation model

The high-resolution model of the wind-driven and thermohaline circulation in the Atlantic Ocean developed in recent years as a “community modeling effort” for the World Ocean Circulation Experiment is examined for the temporal and spatial structure of the deep equatorial current field and its effect on the spreading of North Atlantic Deep Water (NADW). Under seasonally varying wind forcing, the model reveals a system of basin-wide zonal currents of O(5 cm s−1), alternating east-west, and oscillating at an annual period. The current fluctuations are induced by the seasonal cycle of the wind stress in the equatorial Atlantic and show characteristics of long equatorial Rossby waves with westward phase propagation of about 15 cm s−1. The mean flow in the deep western tropical Atlantic is governed by a deep western boundary current (DWBC) with core velocities of more than 10 cm s−1. Only a small fraction of the DWBC branches off at the equator, with correspondingly low mean eastward currents of only about 1 cm s−1. Despite this weak advection along the equator, a well-developed salinity tongue is observed in the model, which is reminiscent of observed property distributions at the upper NADW level. The model evaluation indicates the salinity pattern to be a result of a balance between mean zonal advection and meridional diffusion of salt. The presence of the zonal current oscillations appears to have no significance for the existence of the salinity tongue

The winter monsoon circulation of the northern Arabian Sea and Somali current

The winter monsoon circulation in the northern inflow region of the Somali Current is discussed on the basis of an array of moored acoustic Doppler current profiler and current meter stations deployed during 1995–1996 and a ship survey carried out in January 1998. It is found that the westward inflow into the Somali Current regime occurs essentially south of 11°N and that this inflow bifurcates at the Somali coast, with the southward branch supplying the equatorward Somali Current and the northward one returning into the northwestern Arabian Sea. This northward branch partially supplies a shallow outflow through the Socotra Passage between the African continent and the banks of Socotra and partially feeds into eastward recirculation directly along the southern slopes of Socotra. Underneath this shallow surface flow, southwestward undercurrent flows are observed. Undercurrent inflow from the Gulf of Aden through the Socotra Passage occurs between 100 and 1000 m, with its current core at 700–800 m, and is clearly marked by the Red Sea Water (RSW) salinity maximum. The observations suggest that the maximum RSW inflow out of the Gulf of Aden occurs during the winter monsoon season and uses the Socotra Passage as its main route into the Indian Ocean. Westward undercurrent inflow into the Somali Current regime is also observed south of Socotra, but this flow lacks the RSW salinity maximum. Off the Arabian peninsula, eastward boundary flow is observed in the upper 800 m with a compensating westward flow to the south. The observed circulation pattern is qualitatively compared with recent high-resolution numerical model studies and is found to be in basic agreement

3D cut-cell modelling for high-resolution atmospheric simulations

Owing to the recent, rapid development of computer technology, the resolution of atmospheric numerical models has increased substantially. With the use of next-generation supercomputers, atmospheric simulations using horizontal grid intervals of O(100) m or less will gain popularity. At such high resolution more of the steep gradients in mountainous terrain will be resolved, which may result in large truncation errors in those models using terrain-following coordinates. In this study, a new 3D Cartesian coordinate non-hydrostatic atmospheric model is developed. A cut-cell representation of topography based on finite-volume discretization is combined with a cell-merging approach, in which small cut-cells are merged with neighboring cells either vertically or horizontally. In addition, a block-structured mesh-refinement technique is introduced to achieve a variable resolution on the model grid with the finest resolution occurring close to the terrain surface. The model successfully reproduces a flow over a 3D bell-shaped hill that shows a good agreement with the flow predicted by the linear theory. The ability of the model to simulate flows over steep terrain is demonstrated using a hemisphere-shaped hill where the maximum slope angle is resolved at 71 degrees. The advantage of a locally refined grid around a 3D hill, with cut-cells at the terrain surface, is also demonstrated using the hemisphere-shaped hill. The model reproduces smooth mountain waves propagating over varying grid resolution without introducing large errors associated with the change of mesh resolution. At the same time, the model shows a good scalability on a locally refined grid with the use of OpenMP.Comment: 19 pages, 16 figures. Revised version, accepted for publication in QJRM