Near surface measurements of seismic velocity and Q are useful in a number of situations, for example for when carrying out re-datuming and migration for depth images, or when analysing ground conditions for building. This thesis concentrates on the estimation of surface wave group and phase velocity as well as Q structure through the use of cross correlations of ambient noise recordings. Linearised tomography estimates are made for the British Isles, the Permain Basin of Texas and New Mexico, the Ekofisk Life of Field Seismic (LoFS) array and the Aquistore CO2 storage site. The results correspond well with the known geological structure and/or structure observed in velocity maps by other researchers. For the Ekofisk array a non-linearised non-linear method was also applied and the results estimated by these two methods for the Ekofisk LoFS array are compared. By non-linearised non-linear it is meant that the inversion method is fully non-linear and no linearisation has taken place in the method, this term will be used throughout this thesis for all methods which fall into this category. The tomography results from the two methods had similarities in their general structure but differences in the finer details, and so suggest that the substantial increase in time required for the non-linearised non-linear method is not justified. Linearised tomographic inversion of the Aquistore array was used to determine the potential of using ambient noise tomography for monitoring of CO2 storage sites. It was found that the repeatability of the tomography at the Aquistore site was not good enough to allow ambient noise tomography to be used for monitoring; however, it may be possible to apply the method at other sites. A Q and phase velocity inversion of the Ekofisk array is also presented, with results mostly showing excellent correlation with known geological features. It is shown that the higher frequencies are more sensitive to the effects of sea floor subsidence at the site, while lower frequencies are more sensitive to the effects of faulting. A final near surface method called ambient noise gradiometry was applied, this concentrates on estimating locations of sources of seismic energy within receiver arrays. Ambient noise gradiometry is applied to synthetic and real data for this purpose. It was found that using ambient noise gradiometry allows internal sources of energy to be identified but they produce a bias in the phase velocity tomography result. Two methods of reducing this bias are presented, both of which also provide an estimate of the source term for different sections of time of the recording.

The East African coastal basin of Tanzania, which is characterised by an extensional tectonic style, is located along the passive continental margin of the western Indian Ocean. The present study is concerned with the Mafia Island and the Mafia Channel which together form a subbasin within the north-south elongated coastal sedimentary basin of Tanzania. In the time interval from late Paleozoic to Recent, the passive margin of the region was subjected to a three-fold geological development, namely the Karoo rifting phase (1) which is characterised by extensional tectonics, the Gondwana break-up and opening of the Somali basin (2) which was contemporaneous with the movement of Madagascar off the east African coast in the Mesozoic, and the Cenozoic East African rift system (3). This structural framework made provision to the basin deposition history. The development started with the deposition of the continental, terrigenous, Karoo sequence in the Upper Permian to Lower Jurassic. The Karoo deposition was followed by a series of transgressions and regressions under full marine conditions which started in the Middle Jurassic and continued into the Tertiary. The deposits include marine marls, detrital limestones, fossiliferous shales and calcareous sandstones, reaching in places thicknesses of more than 4000 m of Mesozoic, and more than 6000 m of Cenozoic, sediments. Seismostratigraphic techniques applied in the Mafia Channel and Island identified five deposition sequences separated by regional unconformity surfaces i.e. sequence boundaries. Including the pre-Upper Cretaceous sequence they are; the Upper Cretaceous to Middle Eocene sequence (DS I), the Middle Eocene to Lower Miocene sequence (DS II), the Lower Miocene to Pliocene sequence (DS III), and the Pliocene to Recent sequence (DS IV). In the Mafia Channel up to 6000 m of sediments are present. Their ages range from Middle Eocene to Quaternary. The deposits start with marine shales which are overlain by carbonate rocks of Upper Eocene . These carbonates are in the present investigation regarded as the acoustic basement in the central and northern parts of the study area. On top of Upper Eocene carbonates, deltaic and shallow marine sediments are deposited. Southwards in the Channel, the sequences are located at shallower depths which makes it possible to trace also the Upper Cretaceous sequence with confidence. On the Mafia Island, the deposition on top of the Upper Cretaceous (Campanian) acoustic basement mainly includes deltaic sandstones, followed by intercalations of carbonate and argillaceous rocks. The structural framework reflects the different tectonic regimes which prevailed in the area. Above the acoustic basement structural elements of Mafia Channel and Island are interpreted as originating from the superimposition of the Cenozoic East African rifting event, and from the uplifts of the mainland coast and Mafia Island during Late Eocene time. As a result the central part collapsed and formed an asymmetric sag structure in the channel. These elements are seismically identified and subdivided as (from northwest to southeast), the Dar-es-Salaam Platform Offshore, the Central Mafia Channel and the Mafia Island Rise. These domains are separated by respective NE-SW major faults (MF1, MF2 and MF3) and can be demonstrated in the profiles which run in a NW-SE direction. With respect to petroleum potential, the Mafia Channel and Island indicate a considerable content. Three hydrocarbon plays are introduced, namely; 1) the Neocomian regressive sands of Songosongo play 2) the Upper Eocene limestone and 3) the Upper Oligocene turbidites.

Black shales are sedimentary rocks with a high content of organic carbon, which leads to a dark grayish to black color. Due to their potential to contain oil or gas, black shales are of great interest for the support of the worldwide energy supply. An integrated seismic investigation of the Lower Palaeozoic black shales was carried out at the Danish island Bornholm to locate the shallow-lying Alum Shale layer and its surrounding formations and to characterize its potential as a source rock. Therefore, two seismic experiments at a total of three crossing profiles were carried out in October 2010 and in June 2012 in the southern part of the island. Two different active measurements were conducted with either a weight drop source or a minivibrator. Additionally, the ambient noise field was recorded at the study location over a time interval of about one day, and also a laboratory analysis of borehole samples was carried out. The seismic profiles were positioned as close as possible to two scientific boreholes which were used for comparative purposes. The seismic field data was analyzed with traveltime tomography, surface wave inversion and seismic interferometry to obtain the P-wave and S-wave velocity models of the subsurface. The P-wave velocity models which were determined for all three profiles clearly locate the Alum Shale layer between the Komstad Limestone layer on top and the Læså Sandstone Formation at the base of the models. The black shale layer has P-wave velocities around 3 km/s which are lower compared to the adjacent formations. Due to a very good agreement of the sonic log and the vertical velocity profiles of the two seismic lines, which are directly crossing the borehole where the sonic log was conducted, the reliability of the traveltime tomography is proven. A correlation of the seismic velocities with the content of organic carbon is an important task for the characterization of the reservoir properties of a black shale formation. It is not possible without calibration but in combination with a full 2D tomographic image of the subsurface it gives the subsurface distribution of the organic material. The S-wave model obtained with surface wave inversion of the vibroseis data of one of the profiles images the Alum Shale layer also very well with S-wave velocities around 2 km/s. Although individual 1D velocity models for each of the source positions were determined, the subsurface S-wave velocity distribution is very uniform with a good match between the single models. A really new approach described here is the application of seismic interferometry to a really small study area and a quite short time interval. Also new is the selective procedure of only using time windows with the best crosscorrelation signals to achieve the final interferograms. Due to the small scale of the interferometry even P-wave signals can be observed in the final crosscorrelations. In the laboratory measurements the seismic body waves were recorded for different pressure and temperature stages. Therefore, samples of different depths of the Alum Shale were available from one of the scientific boreholes at the study location. The measured velocities have a high variance with changing pressure or temperature. Recordings with wave propagation both parallel and perpendicular to the bedding of the samples reveal a great amount of anisotropy for the P-wave velocity, whereas the S-wave velocity is almost independent of the wave direction. The calculated velocity ratio is also highly anisotropic with very low values for the perpendicular samples and very high values for the parallel ones. Interestingly, the laboratory velocities of the perpendicular samples are comparable to the velocities of the field experiments indicating that the field measurements are sensitive to wave propagation in vertical direction. The velocity ratio is also calculated with the P-wave and S-wave velocity models of the field experiments. Again, the Alum Shale can be clearly separated from the adjacent formations because it shows overall very low vP/vS ratios around 1.4. The very low velocity ratio indicates the content of gas in the black shale formation. With the combination of all the different methods described here, a comprehensive interpretation of the seismic response of the black shale layer can be made and the hydrocarbon source rock potential can be estimated. / Schwarzschiefer sind Sedimentgesteine, die einen hohen Gehalt an organischem Kohlenstoff aufweisen, was zu einer dunkelgrauen bis schwarzen Färbung führt. Da Schwarzschiefer das Potenzial besitzen, Öl oder Gas zu enthalten und somit zur weltweiten Energieversorgung beitragen könnten, sind sie von großem Interesse. Mit Hilfe der Kombination verschiedener seismischer Messverfahren wurden die Schwarzschiefer des Unteren Paläozoikums auf der dänischen Insel Bornholm untersucht um den oberflächennahen Alaunschiefer und dessen Umgebungsgestein dort zu lokalisieren und sein Potenzial als Muttergestein abzuschätzen. Dafür wurden im Oktober 2010 und im Juni 2012 im südlichen Teil der Insel zwei seismische Experimente auf insgesamt drei sich kreuzenden Profilen durchgeführt. Für zwei aktive seismische Messungen wurden ein Fallgewicht und ein Minivibrator als Quellen genutzt. Zusätzlich wurde im Messgebiet noch das Wellenfeld des umgebenden Rauschens über einen Zeitraum von etwa einem Tag aufgezeichnet. Außerdem wurden Labormessungen an Bohrkernen aus dem Alaunschiefer durchgeführt. Die seismischen Messprofile befanden sich so nah wie möglich an zwei wissenschaftlichen Bohrungen, die für Vergleichszwecke genutzt wurden. Um die P- und S-Wellengeschwindigkeitsmodelle des Untergrundes zu erhalten wurden die seismischen Felddaten mittels Laufzeittomographie, Oberflächenwelleninversion und seismischer Interferometrie ausgewertet. Die P-Wellenmodelle, die für alle drei seismischen Profile erstellt wurden, zeigen den Alaunschiefer zwischen dem Komstad Kalkstein, der den Alaunschiefer überdeckt, und der Læså Sandsteinformation, die die Basis der Modelle bildet. Für die Schwarzschieferschicht ergeben sich mit rund 3 km/s deutlich geringere P-Wellengeschwindigkeiten als für die umgebenden Gesteine. Zwei seismische Profile liegen direkt an einer der Bohrungen, für die verschiedene Bohrloch-Logs durchgeführt wurden. Der Vergleich des Sonic-Logs mit den vertikalen Geschwindigkeitsprofilen beider Modelle am Bohrpunkt zeigt eine sehr gute übereinstimmung aller Geschwindigkeiten. Dies ist ein Indiz für die Plausibilität der durchgeführten Laufzeittomographie. Um die Reservoireigenschaften der Schwarzschieferschicht einordnen zu können, wurde versucht, die seismischen Geschwindigkeiten mit dem Gehalt an organischem Material zu korrelieren. Ohne geeignete Kalibrierung ist diese Korrelation schwierig, kann aber mit Hilfe der Tomographieergebnisse ein zweidimensionales Abbild der Verteilung des organischen Materials im Untergrund liefern. Auch das S-Wellengeschwindigkeitsmodell, welches mit der Oberflächenwelleninversion der Vibroseisdaten erstellt wurde, bildet den Alaunschiefer gut ab. Hierbei zeigen sich S-Wellengeschwindigkeiten um 2 km/s. Obwohl jeweils nur 1D-Modelle für jede Quellposition bestimmt wurden, ergibt sich für die gesamte Untergrundstruktur des untersuchten Profils ein einheitliches Bild der Geschwindigkeiten. Einen sehr neuen Ansatz bildet die Anwendung der seismischen Interferometrie auf ein sehr kleines Untersuchungsgebiet und über einen sehr kurzen Zeitraum. Neu ist außerdem, dass für die Bestimmung der endgültigen Interferogramme nur Zeitfenster der Kreuzkorrelationen ausgewählt werden, in denen die Signalqualität hinreichend gut ist. In den berechneten Kreuzkorrelationen sind sogar P-Wellen enthalten, was auf die geringen Abstände der seismischen Rekorder zurück zu führen ist. Bei den Labormessungen wurden die Raumwellen für verschiedene Drücke und Temperaturen aufgezeichnet. Die Messungen der Geschwindigkeiten sowohl parallel als auch senkrecht zur Schichtung der Proben zeigen eine starke Anisotropie für die P-Welle. Dagegen scheint die S-Wellengeschwindigkeit fast unabhängig von der Ausbreitungsrichtung der Wellen zu sein. Auch das Verhältnis der Geschwindigkeiten weist starke Anisotropie auf. Für die Wellenausbreitung senkrecht zur Schichtung zeigen sich sehr niedrige Werte, die Werte für die Messungen parallel zur Schichtung sind dagegen deutlich erhöht. Ein interessanter Aspekt der aus den Labormessungen resultiert ist, dass die Geschwindigkeit der Messungen senkrecht zur Schichtung mit den Geschwindigkeitswerten der Feldmessungen übereinstimmen. Damit scheinen die Feldmessungen besonders die Ausbreitung der Wellen in vertikaler Richtung zu registrieren. Das Geschwindigkeitsverhältnis wurde auch mit den P- und S-Wellenmodellen der Feldexperimente berechnet. Auch hier hebt sich der Alaunschiefer mit deutlich verringerten Werten um 1.4 vom Umgebungsgestein ab. Solch geringe Werte für das Verhältnis der Geschwindigkeiten deutet auf den Gehalt von Gas im Schwarzschiefer. Mit der Kombination der verschiedenen Methoden ist es möglich, die seismische Antwort der Schwarzschieferschicht umfassend zu beschreiben und Schlussfolgerungen darüber zu ziehen, ob die hier untersuchte Schwarzschieferschicht das Potenzial hat als Kohlenwasserstofflagerstätte zu fungieren.

O presente trabalho de mestrado está inserido nos estudos de refração profunda do Projeto Temático 'Estudos Geofísicos e Modelo Tectônico dos Setores Central e Sudeste da Província Tocantins, Brasil Central'. Nesses estudos foram levantadas três linhas de refração de aproximadamente 300 km de extensão, duas no setor Central da Província Tocantins e uma no setor Sudeste, que é o objeto de estudo deste trabalho. Foram utilizados 111 sismógrafos digitais SGR pertencentes ao programa PASSCAL, instrumentos auxiliares do USGS, e 13 sismógrafos digitais e instrumentos auxiliares do IAG/USP. A linha sísmica teve aproximadamente 300 km de extensão com pontos de registro separados a cada 2,5 km, distribuídos ao longo de estradas principais e secundárias. A cada 50 km, aproximadamente, foi realizada uma explosão, nas explosões dos extremos da linha foram utilizados 1000 kg de explosivo e para a explosão central uma carga de 500 kg. Para a determinação das coordenadas geográficas dos pontos de tiro e de registro, foi utilizado o método diferencial com medidas de GPS. O principal objetivo deste trabalho foi obter como produto final um modelo de velocidades sísmicas contendo as características físicas das principais descontinuidades na crosta terrestre e no manto superior. Para análise e processamento dos dados foram utilizados os pacotes SAC, SU, SEIS. Para a modelagem foram utilizados a teoria do raio e a elaboração de sismogramas sintéticos, do pacote SEIS. Para a elaboração do modelo final foram utilizados os dados das explosões dos pontos extremos e central, tendo em vista que devido a problemas técnicos não foram registrados os sinais das outras 4 explosões. Além disso, as explosões registradas não apresentaram sinais claros em toda a extensão da linha. Devido a tudo isso e considerando as unidades geológicas presentes na região de estudo são sugeridos três modelos de velocidades sísmicas. O primeiro modelo refere-se ao tiro direto (EX31) localizado no extremo sudoeste da linha, sobre a Bacia do Paraná. Para este modelo obteve-se para superfície (0 km) a velocidade inicial de 2 km/s (coberturas); para a profundidade de 0,086 km a velocidade inicial é de 5,15 km/s (basalto); para a profundidade de 0,350 km obteve-se a velocidade inicial de 4,6 km/s (arenito - camada de baixa velocidade); para profundidade de 0,650 km a velocidade inicial é de 5,75 km/s e para profundidade de 4 km obteve-se a velocidade inicial de 6,07 km/s. O segundo modelo refere-se ao tiro reverso (EX34) localizado no centro da linha sobre granitóides do Grupo Araxá. Para este modelo obteve -se para superfície (0 km) a velocidade inicial de 2 km/s; para a profundidade de 0,06 km a velocidade inicial é de 5,69 km/s e para a profundidade de 0,860 km obteve-se a velocidade inicial de 6,25 km/s. Finalmente, o terceiro modelo refere-se ao tiro direto para toda a extensão da linha (300 km). Este modelo foi definido a partir de fases secundárias lidas nos registros e modelos anteriores propostos na literatura. Da superfície até os 4 km iniciais de profundidade este modelo é igual ao primeiro, para uma profundidade de 20 km obteve-se a velocidade inicial de 6,70 km/s e para uma profundidade de 40 km a velocidade é de 8,00 km/s (descontinuidade de MOHO). / This work to fulfil the degree of Master of Sciences is inserted among the deep seismic refraction studies of the Thematic Project 'Geophysical Studies and Tectonic Model of the Tocantins Province Central and Southeast Sectors, Central Brazil'. Three refraction lines, of around 300 km long each, were deployed, two of them in the Central sector and the other in the SE sector, that is subject of the present work. The equipment used in this experiment was composed by 111 SGR digital seismographs belonging to the PASSCAL Program. Complemented with auxiliary instruments from USGS and 13 seismographs belonging to IAG/USP. The space among the recording points was 2.5 km, which were located along main and secondary roads. Every 50 km was fired an explosion with 1000 kg of emulsion in each extreme and 500 kg in the central point. The geographical co-ordinates were determined by using the GPS differential method. The main objective of this work is to obtain as a final product a seismic velocity model with the physical characteristics of the main discontinuities in the crust and upper mantle. The packages SAC, SU and SEIS were used to perform the data analysis and processing. To carry on the modelling were used the ray theory and the synthetic seismograms construction, belonging to the SEIS package Data from the extreme and middle points of the seismic line were used to elaborate the final model, considering that due to technical problems signals from the other four explosions were not recorded. Apart from that, the recorded explosions did not present clear signals all along the extension of the line. Due to these facts, and considering also the geological units present in the studied region, are suggested three seismic velocity models. The first model is referred to the direct shot (EX31), which is localised in the Southwest extreme of the line on the Parana Basin province. In this model we obtained the P wave velocity (VP) of 2 km/sec at the surface, corresponding to the unconsolidated sediments and soil on the top of that basin. At a depth of 86 m we found VP of 5,15 km/sec and at a depth of 350 m the velocity VP of 4,6 km/sec, corresponding to the basalt and sand layers of the Parana Basin. Underlying them, at 650 m of depth we found the basement with VP of 5,75 km/sec and finally at a depth of 4 km there is a layer with VP of 6,07 km/sec, corresponding to a typical upper crust P wave velocity. The second model corresponds to the reverse shot (EX34) that is localised in the middle point of the line on the granitoides of the Araxa Group. For this model we obtained VP of 2 km/sec for the superficial layers, then at a depth of 60 m was obtained V P of 5,69 km/sec and for a depth of 860 m the value of V P is 6,25 km/sec. Finally, the third model belongs to the whole line section (300 km) from the direct shot (EX31). This model was obtained by using the arrivals of secondary phases and the results of models proposed in other works. From the surface down to 4 km of depth this model is similar to the first one. At 20 km of depth there is a layer with VP of 6,70 km/sec, corresponding to the lower crust, with Moho at a depth of 40 km with VP of 8,00 km/sec.

<div>Ensuring the resilience of a state’s transportation network is necessary to guarantee an acceptable quality of life for the people the network serves. A lack of resilience in the wake of a seismic event directly impacts the states’ overall safety and economic vitality. With the recent identification of the Wabash Valley Seismic Zone (WBSV), Department of Transportations (DOTs) like Indiana’s have increased awareness for the vulnerability of their bridge network. The Indiana Department of Transportation (INDOT) has been steadily working to reduce the seismic vulnerability of bridges in the state in particular in the southwest Vincennes District. In the corridor formed by I-69 built in the early 2000s the bridge design is required to consider seismic actions. However, with less recent bridges and those outside the Vincennes District being built without consideration for seismic effects, the potential for vulnerability exists. As such, the objective of this thesis is to develop a robust seismic vulnerability assessment methodology which can assess the overall vulnerability of Indiana’s critical bridge network. </div><div><br></div><div>A representative sample of structures in Indiana’s bridge inventory, which prioritized the higher seismic risk areas, covered the entire state geographically, and ensured robust superstructure details, was chosen. The sample was used to carry a deterministic seismic vulnerability assessment, applicable to all superstructure-substructure combinations. Analysis considerations, such as the calculation of critical capacity measures like moment-curvature and a pushover analysis, are leveraged to accurately account for non-linear effects like force redistribution. This effect is a result of non-simultaneous structural softening in multi-span bridges that maintain piers of varying heights and stiffnesses. These analysis components are incorporated into a dynamic analysis to allow for the more precise identification of vulnerable details in Indiana’s bridge inventory.</div><div><br></div><div>The results of this deterministic seismic assessment procedure are also leveraged to identify trends in the structural response of the sample set. These trends are used to identify limit state thresholds for the development of fragility functions. This conditional probabilistic representation of bridge damage is coupled with the probability of earthquake occurrence to predict the performance of the structure for a given return period. This probabilistic approach alongside a Monte Carlo simulation is applied to assess the vulnerability of linked bridges along key-access corridors throughout the state. With this robust seismic vulnerability methodology, DOTs will have the capability of identifying vulnerable corridors throughout the state allowing for the proactive prioritization of retrofits resulting in the improved seismic performance and resiliency of their transportation network.</div>

No description available.

Seismic waves are used to study the Earth, exploit its hydrocarbon resources, and understand its hazards. Extracting information from seismic waves about the Earth’s subsurface, however, is becoming more challenging as our questions become more complex and our demands for higher resolution increase. This dissertation introduces two new methods that use scattered waves for improving the resolution of subsurface images: natural migration of passive seismic data and convergent full-waveform inversion. In the first part of this dissertation, I describe a method where the recorded seismic data are used to image subsurface heterogeneities like fault planes. This method, denoted as natural migration of backscattered surface waves, provides higher resolution images for near-surface faults that is complementary to surface-wave tomography images. Our proposed method differ from contemporary methods in that it does not (1) require a velocity model of the earth, (2) assumes weak scattering, or (3) have a high computational cost. This method is applied to ambient noise recorded by the US-Array to map regional faults across the American continent. Natural migration can be formulated as a least-squares inversion to furtherer enhance the resolution and the quality of the fault images. This inversion is applied to ambient noise recorded in Long Beach, California to reveal a matrix of shallow subsurface faults. The second part of this dissertation describes a convergent full waveform inversion method for controlled source data. A controlled source excites waves that scatter from subsurface reflectors. The scattered waves are recorded by a large array of geophones. These recorded waves can be inverted for a high-resolution image of the subsurface by FWI, which is typically convergent for transmitted arrivals but often does not converge for deep reflected events. I propose a preconditioning approach that extends the ability of FWI to image deep parts of the velocity model, which significantly improves the chances for finding hydrocarbon deposits.

I address the problem of finding a background model for the estimation of resistivities in the earth from controlled-source electromagnetic (CSEM) data by using seismic data and well logs as constraints. Estimation of resistivities is normally done by trial-and-error, in a process called “inversion”, by finding a model of the earth whose responses match the data to within an acceptable error; what comes out of the inversion is what is put into the model by the geophysicist: it does not come out of the data directly. The premise underlying this thesis is that an earth model can be found that satisfies not only the CSEM data but also the seismic data and any well logs. I present a methodology to determine background resistivities from seismic velocities using rock physics, structural constraints, and depth trends. The physical parameters of the seismic wave equation are different from those in the electromagnetic diffusion equation, so there is no direct link between the governing equations. I therefore use a Bayesian framework to incorporate not only the errors in the data and our limited knowledge of the rock parameters, but also the uncertainty of our chosen and calibrated velocity-to-resistivity transform. To test the methodology I use a well log from the North Sea Harding South oil and gas field to calibrate the transform, and apply it to seismic velocities of the nearby Harding Central oil and gas field. I also use short-offset CSEM inversions to estimate the electric anisotropy and to improve the shallow part of the resistivity model, where there is no well control. Three-dimensional modelling of this resistivity model predicts the acquired CSEM data within the estimated uncertainty. This methodology makes it possible to estimate background resistivities from seismic velocities, well logs, and other available geophysical and geological data. Subsequent CSEM surveys can then focus on finding resistive anomalies relative to this background model; these are, potentially, hydrocarbon-bearing formations.

No description available.

Deep tectonic tremor has been observed in major subduction zones, strike-slip faults, inland faulting systems, and arc-continent collision environments around the Pacific Rim. However, detailed space-time evolution of its source locations remains enigmatic because of difficulties in detecting and locating tremor accurately. In 2011, we installed two dense, small-aperture seismic arrays aiming to detect ambient tremor source beneath southern Central Range in Taiwan. We recorded continuous waveforms for a total of 134 days, including tremor triggered by the great 2011 Mw9.0 Tohoku earthquake. We use the broadband frequency-wavenumber beamforming and the moving-window grid-search methods to compute array parameters for detecting seismic signals. The obtained array parameters closely match both relocated local earthquakes and triggered tremor bursts located by an envelope cross-correlations method, indicating the robustness of our array technique. We identify tremor signals with coherent waveforms and deep incidence angles and detect tremor for 44 days among the 134-day study period. The total duration is 1,481-minute, which is 3-6 times more than that detected by the envelope cross-correlations method. In some cases, we observe rapid tremor migration with a speed at the order of 40-50 km/hour that is similar to the speed of fast tremor migration along-dip on narrow streaks in Japan and Cascadia. Our results suggest that dense array techniques are capable of capturing detailed spatiotemporal evolutions of tremor behaviors in southern Taiwan.