P-Wave Velocity in Precambrian Basement (Pg) and P-Wave Time Terms in and Around Albuquerque Basin, New Mexico

Total Page:16

File Type:pdf, Size:1020Kb

P-Wave Velocity in Precambrian Basement (Pg) and P-Wave Time Terms in and Around Albuquerque Basin, New Mexico P-wave velocity in Precambrian basement (Pg) and P-wave time terms in and around Albuquerque Basin, New Mexico Lawrence H. Jaksha, New Mexico Institute of Mining and Technology (Adjunct faculty, retired), Socorro, New Mexico 87801 Abstract Geology T = a + b + D/V (1) Where: A time-term method of analyzing seismic The Albuquerque Basin is one of a number T is the shot to station travel time. refraction data was used to process 62 of northerly elongate crustal features that a is the delay time associated with time-distance pairs observed on a 13 station in aggregate make up the Rio Grande the receiver. seismic network from 10 shot points in and depression as identified by Bryan (1938) and around Albuquerque Basin, New Mexico. b is the delay time associated with The time terms derived from the analysis described as the Rio Grande rift by Kelley the shot. range from 1.50 seconds at the Albuquerque (1952). The basin measures approximately D/V is the time the signal spends tra- volcanoes down to -.047 seconds at the 160 km long by 50 km in wide, and the surface versing the high-speed refractor Albuquerque Seismological Laboratory. In of the basin floor is roughly 1.6 km above sea (distance divided by velocity). a general way the time terms correlate with level. Parts of the basin probably began to the geologic foundation beneath the site: form during the Laramide, but most of the Delay time is defined by Sheriff (1973) as Large values are calculated for stations and features seen today are the result of Miocene the additional time taken for a wave to follow shots within the basin, intermediate values and younger faulting and uplift. The basin is a trajectory to and along a buried marker over on consolidated sedimentary rocks and linked through bedrock constrictions to the Tertiary intrusions, and the smallest on Pre- that which would have been taken to follow cambrian outcrops. The P-wave (Pg) velocity Española Basin on the north and to the San the same marker considered hypothetically to in bedrock was found to be 6.0 ± .02 km/s. Marcial Basin on the south. The uplifted east- have been at the ground surface. It is synony- This falls approximately midway between ern margin, fault blocks rotated to the east, mous with “time term” in this study. higher velocities reported to the east and consists of the Sandia–Manzano–Los Pinos An examination of equation (1) suggests west of Albuquerque Basin and lower veloc- Mountains. The less clearly defined western that if the coordinates of the shots and stations ities reported to the north and south. These margin is composed of Ladron Peak, Lucero are known, D can be calculated, and that if the numbers will be valuable in future efforts to uplift, and the Puerco fault belt. Deep hydro- shot origin time is known, T can be estimat- accurately locate shallow nearby earthquakes carbon test wells, seismic reflection profiles, ed from the arrival times. Figure 1 shows the and for studies seeking to unravel seismic and gravity/magnetic surveys in the basin complexities deeper in the crust. arrangement of explosions and stations used reveal a highly faulted and irregular basin in this study. There are 13 station time terms, bottom. These studies tally closely with the 10 shot time terms, and 1 velocity [the a, b, Introduction geologic fabric seen directly in the uplifted and V in equation (1)] to be determined. If margins (Kelley 1977). The basin seems trun- The U.S. Geological Survey in cooperation every station recorded every shot, 130 equa- cated at its northern end by the Miocene–late tions like (1) could be written, and there are with Los Alamos National Laboratory Pleistocene Jemez Mountains and the Valles (LANL) and New Mexico Institute of Mining only 24 unknowns so that the system of equa- caldera. Olsen et al. (1986) reported on a tions is over determined and has a solution. and Technology (NMT) operated a 13 station seismic experiment for the Jemez Mountains seismic network around and within the Multiple linear regression (Draper and Smith that analyzed 290 ray paths and found the 1966) was used to estimate the unknowns, Albuquerque Basin for a 66-month period P-wave velocity in bedrock to be 5.86 ± .01 between 1976 and 1981 with the goal of map- their associated standard errors, and the stan- km/s. The southernmost part of the Albu- dard deviation of the solution (Murdock and ping contemporary seismicity. In the course querque Basin is underlain by the northern of this earthquake monitoring, explosions Jaksha 1980). Berry and West (1966) discussed edge of the Socorro Magma Body described the errors associated with the time-term meth- originating from nearby mining operations by Sanford et al. (1977). Singer (1989) studied were routinely recorded, along with special od and stated that the estimates should be the entire crust beneath a seismic network taken as indicators of “goodness of fit” but shots set off for scientific research by the Uni- centered on Socorro and from an analysis versity of New Mexico, the U.S. Department not as accurate fiducial limit calculations. of 265 ray paths reported that the P-wave The time-term method assumes a starting of Defense, and LANL. These data were used velocity in bedrock was best modeled as to estimate some of the seismic parameters model that includes: increasing from 5.71 km/s at the top of the 1. A high velocity layer that is constant of the shallow crust beneath the network. Precambrian section to just under 6.0 km/s This report gives an estimate for the P-wave beneath the study area. at approximately 10 km depth. This study 2. Dip on the refractor surface is small. velocity in the Precambrian basement (Pg) fills the gap between the Valles caldera to and the P-wave time terms (delay times) 3. Curvature on the refractor is small. the north and the contemporary magma 4. The velocity beneath the sites varies associated with each of the recording stations injection proposed to the south of the basin. and shot sites along with the standard errors only with depth. for the computed values. These numbers The method produces relative time terms are important for the accurate location of Method unless a recording station and an explosion earthquakes and for deducing structure at occupy the same site, in which case absolute The time-term method of analyzing seismic deeper levels in the crust. The delay times values are calculated. In this study a station refraction data was proposed by Scheideg- also are related simply to the thickness of at the Occidental Mineral Corporation mine ger and Willmore (1957) and has been used the Phanerozoic section above the basement, (OXM) recorded the first arrival from explo- successfully by many others. The method or, where the site is already located on base- sions at LANL shot points SP5 and SP6, and assumes that the travel time for the first ment, they are related mainly to the amount on the White Sands Missile Range (WSMR) a arriving refracted signal from an explosion of confining pressure in the rocks beneath station at the site of DICE THROW recorded can be partitioned into 3 terms: the explosion or recording station. the shot at DISTANT RUNNER 2. November 2013, Volume 35, Number 4 NEW MEXICO GEOLOGY 95 the permanent seismic stations along with TAOS their elevations are given in Table 2. A SAN RIO ARRIBA description of the instruments in the field JUAN (short period vertical sensors and radio SP4 telemetry) and at the recording site is given EUM by Jaksha et al. (1977). The seismic data were recorded on 16 mm film displayed at Valles 1 cm per second. Two analysts timed the Los Alamos caldera refracted P-wave arrivals independently, Nacimiento Mts. Nacimiento SP6 TSP and if the results varied by more than ± 100 ms the data were rejected. Shot to station Mc KINLEY SP5 distances were computed using DISTAZ, an unpublished computer program from COH SANTA CDN FE the U.S. Geological Survey, Menlo Park, OXM California. The data were then assembled SANDOVAL into matrix algebra format and processed San Pedro- using a TI-89 calculator manufactured by MTL Ortiz porphyry belt Texas Instruments Inc. VOL GNM Residuals for each ray path were calculat- Albuquerque SAN JM MIGUEL ed, and those with values greater than 200 TJQ ms, the estimate for our ability to read the BERNALILLO CIBOLA seismograms, were reviewed, rejected, or ABQ EST both. These data are shown in Figure 2 and CRF Pedernal presented in tabular form in the Appendix. MLM Basin e Hills Albuquerqued Distance constraints were used to look for n ra aberrations associated with arrivals from G TORRANCE o shots too close (direct arrivals) or too far (a VALENCIA Ri deeper refractor) from the stations. It was found that Pg appears as a first arrival in LAD the range 10–150 km in and around the Albuquerque Basin and that its velocity is LPM 6.0 ± .02 km/s. This result is in agreement with a refraction line partially in the basin Plateau Colorado analyzed by Olsen et al. (1979). WTX Estimates for the station and shot time Socorro terms are listed in Table 3. The standard Great Great Plains LINCOLN error for the estimate is given along with Rio Grande rift SOCORRO the number of observations associated with each site. The standard deviation for the entire solution is 103 ms. The estimate for DT our ability to read the arrival times given above is 200 ms, so the solution appears to MRC be well determined.
Recommended publications
  • Velocity-Corrected Area Calculation SCIEX PA 800 Plus Empower
    Velocity-corrected area calculation: SCIEX PA 800 Plus Empower Driver version 1.3 vs. 32 Karat™ Software Firdous Farooqui1, Peter Holper1, Steve Questa1, John D. Walsh2, Handy Yowanto1 1SCIEX, Brea, CA 2Waters Corporation, Milford, MA Since the introduction of commercial capillary electrophoresis (CE) systems over 30 years ago, it has been important to not always use conventional “chromatography thinking” when using CE. This is especially true when processing data, as there are some key differences between electrophoretic and chromatographic data. For instance, in most capillary electrophoresis separations, peak area is not only a function of sample load, but also of an analyte’s velocity past the detection window. In this case, early migrating peaks move past the detection window faster than later migrating peaks. This creates a peak area bias, as any relative difference in an analyte’s migration velocity will lead to an error in peak area determination and relative peak area percent. To help minimize Figure 1: The PA 800 Plus Pharmaceutical Analysis System. this bias, peak areas are normalized by migration velocity. The resulting parameter is commonly referred to as corrected peak The capillary temperature was maintained at 25°C in all area or velocity corrected area. separations. The voltage was applied using reverse polarity. This technical note provides a comparison of velocity corrected The following methods were used with the SCIEX PA 800 Plus area calculations using 32 Karat™ and Empower software. For Empower™ Driver v1.3: both, standard processing methods without manual integration were used to process each result. For 32 Karat™ software, IgG_HR_Conditioning: conditions the capillary Caesar integration1 was turned off.
    [Show full text]
  • The Origins of Velocity Functions
    The Origins of Velocity Functions Thomas M. Humphrey ike any practical, policy-oriented discipline, monetary economics em- ploys useful concepts long after their prototypes and originators are L forgotten. A case in point is the notion of a velocity function relating money’s rate of turnover to its independent determining variables. Most economists recognize Milton Friedman’s influential 1956 version of the function. Written v = Y/M = v(rb, re,1/PdP/dt, w, Y/P, u), it expresses in- come velocity as a function of bond interest rates, equity yields, expected inflation, wealth, real income, and a catch-all taste-and-technology variable that captures the impact of a myriad of influences on velocity, including degree of monetization, spread of banking, proliferation of money substitutes, devel- opment of cash management practices, confidence in the future stability of the economy and the like. Many also are aware of Irving Fisher’s 1911 transactions velocity func- tion, although few realize that it incorporates most of the same variables as Friedman’s.1 On velocity’s interest rate determinant, Fisher writes: “Each per- son regulates his turnover” to avoid “waste of interest” (1963, p. 152). When rates rise, cashholders “will avoid carrying too much” money thus prompting a rise in velocity. On expected inflation, he says: “When...depreciation is anticipated, there is a tendency among owners of money to spend it speedily . the result being to raise prices by increasing the velocity of circulation” (p. 263). And on real income: “The rich have a higher rate of turnover than the poor. They spend money faster, not only absolutely but relatively to the money they keep on hand.
    [Show full text]
  • VELOCITY from Our Establishment in 1957, We Have Become One of the Oldest Exclusive Manufacturers of Commercial Flooring in the United States
    VELOCITY From our establishment in 1957, we have become one of the oldest exclusive manufacturers of commercial flooring in the United States. As one of the largest privately held mills, our FAMILY-OWNERSHIP provides a heritage of proven performance and expansive industry knowledge. Most importantly, our focus has always been on people... ensuring them that our products deliver the highest levels of BEAUTY, PERFORMANCE and DEPENDABILITY. (cover) Velocity Move, quarter turn. (right) Velocity Move with Pop Rojo and Azul, quarter turn. VELOCITY 3 velocity 1814 style 1814 style 1814 style 1814 color 1603 color 1604 color 1605 position direction magnitude style 1814 style 1814 style 1814 color 1607 color 1608 color 1609 reaction move constant style 1814 color 1610 vector Velocity Vector, quarter turn. VELOCITY 5 where to use kinetex Healthcare Fitness Centers kinetex overview Acute care hospitals, medical Health Clubs/Gyms office buildings, urgent care • Cardio Centers clinics, outpatient surgery • Stationary Weight Centers centers, outpatient physical • Dry Locker Room Areas therapy/rehab centers, • Snack Bars outpatient imaging centers, etc. • Offices Kinetex® is an advanced textile composite flooring that combines key attributes of • Cafeteria, dining areas soft-surface floor covering with the long-wearing performance characteristics of • Chapel Retail / Mercantile hard-surface flooring. Created as a unique floor covering alternative to hard-surface Wholesale / Retail merchants • Computer room products, J+J Flooring’s Kinetex encompasses an unprecedented range of • Corridors • Checkout / cash wrap performance attributes for retail, healthcare, education and institutional environments. • Diagnostic imaging suites • Dressing rooms In addition to its human-centered qualities and highly functional design, Kinetex • Dry physical therapy • Sales floor offers a reduced environmental footprint compared to traditional hard-surface options.
    [Show full text]
  • Chapter 3 Motion in Two and Three Dimensions
    Chapter 3 Motion in Two and Three Dimensions 3.1 The Important Stuff 3.1.1 Position In three dimensions, the location of a particle is specified by its location vector, r: r = xi + yj + zk (3.1) If during a time interval ∆t the position vector of the particle changes from r1 to r2, the displacement ∆r for that time interval is ∆r = r1 − r2 (3.2) = (x2 − x1)i +(y2 − y1)j +(z2 − z1)k (3.3) 3.1.2 Velocity If a particle moves through a displacement ∆r in a time interval ∆t then its average velocity for that interval is ∆r ∆x ∆y ∆z v = = i + j + k (3.4) ∆t ∆t ∆t ∆t As before, a more interesting quantity is the instantaneous velocity v, which is the limit of the average velocity when we shrink the time interval ∆t to zero. It is the time derivative of the position vector r: dr v = (3.5) dt d = (xi + yj + zk) (3.6) dt dx dy dz = i + j + k (3.7) dt dt dt can be written: v = vxi + vyj + vzk (3.8) 51 52 CHAPTER 3. MOTION IN TWO AND THREE DIMENSIONS where dx dy dz v = v = v = (3.9) x dt y dt z dt The instantaneous velocity v of a particle is always tangent to the path of the particle. 3.1.3 Acceleration If a particle’s velocity changes by ∆v in a time period ∆t, the average acceleration a for that period is ∆v ∆v ∆v ∆v a = = x i + y j + z k (3.10) ∆t ∆t ∆t ∆t but a much more interesting quantity is the result of shrinking the period ∆t to zero, which gives us the instantaneous acceleration, a.
    [Show full text]
  • Moment of Inertia
    MOMENT OF INERTIA The moment of inertia, also known as the mass moment of inertia, angular mass or rotational inertia, of a rigid body is a quantity that determines the torque needed for a desired angular acceleration about a rotational axis; similar to how mass determines the force needed for a desired acceleration. It depends on the body's mass distribution and the axis chosen, with larger moments requiring more torque to change the body's rotation rate. It is an extensive (additive) property: for a point mass the moment of inertia is simply the mass times the square of the perpendicular distance to the rotation axis. The moment of inertia of a rigid composite system is the sum of the moments of inertia of its component subsystems (all taken about the same axis). Its simplest definition is the second moment of mass with respect to distance from an axis. For bodies constrained to rotate in a plane, only their moment of inertia about an axis perpendicular to the plane, a scalar value, and matters. For bodies free to rotate in three dimensions, their moments can be described by a symmetric 3 × 3 matrix, with a set of mutually perpendicular principal axes for which this matrix is diagonal and torques around the axes act independently of each other. When a body is free to rotate around an axis, torque must be applied to change its angular momentum. The amount of torque needed to cause any given angular acceleration (the rate of change in angular velocity) is proportional to the moment of inertia of the body.
    [Show full text]
  • Rotation: Moment of Inertia and Torque
    Rotation: Moment of Inertia and Torque Every time we push a door open or tighten a bolt using a wrench, we apply a force that results in a rotational motion about a fixed axis. Through experience we learn that where the force is applied and how the force is applied is just as important as how much force is applied when we want to make something rotate. This tutorial discusses the dynamics of an object rotating about a fixed axis and introduces the concepts of torque and moment of inertia. These concepts allows us to get a better understanding of why pushing a door towards its hinges is not very a very effective way to make it open, why using a longer wrench makes it easier to loosen a tight bolt, etc. This module begins by looking at the kinetic energy of rotation and by defining a quantity known as the moment of inertia which is the rotational analog of mass. Then it proceeds to discuss the quantity called torque which is the rotational analog of force and is the physical quantity that is required to changed an object's state of rotational motion. Moment of Inertia Kinetic Energy of Rotation Consider a rigid object rotating about a fixed axis at a certain angular velocity. Since every particle in the object is moving, every particle has kinetic energy. To find the total kinetic energy related to the rotation of the body, the sum of the kinetic energy of every particle due to the rotational motion is taken. The total kinetic energy can be expressed as ..
    [Show full text]
  • Circular Motion Angular Velocity
    PHY131H1F - Class 8 Quiz time… – Angular Notation: it’s all Today, finishing off Chapter 4: Greek to me! d • Circular Motion dt • Rotation θ is an angle, and the S.I. unit of angle is rad. The time derivative of θ is ω. What are the S.I. units of ω ? A. m/s2 B. rad / s C. N/m D. rad E. rad /s2 Last day I asked at the end of class: Quiz time… – Angular Notation: it’s all • You are driving North Highway Greek to me! d 427, on the smoothly curving part that will join to the Westbound 401. v dt Your speedometer is constant at 115 km/hr. Your steering wheel is The time derivative of ω is α. not rotating, but it is turned to the a What are the S.I. units of α ? left to follow the curve of the A. m/s2 highway. Are you accelerating? B. rad / s • ANSWER: YES! Any change in velocity, either C. N/m magnitude or speed, implies you are accelerating. D. rad • If so, in what direction? E. rad /s2 • ANSWER: West. If your speed is constant, acceleration is always perpendicular to the velocity, toward the centre of circular path. Circular Motion r = constant Angular Velocity s and θ both change as the particle moves s = “arc length” θ = “angular position” when θ is measured in radians when ω is measured in rad/s 1 Special case of circular motion: Uniform Circular Motion A carnival has a Ferris wheel where some seats are located halfway between the center Tangential velocity is and the outside rim.
    [Show full text]
  • Lecture 24 Angular Momentum
    LECTURE 24 ANGULAR MOMENTUM Instructor: Kazumi Tolich Lecture 24 2 ¨ Reading chapter 11-6 ¤ Angular momentum n Angular momentum about an axis n Newton’s 2nd law for rotational motion Angular momentum of an rotating object 3 ¨ An object with a moment of inertia of � about an axis rotates with an angular speed of � about the same axis has an angular momentum, �, given by � = �� ¤ This is analogous to linear momentum: � = �� Angular momentum in general 4 ¨ Angular momentum of a point particle about an axis is defined by � � = �� sin � = ��� sin � = �-� = ��. � �- ¤ �⃗: position vector for the particle from the axis. ¤ �: linear momentum of the particle: � = �� �⃗ ¤ � is moment arm, or sometimes called “perpendicular . Axis distance.” �. Quiz: 1 5 ¨ A particle is traveling in straight line path as shown in Case A and Case B. In which case(s) does the blue particle have non-zero angular momentum about the axis indicated by the red cross? A. Only Case A Case A B. Only Case B C. Neither D. Both Case B Quiz: 24-1 answer 6 ¨ Only Case A ¨ For a particle to have angular momentum about an axis, it does not have to be Case A moving in a circle. ¨ The particle can be moving in a straight path. Case B ¨ For it to have a non-zero angular momentum, its line of path is displaced from the axis about which the angular momentum is calculated. ¨ An object moving in a straight line that does not go through the axis of rotation has an angular position that changes with time. So, this object has an angular momentum.
    [Show full text]
  • Chapter 8: Rotational Motion
    TODAY: Start Chapter 8 on Rotation Chapter 8: Rotational Motion Linear speed: distance traveled per unit of time. In rotational motion we have linear speed: depends where we (or an object) is located in the circle. If you ride near the outside of a merry-go-round, do you go faster or slower than if you ride near the middle? It depends on whether “faster” means -a faster linear speed (= speed), ie more distance covered per second, Or - a faster rotational speed (=angular speed, ω), i.e. more rotations or revolutions per second. • Linear speed of a rotating object is greater on the outside, further from the axis (center) Perimeter of a circle=2r •Rotational speed is the same for any point on the object – all parts make the same # of rotations in the same time interval. More on rotational vs tangential speed For motion in a circle, linear speed is often called tangential speed – The faster the ω, the faster the v in the same way v ~ ω. directly proportional to − ω doesn’t depend on where you are on the circle, but v does: v ~ r He’s got twice the linear speed than this guy. Same RPM (ω) for all these people, but different tangential speeds. Clicker Question A carnival has a Ferris wheel where the seats are located halfway between the center and outside rim. Compared with a Ferris wheel with seats on the outside rim, your angular speed while riding on this Ferris wheel would be A) more and your tangential speed less. B) the same and your tangential speed less.
    [Show full text]
  • Physics 3550, Fall 2012 Two Body, Central-Force Problem Relevant Sections in Text: §8.1 – 8.7
    Two Body, Central-Force Problem. Physics 3550, Fall 2012 Two Body, Central-Force Problem Relevant Sections in Text: x8.1 { 8.7 Two Body, Central-Force Problem { Introduction. I have already mentioned the two body central force problem several times. This is, of course, an important dynamical system since it represents in many ways the most fundamental kind of interaction between two bodies. For example, this interaction could be gravitational { relevant in astrophysics, or the interaction could be electromagnetic { relevant in atomic physics. There are other possibilities, too. For example, a simple model of strong interactions involves two-body central forces. Here we shall begin a systematic study of this dynamical system. As we shall see, the conservation laws admitted by this system allow for a complete determination of the motion. Many of the topics we have been discussing in previous lectures come into play here. While this problem is very instructive and physically quite important, it is worth keeping in mind that the complete solvability of this system makes it an exceptional type of dynamical system. We cannot solve for the motion of a generic system as we do for the two body problem. The two body problem involves a pair of particles with masses m1 and m2 described by a Lagrangian of the form: 1 2 1 2 L = m ~r_ + m ~r_ − V (j~r − ~r j): 2 1 1 2 2 2 1 2 Reflecting the fact that it describes a closed, Newtonian system, this Lagrangian is in- variant under spatial translations, time translations, rotations, and boosts.* Thus we will have conservation of total energy, total momentum and total angular momentum for this system.
    [Show full text]
  • 1 Circular Motion by Prof. Massimiliano Galeazzi, University of Miami If You Spin a Bucket Full of Water in a Vertical Circl
    Circular Motion By Prof. Massimiliano Galeazzi, University of Miami If you spin a bucket full of water in a vertical circle fast enough, the water will stay in the bucket. How fast do you need to spin it to make sure you don’t get showered? What is the force that makes your car turn when you turn the steering wheel? How big is that force? These are some of the questions we will answer in this chapter. The general topic is that of circular motion, that is, objects moving on a circular trajectory. Along the way we will learn about ferris wheels, cars turning without going off the road, etc. But first, we need to improve our math understanding by reviewing the product between vectors. --------------------------------------------------------------------------------------------------------------------------------------------------- MATH INSERT In chapter 1 we have introduced the concept of vectors and found how vectors can be added and or subtracted. We also mentioned the fact that vectors can be multiplied with each other in two different ways, what are called the scalar (or dot) product and the vector (or cross) product. In this chapter we will make use, for the first time, of the dot-product, so it is appropriate to take a look at it before we proceed. 1. The scalar (dot) product between vectors FIG. 1 Let’s consider two vectors, = + + and = + + . The scalar product between and , denoted (“A dot B”, hence the퐴⃗ term퐴푥 횤“̂ dot퐴 product푦횥̂ 퐴푧”),푘� is a 퐵�scalar⃗ 퐵 푥qu횤̂ antity퐵푦횥̂ equal퐵푧푘� to the product between the magnitudes퐴⃗ 퐵�⃗ of the two vectors, times the cosine of the angle between the two: 퐴⃗ ∙ 퐵�⃗ = cos , 휙 [1] where 퐴the⃗ ∙ 퐵� ⃗angle�퐴⃗ �� 퐵�is⃗� measured휙 by drawing the two vectors with the tails in the same point (see Fig.
    [Show full text]
  • 4. Central Forces
    4. Central Forces In this section we will study the three-dimensional motion of a particle in a central force potential. Such a system obeys the equation of motion mx¨ = V (r)(4.1) r where the potential depends only on r = x .Sincebothgravitationalandelectrostatic | | forces are of this form, solutions to this equation contain some of the most important results in classical physics. Our first line of attack in solving (4.1)istouseangularmomentum.Recallthatthis is defined as L = mx x˙ ⇥ We already saw in Section 2.2.2 that angular momentum is conserved in a central potential. The proof is straightforward: dL = mx x¨ = x V =0 dt ⇥ − ⇥r where the final equality follows because V is parallel to x. r The conservation of angular momentum has an important consequence: all motion takes place in a plane. This follows because L is a fixed, unchanging vector which, by construction, obeys L x =0 · So the position of the particle always lies in a plane perpendicular to L.Bythesame argument, L x˙ =0sothevelocityoftheparticlealsoliesinthesameplane.Inthis · way the three-dimensional dynamics is reduced to dynamics on a plane. 4.1 Polar Coordinates in the Plane We’ve learned that the motion lies in a plane. It will turn out to be much easier if we work with polar coordinates on the plane rather than Cartesian coordinates. For this reason, we take a brief detour to explain some relevant aspects of polar coordinates. To start, we rotate our coordinate system so that the angular momentum points in the z-direction and all motion takes place in the (x, y)plane.Wethendefinetheusual polar coordinates x = r cos ✓, y= r sin ✓ –48– Our goal is to express both the velocity and acceleration y ^ ^ θ r in polar coordinates.
    [Show full text]