Improving a convection-permitting model simulation of a cold air outbreak

A convection-permitting local-area model was used to simulate a cold air outbreak crossing from the Norwegian Sea into the Atlantic Ocean near Scotland. A control model run based on an operational configuration of the Met Office UKV high-resolution (1.5 km grid spacing) NWP model was compared to satellite, aircraft and radar data. While the control model captured the large-scale features of the synoptic situation, it was not able to reproduce the shallow (<1.5 km) stratiform layer to the north of the open cellular convection. Liquid water paths were found to be too low in both the stratiform and convective cloud regions. Sensitivity analyses including a modified boundary-layer diagnosis to generate a more well-mixed boundary layer and inhibition of ice formation to lower temperatures improved cloud morphology and comparisons with observational data. Copyright © 2013 Royal Meteorological Society and British Crown Copyright, the Met Office

An evaluation of ocean color model estimates of marine primary productivity in coastal and pelagic regions across the globe

Nearly half of the earth's photosynthetically fixed carbon derives from the oceans. To determine global and region specific rates, we rely on models that estimate marine net primary productivity (NPP) thus it is essential that these models are evaluated to determine their accuracy. Here we assessed the skill of 21 ocean color models by comparing their estimates of depth-integrated NPP to 1156 in situ C-14 measurements encompassing ten marine regions including the Sargasso Sea, pelagic North Atlantic, coastal Northeast Atlantic, Black Sea, Mediterranean Sea, Arabian Sea, subtropical North Pacific, Ross Sea, West Antarctic Peninsula, and the Antarctic Polar Frontal Zone. Average model skill, as determined by root-mean square difference calculations, was lowest in the Black and Mediterranean Seas, highest in the pelagic North Atlantic and the Antarctic Polar Frontal Zone, and intermediate in the other six regions. The maximum fraction of model skill that may be attributable to uncertainties in both the input variables and in situ NPP measurements was nearly 72%. on average, the simplest depth/wavelength integrated models performed no worse than the more complex depth/wavelength resolved models. Ocean color models were not highly challenged in extreme conditions of surface chlorophyll-a and sea surface temperature, nor in high-nitrate low-chlorophyll waters. Water column depth was the primary influence on ocean color model performance such that average skill was significantly higher at depths greater than 250 m, suggesting that ocean color models are more challenged in Case-2 waters (coastal) than in Case-1 (pelagic) waters. Given that in situ chlorophyll-a data was used as input data, algorithm improvement is required to eliminate the poor performance of ocean color NPP models in Case-2 waters that are close to coastlines. Finally, ocean color chlorophyll-a algorithms are challenged by optically complex Case-2 waters, thus using satellite-derived chlorophyll-a to estimate NPP in coastal areas would likely further reduce the skill of ocean color models.

An Ensemble Kalman Filter multi-tracer assimilation: Determining uncertain ocean model parameters for improved climate-carbon cycle projections

An Ensemble Kalman Filter is applied to assimilate observed tracer fields in various combinations in the Bern3D ocean model. Each tracer combination yields a set of optimal transport parameter values that are used in projections with prescribed CO2 stabilization pathways. The assimilation of temperature and salinity fields yields a too vigorous ventilation of the thermocline and the deep ocean, whereas the inclusion of CFC-11 and radiocarbon improves the representation of physical and biogeochemical tracers and of ventilation time scales. Projected peak uptake rates and cumulative uptake of CO2 by the ocean are around 20% lower for the parameters determined with CFC-11 and radiocarbon as additional target compared to those with salinity and temperature only. Higher surface temperature changes are simulated in the Greenland–Norwegian–Iceland Sea and in the Southern Ocean when CFC-11 is included in the Ensemble Kalman model tuning. These findings highlights the importance of ocean transport calibration for the design of near-term and long-term CO2 emission mitigation strategies and for climate projections.

Model simulation output data (CCSM3) related to uplift of African topography and the intensification of upwelling in the Benguela region

The Benguela Current, located off the west coast of southern Africa, is tied to a highly productive upwelling system**1. Over the past 12 million years, the current has cooled, and upwelling has intensified**2, 3, 4. These changes have been variously linked to atmospheric and oceanic changes associated with the glaciation of Antarctica and global cooling**5, the closure of the Central American Seaway**1, 6 or the further restriction of the Indonesian Seaway**3. The upwelling intensification also occurred during a period of substantial uplift of the African continent**7, 8. Here we use a coupled ocean-atmosphere general circulation model to test the effect of African uplift on Benguela upwelling. In our simulations, uplift in the East African Rift system and in southern and southwestern Africa induces an intensification of coastal low-level winds, which leads to increased oceanic upwelling of cool subsurface waters. We compare the effect of African uplift with the simulated impact of the Central American Seaway closure9, Indonesian Throughflow restriction10 and Antarctic glaciation**11, and find that African uplift has at least an equally strong influence as each of the three other factors. We therefore conclude that African uplift was an important factor in driving the cooling and strengthening of the Benguela Current and coastal upwelling during the late Miocene and Pliocene epochs.

On the effects of circulation, sediment resuspension and biological incorporation by diatoms in an ocean model of aluminium, link to model data

The distribution of dissolved aluminium in the West Atlantic Ocean shows a mirror image with that of dissolved silicic acid, hinting at intricate interactions between the ocean cycling of Al and Si. The marine biogeochemistry of Al is of interest because of its potential impact on diatom opal remineralisation, hence Si availability. Furthermore, the dissolved Al concentration at the surface ocean has been used as a tracer for dust input, dust being the most important source of the bio-essential trace element iron to the ocean. Previously, the dissolved concentration of Al was simulated reasonably well with only a dust source, and scavenging by adsorption on settling biogenic debris as the only removal process. Here we explore the impacts of (i) a sediment source of Al in the Northern Hemisphere (especially north of ~ 40° N), (ii) the imposed velocity field, and (iii) biological incorporation of Al on the modelled Al distribution in the ocean. The sediment source clearly improves the model results, and using a different velocity field shows the importance of advection on the simulated Al distribution. Biological incorporation appears to be a potentially important removal process. However, conclusive independent data to constrain the Al / Si incorporation ratio by growing diatoms are missing. Therefore, this study does not provide a definitive answer to the question of the relative importance of Al removal by incorporation compared to removal by adsorptive scavenging.

Is FTA/EPA effective for a developing country to attract FDI? : simulation analysis based on an extended knowledge-capital model

To prepare an answer to the question of how a developing country can attract FDI, this paper explored the factors and policies that may help bring FDI into a developing country by utilizing an extended version of the knowledge-capital model. With a special focus on the effects of FTAs/EPAs between market countries and developing countries, simulations with the model revealed the following: (1) Although FTA/EPA generally ends to increase FDI to a developing country, the possibility of improving welfare through increased demand for skilled and unskilled labor becomes higher as the size of the country declines; (2) Because the additional implementation of cost-saving policies to reduce firm-type/trade-link specific fixed costs ends to depreciate the price of skilled labor by saving its input, a developing country, which is extremely scarce in skilled labor, is better off avoiding the additional option; (3) If a country hopes to enjoy larger welfare gains with EPA, efforts to increase skilled labor in the country, such as investing in education, may be beneficial.

Simulation analysis of the EU ELV/RoHS directives based on an applied general equilibrium model with Melitz-type trade specification

This paper explores the potential usefulness of an AGE model with the Melitz-type trade specification to assess economic effects of technical regulations, taking the case of the EU ELV/RoHS directives as an example. Simulation experiments reveal that: (1) raising the fixed exporting cost to make sales in the EU market brings results that exports of the targeted commodities (motor vehicles and parts for ELV and electronic equipment for RoHS) to the EU from outside regions/countries expand while the domestic trade in the EU shrinks when the importer's preference for variety (PfV) is not strong; (2) if the PfV is not strong, policy changes that may bring reduction in the number of firms enable survived producers with high productivity to expand production to be large-scale mass producers fully enjoying the fruit of economies of scale; and (3) When the strength of the importer's PfV is changed from zero to unity, there is the value that totally changes simulation results and their interpretations.

Construction of an Intraplate Fault System Model of South Australia, and Simulation Tool for the iSERVO Institute Seed Project

To foster ongoing international cooperation beyond ACES (APEC Cooperation for Earthquake Simulation) on the simulation of solid earth phenomena, agreement was reached to work towards establishment of a frontier international research institute for simulating the solid earth: iSERVO = International Solid Earth Research Virtual Observatory institute (http://www.iservo.edu.au). This paper outlines a key Australian contribution towards the iSERVO institute seed project, this is the construction of: (1) a typical intraplate fault system model using practical fault system data of South Australia (i.e., SA interacting fault model), which includes data management and editing, geometrical modeling and mesh generation; and (2) a finite-element based software tool, which is built on our long-term and ongoing effort to develop the R-minimum strategy based finite-element computational algorithm and software tool for modelling three-dimensional nonlinear frictional contact behavior between multiple deformable bodies with the arbitrarily-shaped contact element strategy. A numerical simulation of the SA fault system is carried out using this software tool to demonstrate its capability and our efforts towards seeding the iSERVO Institute.

Holocene experiment with coupled atmosphere-ocean-model ECHAM5/MPI-OM

Changes in the Earth's orbit lead to changes in the seasonal and meridional distribution of insolation. We quantify the influence of orbitally induced changes on the seasonal temperature cycle in a transient simulation of the last 6000 years - from the mid-Holocene to today - using a coupled atmosphere-ocean general circulation model (ECHAM5/MPI-OM) including a land surface model (JSBACH). The seasonal temperature cycle responds directly to the insolation changes almost everywhere. In the Northern Hemisphere, its amplitude decreases according to an increase in winter insolation and a decrease in summer insolation. In the Southern Hemisphere, the opposite is true. Over the Arctic Ocean, decreasing summer insolation leads to an increase in sea-ice cover. The insulating effect of sea ice between the ocean and the atmosphere leads to decreasing heat flux and favors more "continental" conditions over the Arctic Ocean in winter, resulting in strongly decreasing temperatures. Consequently, there are two competing effects: the direct response to insolation changes and a sea-ice insulation effect. The sea-ice insulation effect is stronger, and thus an increase in the amplitude of the seasonal temperature cycle over the Arctic Ocean occurs. This increase is strongest over the Barents Shelf and influences the temperature response over northern Europe. We compare our modeled seasonal temperatures over Europe to paleo reconstructions. We find better agreements in winter temperatures than in summer temperatures and better agreements in northern Europe than in southern Europe, since the model does not reproduce the southern European Holocene summer cooling inferred from the paleo reconstructions. The temperature reconstructions for northern Europe support the notion of the influence of the sea-ice insulation effect on the evolution of the seasonal temperature cycle.

Model simulation output data (CCSM3) related to uplift of Africa and Himalaya/Tibet Plateau and the influence on tropical African vegetation

Hominid evolution in the late Miocene has long been hypothesized to be linked to the retreat of the tropical rainforest in Africa. One cause for the climatic and vegetation change often considered was uplift of Africa, but also uplift of the Himalaya and the Tibetan Plateau was suggested to have impacted rainfall distribution over Africa. Recent proxy data suggest that in East Africa open grassland habitats were available to the common ancestors of hominins and apes long before their divergence and do not find evidence for a closed rainforest in the late Miocene. We used the coupled global general circulation model CCSM3 including an interactively coupled dynamic vegetation module to investigate the impact of topography on African hydro-climate and vegetation. We performed sensitivity experiments altering elevations of the Himalaya and the Tibetan Plateau as well as of East and Southern Africa. The simulations confirm the dominant impact of African topography for climate and vegetation development of the African tropics. Only a weak influence of prescribed Asian uplift on African climate could be detected. The model simulations show that rainforest coverage of Central Africa is strongly determined by the presence of elevated African topography. In East Africa, despite wetter conditions with lowered African topography, the conditions were not favorable enough to maintain a closed rainforest. A discussion of the results with respect to other model studies indicates a minor importance of vegetation-atmosphere or ocean-atmosphere feedbacks and a large dependence of the simulated vegetation response on the land surface/vegetation model.

Comparison between detailed model simulation and artificial neural network for forecasting building energy consumption

There are several ways to attempt to model a building and its heat gains from external sources as well as internal ones in order to evaluate a proper operation, audit retrofit actions, and forecast energy consumption. Different techniques, varying from simple regression to models that are based on physical principles, can be used for simulation. A frequent hypothesis for all these models is that the input variables should be based on realistic data when they are available, otherwise the evaluation of energy consumption might be highly under or over estimated. In this paper, a comparison is made between a simple model based on artificial neural network (ANN) and a model that is based on physical principles (EnergyPlus) as an auditing and predicting tool in order to forecast building energy consumption. The Administration Building of the University of Sao Paulo is used as a case study. The building energy consumption profiles are collected as well as the campus meteorological data. Results show that both models are suitable for energy consumption forecast. Additionally, a parametric analysis is carried out for the considered building on EnergyPlus in order to evaluate the influence of several parameters such as the building profile occupation and weather data on such forecasting. (C) 2008 Elsevier B.V. All rights reserved.

Satellite observation and climate system model simulation of the St. Lawrence Island polynya

The St. Lawrence Island polynya (SLIP) is a commonly occurring winter phenomenon in the Bering Sea, in which dense saline water produced during new ice formation is thought to flow northward through the Bering Strait to help maintain the Arctic Ocean halocline. Winter darkness and inclement weather conditions have made continuous in situ and remote observation of this polynya difficult. However, imagery acquired from the European Space Agency ERS-1 Synthetic Aperture Radar (SAR) has allowed observation of the St. Lawrence Island polynya using both the imagery and derived ice displacement products. With the development of ARCSyM, a high resolution regional model of the Arctic atmosphere/sea ice system, simulation of the SLIP in a climate model is now possible. Intercomparisons between remotely sensed products and simulations can lead to additional insight into the SLIP formation process. Low resolution SAR, SSM/I and AVHRR infrared imagery for the St. Lawrence Island region are compared with the results of a model simulation for the period of 24-27 February 1992. The imagery illustrates a polynya event (polynya opening). With the northerly winds strong and consistent over several days, the coupled model captures the SLIP event with moderate accuracy. However, the introduction of a stability dependent atmosphere-ice drag coefficient, which allows feedbacks between atmospheric stability, open water, and air-ice drag, produces a more accurate simulation of the SLIP in comparison to satellite imagery. Model experiments show that the polynya event is forced primarily by changes in atmospheric circulation followed by persistent favorable conditions: ocean surface currents are found to have a small but positive impact on the simulation which is enhanced when wind forcing is weak or variable.

An alternative plate tectonic model for the Palaeozoic-Early Mesozoic Palaeotethyan evolution of Southeast Asia (northern Thailand-Myanmar)

An alternative model for the geodynamic evolution of Southeast Asia is proposed and inserted in a modern plate tectonic model. The reconstruction methodology is based on dynamic plate boundaries, constrained by data such as spreading rates and subduction velocities; in this way it differs from classical continental drift models proposed so far. The different interpretations about the location of the Palaeotethys suture in Thailand are revised, the Tertiary Mae Yuam fault is seen as the emplacement of the suture. East of the suture we identify an Indochina derived terrane for which we keep the name Shan-Thai, formerly used to identify the Cimmerian block present in Southeast Asia, now called Sibumasu. This nomenclatural choice was made on the basis of the geographic location of the terrane (Eastern Shan States in Burma and Central Thailand) and in order not to introduce new confusing terminology. The closure of the Eastern Palaeotethys is related to a southward subduction of the ocean, that triggered the Eastern Neotethys to open as a back-arc, due to the presence of Late Carboniferous-Early Permian arc magmatism in Mergui (Burma) and in the Lhasa block (South Tibet), and to the absence of arc magmatism of the same age East of the suture. In order to explain the presence of Carboniferous-Early Permian and Permo-Triassic volcanic arcs in Cambodia, Upper Triassic magmatism in Eastern Vietnam and Lower Permian-Middle Permian arc volcanites in Western Sumatra, we introduce the Orang Laut terranes concept. These terranes were detached from Indochina and South China during back-arc opening of the Poko-Song Ma system, due to the westward subduction of the Palaeopacific. This also explains the location of the Cathaysian West Sumatra block to the West of the Cimmerian Sibumasu block.