hierarchical segmentation of the malawi rift the influence of inherited lithospheric heterogeneity and kinematics in the evolution of continental rifts we used detailed analysis of shuttle radar topography mission digital elevation model and observations from aeromagnetic data to examine the influence of inherited lithospheric heterogeneity and kinematics in the segmentation of largely amagmatic continental rifts we focused on the cenozoic malawi rift which represents the southern extension of the western branch of the east african rift system this north trending rift traverses precambrian and paleozoic mesozoic structures of different orientations we found that the rift can be hierarchically divided into first order and second order segments in the first order segmentation we divided the rift into northern central and southern sections in its northern section the rift follows paleoproterozoic and neoproterozoic terrains with structural grain that favored the localization of extension within well developed border faults the central section occurs within mesoproterozoic neoproterozoic terrain with regional structures oblique to the rift extent we propose that the lack of inherited lithospheric heterogeneity favoring extension localization resulted in the development of the rift in this section as a shallow graben with undeveloped border faults in the southern section mesoproterozoic neoproterozoic rocks were reactivated and developed the border faults in the second order segmentation only observed in the northern section we divided the section into five segments that approximate four half grabens asymmetrical grabens with alternating polarities the change of polarity coincides with flip over full grabens occurring within overlap zones associated with 150 km long alternating border faults segments the inherited lithospheric heterogeneity played the major role in facilitating the segmentation of the malawi rift during its opening resulting from extension

control of pre rift lithospheric structure on the architecture and evolution of continental rifts insights from the main ethiopian rift east africa we investigate the along axis variations in architecture segmentation and evolution of the main ethiopian rift mer east africa and relate these characteristics to the regional geology lithospheric structure and surface processes we first illustrate significant along axis variations in basin architecture through analysis of simplified geological cross sections in different rift sectors we then integrate this information with a new analysis of ethiopian topography and hydrography to illustrate how rift architecture basin symmetry asymmetry is reflected in the margin topography and has been likely amplified by a positive feedback between tectonics flexural uplift and surface processes fluvial erosion unloading this analysis shows that 70 of the 500 km long mer is asymmetric with most of the asymmetric rift sectors being characterized by a master fault system on the eastern margin we finally relate rift architecture and segmentation to the regional geology and geophysical constraints on the lithosphere we provide strong evidence that rift architecture is controlled by the contrasting nature of the lithosphere beneath the homogeneous strong somalian plateau and the weaker more heterogeneous ethiopian plateau differences originating from the presence of pre rift zones of weakness on the ethiopian plateau and likely amplified by surface processes the data provided by this integrated analysis suggest that asymmetric rifts may directly progress to focused axial tectonic magmatic activity without transitioning into a symmetric rifting stage these observations have important implications for the asymmetry of continental rifts and conjugate passive margins worldwide

ice pumps and their rates an ice pump is a heat engine driven by the change of freezing point with pressure which will melt ice at depth in the ocean and deposit it at a shallower location it is self starting calculations of the maximum magnitude of this effect are made which show good agreement with field data available for sea and lake ice the discussion is applied to the general case of a moving pack ice sheet with a well mixed surface layer and to floating ice shelves the rate of melt from an 11 m deep pressure ridge keel due to ice pumping is estimated as 26 cm year and that from the front of the ross ice shelf at mcmurdo sound antarctica is estimated as 5 m year for the level of water movement noted in the authors field observations far from the ice front pumping between shelf areas of different thickness will still occur with tidal motion providing the necessary water exchange but its magnitude is now limited by the ability to remove the potentially stable layer of melt water out of the system it is important to realize that the pumping does not depend on the availability of sensible heat in the water column and its effects are additional to any melting caused by the advection of warmer water to the ice water interface

antarctic ice sheet melting in the southeast pacific the first oceanographic measurements across a deep channel beneath the calving front of pine island glacier reveal a sub ice circulation driven by basal melting of 10 12 m yr 1 a salt box model described here gives a melt rate similar to that of ice balance and numerical models 5 50 times higher than averages for the george vi and ross ice shelves melting is fueled by relatively warm circumpolar deep water that floods the deep floor of the amundsen and bellingshausen sea continental shelves reaching the deep draft of this floating glacier a revised melt rate for ice shelves in the southeast pacific sector raises circumpolar ice shelf melting to 756 gt yr 1 given prior estimates of surface accumulation and iceberg calving this suggests that the antarctic ice sheet is currently losing mass to the ocean

on the dynamics of β plumes abstract the dynamics of zonal recirculating flows plumes driven by a source of potential vorticity or a sink of mass located at the eastern boundary of an ocean basin are investigated using analytical solutions of the barotropic linear steady quasigeostrophic equation with bottom and lateral friction by scaling the ratio of the strength of the recirculation to that of the sink by that from the inviscid solution the regimes are identified in which friction becomes a dominant factor the primary new finding of this study is that the recirculating flow component disappears due to frictional effects when the meridional extent of the sink becomes small unlike for the zonal extent of the sink which affects the recirculating component only if it is on the order of the frictional boundary layer scale the deviation of the recirculating flow strength from that given by the inviscid solution is apparent for meridional sink scales which are much larger than the frictional boundary layer scale also lo

on the role of coastal troughs in the circulation of warm circumpolar deep water on antarctic shelves abstractoceanic exchanges across the continental shelves of antarctica play an important role in biological systems and the mass balance of ice sheets the focus of this study is on the mechanisms responsible for the circulation of warm circumpolar deep water cdw within troughs running perpendicular to the continental shelf this is examined using process oriented numerical experiments with an eddy resolving 1 km 3d ocean model that includes a static and thermodynamically active ice shelf three mechanisms that create a significant onshore flow within the trough are identified 1 a deep onshore flow driven by the melt of the ice shelf 2 interaction between the longshore mean flow and the trough and 3 interaction between a rossby wave along the shelf break and the trough in each case the onshore flow is sufficient to maintain the warm temperatures underneath the ice shelf and basal melt rates of o 1 m yr 1 the third mechanism in particular reproduces several features revealed by moorings from m

controls on early rift geometry new perspectives from the bilila mtakataka fault malawi we use the 110 km long bilila mtakataka fault in the amagmatic southern east african rift malawi to investigate the controls on early rift geometry at the scale of a major border fault morphological variations along the 14 8 m high scarp define six 10 40 km long segments which are either foliation parallel or oblique to both foliation and the current regional extension direction as the scarp is neither consistently parallel to foliation nor well oriented for the current regional extension direction we suggest the segmented surface expression is related to the local reactivation of well oriented weak shallow fabrics above a broadly continuous structure at depth using a geometrical model the geometry of the best fitting subsurface structure is consistent with the local strain field from recent seismicity in conclusion within this early rift pre existing weaknesses only locally control border fault geometry at subsurface

ice shelf grounding zone structure from icesat laser altimetry 1 we present a technique for investigating the grounding zone gz of antarctic ice shelves using laser altimetry from the ice cloud and land elevation satellite icesat most surface height variability in the gz is easily resolved by the icesat laser s 65 m footprint and 172 m along track spacing comparisons of repeated tracks sampled at different phases of the ocean tide identify the landward and seaward limits of tide forced ice flexure providing gz location and width information for each track using icesat data in the institute ice stream region of southern ronne ice shelf we demonstrate that the location of the gz based on feature identification in satellite imagery or digital elevation models may be in error by several km our results show that icesat will contribute significantly to improving knowledge of gz structure and to studies requiring accurate gz locations e g ice mass balance calculations and ice sheet ocean modeling

terrestrial heat flow in the malawi rifted zone east africa implications for tectono thermal inheritance in continental rift basins abstract rift basins show the widest range of heat flow values compared to passive continental margins and intracratonic sag basins elevated heat flow is often used to infer asthenospheric upwelling or subsurface magmatic activity for magma assisted rifts despite the lack of surface volcanism for magma poor rifts it is possible that magmatic bodies will be present at lower and mid crustal levels but yet to breach the surface thus necessitating consideration of heat flow from tectono thermal inheritance in addition to more recent magmatism we use aeromagnetic data to investigate the terrestrial heat flow distribution in the malawi rifted zone mrz in east africa defined as the neogene malawi rift and surrounding permo triassic karoo rift basins we use the two dimensional power density spectrum technique to estimate the curie point depth cpd geothermal gradient and heat flow beneath the mrz our results reveal predominance of shallow cpds 18 20 km high geothermal gradients 29 32 c km and elevated heat flow 70 82 mw m 2 within the karoo rift basins along the malawi rift geothermal gradient 25 27 c km and heat flow 60 66 mw m 2 are generally lower except pronounced high heat flow 70 82 mw m 2 at the rungwe volcanic province rvp in the north and at the central rift segment where the rift overprints karoo rift basins the surrounding cratonic blocks show deeper cpds 24 27 km lower geothermal gradients 22 24 c km and lower heat flow 53 63 mw m 2 we interpret that the elevated heat flow in the karoo rift basins are related to residual heat flow from the permo triassic rifting phase replenished by later jurassic cretaceous diking events apart from the thermal anomaly beneath the rvp which is due to neogene magmatism other important targets for geothermal explorations in the mrz are the residual thermal anomalies in the karoo rift basins we infer that in areas of active magma poor rifting pronounced heat flow may not only occur in regions of asthenospheric upwelling but may concentrate in the ancient magmatic rift segments our findings have implications for crustal stretching tectono thermal inheritance and geothermal energy potentials in east africa and other segments of magma poor rifting

a 1 3 million year record of synchronous faulting in the hangingwall and border fault of a half graben in the malawi nyasa rift abstract this paper analyzes throw depth t z profiles from a high resolution 2d reflection seismic grid in the central basin of lake malawi to investigate whether evidence exists 1 for migration of faulting away from the border fault of the half graben and 2 that faults in the hangingwall lengthened over the last 1 3 million years we use the high precision age model from a 2005 scientific drilling project in our study area to constrain the ages of our seismic horizons and examine a fault array and two individual faults within the hangingwall of the central basin border fault we account for climatic and sedimentological controls on stratal growth with a lake level curve that accompanies the age model a comparison of our hangingwall t z profiles with published throw distance t x profiles for the border fault shows synchronous faulting over the last 1 3 m y rather than basinward migration of faulting furthermore we find no evidence for significant propagation of the tips of the hangingwall faults in the last 1 3 m y and conclude that the lack of basinward migration of faulting is a consequence of strain localization on faults established at an early stage in basin development

simulation of subice shelf melt rates in a general circulation model velocity dependent transfer and the role of friction two parameterizations of turbulent boundary layer processes at the interface between an ice shelf and the ocean beneath are investigated in terms of their impact on simulated melt rates and feedbacks the parameterizations differ in the transfer coefficients for heat and freshwater fluxes in their simplest form they are assumed constant and hence are independent of the velocity of ocean currents at the ice shelf base an augmented melt rate parameterization accounts for frictional turbulence via transfer coefficients that do depend on boundary layer current velocities via a drag law in simulations with both parameterizations for idealized as well as realistic cavity geometries under pine island ice shelf west antarctica significant differences in melt rate patterns between the velocity independent and velocity dependent formulations are found while patterns are strongly correlated to those of thermal forcing for velocity independent transfer coefficients melting in the case of velocity dependent coefficients is collocated with regions of high boundary layer currents in particular where rapid plume outflow occurs both positive and negative feedbacks between melt rates boundary layer temperature velocities and buoyancy fluxes are identified melt rates are found to increase with increasing drag coefficient inline image in agreement with plume model simulations but optimal values of cd inferred from plume models are not easily transferable uncertainties therefore remain both regarding simulated melt rate spatial distributions and magnitudes

the role of coseismic coulomb stress changes in shaping the hard link between normal fault segments the mechanism and evolution of fault linkage is important in the growth and development of large faults here we investigate the role of coseismic stress changes in shaping the hard links between parallel normal fault segments or faults by comparing numerical models of the coulomb stress change from simulated earthquakes on two en echelon fault segments to natural observations of hard linked fault geometry we consider three simplified linking fault geometries 1 fault bend 2 breached relay ramp and 3 strike slip transform fault we consider scenarios where either one or both segments rupture and vary the distance between segment tips fault bends and breached relay ramps are favoured where segments underlap or when the strike perpendicular distance between overlapping segments is less than 20 of their total length matching all documented examples transform fault linkage geometries are preferred when overlapping segments are laterally offset at larger distances few transform faults exist in continental extensional settings and our model suggests that propagating faults or fault segments may first link through fault bends or breached ramps before reaching sufficient overlap for a transform fault to develop our results suggest that coulomb stresses arising from multi segment ruptures or repeated earthquakes are consistent with natural observations of the geometry of hard links between parallel normal fault segments

analogue modelling of orthogonal and oblique rifting abstract the three dimensional progressive development of orthogonal and oblique rift systems has been simulated using scaled analogue sandbox models the orientation of the zone of rifting where angle between the rift axis and the extension direction was varied systematically from 90 orthogonal rifting to 75 60 and 45 oblique rifting completed models were impregnated and sectioned both vertically and horizontally orthogonal rift models are characterized by long straight rift border faults and intra rift faults perpendicular to the extension direction oblique rift models are characterized by shorter segmented rift border faults and intrarift faults in the oblique models the rift border faults formed en echelon arrays parallel to the inderlying zone of basement stretching whereas the intra rift faults formed at a high angle to the extension vector increased extension in the moderately oblique and oblique rift systems 60 and 45 caused the intra rift faults to rotate towards parallelism with the rift border faults no strike slip or oblique slip transfer faults were found in any experiments rift border faults and intra rift faults in the models are segmented along strike with relay ramp structures formed between overlapping fault tips along the strike of the rift variations in half graben polarities and offsets of depocentres are accomplished by accommodation zones formed by interlocking conjugate fault arrays the experimental fault patterns compare well with the fault architectures found in those modern or relatively young rift systems where well defined basement controlled strike slip or oblique slip transfer faults are absent an accommodation zone model consisting of interlocking conjugate extension faults is proposed for the along strike switching of half graben asymmetries and offsetting of depocentres in rift systems the fault geometries in the analogue models can be used as templates for recognition of fault styles in both orthogonal and in oblique rift environments

evidence from high resolution topography for multiple earthquakes on high slip to length fault scarps the bilila mtakataka fault malawi geomorphological features such as fault scarps and stream knickpoints are indicators of recent fault activity determining whether these features formed during a single earthquake or over multiple earthquakes cycles has important implications for the interpretation of the size and frequency of past events here we focus on the bilila mtakataka fault malawi where the 20 m high fault scarps exceed the height expected from a single earthquake rupture we use a high resolution digital elevation model 1 m to identify complexity in the fault scarp and knickpoints in river profiles of 39 selected scarp profiles 20 showed evidence of either multi scarps or composite scarps and of the seven selected river and stream profiles five showed evidence for multiple knickpoints a near uniform distribution of vertical offsets on the sub scarps suggests they were formed by separate earthquakes these independent methods agree that at least two earthquakes have occurred with an average vertical offset per event of 10 and 12 m this contrasts earlier studies which proposed that this scarp formed during a single event and demonstrates the importance of high resolution topographic data for understanding tectonic geomorphology we use a one dimensional diffusion model of scarp degradation to demonstrate how fault splays form multi scarps and estimate the diffusion age t of the bilila mtakataka fault scarp to be 48 m25 m2 corresponding to 6400 4000 years since formation we calculate that a continuous rupture would equate to a mw 7 8 0 3 earthquake greater than the largest seismic event previously recorded in east africa

lithospheric structure of the malawi rift implications for magma poor rifting processes this project was supported by the national science foundation continental dynamics nsf cd grant ear 1255233 and office of international and integrative activities oiia grant iia 1358150 this project is also supported by the nsf earthcube integration grant 1740627 we thank the computational infrastructure for geodynamics geodynamics org which is funded by the national science foundation under award ear 0949446 and ear 1550901 for supporting the development of advanced solver for problems in earth s convection

estimation of ice shelf melt rate in the presence of a thermohaline staircase abstractdiffusive convection favorable thermohaline staircases are observed directly beneath george vi ice shelf antarctica a thermohaline staircase is one of the most pronounced manifestations of double diffusive convection cooling and freshening of the ocean by melting ice produces cool freshwater above the warmer saltier water the water mass distribution favorable to a type of double diffusive convection known as diffusive convection while the vertical distribution of water masses can be susceptible to diffusive convection none of the observations beneath ice shelves so far have shown signals of this process and its effect on melting ice shelves is uncertain the melt rate of ice shelves is commonly estimated using a parameterization based on a three equation model which assumes a fully developed unstratified turbulent flow over hydraulically smooth surfaces these prerequisites are clearly not met in the presence of a thermohaline staircase the basal melt rate is estimated by applying an ex

structural evolution of lake malaŵi africa over 2 000 km of seismic reflection profiles have been collected in lake mala i nyasa the southernmost of the east african rift lakes these studies reported here reveal a severely faulted lake floor bounded by rift structures of pliocene or greater age faults that penetrate the uppermost sedimentary units and numerous earthquakes1 attest to continuing activity in the rift inter basinal differences in fault spacing orientation and sedimentary thicknesses have led to a subdivision of lake mala i into four discrete structural provinces these provinces only now are linked by faults attributed to the most recent episodes of extension an evolutionary sequence of 1 initial block faulting and subsidence 2 fragmentation of fault blocks 3 uplift and 4 renewed subsidence and rotation of fault blocks is proposed for the central part of the mala i rift

oblique extensional tectonics in the malawi rift africa the malawi rift is the southern segment of the western branch of the east african rift system there has been little agreement as to whether this rift is opening perpendicularly or obliquely to its north south axis and little agreement on a mechanism for the formation of the asymmetric basins that make up the rift our analysis combines landsat image interpretation and field measurement of the latest slip direction on faults to produce a kinematic model of rift evolution the fault slip data are consistently compatible with a late cenozoic northwest southeast directed divergence between the nubian and somalian blocks plates across the rift the northwest striking livingstone fault at the northern limit of the malawi rift and the northwest striking zambezi fault zone to the south serve as dextral intracontinental transform zones between left stepping segments of the western branch of the east african rift system within the rift second order northwest striking transfer faults separate adjacent asymmetric basins the transform and transfer faults form a pattern that is similar to overlapping right stepping strike slip faults with a releasing geometry fault surfaces with two sets of striations indicate a progression from predominantly strike slip toward oblique normal faulting we hypothesize that during initial rift formation most northwest striking faults were dextral and that fault interaction between these structures resulted in curved stress trajectories and a pull apart component arcuate north to northeast striking border faults and subsiding basins then formed between these northwest striking faults in response to the right stepping dextral shear the maximum distance over which fault interaction is significant may control basin dimensions as opening of the rift proceeds slip partitioning and gravity failure result in increasing normal slip perpendicular to major escarpments and decreasing strike slip parallel to them

impact of tide topography interactions on basal melting of larsen c ice shelf antarctica 1 basal melting of ice shelves around antarctica contributes to formation of antarctic bottom water and can affect global sea level by altering the offshore flow of grounded ice streams and glaciers tides influence ice shelf basal melt rate wb by contributing to ocean mixing and mean circulation as well as thermohaline exchanges with the ice shelf we use a three dimensional ocean model thermodynamically coupled to a nonevolving ice shelf to investigate the relationship between topography tides andwb for larsen c ice shelf lcis in the northwestern weddell sea antarctica using our best estimates of ice shelf thickness and seabed topography we find that the largest modeled lcis melt rates occur in the northeast where our model predicts strong diurnal tidal currents 0 4 m s 1 this distribution is significantly different from models with no tidal forcing which predict largest melt rates along the deep grounding lines we compare several model runs to explore melt rate sensitivity to geometry initial ocean potential temperature 0 thermodynamic parameterizations of heat and freshwater ice ocean exchange and tidal forcing the resulting range of lcis averagedwb is 0 11 0 44 m a 1 the spatial distribution of wb is very sensitive to model geometry and thermodynamic parameterization while the overall magnitude of wb is influenced by 0 these sensitivities in wbpredictions reinforce a need for high resolution maps of ice draft and sub ice shelf seabed topography together with ocean temperature measurements at the ice shelf front to improve representation of ice shelves in coupled climate system models

melting of ice shelves and the mass balance of antarctica we calculate the present ice budget for antarctica from measurements of accumulation minus iceberg calving run off and in situ melting beneath the floating ice shelves the resulting negative mass balance of 469 gt year 1 differs substantially from other recent estimates but some components are subject to high temporal variability and budget uncertainties of 20 50 annual accumulation from an earlier review is adjusted to include the antarctic peninsula for a total of 2144 gt year 1 an iceberg production rate of 2016 gt year 1 is obtained from the volume of large icebergs calculated from satellite images since 1978 and from the results of an international iceberg census project ice shelf melting of 544 gt year 1 is derived from physical and geochemical observations of meltwater outflow glaciological field studies and modeling of the sub ice ocean circulation the highest melt rates occur near ice fronts and deep within sub ice cavities run off from the ice sheet surface and from beneath the grounded ice is taken to be 53 gt year 1 less than half of the negative mass balance need come from the grounded ice to account for the unattributed 0 45 mm year 1 in the ipcc best estimate of the recent global sea level rise

an exploration of the role of transient eddies in determining the transport of a zonally reentrant current the meridional ekman transport in a zonally reentrant channel may be balanced by diabatic circulations standing eddies associated with topography or by lagrangian mean eddy mass fluxes a simple model is used to explore the interaction between these mechanisms a key assumption of this study is that diabatic forcing in the poleward edge of the channel acts to create lighter fluid as is the case with net freshwater fluxes into the southern ocean for weak wind forcing or strong diabatic constraint a simple scaling argument accurately predicts the level of baroclinic shear however given our understanding of the relative magnitudes of ekman flux and deep upwelling this is not the appropriate parameter range for the antarctic circumpolar current with stronger wind stresses eddies are prominent with baroclinic instability initially developing in the vicinity of large topography arguments have been advanced by a number of authors that baroclinic instability should limit the velocity shear leading to a stiff upper limit on the transport of the current however in the simulations presented here baroclinic instability is largely confined to the region of topographic highs and the approach to a current that is independent of the wind stress occurs gradually several recent parameterizations of transient eddy fluxes do not reproduce key features of the observed behavior

high resolution electrical resistivity and aeromagnetic imaging reveal the causative fault of the 2009 mw 6 0 karonga malawi earthquake

numerical modeling of ocean ice interactions under pine island bay s ice shelf 1 a two dimensional numerical model is used to simulate the dynamics of buoyant meltwater rich plumes flowing beneath the ice shelf occupying much of pine island bay west antarctica recent studies have shown that this ice shelf along with all others fringing the amundsen sea is thinning rapidly in the model both the coriolis effect and subshelf topography are important in controlling plume dynamics and the spatial distribution of ice melt melt is concentrated in a narrow zone within 20 km of the grounding line where steep subshelf slopes and access to warm ambient water allow melt rates to exceed 100 m yr 1 the plume generated by entrainment of ambient water into the meltwater in these areas is guided by the topography of the ice shelf underside and exits the ice shelf at three distinct outflow locations melt rates generated along the course of the plume are higher approximately 2 5 than rates elsewhere the model suggests that the observed ice shelf thinning rates could have resulted from a hypothetical instantaneous 0 25 c warming of the ambient water entrained by the plume a context for this value is provided by the 40 year warming trend documented by jacobs et al 2002 for circumpolar deep water in the nearby ross sea

the role of pre existing precambrian structures in rift evolution the albertine and rhino grabens uganda abstract we integrated shuttle radar topography mission srtm digital elevation models dem airborne magnetic radiometric and three dimensional full tensor gravity gradiometry 3d ftg data to investigate the role of precambrian structures in the evolution of the largely amagmatic miocene recent aged albertine and rhino grabens in uganda these grabens represent the northern segment of the western branch of the east african rift system ears the two ne trending grabens are connected by a right stepping transfer zone and they extend within the archean paleoproterozoic northeast congo block which represents the northeastern extension of the congo craton our results show the following and highlight the importance of pre existing structures in the evolution of continental rift systems 1 the ne extent of the albertine full graben is controlled by ne trending precambrian fabric and the graben terminates at its northeastern end when it encounters a multiply folded precambrian basement terrain with poorly developed nw trending structural grain additionally the northeastern termination of the albertine graben coincides with the presence of nw trending right stepping high density bodies within the precambrian terrain 2 the transfer zone between the albertine and rhino grabens is controlled by ne trending precambrian structures which might have facilitated the development of relay ramp faults 3 strain localization within the better developed southeastern border fault of the rhino half graben is facilitated by the presence of precambrian structures better aligned in a ne direction in the southeastern part of the basin compared to its northwestern part 4 further to the northeast the rhino graben is segmented and transitions into a narrower ene trending half graben with a better developed border fault on its northwestern side this segmentation coincides with the presence of n trending precambrian structures 5 the rhino graben terminates farther northeast against the nw trending precambrian aswa shear zone a prominent structure with complex but generally nw trending fabric

modelling the response of ice shelf basal melting to different ocean cavity environmental regimes we present simulation results from a version of the regional ocean modeling system modified for ice shelf ocean interaction including the parameterisation of basal melting by molecular diffusion alone simulations investigate the differences in melting for an idealised ice shelf experiencing a range of cold to hot ocean cavity conditions both the pattern of melt and the location of maximum melt shift due to changes in the buoyancy driven circulation in a different way to previous studies tidal forcing increases both the circulation strength and melting with the strongest impact on the cold cavity case our results highlight the importance of including a complete melt parameterisation and tidal forcing in response to the 2 4 c ocean warming initially applied to a cold cavity ice shelf we find that melting will increase by about an order of magnitude 24 with tides and 41 without tides

investigation of land ice ocean interaction with a fully coupled ice ocean model 1 model description and behavior 1 antarctic ice shelves interact closely with the ocean cavities beneath them with ice shelf geometry influencing ocean cavity circulation and heat from the ocean driving changes in the ice shelves as well as the grounded ice streams that feed them we present a new coupled model of an ice stream ice shelf ocean system that is used to study this interaction the model is capable of representing a moving grounding line and dynamically responding ocean circulation within the ice shelf cavity idealized experiments designed to investigate the response of the coupled system to instantaneous increases in ocean temperature show ice ocean system responses on multiple timescales melt rates and ice shelf basal slopes near the grounding line adjust in 1 2 years and downstream advection of the resulting ice shelf thinning takes place on decadal timescales retreat of the grounding line and adjustment of grounded ice takes place on a much longer timescale and the system takes several centuries to reach a new steady state during this slow retreat and in the absence of either an upward or downward sloping bed or long term trends in ocean heat content the ice shelf and melt rates maintain a characteristic pattern relative to the grounding line

large scale oceanographic constraints on the distribution of melting and freezing under ice shelves previous studies suggest that ice shelves experience asymmetric melting and freezing topography may constrain oceanic circulation and thus basal melt freeze patterns through its influence on the potential vorticity pv field however melting and freezing induce a local circulation that may modify locations of heat transport to the ice shelf this paper investigates the influence of buoyancy fluxes on locations of melting and freezing under different bathymetric conditions an idealized set of numerical simulations the decoupled simulations employs spatially and temporally fixed diapycnal fluxes these experiments in combination with scaling considerations indicate that while flow in the interior is governed by large scale topographic gradients recirculation plumes dominate near buoyancy fluxes thermodynamically decoupled models are then compared to those in which ice ocean heat and freshwater fluxes are driven by the interior flow the coupled simulations near the southern boundary strong cyclonic flow forced by melt induced upwelling drives inflow and melting to the east recirculation is less evident in the upper water column as shoaling of meltwater freshened layers dissipates the dynamic influence of buoyancy forcing yet freezing remains intensified in the west in coupled simulations the flow throughout the cavity is relatively insensitive to bathymetry stratification the slope of the ice shelf and strong meridionally distributed buoyancy fluxes weaken its influence

modeling thermodynamic ice ocean interactions at the base of an ice shelf abstract models of ocean circulation beneath ice shelves are driven primarily by the heat and freshwater fluxes that are associated with phase changes at the ice ocean boundary their behavior is therefore closely linked to the mathematical description of the interaction between ice and ocean that is included in the code an hierarchy of formulations that could be used to describe this interaction is presented the main difference between them is the treatment of turbulent transfer within the oceanic boundary layer the computed response to various levels of thermal driving and turbulent agitation in the mixed layer is discussed as is the effect of various treatments of the conductive heat flux into the ice shelf the performance of the different formulations that have been used in models of sub ice shelf circulation is assessed in comparison with observations of the turbulent heat flux beneath sea ice formulations that include an explicit parameterization of the oceanic boundary layer give results that

observation and parameterization of ablation at the base of ronne ice shelf antarctica parameterizations of turbulent transfer through the oceanic boundary layer beneath an ice shelf are tested using direct measurements of basal ablation observations were made in the southwestern part of ronne ice shelf about 500 km from open water the mean basal ablation rate was measured over a month long and a year long period using phase sensitive radar to record the thinning of the ice shelf ocean temperatures were observed within about 25 m of the ice shelf base over the period of the radar observations while the tidally dominated ocean currents were estimated from tidal analysis of collocated current observations from an earlier period ablation rates derived using these ocean data and a number of bulk parameterizations of turbulent transfer within the boundary layer are compared with the direct measurements the ablation rates derived using a parameterization that explicitly includes the impact of ocean currents on the turbulent transfer of heat and salt match the observations to within 40 with suitable tuning of the drag coefficient the mismatch can be reduced below the level of the observational errors equally good agreement can be obtained with two slightly simpler current dependent parameterizations that use constant turbulent transfer coefficients andtheoptimalvaluesforthecoefficientsatthis particularlocationon ronneice shelfaregiven

sill influenced exchange flows in ice shelf cavities abstractbathymetric sills are important features in the ocean filled cavities beneath a few fast retreating ice shelves in west antarctica and northern greenland the sills can be high enough to ob

role of weak zone orientation in continental lithosphere extension 1 lithosphere extension and continental breakup axes are often sub parallel to orogenic belts and suture zones in an attempt to understand the relation between weak zone orientation and extension direction a 3 d numerical model is developed with which the following aspects are studied the relation between weak zone axis and axis of maximum crustal deformation the amount and distribution of crustal thinning and subsequent implications for rift a symmetry the results suggest that upon oblique extension rifts develop within the weak zone that individually cross the inherited structure while as a group they follow the weak trend this results in an alternating rift asymmetry the pattern of crustal thinning topography and thermal structure are not symmetric around the rift axes of the major rift zones and change along axis over such a rift this is a possible explanation for alternating asymmetric rift structures such as observed in the east african rift system

tectonic inheritance and continental rift architecture numerical and analogue models of the east african rift system 1 the western branch of the east african rift is composed of an arcuate succession of elongate asymmetric basins which differ in terms of interaction geometry fault architecture and kinematics and patterns of uplift subsidence and erosion sedimentation the basins are located within proterozoic mobile belts at the edge of the strong tanzanian craton surface geology suggests that the geometry of these weak zones is an important parameter in controlling rift development and architecture although other processes have been proposed in this study we use lithosphere scale numerical models and crustal scale analogue experiments to shed light on the relations between preexisting structures and rift architecture results illustrate that on a regional scale rift localization within the mobile belts at the curved craton s western border results in an arcuate rift system which implies that under a constant extensional stress field part of the western branch experienced orthogonal extension and part oblique extension largest depocenters are predicted to form mostly orthogonal to the extension direction and smaller depocenters will form along the oblique parts of the rift the varying extension direction along the rift zone furthermore results in lengthwise varying rift asymmetry segmentation characteristics and border fault architecture trend length and kinematics analogue models predict that discrete upper crustal fabrics may influence the location of accommodation zones and control the architecture of extension related faults at a local scale models support that fabric reactivation is responsible for the oblique slip kinematics on faults and for the development of z shaped or arcuate normal faults typically documented in nature

crustal structure beneath the malawi and luangwa rift zones and adjacent areas from ambient noise tomography abstract crustal shear wave velocity structure beneath the malawi and luangwa rift zones mrz and lrz respectively and adjacent regions in southern africa is imaged using fundamental mode rayleigh waves recorded by 31 safari seismic arrays for african rift initiation stations dispersion measurements estimated from empirical green s functions are used to construct 2 d phase velocity maps for periods between 5 and 28 s the resulting rayleigh wave phase velocities demonstrate significant lateral variations and are in general agreement with known geological features and tectonic units within the study area subsequently we invert rayleigh wave phase velocity dispersion curves to construct a 3 d shear wave velocity model beneath the mrz and lrz low velocity anomalies are found in the upper most crust probably reflecting the sedimentary cover the mid crust of the mrz is characterized by an 3 7 low velocity anomaly which cannot be adequately explained by higher than normal temperatures alone instead other factors such as magmatic intrusion partial melting and fluid filled deep crustal faults might also play a role thinning of the crust of a few kilometers beneath the rifts is revealed by the inversion a compilation of crustal thicknesses and velocities beneath the world s major continental rifts suggests that both the mrz and lrz are in the category of rifts beneath which the crust has not been sufficiently thinned to produce widespread syn rifting volcanisms

thermohaline circulation and interaction between ice shelf cavities and the adjacent open ocean the circulation system in an ice shelf cavity is driven by buoyancy fluxes due to melting and freezing of ice and horizontal pressure gradients at the interface between the cavity and the open ocean hence the inflow and outflow pattern and the hydrography in the open ocean influence the general hydrographic condition in the cavity which at least provides the potential for melting and freezing processes applying a three dimensional ocean general circulation model to an idealized ice shelf cavity geometry coupled with an open ocean at a topographic ice shelf barrier we found an important parameter controlling the interaction between these two systems idealized studies for different ice shelf and sea bottom topographies and forcing mechanisms for the open ocean show that the ice shelf edge represents a natural barrier for barotropic interaction because of the sudden decrease in water column thickness since the water column thickness and the coriolis force determine the characteristics for geostrophic flow separated circulation systems arise for the open ocean and the ice shelf cavity only in areas where constant water column thickness and from the oceanographic point of view constant f h contours can be observed across the barrier an increased barotropic current can surmount the ice edge and ventilate the water mass beneath the ice shelf this is only the case at lateral sloping sidewalls or at deep depressions which can be found for example in the southern weddell sea in all other cases the circulation in the ice shelf cavity is closed and almost unaffected by the hydrography outside the barrier

recent dramatic thinning of largest west antarctic ice stream triggered by oceans a growing body of observational data suggests that pine island glacier pig is changing on decadal or shorter timescales these changes may have far reaching consequences for the future of the west antarctic ice sheet wais and global sea levels because of pig s role as the ice sheet s primary drainage portal we test the hypothesis that these changes are triggered by the adjoining ocean specifically we employ an advanced numerical ice flow model to simulate the effects of perturbations at the grounding line on pig s dynamics the speed at which these changes are propagated upstream implies a tight coupling between ice sheet interior and surrounding ocean

oceanographic conditions south of berkner island beneath filchner ronne ice shelf antarctica we have made oceanographic measurements at two sites beneath the southern filchner ronne ice shelf hot water drilled access holes were made during january 1999 allowing conductivity temperature depth ctd profiling and the deployment of instrument moorings the ctd profiles show that the entire water column is below the surface freezing point we estimate the summer flux of water between the two sites to be 2 106 m3 s 1 the summer potential temperature salinity properties of the water column suggest that this flow is part of a recirculation in the deepest part of the subice shelf cavity and the filchner depression the recirculation is driven by a combination of the melting of deep basal ice and the freezing that results from the depressurization of the cold buoyant water as it ascends the ice shelf base the source of the water was high salinity shelf water hssw produced in the ronne depression this is the water that provides the external heat necessary for the strong melting at the deep grounding lines in the vicinity of foundation ice stream instruments moored at the drill sites show that during the winter hssw formed on the berkner shelf flows beneath the ice shelf and largely displaces the recirculating water from the two sites this provides an externally driven through flow that is warmer nearer the surface freezing point and slower than the internal recirculation and which is low enough in density to escape the filchner depression

seismic anisotropy and mantle dynamics beneath the malawi rift zone east africa sks skks and pks splitting parameters measured at 34 seismic stations that we deployed in the vicinity of the cenozoic malawi rift zone mrz of the east african rift system demonstrate systematic spatial variations with an average splitting time of 1 0 0 3 s the overall ne sw fast orientations are consistent with absolute plate motion apm models of the african plate constructed under the assumption of no net rotation of the global lithosphere and are inconsistent with predicted apm directions from models employing a fixed hotspot reference frame they also depart considerably from the trend of most of the major tectonic features these observations together with the results of anisotropy depth estimation using the spatial coherency of the splitting parameters suggest a mostly asthenospheric origin of the observed azimuthal anisotropy the single layered anisotropy observed at 30 and two layered anisotropy observed at four of the 34 stations can be explained by apm related simple shear within the rheologically transitional layer between the lithosphere and asthenosphere as well as by the horizontal deflection of asthenospheric flow along the southern and western edges of a continental block with relatively thick lithosphere revealed by previous seismic tomography and receiver function investigations this first regional scale shear wave splitting investigation of the mrz suggests the absence of rifting related active mantle upwelling or small scale mantle convection and supports a passive rifting process for the mrz

the effects of rotation and ice shelf topography on frazil laden ice shelf water plumes a model of the dynamics and thermodynamics of a plume of meltwater at the base of an ice shelf is presented such ice shelf water plumes may become supercooled and deposit marine ice if they rise because of the pressure decrease in the in situ freezing temperature so the model incorporates both melting and freezing at the ice shelf base and a multiple size class model of frazil ice dynamics and deposition the plume is considered in two horizontal dimensions so the influence of coriolis forces is incorporated for the first time it is found that rotation is extremely influential with simulated plumes flowing in near geostrophy because of the low friction at a smooth ice shelf base as a result an ice shelf water plume will only rise and become supercooled and thus deposit marine ice if it is constrained to flow upslope by topography this result agrees with the observed distribution of marine ice under filchner ronne ice shelf antarctica in addition it is found that the model only produces reasonable marine ice formation rates when an accurate ice shelf draft is used implying that the characteristics of real ice shelf water plumes can only be captured using models with both rotation and a realistic topography

growth of a normal fault system observations from the lake malawi basin of the east african rift abstract we have studied the growth history of the usisya normal fault system which bounds the west side of the lake malawi basin one of the largest rift basins in the east african rift system the lake malawi rift basin is a complex of intrabasins and intrabasin highs formed by three major fault segments each about 100 km in length the basin has been actively subsiding since the late miocene and has a maximum depth of about 3 km unlike the majority of previous studies of fault systems we are able to define the temporal evolution of this fault system using the patterns of sediment infill this is due almost exclusively to the availability of a high density seismic reflection network in combination with several clearly identifiable temporal marker beds an analysis of the seismic reflection data reveals that growth of the basin bounding faults occurred in the following sequence 1 in the early stages of rifting starting about 8 6 ma the northern fault was the most active 2 then extension shifted to the southernmost fault segment and lasted until about 2 5 ma 3 during the interval between about 2 3 ma and 1 6 ma the central segment was most active prior to the last interval the central segment accrued the least total displacement of the three segments this contradicts the common notion that the location of maximum displacement remains at a fixed location from fault inception or that the largest fault in any population of faults always maintained the highest displacement rate consistent with observations on smaller faults there is a marked increase in the displacement gradient on individual faults in the regions of overlapping segments although there is an observable asymmetry in individual segment displacement profiles there is a clear evolution toward a flattened bell shaped total displacement profile for the fault system that is consistent with the shape and scaling relationships of displacement vs length that are observed on a wide range of individual normal faults this suggests that the irregularities in the shape and scaling relationships observed in complex fault systems will eventually smooth out moreover the observed growth pattern suggests that the profile of the fault system will progress toward that of isolated faults provided the system is allowed sufficient time to evolve this in turn describes a process for maintaining a self similar scaling observed in large populations of faults which span several orders of magnitude in fault length

turbulent heat exchange between water and ice at an evolving ice water interface we conduct laboratory experiments on the time evolution of an ice layer cooled from below and subjected to a turbulent shear flow of warm water from above our study is motivated by observations of warm water intrusion into the ocean cavity under antarctic ice shelves accelerating the melting of their basal surfaces the strength of the applied turbulent shear flow in our experiments is represented in terms of its reynolds number re which is varied over the range 2 0 times 10 3 leqslant re leqslant 1 0 times 10 4 depending on the water temperature partial transient melting of the ice occurs at the lower end of this range of re and complete transient melting of the ice occurs at the higher end following these episodes of transient melting the ice reforms at a rate that is independent of re we fit our experimental measurements of ice thickness and temperature to a one dimensional model for the evolution of the ice thickness in which the turbulent heat transfer is parameterized in terms of the friction velocity of the shear flow applying our model to field measurements at a site under the antarctic pine island glacier ice shelf yields a predicted melt rate that exceeds present day observations

distribution of marine ice beneath the amery ice shelf we present a map of the marine ice accreted to the base of the amery ice shelf ais east antarctica this map is obtained by converting a digital elevation model dem of the ais generated from satellite radar altimeter data to an ice thickness map assuming hydrostatic equilibrium and subtracting from that a second ice thickness map derived from airborne radio echo sounding res measurements the res signal does not penetrate the marine ice so the measurement is only to the meteoric marine ice boundary and therefore the difference between the two maps is the marine ice thickness the marine ice is up to 190 m thick and accounts for about 9 of the shelf volume it is concentrated in the northwest of the shelf a result of the clockwise ocean circulation in the cavity below

computing marine ice thickness at an ice shelf base the freezing of sea water to the base of an ice shelf can give rise to large patches of accumulated ice a phenomenon known as marine ice in this study a numerical method is presented for calculating the thickness of the marine ice layer using an ice shelf ocean model the present day modeling paradigm of ice shelf ocean interaction usually involves the fixed specification of the ice shelf geometry while the ocean circulation in the cavity beneath the ice shelf evolves freely this approach relies on several assumptions such as steady state ice shelf thickness and ice shelf flow fields in order to make reasonable quantitative estimates of the thermodynamic exchange processes occurring at the ice shelf base this paper discusses the impact of these and other assumptions on the estimation of the thickness of the marine ice layer model simulation results are presented for an idealized ice shelf ocean configuration as a demonstration of the feasibility of the numerical method a sensitivity analysis is given so as to quantify the relative uncertainty in the marine ice thickness that arises from uncertainties in the model input parameters these being principally the ice shelf flow field the basal accumulation rate and the ice shelf thickness field

buoyancy driven circulation in an ocean basin with isopycnals intersecting the sloping boundary abstract the dynamics that govern the spreading of a convectively formed water mass in an ocean with sloping boundaries are examined using an isopycnal model that permits the interface between the layers to intersect the sloping boundaries the simulations presented here use a two layer configuration to demonstrate some of the pronounced differences in a baroclinically forced flow between the response in a basin with a flat bottom and vertical walls and a more realistic basin bounded by a sloping bottom each layer has a directly forced signal that propagates away from the forcing along the potential vorticity pv contours of that layer paired opposed boundary currents are generated by refracted topographic rossby waves rather than kelvin waves it is impossible to decompose the flow into globally independent baroclinic and barotropic modes topography causes the barotropic i e depth averaged response to buoyancy forcing to be just as strong as the baroclinic response because layer pv contours di

evidence for basal marine ice in the filchner ronne ice shelf the filchner ronne ice shelf which drains most of the marine based portions of the west antarctic ice sheet is the largest ice shelf on earth by volume the origin and properties of the ice that constitutes this shelf are poorly understood because a strong reflecting interface within the ice and the diffuse nature of the ice ocean interface make seismic and radio echo sounding data difficult to interpret1 2 ice in the upper part of the shelf is of meteoric origin but it has been proposed2 5 that a basal layer of saline ice accumulates from below here we present the results of an analysis of the physical and chemical characteristics of an ice core drilled almost to the bottom of the ronne ice shelf we observe a change in ice properties at about 150 m depth which we ascribe to a change from meteoric ice to basal marine ice the basal ice is very different from sea ice formed at the ocean surface and we propose a formation mechanism in which ice platelets in the water column accrete to the bottom of the ice shelf

on the connection between the mediterranean outflow and the azores current as the salty and dense mediteranean overflow exits the strait of gibraltar and descends rapidly in the gulf of cadiz it entrains the fresher overlying subtropical atlantic water a minimal model is put forth in this study to show that the entrainment process associated with the mediterranean outflow in the gulf of cadiz can impact the upper ocean circulation in the subtropical north atlantic ocean and can be a fundamental factor in the establishment of the azores current two key simplifications are applied in the interest of producing an economical model that captures the dominant effects the first is to recognize that in a vertically asymmetric twolayer system a relatively shallow upper layer can be dynamically approximated as a single layer reduced gravity controlled barotropic system and the second is to apply quasigeostrophic dynamics such that the volume flux divergence effect associated with the entrainment is represented as a source of potential vorticity two sets of computations are presented within the 1 layer framework a primitive equation based computation which includes the divergent flow effects is first compared with the equivalent quasigeostrophic formulation the upper ocean cyclonic eddy generated by the loss of mass over a localized area elongates westward under the influence of the b effect until the flow encounters the western boundary in the steady state the circulation pattern consists of bidirectional zonal flows with a limited meridional extent eastward to the south of the sink and westward to the north of the sink the localized sink drives a horizontal circulation in the interior ocean whose strength is approximately an order of magnitude greater than the sink s strength it is demonstrated that the induced circulation in the far field from a localized sink is insensitive to the neglect of the divergent flow component a set of parameter sensitivity experiments is then undertaken with the quasigeostrophic model for an idealized midlatitude circulation driven both by wind forcing and thermohaline flow through the open southern and northern boundaries when a sink near the eastern boundary is superimposed on the idealized midlatitude circulation it is shown to alter significantly the upper ocean flow and induce an eastward zonal current which resembles the azores current in location and transport this mechanism also generates a westward current to the north of the sink location which could be associated with the azores countercurrent an extensive series of sensitivity experiments is conducted to determine the response of this current system to changes in the boundary layer processes sink strength sink distribution model resolution and wind forcing the magnitude of the current transports is found to be sensitive to the sink intensity and to its distance from the coastline

normal vs strike slip faulting during rift development in east africa the malawi rift kinematic analysis of neogene and quaternary faults demonstrates that the direction of extension in the malawi rift rotated from east northeast to southeast rift development commenced with the formation of half grabens bounded by northwest north and northeast striking normal faults owing to slightly oblique rifting the northwest striking faults in the northernmost rift segment show a small dextral oblique slip component whereas north and northeast oriented faults in the central part of the rift display a sinistral oblique slip component this first event resulted in block faulting and basin subsidence which is largely responsible for the present day basin morphology of lake malawi a major change in fault kinematics occurred because of rotation of the extension direction and permutation of the principal compressive and intermediate axes the structural pattern inherited from the first rifting phase was no longer suitably oriented to accommodate extensional deformation and strike slip faulting assumed a major role the strike slip regime amplified uplift of basement ridges within the rift in regions of local transpression but it also created alluvial basins because of local transtension this new kinematic style is compatible with the recent seismicity older faults that show mainly the first deformational increment are restricted to the outermost parts of the rift toward the center the faults depict an increase in strike slip components of movement suggesting deformation propagation toward the rift center which results in a narrowing of the active rift environments with time

is the south pacific helium 3 plume dynamically active it is suggested that the hydrothermal vents of the south pacific rise produce a beta governed circulation at mid depth and that perhaps the associated plume of excess3he lupton and craig 1 points westward because of the dynamics of this circulation rather than as a passive tracer

adaptation of an isopycnic coordinate ocean model for the study of circulation beneath ice shelves abstract much of the antarctic coastline comprises large floating ice shelves beneath which waters from the open ocean circulate the interaction of the seawater with the base of these ice shelves has a bearing both on the rate at which antarctic bottom water is formed and on the mass balance of the ice sheet an isopycnic coordinate ocean general circulation model has been modified so as to allow the incorporation of a floating ice shelf as an upper boundary to the model domain the modified code admits the introduction of an arbitrary surface pressure field and includes new algorithms for the diagnosis of entrainment into and detrainment from the surface mixed layer special care is needed in handling the cases where the mixed layer and isopycnic interior layers interact with surface and basal topography the modified model is described in detail and then applied to an idealized ice shelf ocean geometry simple tests with zero surface buoyancy forcing indicate that the introduction of the static s

the role of inherited crustal and lithospheric architecture during the evolution of the red sea insights from three dimensional analogue experiments abstract the red sea represents an incipient ocean basin that is in transition from continental rifting to ocean initiation in this paper we assess how pre existing mantle and crustal heterogeneities within the african and arabian lithospheres influenced rift development and propagation using three dimensional 3d brittle ductile analogue laboratory experiments with a rotational boundary condition here we investigate the effect of three different configurations i pre existing rheological heterogeneities in the crust ii in the lithospheric mantle and iii in both the experimental outcomes show that a linear weak zone in the lithospheric mantle layer exerts a first order control on the overall mechanical behaviour of extending lithosphere rift architecture and style of rifting crustal heterogeneities control second order features such as the trends of the main boundary escarpments along the rift axis our experiments reproduce the geological history of the red sea and provide key insights on its deformation history in space and time we provide for the first time 3d dynamic constraints on how pre existing weak zones in the lithosphere have controlled the temporal evolution of the red sea and determined its distinctive present day morphology the results find that certain morphological elements of the red sea can only be explained by the combined effects of crustal and mantle lithosphere heterogeneities when applied to other natural examples our findings provide significant new insights on the role of crustal and lithospheric heterogeneities during continental rifting

breaking up the hanging wall of a rift border fault the 2009 karonga earthquakes malawi 1 the southern east african rift has an unusually large seismogenic thickness 35 40 km which is responsible for wide tilted basins and extremely long faults with the potential for m7 8 normal faulting earthquakes from 6 19 december 2009 a shallow earthquake sequence four of mw 5 5 hit the karonga region of northern lake malawi the location is 50 km west of the rift bounding livingstone fault within the hanging wall we used seismology and insar to obtain source parameters and combined this with information on rift structure from geomorphology and seismic profiles the deformation is consistent with rupture of a shallow west dipping fault with no evidence for the involvement of magmatic fluids although the livingstone fault dominates local geomorphology the karonga earthquakes demonstrate that the hanging wall block is actively breaking up reflecting temporal and spatial migration of activity or the release of stresses within it

redox states of archean surficial environments the importance of h2 g instead of o2 g for weathering reactions abstract redox states of the archean eon have been constrained by various lines of evidence including atmospheric photochemical and ecological models mass independent fractionations of sulfur isotopes fe depletion of paleosols and preservation of diagnostic detrital minerals although these lines of evidence present seemingly consistent upper limits on po2 g they are conceptually contradictory about the redox state of archean surficial environments atmospheric photochemical and ecological modeling studies suggest weakly reducing environments under redox states represented by moderate h2 g levels however current interpretations of fe depletion in paleosols and the preservation of detrital minerals are based on low o2 g levels at which the reducing detrital minerals are thermodynamically unstable and survive because of slow kinetics of oxidative weathering in this study we show that under the redox state indicated by the archean ph2 g range fe2 and reducing fe ii minerals are actually thermodynamically stable and have no tendency to be oxidized we emphasize that ph2 g values are orders of magnitude higher than po2 g in the archean atmosphere and that h2 g and o2 g are not in equilibrium the redox states of archean surface environments behave as though they were controlled by the more abundant h2 instead of the very low o2 as in modern anoxic basins weathering in this case should have involved non redox acidic dissolution of fe ii species or reductive reaction of fe iii species fe ii depleted paleosols and the preservation of relatively reduced detrital minerals are natural consequences of their thermodynamic stabilities in the archean eon s reducing environments rather than slow kinetics of oxidizing reactions after the appearance of oxygenic photosynthesis probably in the middle late archean locally oxygenated environments could have existed while the atmosphere as a whole remained anoxic the profile of redox states on the archean surface seems to be a reverse analogue to the modern earth although oxidizing dissolution of transition metals could happen in o2 oases quick reduction of oxyanions by abundant reductants such as aqueous fe2 fe ii minerals and h2 might have restricted riverine transport of oxyanions and potentially complicate the interpretation of signals of o2 whiffs in marine sediments

active fault scarps in southern malawi and their implications for the distribution of strain in incipient continental rifts the distribution of deformation during the early stages of continental rifting is an important constraint on our understanding of continental breakup incipient rifting in east africa has been considered to be dominated by slip along rift border faults with a subsequent transition to focussed extension on axial segments in thinned crust and or with active magmatism here we study high resolution satellite data of the zomba graben in southern malawi an amagmatic rift whose topography is dominated by the west dipping zomba fault we document evidence for five sub parallel fault scarps between 13 and 51 km long spaced 10 15 km apart the scarps consist of up to five segments between 4 18 km long separated by minima in scarp height and river knickpoints the maximum height of each fault scarp ranges from 9 5 4 2 m to 35 3 14 6 m with the highest scarp measured on the intrabasin chingale step fault we estimate that the scarps were formed by multiple earthquakes of up to mw7 1 and represent a previously unrecognized seismic hazard our calculations show that 55 24 of extensional strain is accommodated across intrabasin faults within the 50 km wide rift this demonstrates that a significant proportion of displacement can occur on intrabasin faults during early stage rifting even in thick continental lithosphere with no evidence for magmatic fluids

4 d evolution of rift systems insights from scaled physical models the four dimensional 4 d evolution of brittle fault systems in orthogonal oblique and offset rift systems has been simulated by scaled sandbox models using dry cohesionless fine grained quartz sand extensional deformation in the models was controlled by the orientation and geometry of a zone of stretching at the base of the model the results of these analog model studies are compared with natural examples of rift fault systems rift basins produced by orthogonal and oblique rifting are defined by segmented border fault systems parallel to the rift axes and by intrarift fault systems that are subperpendicular to the extension direction segmentation of the rift margin increases with increase in obliquity of the rift axis resulting in a consequent increase in displacement on intrarift fault systems offset rift models are characterized by highly segmented border faults and offset subbasins in the rift zone along strike displacement transfer in the model rifts occurred as a result of formation of two types of accommodation zones high relief extension parallel accommodation zones typically are found in 60 degrees rifts and above left steps in offset rift systems changes in fault polarities in these accommodation zones were achieved by interlocking arrays of conjugate extensional faults the second type of accommodation zone was generally oblique to the extension direction and consisted of conjugate fault arrays having rotated tips that bounded a low relief oblique slip zone or grabens these typically are found in highly oblique rift systems 45 degrees and above right steps in offset rift models

tectonic controls on rift basin morphology evolution of the northern malawi nyasa rift radiometric k ar and ar 40 ar 39 age determinations of volcanic and volcaniclastic rocks combined with structural gravity and seismic reflection data are used to constrain the age of sedimentary strata contained within the seismically and volcanically active northern malawi nyasa rift and to characterize changes in basin and flank morphologies with time faulting and volcanism within the tukuyu karonga basin began at approximately 8 6 ma when sediments were deposited in abroad initially asymmetric lake basin bounded on its northeastern side by a border fault system with minor topographic relief extensions primarily by a slip along the border fault and subsequent regional isostatic compensation led to the development of a 5 km deep basin bounded by broad uplifted flanks along the low relief basin margin opposite border fault younger stratigraphic sequences commonly onlap older wedge shaped sequences although their internal geometry is often progradational intrabasinal faulting flankuplift and basaltic and felsic volcanism from centers at the northern end of the basin became more important at about 2 5 ma when cross rift transfer faults developed to link the tukuyu karonga basin to the rukwa basin local uplift and volcanic construction at the northern end of the basin led to a southeastward shift in the basin s depocenter sequence boundaries are commonly erosional along this low relief hanging wall margin and conformable in the deep lake basin the geometry of stratigraphic sequences and the distribution of the erosion indicate that horizontal and vertical crustal movements both across and along the length of the rift basin led to changes in levels of the lake irrespective of paleoclimatic fluctuations