Quantifying rates of landscape evolution and tectonic processes by thermochronology and numerical modeling of crustal heat transport using PECUBE

PECUBE is a three-dimensional thermal-kinematic code capable of solving the heat production-diffusion-advection equation under a temporally varying surface boundary condition. It was initially developed to assess the effects of time-varying surface topography (relief) on low-temperature thermochronological datasets. Thermochronometric ages are predicted by tracking the time-temperature histories of rock-particles ending up at the surface and by combining these with various age-prediction models. In the decade since its inception, the PECUBE code has been under continuous development as its use became wider and addressed different tectonic-geomorphic problems. This paper describes several major recent improvements in the code, including its integration with an inverse-modeling package based on the Neighborhood Algorithm, the incorporation of fault-controlled kinematics, several different ways to address topographic and drainage change through time, the ability to predict subsurface (tunnel or borehole) data, prediction of detrital thermochronology data and a method to compare these with observations, and the coupling with landscape-evolution (or surface-process) models. Each new development is described together with one or several applications, so that the reader and potential user can clearly assess and make use of the capabilities of PECUBE. We end with describing some developments that are currently underway or should take place in the foreseeable future. (C) 2012 Elsevier B.V. All rights reserved.

Western Préalpes Médianes Romandes: Timing and structure. A review.

Between the original position and their present day location as klippen, the Prealpes Medianes underwent a complex history of paleotectonics and alpine tectonics. Due to the opening of the Piemont ocean the Brianconnais sedimentation realm of the Prealpes Medianes evolved as a rim basin of the northern passive margin during Jurassic to Eocene times. Different paleotectonic features (normal faults, synsedimentary growth structures, inversion structures) developed and were active above a basal detachment in evaporitic layers. The tectonic movements were a consequence of thermal events in the crust. Isolated from the Iberic continent at the end of the Late Cretaceous, the Brianconnais exotic terrain was incorporated into the accretionary prism of the closing Piemont ocean and the incipient alpine orogeny during the Lutetian-Bartonian. The Prealpes Medianes were detached from their homeland during the Bartonian-Priabonian and were transported onto the foreland. The tectonic style is one of a thin-skinned foreland fold and thrust belt. Fault associated fold development above a main decollement, together with internal deformation, represent the Prealpes Medianes main structural features. The very low-grade metamorphic conditions have their origin in the heat flux induced by tectonic burial by overriding nappes in the accretionary prism. After having been transported on top of the developing Helvetic nappes the Prealpes were emplaced in their present day position in front of the Alpine mountain belt during Oligocene times. Post-emplacement and out of sequence thrusting, possibly younger than Oligocene, is observed and can be related to thrusting in the sedimentary substratum and the basement.

Stability analysis of the 2007 Chehalis lake landslide based on long-range terrestrial photogrammetry and airborne LIDAR data

On December 4th 2007, a 3-Mm3 landslide occurred along the northwestern shore of Chehalis Lake. The initiation zone is located at the intersection of the main valley slope and the northern sidewall of a prominent gully. The slope failure caused a displacement wave that ran up to 38 m on the opposite shore of the lake. The landslide is temporally associated with a rain-on-snow meteorological event which is thought to have triggered it. This paper describes the Chehalis Lake landslide and presents a comparison of discontinuity orientation datasets obtained using three techniques: field measurements, terrestrial photogrammetric 3D models and an airborne LiDAR digital elevation model to describe the orientation and characteristics of the five discontinuity sets present. The discontinuity orientation data are used to perform kinematic, surface wedge limit equilibrium and three-dimensional distinct element analyses. The kinematic and surface wedge analyses suggest that the location of the slope failure (intersection of the valley slope and a gully wall) has facilitated the development of the unstable rock mass which initiated as a planar sliding failure. Results from the three-dimensional distinct element analyses suggest that the presence, orientation and high persistence of a discontinuity set dipping obliquely to the slope were critical to the development of the landslide and led to a failure mechanism dominated by planar sliding. The three-dimensional distinct element modelling also suggests that the presence of a steeply dipping discontinuity set striking perpendicular to the slope and associated with a fault exerted a significant control on the volume and extent of the failed rock mass but not on the overall stability of the slope.

40Ar/39Ar dating of Penninic Front tectonic displacement (W Alps) during the Lower Oligocene (31-34 Ma)

Direct absolute dating of the Penninic Frontal Thrust tectonic motion is achieved using the Ar-40/Ar-39 technique in the Pelvoux Crystalline Massif (Western Alps). The dated phengites were formed syn-kinematically in shear zones. They underline the brittle-ductile stretching lineation, pressure-shadow fibres and slickensides consistent with underthrusting of the European continental slab below the propagating Penninic Thrust. Chlorite-phengite thermobarometry yields 10-15 km and T similar to 280 degrees C, while Ar-40/Ar-39 phengite ages mainly range between 34 and 30 Ma, with one younger age at 27 Ma. This Early Oligocene age range matches a major tectonic rearrangement of the Alpine chain. Preservation of prograde Ar-40/Ar-39 ages is ascribed to passive exhumation of the Pelvoux shear zone network, sandwiched between more external thrusts and the Penninic Front reactivated as an E-dipping detachment fault. Partial resetting in the Low Temperature part of argon spectra below 24 Ma is ascribed to brittle deformation and alteration of phengites.

Imaging molasse and quaternary sediments in Lake Geneva, Switzerland, with 3-D high-resolution seismic reflection methods: A case study

A high-resolution three-dimensional (3-D) seismic reflection survey was conducted in Lake Geneva, near the city of Lausanne, Switzerland, as part of a project for developing such seismic techniques. Using a single 48-channel streamer, the 3-D site with an area of 1200 m x 600 m was surveyed in 10 days. A variety of complex geologic structures (e.g. thrusts, folds, channel-fill) up to similar to150 m below the water bottom were obtained with a 15 in.(3) water gun. The 3-D data allowed the construction of an accurate velocity model and the distinction of five major seismic facies within the Lower Freshwater Molasse (Aquitanian) and the Quaternary sedimentary units. Additionally, the Plateau Molasse (PM) and Subalpine Molasse (SM) erosional surface, "La Paudeze" thrust fault (PM-SM boundary) and the thickness of Quaternary sediments were accurately delineated in 3-D.

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.

Accreted fragments of the Late Cretaceous Caribbean-Colombian Plateau in Ecuador

The eastern part of the Western Cordillera of Ecuador includes fragments of an Early Cretaceous ( approximate to 123 Ma) oceanic plateau accreted around 85-80 Ma (San Juan unit). West of this unit and in fault contact with it, another oceanic plateau sequence (Guaranda unit) is marked by the occurrence of picrites, ankaramites, basalts, dolerites and shallow level gabbros. A comparable unit is also exposed in northwestern coastal Ecuador (Pedernales unit). Picrites have LREE-depleted patterns, high epsilonNd(i) and very low Pb isotopic ratios, suggesting that they were derived from an extremely depleted source. In contrast, the ankaramites and Mg-rich basalts are LREE-enriched and have radiogenic Pb isotopic compositions similar to the Galapagos HIMU component; their epsilonNd(i) are slightly lower than those of the picrites. Basalts, dolerites and gabbros differ from the picrites and ankaramites by flat rare earth element (REE) patterns and lower epsilonNd; their Pb isotopic compositions are intermediate between those of the picrites and ankaramites. The ankaramites, Mg-rich basalts, and picrites differ from the lavas from the San Juan-Multitud Unit by higher Pb ratios and lower epsilonNd(i). The Ecuadorian and Gorgona 88-86 Ma picrites are geochemically similar. The Ecuadorian ankaramites and Mg-rich basalts share with the 92-86 Ma Mg-rich basalts of the Caribbean-Colombian Oceanic Plateau (CCOP) similar trace element and Nd and Pb isotopic chemistry. This suggests that the Pedernales and Guaranda units belong to the Late Cretaceous CCOR The geochemical diversity of the Guaranda and Pedernales rocks illustrates the heterogeneity of the CCOP plume source and suggests a multi-stage model for the emplacement of these rocks. Stratigraphic and geological relations strongly suggest that the Guaranda unit was accreted in the late Maastrichtian (approximate to 68-65 Ma). (C) 2002 Elsevier Science B.V. All rights reserved.

Stratigraphic significance of siliceous microfossils collected during Nautiperc dives (Off Peru, 5°-6° S)

The geological evolution of the northern Peru convergent margin can be traced using samples collected during deep-sea dives of the submersible Nautile. In the Paita area (5 degrees-6 degrees S), the sedimentary sequence was intensively sampled along the main scarp of the middle slope area. It consists of Upper Miocene (7-9 Ma) to Pleistocene siltstone, sandstone and rare dolostone. The age distribution of these samples is the basis for a new geologic interpretation of the multichannel seismic line CDP3. Siliceous microfossils (both diatoms and radiolarians) show influence of both cold and temperate waters (local species mixed with upwelling ones). Diatom assemblages studied from the NP1-13 and NP1-15 dives bear a strong resemblance to assemblages from the Pisco Formation of southern Peru. Micropaleontological data from siliceous microfossils, provide evidence for two main unconformities, one is at the base of the Quaternary sequence and the other corresponds to a hiatus of 1 Myr, separating the Upper Miocene (7-8 Ma) sediments from uppermost Miocene (5-6 Ma) sediments. During the past 400 kyr, a wide rollover fold developed in the middle slope area associated with a major seaward dipping detachment fault. A catastrophic debris avalanche occurred as the result of an oversteepening of the landward flank of the rollover fold. The gravity failure of the slope, recognized by SeaBEAM and hydrosweep mapping, displaced enough material to produce a destructive tsunami which occurred 13.8 +/- 2.7 kyr ago.

Exhumation history of eastern Ladakh revealed by 40Ar/39Ar and fission-track ages: the Indus River-Tso Morari transect, NW Himalaya

Fission-track and (40)Ar/(39)Ar ages place time constraints on the exhumation of the North Himalayan nappe stack, the Indus Suture Zone and Molasse, and the Transhimalayan Batholith in eastern Ladakh (NW India). Results from this and previous studies on a north-south transect passing near Tso Morari Lake suggest that the SW-directed North Himalayan nappe stack (comprising the Mata, Tetraogal and Tso Morari nappes) was emplaced and metamorphosed by c. 50-45 Ma, and exhumed to moderately shallow depths (c. 10 km) by c. 45-40 Ma. From the mid-Eocene to the present, exhumation continued at a steady and slow rate except for the root zone of the Tso Morari nappe, which cooled faster than the rest of the nappe stack. Rapid cooling occurred at c. 20 Ma and is linked to brittle deformation along the normal Ribil-Zildat Fault concomitant with extrusion of the Crystalline nappe in the south. Data from the Indus Molasse suggest that sediments were still being deposited during the Miocene.

The hypogene iron oxide copper-gold mineralization in the Mantoverde District, Northern Chile

The Mantoverde iron oxide copper-gold (IOCC) district, northern Chile, is known for its Cu production from supergene ores. Recently, exploration outlined an additional hypogene ore resource of 440 Mt with 0.56 percent Cu, and 0.12 g/t An. The hypogene sulfide mineralization occurs mainly as chalcopyrite and pyrite, typically in specularite or magnetite-cemented breccias and associated stockworks. The host rocks underwent variably intense K feldspar alteration, chloritization, sericitization, silicification, and/or carbonatization. A district scale Na(-Ca) alteration is absent. The IOCC mineralization in the district shows a strong tectonic control by northwest- to north-northwest-trending brittle structures. Large Cu sulfide-rich veins or Cu sulfide-cemented breccias are absent. Therefore, head grades of 4 percent Cu are an exception. There is a positive correlation between Cu and An grades. Gold is probably contained mostly in chalcopyrite and pyrite. Elevated concentrations of light rare-earth elements (LREE) occur locally but are attributed to redistribution of LREE within the deposits rather than to derivation from external sources. The Cu-Au ores in the Mantoverde district are low in and have relatively low contents in heavy metals that are potentially hazardous to the environment, such as As (avg 14 ppm), Hg (<5 ppm), or Cd (<0.2 ppm). The sulfur isotope ratios of chalcopyrite from the IOCC deposits lie between -5.6 and 8.9 per mil delta(34)S(VCDT). They show systematic variations within the district, which are interpreted to reflect relative distance to inferred fluid conduits and the level of deposition within the hydrothermal system. Most initial (87)Sr/(86)Sr values of altered volcanic rocks and hydrothermal calcite from the Mantoverde district are between 0.7031 and 0.7060 and are similar to those of the igneous rocks of the region. Lead isotope ratios of chalcopyrite are consistent with Pb (and by inference Cu) derived from Early Cretaceous magmatism. The sulfur, strontium, and lead isotope data of chalcopyrite, calcite gangue, or altered host rocks, respectively, are compatible with a genetic model that involves cooling of metal and sulfur-bearing magmatic-hydrothermal fluids that mix with meteoric waters or seawater at relatively shallow crustal levels. An additional exotic sulfur input is likely, though not required, for the copper mineralization. Apart from the IOCC. deposits, there are a number of smaller magnetite(-apatite) bodies in the district. These are geologically similar to the Cu-Au-bearing magnetite bodies, but are related to splays of the north-south-trending Atacama fault zone and differ in alteration and texture.

The Tertiary structural and thermal evolution of the central Alps - Compressional and extensional structures in an orogenic belt

The Western Alpine Are has been created during the Cretaceous and the Tertiary orogenies. The interference patterns of the Tertiary structures suggest their formation during continental collision of the European and the Adriatic Plates, with an accompanying anticlockwise rotation of the Adriatic indenter. Extensional structures are mainly related to ductile deformation by simple shear. These structures developed at a deep tectonic level, in granitic crustal rocks, at depths in excess of 10 km. In the early Palaeogene period of the Tertiary Orogeny, the main Tertiary nappe emplacement resulted from a NW-thrusting of the Austroalpine, Penninic and Helvetic nappes. Heating of the deep zone of the Upper Cretaceous and Tertiary nappe stack by geothermal heat flow is responsible for the Tertiary regional metamorphism, reaching amphibolite-facies conditions in the Lepontine Gneiss Dome (geothermal gradient 25 degrees C/ km). The Tertiary thrusting occurred mainly during prograde metamorphic conditions with creation of a penetrative NW-SE-oriented stretching lineation, X(1) (finite extension), parallel to the direction of simple shear. Earliest cooling after the culmination of the Tertiary metamorphism, some 38 Ma ago, is recorded by the cooling curves of the Monte Rosa and Mischabel nappes to the west and the Suretta Nappe to the east of the Lepontine Gneiss Dome. The onset of dextral transpression, with a strong extension parallel to the mountain belt, and the oldest S-vergent `'backfolding'' took place some 35 to 30 Ma ago during retrograde amphibolite-facies conditions and before the intrusion of the Oligocene dikes north of the Periadriatic Line. The main updoming of the Lepontine Gneiss Dome started some 32-30 Ma ago with the intrusion of the Bergell tonalites and granodiorites, concomitant with S-vergent backfolding and backthrusting and dextral strike-slip movements along the Tonale and Canavese Lines (Argand's Insubric phase). Subsequently, the center of main updoming migrated slowly to the west, reaching the Simplon region some 20 Ma ago. This was contemporaneous with the westward migration of the Adriatic indenter. Between 20 Ma and the present, the Western Aar Massif-Toce culmination was the center of strong uplift. The youngest S-vergent backfolds, the Glishorn anticline and the Berisal syncline fold the 12 Ma Rb/Sr biotite isochron and are cut by the 11 Ma old Rhone-Simplon Line. The discrete Rhone-Simplon Line represents a late retrograde manifestation in the preexisting ductile Simplon Shear Zone. This fault zone is still active today. The Oligocene-Neogene dextral transpression and extension in the Simplon area were concurrent with thrusting to the northwest of the Helvetic nappes, the Prealpes (35-15 Ma) and with the Jura thin-skinned thrust (11-3 Ma). It was also contemporaneous with thrusting to the south of the Bergamasc (> 35-5 Ma) and Milan thrusts (16-5 Ma).

Jurassic sedimentary record and tectonic evolution of the northwestern corner of Africa

The history of the opening seaway from the westernmost Tethys to the Central Atlantic is traced by the analysis of the sedimentary facies development in the external Rif basin of Northern Morocco and the geological and seismic data from the Moroccan Atlantic continental margin. In the Rif basin, after the early Sinemurian, sedimentary facies dated by ammonites, foraminifers and brachiopods, indicate a progression of rapid subsidence resulting from extensional tectonic (tilted blocks, escarpment fault breccias, neptunian dykes etc.) from the N and NE to the S and SW. From the Toarcian to the Bajocian, deltas progress from the W and SW into the `'Rides sud-rifaines'' realm. From the late Bathonian to the Oxfordian, deep-sea fans develop in the external Rif. During the same period, deltaic sediments fill in the Middle Atlas basin of Eastern Morocco and progress into the external Rif. The top of the Jurassic is characterised by carbonate deposits. At the northwestern corner of Africa, the subsidence of the sedimentary basins by rifting is initiated in the late Triassic; however, at the Mazagan transect of the Atlantic continental margin, the tectonic pattern characteristic of a passive continental margin appears clearly only in the early Jurassic. At the foot of the Mazagan escarpment, the sedimentary record shows a foundering of the first bloc during early to middle Lias. A thermal uplift phase is indicated by emersion of the African margin shoulder in late Liassic, and thermal relaxation starts in the middle Jurassic. The morphology of this transect, compared with the conjugate side of the American continent is most easily explained by the uniform sense simple shear model.

Synorogenic extension in the Tethyan Himalaya documented by structural studies and the Kubler index, Lachung La area, NW India

Tectonic observations in the Tethyan Himalaya reveal an important extensional event that succeeds the emplacement of SW-verging nappes. A major thrust, called the Kum Tso Thrust, has been backfolded and reactivated by normal faulting associated with this event. Measurements of the Kubler index, coupled with characterization of clay-size paragenesis show the effect of normal faulting on the regional metamorphic zonation and indicate that important extension zones, like the Sarchu-Lachung La Normal Fault Zone (SLFZ), exist within the Tethyan Himalaya. Diagenetic limestones from within the SLFZ are characterized by the occurrence of mixed-layered clay phases, kaolinite and an illite with a 001 peak >0.4 Delta degrees2 theta. This zone is bordered by two anchizonal-to-epizonal zones, where illite peaks become narrower. Further to the NE the successive appearance of biotite, chloritoid, garnet and garnet-staurolite-kyanite assemblapes testifies to an increase in metamorphic grade. The cataclastic samples from the normal faults contain kaolinite, smectite and a `broad' illite, indicating that extension occurs under diagenetic conditions.

High-resolution seismic images of potentially seismogenic structures beneath the northwest Canterbury Plains, New Zealand

The transpressional boundary between the Australian and Pacific plates in the central South Island of New Zealand comprises the Alpine Fault and a broad region of distributed strain concentrated in the Southern Alps but encompassing regions further to the east, including the northwest Canterbury Plains. Low to moderate levels of seismicity (e. g., 2 > M 5 events since 1974 and 2 > M 4.0 in 2009) and Holocene sediments offset or disrupted along rare exposed active fault segments are evidence for ongoing tectonism in the northwest plains, the surface topography of which is remarkably flat and even. Because the geology underlying the late Quaternary alluvial fan deposits that carpet most of the plains is not established, the detailed tectonic evolution of this region and the potential for larger earthquakes is only poorly understood. To address these issues, we have processed and interpreted high-resolution (2.5 m subsurface sampling interval) seismic data acquired along lines strategically located relative to extensive rock exposures to the north, west, and southwest and rare exposures to the east. Geological information provided by these rock exposures offer important constraints on the interpretation of the seismic data. The processed seismic reflection sections image a variably thick layer of generally undisturbed younger (i.e., < 24 ka) Quaternary alluvial sediments unconformably overlying an older (> 59 ka) Quaternary sedimentary sequence that shows evidence of moderate faulting and folding during and subsequent to deposition. These Quaternary units are in unconformable contact with Late Cretaceous-Tertiary interbedded sedimentary and volcanic rocks that are highly faulted, folded, and tilted. The lowest imaged unit is largely reflection-free Permian Triassic basement rocks. Quaternary-age deformation has affected all the rocks underlying the younger alluvial sediments, and there is evidence for ongoing deformation. Eight primary and numerous secondary faults as well as a major anticlinal fold are revealed on the seismic sections. Folded sedimentary and volcanic units are observed in the hanging walls and footwalls of most faults. Five of the primary faults represent plausible extensions of mapped faults, three of which are active. The major anticlinal fold is the probable continuation of known active structure. A magnitude 7.1 earthquake occurred on 4 September 2010 near the southeastern edge of our study area. This predominantly right-lateral strike-slip event and numerous aftershocks (ten with magnitudes >= 5 within one week of the main event) highlight the primary message of our paper: that the generally flat and topographically featureless Canterbury Plains is underlain by a network of active faults that have the potential to generate significant earthquakes.

Towards robust network based complex systems : from evolutionary cellular automata to biological models

Abstract Sitting between your past and your future doesn't mean you are in the present. Dakota Skye Complex systems science is an interdisciplinary field grouping under the same umbrella dynamical phenomena from social, natural or mathematical sciences. The emergence of a higher order organization or behavior, transcending that expected of the linear addition of the parts, is a key factor shared by all these systems. Most complex systems can be modeled as networks that represent the interactions amongst the system's components. In addition to the actual nature of the part's interactions, the intrinsic topological structure of underlying network is believed to play a crucial role in the remarkable emergent behaviors exhibited by the systems. Moreover, the topology is also a key a factor to explain the extraordinary flexibility and resilience to perturbations when applied to transmission and diffusion phenomena. In this work, we study the effect of different network structures on the performance and on the fault tolerance of systems in two different contexts. In the first part, we study cellular automata, which are a simple paradigm for distributed computation. Cellular automata are made of basic Boolean computational units, the cells; relying on simple rules and information from- the surrounding cells to perform a global task. The limited visibility of the cells can be modeled as a network, where interactions amongst cells are governed by an underlying structure, usually a regular one. In order to increase the performance of cellular automata, we chose to change its topology. We applied computational principles inspired by Darwinian evolution, called evolutionary algorithms, to alter the system's topological structure starting from either a regular or a random one. The outcome is remarkable, as the resulting topologies find themselves sharing properties of both regular and random network, and display similitudes Watts-Strogtz's small-world network found in social systems. Moreover, the performance and tolerance to probabilistic faults of our small-world like cellular automata surpasses that of regular ones. In the second part, we use the context of biological genetic regulatory networks and, in particular, Kauffman's random Boolean networks model. In some ways, this model is close to cellular automata, although is not expected to perform any task. Instead, it simulates the time-evolution of genetic regulation within living organisms under strict conditions. The original model, though very attractive by it's simplicity, suffered from important shortcomings unveiled by the recent advances in genetics and biology. We propose to use these new discoveries to improve the original model. Firstly, we have used artificial topologies believed to be closer to that of gene regulatory networks. We have also studied actual biological organisms, and used parts of their genetic regulatory networks in our models. Secondly, we have addressed the improbable full synchronicity of the event taking place on. Boolean networks and proposed a more biologically plausible cascading scheme. Finally, we tackled the actual Boolean functions of the model, i.e. the specifics of how genes activate according to the activity of upstream genes, and presented a new update function that takes into account the actual promoting and repressing effects of one gene on another. Our improved models demonstrate the expected, biologically sound, behavior of previous GRN model, yet with superior resistance to perturbations. We believe they are one step closer to the biological reality.